Sequence and comparative analysis of the chicken genome provide unique perspectives on vertebrate evolution.

We present here a draft genome sequence of the red jungle fowl, Gallus gallus. Because the chicken is a modern descendant of the dinosaurs and the first non-mammalian amniote to have its genome sequenced, the draft sequence of its genome--composed of approximately one billion base pairs of sequence and an estimated 20,000-23,000 genes--provides a new perspective on vertebrate genome evolution, while also improving the annotation of mammalian genomes. For example, the evolutionary distance between chicken and human provides high specificity in detecting functional elements, both non-coding and coding. Notably, many conserved non-coding sequences are far from genes and cannot be assigned to defined functional classes. In coding regions the evolutionary dynamics of protein domains and orthologous groups illustrate processes that distinguish the lineages leading to birds and mammals. The distinctive properties of avian microchromosomes, together with the inferred patterns of conserved synteny, provide additional insights into vertebrate chromosome architecture.

GSVA: gene set variation analysis for microarray and RNA-seq data

Gene set enrichment (GSE) analysis is a popular framework for condensing information from gene expression profiles into a pathway or signature summary. The strengths of this approach over single gene analysis include noise and dimension reduction, as well as greater biological interpretability. As molecular profiling experiments move beyond simple case-control studies, robust and flexible GSE methodologies are needed that can model pathway activity within highly heterogeneous data sets. To address this challenge, we introduce Gene Set Variation Analysis (GSVA), a GSE method that estimates variation of pathway activity over a sample population in an unsupervised manner. We demonstrate the robustness of GSVA in a comparison with current state of the art sample-wise enrichment methods. Further, we provide examples of its utility in differential pathway activity and survival analysis. Lastly, we show how GSVA works analogously with data from both microarray and RNA-seq experiments. GSVA provides increased power to detect subtle pathway activity changes over a sample population in comparison to corresponding methods. While GSE methods are generally regarded as end points of a bioinformatic analysis, GSVA constitutes a starting point to build pathway-centric models of biology. Moreover, GSVA contributes to the current need of GSE methods for RNA-seq data. GSVA is an open source software package for R which forms part of the Bioconductor project and can be downloaded at http://www.bioconductor.org.

Molecular analysis of a novel gene cluster encoding an insect toxin in plant-associated strains of Pseudomonas fluorescens

Pseudomonas fluorescens CHA0 and the related strain Pf-5 are well-characterized representatives of rhizosphere bacteria that have the capacity to protect crop plants from fungal root diseases, mainly by releasing a variety of exoproducts that are toxic to plant pathogenic fungi. Here, we report that the two plant-beneficial pseudomonads also exhibit potent insecticidal activity. Anti-insect activity is linked to a novel genomic locus encoding a large protein toxin termed Fit (for P. fluorescensinsecticidal toxin) that is related to the insect toxin Mcf (Makes caterpillars floppy) of the entomopathogen Photorhabdus luminescens, a mutualist of insect-invading nematodes. When injected into the haemocoel, even low doses of P. fluorescens CHA0 or Pf-5 killed larvae of the tobacco hornworm Manduca sexta and the greater wax moth Galleria mellonella. In contrast, mutants of CHA0 or Pf-5 with deletions in the Fit toxin gene were significantly less virulent to the larvae. When expressed from an inducible promoter in a non-toxic Escherichia coli host, the Fit toxin gene was sufficient to render the bacterium toxic to both insect hosts. Our findings establish the Fit gene products of P. fluorescens CHA0 and Pf-5 as potent insect toxins that define previously unappreciated anti-insect properties of these plant-colonizing bacteria

Characterization and Control Strategies of an Integrated Chemical-Biological System for the Remediation of Toxic Pollutants in Wastewater: A case of study

In a previous work, a hybrid system consisting of an advanced oxidation process (AOP) named Photo-Fenton (Ph-F) and a fixed bed biological treatment operating as a sequencing batch biofilm reactor (SBBR) was started-up and optimized to treat 200 mg·L-1 of 4-chlorophenol (4-CP) as a model compound. In this work, studies of reactor stability and control as well as microbial population determination by molecular biology techniques were carried out to further characterize and control the biological reactor. Results revealed that the integrated system was flexible and even able to overcome toxic shock loads. Oxygen uptake rate (OUR) in situ was shown to be a valid tool to control the SBBR operation, to detect toxic conditions to the biomass, and to assess the recovery of performance. A microbial characterization by 16S rDNA sequence analysis reveals that the biological population was varied, although about 30% of the bacteria belonged to the Wautersia genus.

Sequencing and promoter analysis of the nifENXorf3orf5fdxAnifQ operon from Azospirillum brasilense Sp7

A 40-kb DNA region containing the major cluster of nif genes has been isolated from the Azospirillum brasilense Sp7 genome. In this region three nif operons have been identified: nifHDKorf1Y, nifENXorf3orf5fdxAnifQ and orf2nifUSVorf4. The operons containing nifENX and nifUSV genes are separated from the structural nifHDKorf1Y operon by about 5 kb and 10 kb, respectively. The present study shows the sequence analysis of the 6045-bp DNA region containing the nifENX genes. The deduced amino acid sequences from the open reading frames were compared to the nif gene products of other diazotrophic bacteria and indicate the presence of seven ORFs, all reading in the same direction as that of the nifHDKorf1Y operon. Consensus sigma54 and NifA-binding sites are present only in the promoter region upstream of the nifE gene. This promoter is activated by NifA protein and is approximately two-times less active than the nifH promoter, as indicated by the ß-galactosidase assays. This result suggests the differential expression of the nif genes and their respective products in Azospirillum.

Isolation and analysis of a corprinus cinereus DNA fragment showing homology to a fimbrial cDNA from ustilago violacea

Surface proteinaceous fibrils, termed fimbriae, were first identified on gram negative bacteria in the 1940s. Fungal fimbriae, discovered some 25 years later, are found on members of all fungal classes. In the present study, polyclonal antiserum raised against the fimbrial proteins of U. vio/acea were used in order to identify antigenically related proteins from Coprinus cinereus and Schizophy//um commune. Two polypeptides with molecular masses of 37 and 39 kDa from C. cinereus were observed and confirm earlier results. A single previously unidentified 50 kDa polypeptide in S. commune crossreacted with the antiserum. The 50 kDa protein was found to consist of 3 isoforms with isoelectric points ranging from 5.6 to 5.8. A fimbrial cDNA derived from U. vio/acea was used to identify DNA restriction fragments from C. cinereus and S. commune showing homology to the fimbrial transcript of U. vio/acea. Heterologous hybridization with this cDNA was used in order to screen a C. cinereus genomic DNA library. A single clone, A2-3A, with a 14 kbp insert showed strong homology to the pfim3-1 cDNA. The region of homology, a 700 bp Xba I fragment, was subcloned into pUG19. This plasmid was refered to as pXX8. DNA sequence determinations of pXX8 and adjacent fragments from A2-3A suggested that the cloned DNA was a portion of the rONA repeat encoding the small subunit rRNA. DNA sequence analysis of pfim3-1 yielded an incomplete open reading frame. The predicted amino acid sequence codes for a 206 amino acid, 22 kDa polypeptide which contains a domain similar to a transmembrane domain from rat leukocyte antigen, GDS3. As well, an untranslated 576 nucleotide domain showed 81 % homology to pXX8 and 830/0 homology to the 188 rRNA sequence of Ustilago maydis. This sequence was found adjacent to a region of adenine-thymine base pairs presumed to represent the polyadenylation sequence of the fimbrial transcript. The size and extent of homology is sufficient to account for the hybridization of pfim3-1 to rDNA. It is suggested that this domain represents a completely novel regulatory domain within eukaryotes that may enable the observed rapid regeneration of fimbriae in U. violacea.

Chemical and biological characterization of four new linear cationic α-helical peptides from the venoms of two solitary eumenine wasps

Four novel peptides were isolated from the venoms of the solitary eumenine wasps Eumenes rubrofemoratus and Eumenes fraterculus. Their sequences were determined by MALDI-TOF/TOF (matrix assisted laser desorption/ionization time-of-flight mass spectrometry) analysis, Edman degradation and solid-phase synthesis. Two of them, eumenitin-R (LNLKGLIKKVASLLN) and eumenitin-F (LNLKGLFKKVASLLT), are highly homologous to eumenitin, an antimicrobial peptide from a solitary eumenine wasp, whereas the other two, EMP-ER (FDIMGLIKKVAGAL-NH 2) and EMP-EF (FDVMGIIKKIAGAL-NH 2), are similar to eumenine mastoparan-AF (EMP-AF), a mast cell degranulating peptide from a solitary eumenine wasp. These sequences have the characteristic features of linear cationic cytolytic peptides; rich in hydrophobic and basic amino acids with no disulfide bond, and accordingly, they can be predicted to adopt an amphipathic α-helix secondary structure. In fact, the CD (circular dichroism) spectra of these peptides showed significant α-helical conformation content in the presence of TFE (trifluoroethanol), SDS (sodium dodecylsulfate) and asolectin vesicles. In the biological evaluation, all the peptides exhibited a significant broad-spectrum antimicrobial activity, and moderate mast cell degranulation and leishmanicidal activities, but showed virtually no hemolytic activity. © 2011 Elsevier Ltd.

Morphological, molecular and cross-breeding analysis of geographic populations of coconut-mite associated predatory mites identified as Neoseiulus baraki: evidence for cryptic species?

Surveys were conducted in Brazil, Benin and Tanzania to collect predatory mites as candidates for control of the coconut mite Aceria guerreronis Keifer, a serious pest of coconut fruits. At all locations surveyed, one of the most dominant predators on infested coconut fruits was identified as Neoseiulus baraki Athias-Henriot, based on morphological similarity with regard to taxonomically relevant characters. However, scrutiny of our own and published descriptions suggests that consistent morphological differences may exist between the Benin population and those from the other geographic origins. In this study, we combined three methods to assess whether these populations belong to one species or a few distinct, yet closely related species. First, multivariate analysis of 32 morphological characters showed that the Benin population differed from the other three populations. Second, DNA sequence analysis based on the mitochondrial cytochrome oxidase subunit I (COI) showed the same difference between these populations. Third, cross-breeding between populations was unsuccessful in all combinations. These data provide evidence for the existence of cryptic species. Subsequent morphological research showed that the Benin population can be distinguished from the others by a new character (not included in the multivariate analysis), viz. the number of teeth on the fixed digit of the female chelicera.

Decreasing the number of false positives in sequence classification

Abstract Background A large number of probabilistic models used in sequence analysis assign non-zero probability values to most input sequences. To decide when a given probability is sufficient the most common way is bayesian binary classification, where the probability of the model characterizing the sequence family of interest is compared to that of an alternative probability model. We can use as alternative model a null model. This is the scoring technique used by sequence analysis tools such as HMMER, SAM and INFERNAL. The most prevalent null models are position-independent residue distributions that include: the uniform distribution, genomic distribution, family-specific distribution and the target sequence distribution. This paper presents a study to evaluate the impact of the choice of a null model in the final result of classifications. In particular, we are interested in minimizing the number of false predictions in a classification. This is a crucial issue to reduce costs of biological validation. Results For all the tests, the target null model presented the lowest number of false positives, when using random sequences as a test. The study was performed in DNA sequences using GC content as the measure of content bias, but the results should be valid also for protein sequences. To broaden the application of the results, the study was performed using randomly generated sequences. Previous studies were performed on aminoacid sequences, using only one probabilistic model (HMM) and on a specific benchmark, and lack more general conclusions about the performance of null models. Finally, a benchmark test with P. falciparum confirmed these results. Conclusions Of the evaluated models the best suited for classification are the uniform model and the target model. However, the use of the uniform model presents a GC bias that can cause more false positives for candidate sequences with extreme compositional bias, a characteristic not described in previous studies. In these cases the target model is more dependable for biological validation due to its higher specificity.

Detailed analysis of sequence changes occurring during vlsE antigenic variation in the mouse model of Borrelia burgdorferi infection.

Lyme disease Borrelia can infect humans and animals for months to years, despite the presence of an active host immune response. The vls antigenic variation system, which expresses the surface-exposed lipoprotein VlsE, plays a major role in B. burgdorferi immune evasion. Gene conversion between vls silent cassettes and the vlsE expression site occurs at high frequency during mammalian infection, resulting in sequence variation in the VlsE product. In this study, we examined vlsE sequence variation in B. burgdorferi B31 during mouse infection by analyzing 1,399 clones isolated from bladder, heart, joint, ear, and skin tissues of mice infected for 4 to 365 days. The median number of codon changes increased progressively in C3H/HeN mice from 4 to 28 days post infection, and no clones retained the parental vlsE sequence at 28 days. In contrast, the decrease in the number of clones with the parental vlsE sequence and the increase in the number of sequence changes occurred more gradually in severe combined immunodeficiency (SCID) mice. Clones containing a stop codon were isolated, indicating that continuous expression of full-length VlsE is not required for survival in vivo; also, these clones continued to undergo vlsE recombination. Analysis of clones with apparent single recombination events indicated that recombinations into vlsE are nonselective with regard to the silent cassette utilized, as well as the length and location of the recombination event. Sequence changes as small as one base pair were common. Fifteen percent of recovered vlsE variants contained "template-independent" sequence changes, which clustered in the variable regions of vlsE. We hypothesize that the increased frequency and complexity of vlsE sequence changes observed in clones recovered from immunocompetent mice (as compared with SCID mice) is due to rapid clearance of relatively invariant clones by variable region-specific anti-VlsE antibody responses.

Transcriptional regulation and functional analysis of the human Wilms' tumor 1 gene

The Wilms' tumor 1 gene (WT1) encodes a zinc-finger transcription factor and is expressed in urogenital, hematopoietic and other tissues. It is expressed in a temporal and spatial manner in both embryonic and adult stages. To obtain a better understanding of the biological function of WT1, we studied two aspects of WT1 regulation: one is the identification of tissue-specific cis-regulatory elements that regulate its expression, the other is the downstream genes which are modulated by WT1.^ My studies indicate that in addition to the promoter, other regulatory elements are required for the tissue specific expression of this gene. A 259-bp hematopoietic specific enhancer in intron 3 of the WT1 gene increased the transcriptional activity of the WT1 promoter by 8- to 10-fold in K562 and HL60 cells. Sequence analysis revealed both GATA and c-Myb motifs in the enhancer fragment. Mutation of the GATA motif decreased the enhancer activity by 60% in K562 cells. Electrophoretic mobility shift assays showed that both GATA-1 and GATA-2 proteins in K562 nuclear extracts bind to this motif. Cotransfection of the enhancer containing reporter construct with a GATA-1 or GATA-2 expression vector showed that both GATA-1 and GATA-2 transactivated this enhancer, increasing the CAT reporter activity 10-15 fold and 5-fold respectively. Similar analysis of the c-Myb motif by cotransfection with the enhancer CAT reporter construct and a c-Myb expression vector showed that c-Myb transactivated the enhancer by 5-fold. A DNase I-hypersensitive site has been identified in the 258 bp enhancer region. These data suggest that GATA-1 and c-Myb are responsible for the activity of this enhancer in hematopoietic cells and may bind to the enhancer in vivo. In the process of searching for cis-regulatory elements in transgenic mice, we have identified a 1.0 kb fragment that is 50 kb downstream from the promoter and is required for the central nervous system expression of WT1.^ In the search for downstream target genes of WT1, we noted that the proto-oncogene N-myc is coexpressed with the tumor suppressor gene WT1 in the developing kidney and is overexpressed in many Wilms' tumors. Sequence analysis revealed eleven consensus WT1 binding sites located in the 1 kb mouse N-myc promoter. We further showed that the N-myc promoter was down-regulated by WT1 in transient transfection assays. Electrophoretic mobility shift assays showed that oligonucleotides containing the WT1 motifs could bind WT1 protein. Furthermore, a Denys-Drash syndrome mutant of WT1, R394W, that has a mutation in the DNA binding domain, failed to repress the N-myc promoter. This suggests that the repression of the N-myc promoter is mediated by DNA binding of WT1. This finding helps to elucidate the relationship of WT1 and N-myc in tumorigenesis and renal development. ^

Simultaneous analysis of nuclear and mitochondrial DNA, mRNA and miRNA from backspatter from inside parts of firearms generated by shots at "triple contrast" doped ballistic models

When a firearm projectile hits a biological target a spray of biological material (e.g., blood and tissue fragments) can be propelled from the entrance wound back towards the firearm. This phenomenon has become known as "backspatter" and if caused by contact shots or shots from short distances traces of backspatter may reach, consolidate on, and be recovered from, the inside surfaces of the firearm. Thus, a comprehensive investigation of firearm-related crimes must not only comprise of wound ballistic assessment but also backspatter analysis, and may even take into account potential correlations between these emergences. The aim of the present study was to evaluate and expand the applicability of the "triple contrast" method by probing its compatibility with forensic analysis of nuclear and mitochondrial DNA and the simultaneous investigation of co-extracted mRNA and miRNA from backspatter collected from internal components of different types of firearms after experimental shootings. We demonstrate that "triple contrast" stained biological samples collected from the inside surfaces of firearms are amenable to forensic co-analysis of DNA and RNA and permit sequence analysis of the entire mtDNA displacement-loop, even for "low template" DNA amounts that preclude standard short tandem repeat DNA analysis. Our findings underscore the "triple contrast" method's usefulness as a research tool in experimental forensic ballistics.

Identification and analysis of the plant peroxisomal targeting signal 1 receptor NtPEX5

Protein translocation into peroxisomes takes place via recognition of a peroxisomal targeting signal present at either the extreme C termini (PTS1) or N termini (PTS2) of matrix proteins. In mammals and yeast, the peroxisomal targeting signal receptor, Pex5p, recognizes the PTS1 consisting of -SKL or variants thereof. Although many plant peroxisomal matrix proteins are transported through the PTS1 pathway, little is known about the PTS1 receptor or any other peroxisome assembly protein from plants. We cloned tobacco (Nicotiana tabacum) cDNAs encoding Pex5p (NtPEX5) based on the protein’s interaction with a PTS1-containing protein in the yeast two-hybrid system. Nucleotide sequence analysis revealed that the tobacco Pex5p contains seven tetratricopeptide repeats and that NtPEX5 shares greater sequence similarity with its homolog from humans than from yeast. Expression of NtPEX5 fusion proteins, consisting of the N-terminal part of yeast Pex5p and the C-terminal region of NtPEX5, in a Saccharomyces cerevisiae pex5 mutant restored protein translocation into peroxisomes. These experiments confirmed the identity of the tobacco protein as a PTS1 receptor and indicated that components of the peroxisomal translocation apparatus are conserved functionally. Two-hybrid assays showed that NtPEX5 interacts with a wide range of PTS1 variants that also interact with the human Pex5p. Interestingly, the C-terminal residues of some of these peptides deviated from the established plant PTS1 consensus sequence. We conclude that there are significant sequence and functional similarities between the plant and human Pex5ps.

Phylogenetic analysis of “Volvocacae” for comparative genetic studies

Sequence analysis based on multiple isolates representing essentially all genera and species of the classic family Volvocaeae has clarified their phylogenetic relationships. Cloned internal transcribed spacer sequences (ITS-1 and ITS-2, flanking the 5.8S gene of the nuclear ribosomal gene cistrons) were aligned, guided by ITS transcript secondary structural features, and subjected to parsimony and neighbor joining distance analysis. Results confirm the notion of a single common ancestor, and Chlamydomonas reinharditii alone among all sequenced green unicells is most similar. Interbreeding isolates were nearest neighbors on the evolutionary tree in all cases. Some taxa, at whatever level, prove to be clades by sequence comparisons, but others provide striking exceptions. The morphological species Pandorina morum, known to be widespread and diverse in mating pairs, was found to encompass all of the isolates of the four species of Volvulina. Platydorina appears to have originated early and not to fall within the genus Eudorina, with which it can sometimes be confused by morphology. The four species of Pleodorina appear variously associated with Eudorina examples. Although the species of Volvox are each clades, the genus Volvox is not. The conclusions confirm and extend prior, more limited, studies on nuclear SSU and LSU rDNA genes and plastid-encoded rbcL and atpB. The phylogenetic tree suggests which classical taxonomic characters are most misleading and provides a framework for molecular studies of the cell cycle-related and other alterations that have engendered diversity in both vegetative and sexual colony patterns in this classical family.

Sequence of the FRA3B common fragile region: Implications for the mechanism of FHIT deletion

The hypothesis that chromosomal fragile sites may be “weak links” that result in hot spots for cancer-specific chromosome rearrangements was supported by the discovery that numerous cancer cell homozygous deletions and a familial translocation map within the FHIT gene, which encompasses the common fragile site, FRA3B. Sequence analysis of 276 kb of the FRA3B/FHIT locus and 22 associated cancer cell deletion endpoints shows that this locus is a frequent target of homologous recombination between long interspersed nuclear element sequences resulting in FHIT gene internal deletions, probably as a result of carcinogen-induced damage at FRA3B fragile sites.