Navigating the labyrinth: a guide to sequence-based, community ecology of arbuscular mycorrhizal fungi

Data generated from next generation sequencing (NGS) will soon comprise the majority of information about arbuscular mycorrhizal fungal (AMF) communities. Although these approaches give deeper insight, analysing NGS data involves decisions that can significantly affect results and conclusions. This is particularly true for AMF community studies, because much remains to be known about their basic biology and genetics. During a workshop in 2013, representatives from seven research groups using NGS for AMF community ecology gathered to discuss common challenges and directions for future research. Our goal was to improve the quality and accessibility of NGS data for the AMF research community. Discussions spanned sampling design, sample preservation, sequencing, bioinformatics and data archiving. With concrete examples we demonstrated how different approaches can significantly alter analysis outcomes. Failure to consider the consequences of these decisions may compound bias introduced at each step along the workflow. The products of these discussions have been summarized in this paper in order to serve as a guide for any researcher undertaking NGS sequencing of AMF communities.

Transcriptome analysis for discovering candidate genes involve in embryogenesis in coconut (Cocos nucifera L.) through 454 pyrosequencing

Coconut, Cocos nucifera L. is a major plantation crop, which ensures income for millions of people in the tropical region. Detailed molecular studies on zygotic embryo development would provide valuable clues for the identification of molecular markers to improve somatic embryogenesis. Since there is no ongoing genome project for this species, coconut expressed sequence tags (EST) would be an interesting technique to identify important coconut embryo specific genes as well as other functional genes in different biochemical pathways. The goal of this study was to analyse the ESTs by examining the transcriptome data of the different embryo tissue types together with one somatic tissue. Here, four cDNA libraries from immature embryo, mature embryo, microspore derived embryo and mature leaves were constructed. cDNA was sequenced by the Roche-454 GS-FLX system and assembled into 32621 putative unigenes and 155017 singletons. Of these unigenes, 18651 had significant sequence similarities to non-redundant protein database, from which 16153 were assigned to one or more gene ontology categories. Homologue genes, which are responsible for embryo development such as chitinase, beta-1,3-glucanase, ATP synthase CF0 subunit, thaumatin-like protein and metallothionein-like protein were identified among the embryo EST collection. Of the unigenes, 6694 were mapped into 139 KEGG pathways including carbohydrate metabolism, energy metabolism, lipid metabolism, amino acid metabolism and nucleotide metabolism. This collection of 454-derived EST data generated from different tissue types provides a significant resource for genome wide studies and gene discovery of coconut, a non-model species.

Proteomic analysis of the medicinal plant Artemisia annua: data from leaf and trichome extracts

This article contains raw and processed data related to research published by Bryant et al. [1]. Data was obtained by MS-based proteomics, analysing trichome-enriched, trichome-depleted and whole leaf samples taken from the medicinal plant Artemisia annua and searching the acquired MS/MS data against a recently published contig database [2] and other genomic and proteomic sequence databases for comparison. The processed data shows that an order-of-magnitude more proteins have been identified from trichome-enriched Artemisia annua samples in comparison to previously published data. Proteins known to have a role in the biosynthesis of artemisinin and other highly abundant proteins were found which imply additional enzymatically driven processes occurring within the trichomes that are significant for the biosynthesis of artemisinin.

Light modulates the melanophore response to alpha-MSH in Xenopus laevis: An analysis of the signal transduction crosstalk mechanisms involved

Melanin granule (melanosome) dispersion within Xenopus laevis melanophores is evoked either by light or alpha-MSH. We have previously demonstrated that the initial biochemical steps of light and alpha-MSH signaling are distinct, since the increase in cAMP observed in response to alpha-MSH was not seen after light exposure. cAMP concentrations in response to alpha-MSH were significantly lower in cells pre-exposed to light as compared to the levels in dark-adapted melanophores. Here we demonstrate the presence of an adenylyl cyclase (AC) in the Xenopus melanophore, similar to the mammalian type IX which is inhibited by Ca(2+)-calmodulin-activated phosphatase. This finding supports the hypothesis that the cyclase could be negatively modulated by a light-promoted Ca(2+) increase. In fact, the activity of calcineurin PP2B phosphatase was increased by light, which could result in AC IX inhibition, thus decreasing the response to alpha-MSH. St-Ht31, a disrupting agent of protein kinase A (PKA)-anchoring kinase A protein (AKAP) complex totally blocked the melanosome dispersing response to alpha-MSH, but did not impair the photo-response in Xenopus melanophores. Sequence comparison of a melanophore AKAP partial clone with GenBank sequences showed that the anchoring protein was a gravin-like adaptor previously sequenced from Xenopus non-pigmentary tissues. Co-immunoprecipitation of Xenopus AKAP and the catalytic subunit of PKA demonstrated that PKA is associated with AKAP and it is released in the presence of alpha-MSH. We conclude that in X laevis melanophores, AKAP12 (gravin-like) contains a site for binding the inactive PKA thus compartmentalizing PKA signaling and also possesses binding sites for PKC. Light diminishes alpha-MSH-induced increase of cAMP by increasing calcineurin (PP2B) activity, which in turn inhibits adenylyl cyclase type IX, and/or by activating PKC, which phosphorylates the gravin-like molecule, thus destabilizing its binding to the cell membrane. (C) 2009 Elsevier Inc. All rights reserved.

Functionally conserved cis-regulatory elements of COL18A1 identified through zebrafish transgenesis

Type XVIII collagen is a component of basement membranes, and expressed prominently in the eye, blood vessels, liver, and the central nervous system. Homozygous mutations in COL18A1 lead to Knobloch Syndrome, characterized by ocular defects and occipital encephalocele. However, relatively little has been described on the role of type XVIII collagen in development, and nothing is known about the regulation of its tissue-specific expression pattern. We have used zebrafish transgenesis to identify and characterize cis-regulatory sequences controlling expression of the human gene. Candidate enhancers were selected from non-coding sequence associated with COL18A1 based on sequence conservation among mammals. Although these displayed no overt conservation with orthologous zebrafish sequences, four regions nonetheless acted as tissue-specific transcriptional enhancers in the zebrafish embryo, and together recapitulated the major aspects of col18a1 expression. Additional post-hoc computational analysis on positive enhancer sequences revealed alignments between mammalian and teleost sequences, which we hypothesize predict the corresponding zebrafish enhancers; for one of these, we demonstrate functional overlap with the orthologous human enhancer sequence. Our results provide important insight into the biological function and regulation of COL18A1, and point to additional sequences that may contribute to complex diseases involving COL18A1. More generally, we show that combining functional data with targeted analyses for phylogenetic conservation can reveal conserved cis-regulatory elements in the large number of cases where computational alignment alone falls short. (C) 2009 Elsevier Inc. All rights reserved.

Sequence and function of lysosomal and digestive cathepsin D-like proteinases of Musca domestica midgut

Musca domestica larvae display in anterior and middle midgut contents, a proteolytic activity with pH optimum of 3.0-3.5 and kinetic properties like cathepsin D. Three cDNAs coding for preprocathepsin D-like proteinases (ppCAD 1, ppCAD 2, ppCAD 3) were cloned from a M. domestica midgut cDNA library. The coded protein sequences included the signal peptide, propeptide and mature enzyme that has all conserved catalytic and substrate binding residues found in bovine lysosomal cathepsin D. Nevertheless, ppCAD 2 and ppCAD 3 lack the characteristic proline loop and glycosylation sites. A comparison among the sequences of cathepsin D-like enzymes from some vertebrates and those found in M. domestica and in the genomes of Aedes aegypti, Drosophila melanogaster, Tribolium castaneum, and Bombyx mori showed that only flies have enzymes lacking the proline loop (as defined by the motif: DxPxPx(G/A)P), thus resembling vertebrate pepsin. ppCAD 3 should correspond to the digestive cathepsin D-like proteinase (CAD) found in enzyme assays because: (1) it seems to be the most expressed CAD, based on the frequency of ESTs found. (2) The mRNA for CAD 3 is expressed only in the anterior and proximal middle midgut. (3) Recombinant procathepsin D-like proteinase (pCAD 3), after auto-activation has a pH optimum of 2.5-3.0 that is close to the luminal pH of M. domestica midgut. (4) Immunoblots of proteins from different tissues revealed with anti-pCAD 3 serum were positive only in samples of anterior and middle midgut tissue and contents. (5) CAD 3 is localized with immunogold inside secretory vesicles and around microvilli in anterior and middle midguit cells. The data support the view that on adapting to deal with a bacteria-rich food in an acid midgut region, M. domestica digestive CAD resulted from the same archetypical gene as the intracellular cathepsin D, paralleling what happened with vertebrates. The lack of the proline loop may be somehow associated with the extracellular role of both pepsin and digestive CAD 3. (C) 2009 Elsevier Ltd. All rights reserved.

The Genetics of Alzheimer`s Disease in Brazil: 10 Years of Analysis in a Unique Population

Alzheimer`s Disease (AD) is the most common type of dementia among the elderly, with devastating consequences for the patient, their relatives, and caregivers. More than 300 genetic polymorphisms have been involved with AD, demonstrating that this condition is polygenic and with a complex pattern of inheritance. This paper aims to report and compare the results of AD genetics studies in case-control and familial analysis performed in Brazil since our first publication, 10 years ago. They include the following genes/markers: Apolipoprotein E (APOE), 5-hidroxytryptamine transporter length polymorphic region (5-HTTLPR), brain-derived neurotrophin factor (BDNF), monoamine oxidase A (MAO-A), and two simple-sequence tandem repeat polymorphisms (DXS1047 and D10S1423). Previously unpublished data of the interleukin-1 alpha (IL-1 alpha) and interleukin-1 beta (IL-1 beta) genes are reported here briefly. Results from others Brazilian studies with AD patients are also reported at this short review. Four local families studied with various markers at the chromosome 21, 19, 14, and 1 are briefly reported for the first time. The importance of studying DNA samples from Brazil is highlighted because of the uniqueness of its population, which presents both intense ethnical miscegenation, mainly at the east coast, but also clusters with high inbreeding rates in rural areas at the countryside. We discuss the current stage of extending these studies using high-throughput methods of large-scale genotyping, such as single nucleotide polymorphism microarrays, associated with bioinformatics tools that allow the analysis of such extensive number of genetics variables, with different levels of penetrance. There is still a long way between the huge amount of data gathered so far and the actual application toward the full understanding of AD, but the final goal is to develop precise tools for diagnosis and prognosis, creating new strategies for better treatments based on genetic profile.

Croton laceratoglandulosus (Euphorbiaceae s.s.), a new glandular-stipulate species from Brazil and Bolivia, and its systematic position based on molecular analysis

Croton laceratoglandulosus, a new species from the dry forests of the Brazilian states of Piaui, Ceara, Bahia and Minas Gerais, and the Bolivian department of Santa Cruz, is described and illustrated here. Molecular sequence data demonstrate that it is most closely related to the taxa of Croton section Cascarilla, and not to sections Medea or Barhamia, which also have glandular calyces and laciniate stipules. (C) 2008 The Linnean Society of London, Botanical Journal of the Linnean Society, 2008, 158, 493-498.

Coding behavioural data for cladistic analysis: using dynamic homology without parsimony

Many of the controversies around the concept of homology rest on the subjectivity inherent to primary homology propositions. Dynamic homology partially solves this problem, but there has been up to now scant application of it outside of the molecular domain. This is probably because morphological and behavioural characters are rich in properties, connections and qualities, so that there is less space for conflicting character delimitations. Here we present a new method for the direct optimization of behavioural data, a method that relies on the richness of this database to delimit the characters, and on dynamic procedures to establish character state identity. We use between-species congruence in the data matrix and topological stability to choose the best cladogram. We test the methodology using sequences of predatory behaviour in a group of spiders that evolved the highly modified predatory technique of spitting glue onto prey. The cladogram recovered is fully compatible with previous analyses in the literature, and thus the method seems consistent. Besides the advantage of enhanced objectivity in character proposition, the new procedure allows the use of complex, context-dependent behavioural characters in an evolutionary framework, an important step towards the practical integration of the evolutionary and ecological perspectives on diversity. (C) The Willi Hennig Society 2010.

Genomic Analysis of Wild Tomato Introgressions Determining Metabolism- and Yield-Associated Traits

With the aim of determining the genetic basis of metabolic regulation in tomato fruit, we constructed a detailed physical map of genomic regions spanning previously described metabolic quantitative trait loci of a Solanum pennellii introgression line population. Two genomic libraries from S. pennellii were screened with 104 colocated markers from five selected genomic regions, and a total of 614 bacterial artificial chromosome (BAC)/cosmids were identified as seed clones. Integration of sequence data with the genetic and physical maps of Solanum lycopersicum facilitated the anchoring of 374 of these BAC/cosmid clones. The analysis of this information resulted in a genome-wide map of a nondomesticated plant species and covers 10% of the physical distance of the selected regions corresponding to approximately 1% of the wild tomato genome. Comparative analyses revealed that S. pennellii and domesticated tomato genomes can be considered as largely colinear. A total of 1,238,705 bp from both BAC/cosmid ends and nine large insert clones were sequenced, annotated, and functionally categorized. The sequence data allowed the evaluation of the level of polymorphism between the wild and cultivated tomato species. An exhaustive microsynteny analysis allowed us to estimate the divergence date of S. pennellii and S. lycopersicum at 2.7 million years ago. The combined results serve as a reference for comparative studies both at the macrosyntenic and microsyntenic levels. They also provide a valuable tool for fine-mapping of quantitative trait loci in tomato. Furthermore, they will contribute to a deeper understanding of the regulatory factors underpinning metabolism and hence defining crop chemical composition.

The R2 mobile element of Rhynchosciara americana: Molecular, cytological and dynamic aspects

Ribosomal RNA genes are encoded by large units clustered (18S, 5S, and 28S) in the nucleolar organizer region in several organisms. Sometimes additional insertions are present in the coding region for the 28S rDNA. These insertions are specific non-long terminal repeat retrotransposons that have very restricted integration targets within the genome. The retrotransposon present in the genome of Rhynchosciara americana, RaR2, was isolated by the screening of a genomic library. Sequence analysis showed the presence of conserved regions, such as a reverse transcriptase domain and a zinc finger motif in the amino terminal region. The insertion site was highly conserved in R. americana and a phylogenetic analysis showed that this element belongs to the R2 clade. The chromosomal localization confirmed that the RaR2 mobile element was inserted into a specific site in the rDNA gene. The expression level of RaR2 in salivary glands during larval development was determined by quantitative RT-PCR, and the increase of relative expression in the 3P of the fourth instar larval could be related to intense gene activity characteristic of this stage. 5`-Truncated elements were identified in different DNA samples. Additionally, in three other Rhynchosciara species, the R2 element was present as a full-length element.

Structural and phylogenetic relationship of ORF 31 from the Anticarsia gemmatalis MNPV to poly (ADP-ribose) polymerases (PARP)

ORF 31 is a unique baculovirus gene in the genome of Anticarsia gemmatalis multiple nucleopolyhedrovirus isolate 2D (AgMNPV-2D). It encodes a putative polypeptide of 369 aa homologous to poly (ADP-ribose) polymerase (PARP) found in the genomes of several organisms. Moreover, we found a phylogenetic association with Group I PARP proteins and a 3D homology model of its conserved PARP C-terminal catalytic domain indicating that had almost an exact spatial superimposition of < 1 angstrom with other PARP available structures. The 5` end of ORF 31 mRNA was located at the first nucleotide of a CATT motif at position -27. Using real-time PCR we detected transcripts at 3 h post-infection (p.i.) increasing until 24 h p.i., which coincides with the onset of DNA replication, suggestive of a possible role in DNA metabolism.

Analysis of the genome of Spodoptera frugiperda nucleopolyhedrovirus (SfMNPV-19) and of the high genomic heterogeneity in group II nucleopolyhedroviruses

The genome of the most virulent among 22 Brazilian geographical isolates of Spodoptera frugiperda nucleopolyhedrovirus, isolate 19 (SfMNPV-1 9), was completely sequenced and shown to comprise 132 565 bp and 141 open reading frames (ORFs). A total of 11 ORFs with no homology to genes in the GenBank database were found. Of those, four had typical baculovirus; promoter motifs and polyadenylation sites. Computer-simulated restriction enzyme cleavage patterns of SfMNPV-1 9 were compared with published physical maps of other SfMNPV isolates. Differences were observed in terms of the restriction profiles and genome size. Comparison of SfMNPV-1 9 with the sequence of the SfMNPV isolate 3AP2 indicated that they differed due to a 1427 bp deletion, as well as by a series of smaller deletions and point mutations. The majority of genes of SfMNPV-1 9 were conserved in the closely related Spodoptera exigua NPV (SeMNPV) and Agrotis segetum NPV (AgseMNPV-A), but a few regions experienced major changes and rearrangements. Synthenic maps for the genomes of group 11 NPVs revealed that gene collinearity was observed only within certain clusters. Analysis of the dynamics of gene gain and loss along the phylogenetic tree of the NPVs showed that group 11 had only five defining genes and supported the hypothesis that these viruses form ten highly divergent ancient lineages. Crucially, more than 60% of the gene gain events followed a power-law relation to genetic distance among baculoviruses, indicative of temporal organization in the gene accretion process.

Molecular analysis of a mutant Anticarsia gemmatalis multiple nucleopolyhedrovirus (AgMNPV) shows an interruption of an inhibitor of apoptosis gene (iap-3) by a new class-II piggyBac-related insect transposon

A new piggyBac-related transposable element (TE) was found in the genome of a mutant Anticarsia gemmatalis multiple nucleopolyhedrovirus interrupting an inhibitor of apoptosis gene. This mutant virus induces apoptosis upon infection of an Anticarsia gemmatalis cell line, but not in a Trichoplusia ni cell line. The sequence of the new TE (which was named IDT for iap disruptor transposon) has 2531 bp with two DNA sequences flanking a putative Transposase (Tpase) ORF of 1719 bp coding for a protein with 572 amino acids. These structural features are similar to the piggyBac TE, also reported for the first time in the genome of a baculovirus. We have also isolated variants of this new TE from different lepidopteran insect cells and compared their Tpase sequences.

High genetic diversity in field isolates of Trypanosoma theileri assessed by analysis of cathepsin L-like sequences disclosed multiple and new genotypes infecting cattle in Thailand

In this study, we describe the first survey in Thailand of Trypanosoma theileri, a widespread and prevalent parasite of cattle that is transmitted by tabanid flies. Investigation of 210 bovine blood samples of Thai cattle from six farms by hematocrit centrifuge technique (HCT) revealed 14 samples with trypanosomes morphologically compatible to T. theileri. Additional animals were positive for T. theileri by PCR based on the Cathepsin L-like sequence (TthCATL-PCR) despite negative by HCT, indicating cryptic infections. Results revealed a prevalence of 26 +/- 15% (95% CI) of T. theileri infection. Additionally, 12 samples positive for T. theileri were detected in cattle from other 11 farms. From a total of 30 blood samples positive by HCT and/or PCR from 17 farms, seven were characterized to evaluate the genetic polymorphism of T. theileri through sequence analysis of PCR-amplified CATL DNA sequences. All CATL sequences of T. theileri from Thai cattle clustered with sequences of the previously described phylogenetic lineages TthI and TthII, supporting only two major lineages of T. theileri in cattle around the world. However, 11 of the 29 CATL sequences analyzed showed to be different, disclosing an unexpectedly large polymorphic genetic repertoire, with multiple genotypes of T. theileri not previously described in other countries circulating in Thai cattle. (C) 2011 Elsevier B.V. All rights reserved.