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.

Towards the main sequence: detailed analysis of weak line and post-T Tauri stars

A detailed study was performed for a sample of low-mass pre-main-sequence (PMS) stars, previously identified as weak-line T Tauri stars, which are compared to members of the Tucanae and Horologium Associations. Aiming to verify if there is any pattern of abundances when comparing the young stars at different phases, we selected objects in the range from 1 to 100 Myr, which covers most of PMS evolution. High-resolution optical spectra were acquired at European Southern Observatory and Observatorio do Pico dos Dias. The stellar fundamental parameters effective temperature and gravity were calculated by excitation and ionization equilibria of iron absorption lines. Chemical abundances were obtained via equivalent width calculations and spectral synthesis for 44 per cent of the sample, which shows metallicities within 0.5 dex solar. A classification was developed based on equivalent width of Li I 6708 angstrom and Ha lines and spectral types of the studied stars. This classification allowed a separation of the sample into categories that correspond to different evolutive stages in the PMS. The position of these stars in the Hertzsprung-Russell diagram was also inspected in order to estimate their ages and masses. Among the studied objects, it was verified that our sample actually contains seven weak-line T Tauri stars, three are Classical T Tauri, 12 are Fe/Ge PMS stars and 21 are post-T Tauri or young main-sequence stars. An estimation of circumstellar luminosity was obtained using a disc model to reproduce the observed spectral energy distribution. Most of the stars show low levels of circumstellar emission, corresponding to less than 30 per cent of the total emission.

Generic relationships and dating of lineages in Winteraceae based on nuclear (ITS) and plastid (rpS16 and psbA-trnH) sequence data

Phylogenetic analyses of representative species from the five genera of Winteraceae (Drimys, Pseudowintera, Takhtajania, Tasmannia, and Zygogynum s.l.) were performed using ITS nuclear sequences and a combined data-set of ITS + psbA-trnH + rpS16 sequences (sampling of 30 and 15 species, respectively). Indel informativity using simple gap coding or gaps as a fifth character was examined in both data-sets. Parsimony and Bayesian analyses support the monophyly of Drimys, Tasmannia, and Zygogynum s.l., but do not support the monophyly of Belliolum, Zygogynum s.s., and Bubbia. Within Drimys, the combined data-set recovers two subclades. Divergence time estimates suggest that the splitting between Drimys and its sister clade (Pseudowintera + Zygogynum s.l.) occurred around the end of the Cretaceous; in contrast, the divergence between the two subclades within Drimys is more recent (15.5-18.5 MY) and coincides in time with the Andean uplift. Estimates suggest that the earliest divergences within Winteraceae could have predated the first events of Gondwana fragmentation. (C) 2009 Elsevier Inc. All rights reserved.

Relation of pretreatment sequence diversity in NS5A region of HCV genotype 1 with immune response between pegylated-INF/ribavirin therapy outcomes

Hepatitis C virus (HCV) infection frequently persists despite substantial virus-specific immune responses and the combination of pegylated interferon (INF)-alpha and ribavirin therapy. Major histocompatibility complex class I restricted CD8+ T cells are responsible for the control of viraemia in HCV infection, and several studies suggest protection against viral infection associated with specific HLAs. The reason for low rates of sustained viral response (SVR) in HCV patients remains unknown. Escape mutations in response to cytotoxic T lymphocyte are widely described; however, its influence in the treatment outcome is ill understood. Here, we investigate the differences in CD8 epitopes frequencies from the Los Alamos database between groups of patients that showed distinct response to pegylated alpha-INF with ribavirin therapy and test evidence of natural selection on the virus in those who failed treatment, using five maximum likelihood evolutionary models from PAML package. The group of sustained virological responders showed three epitopes with frequencies higher than Non-responders group, all had statistical support, and we observed evidence of selection pressure in the last group. No escape mutation was observed. Interestingly, the epitope VLSDFKTWL was 100% conserved in SVR group. These results suggest that the response to treatment can be explained by the increase in immune pressure, induced by interferon therapy, and the presence of those epitopes may represent an important factor in determining the outcome of therapy.

Mobius Sequence in a Girl and Arthrogryposis in her Half-Brother: Distinct Phenotypes Caused By Prenatal Injuries

Mobius sequence is a congenital facial and abducens nerve palsy, frequently associated to abnormalities of extremities. Arthrogryposis multiplex congenital is defined as a congenital fixation of multiple joints seldom of neurogenic origin. Both sequences must have a genetic origin, but usually are sporadic cases related to environmental factors such as drugs exposition and maternal trauma. A 5-year-old girl and a 1-year-old boy were born with Mobius sequence and arthrogryposis multiplex congenital, respectively. During pregnancies, the mother had vaginal bleeding at 7 weeks and used crack (free-based cocaine) in the first trimester, respectively. The girl also has equinovarus talipes and autistic behavior. The boy has arthrogryposis with flexion contractures of the feet and knees. A vascular disruption, due to hemorrhage and cocaine exposure, causing a transient ischemic insult to embryos in a critical period of development may be responsible for distinct phenotypes in these cases.

The DNA puff 4C expresses a salivary secretion protein in Trichosia pubescens (Diptera; Sciaridae)

DNA puffs are genomic regions of polytene chromosomes that undergo developmentally controlled DNA amplification and transcription in salivary glands of sciarid flies. Here, we tested the hypothesis that DNA puff genes code for salivary proteins in Trichosia pubescens. To do that, we generated antibodies against saliva and immunoscreened a cDNA library made from salivary glands. We isolated clones corresponding to DNA puff regions, including clone D-50 that contained the entire coding sequence of the previously isolated C4B1 gene from puff 4C. Indeed, we showed that puff 4C is a DNA puff region detecting its local transcription and its extra rounds of DNA incorporation compared to neighboring regions. We further confirmed D-50 clone identity in Western blots reacted with the anti-saliva anitiserum. We detected a recombinant protein expressed by this clone that had the expected size for a full-length product of the gene. We end with a discussion of the relationship between DNA puff genes and their products.

Bothrops insularis venomics: A proteomic analysis supported by transcriptomic-generated sequence data

A joint transcriptomic and proteomic approach employing two-dimensional electrophoresis, liquid chromatography and mass spectrometry was carried out to identify peptides and proteins expressed by the venom gland of the snake Bothrops insularis, an endemic species of Queimada Grande Island, Brazil. Four protein families were mainly represented in processed spots, namely metalloproteinase, serine proteinase, phospholipase A(2) and lectin. Other represented families were growth factors, the developmental protein G10, a disintegrin and putative novel bradykinin-potentiating peptides. The enzymes were present in several isoforms. Most of the experimental data agreed with predicted values for isoelectric point and M(r) of proteins found in the transcriptome of the venom gland. The results also support the existence of posttranslational modifications and of proteolytic processing of precursor molecules which could lead to diverse multifunctional proteins. This study provides a preliminary reference map for proteins and peptides present in Bothrops insularis whole venom establishing the basis for comparative studies of other venom proteomes which could help the search for new drugs and the improvement of venom therapeutics. Altogether, our data point to the influence of transcriptional and post-translational events on the final venom composition and stress the need for a multivariate approach to snake venomics studies. (c) 2009 Elsevier B.V. All rights reserved.

Mitochondrial sequence data and Dipsacales phylogeny: Mixed models, partitioned Bayesian analyses, and model selection

Phylogenetic analyses of chloroplast DNA sequences, morphology, and combined data have provided consistent support for many of the major branches within the angiosperm, clade Dipsacales. Here we use sequences from three mitochondrial loci to test the existing broad scale phylogeny and in an attempt to resolve several relationships that have remained uncertain. Parsimony, maximum likelihood, and Bayesian analyses of a combined mitochondrial data set recover trees broadly consistent with previous studies, although resolution and support are lower than in the largest chloroplast analyses. Combining chloroplast and mitochondrial data results in a generally well-resolved and very strongly supported topology but the previously recognized problem areas remain. To investigate why these relationships have been difficult to resolve we conducted a series of experiments using different data partitions and heterogeneous substitution models. Usually more complex modeling schemes are favored regardless of the partitions recognized but model choice had little effect on topology or support values. In contrast there are consistent but weakly supported differences in the topologies recovered from coding and non-coding matrices. These conflicts directly correspond to relationships that were poorly resolved in analyses of the full combined chloroplast-mitochondrial data set. We suggest incongruent signal has contributed to our inability to confidently resolve these problem areas. (c) 2007 Elsevier Inc. All rights reserved.

Transporters involved in glucose and water absorption in the Dysdercus peruvianus (Hemiptera: Pyrrhocoridae) anterior midgut

Little is known about insect intestinal sugar absorption, in spite of the recent findings, and even less has been published regarding water absorption. The aim of this study was to shed light on putative transporters of water and glucose in the insect midgut Glucose and water absorptions by the anterior ventriculus of Dysdercus peruvianus midgut were determined by feeding the insects with a glucose and a non-absorbable dye solution, followed by periodical dissection of insects and analysis of ventricular contents. Glucose absorption decreases glucose/dye ratios and water absorption increases dye concentrations. Water and glucose transports are activated (water 50%, glucose 33%) by 50 mM K(2)SO(4) and are inhibited (water 46%, glucose 82%) by 0.2 mM phloretin, the inhibitor of the facilitative hexose transporter (GLUT) or are inhibited (water 45%, glucose 35%) by 0.1 mM phlorizin, the inhibitor of the Na(+)-glucose cotransporter (SGLT). The results also showed that the putative SGLT transports about two times more water relative to glucose than the putative GLUT. These results mean that D. peruvianus uses a GLUT-like transporter and an SGLT-like transporter (with K(+) instead of Na(+)) to absorb dietary glucose and water. A cDNA library from D. peruvianus midgut was screened and we found one sequence homologous to GLUT1, named DpGLUT, and another to a sodium/solute symporter, named DpSGLT. Semi-quantitative RT-PCR studies revealed that DpGLUT and DpSGLTs mRNA were expressed in the anterior midgut, where glucose and water are absorbed, but not in fat body, salivary gland and Malpighian tubules. This is the first report showing the involvement of putative GLUT and SGLT in both water and glucose midgut absorption in insects. (C) 2010 Elsevier Inc. All rights reserved.

Expression of alpha-synuclein is increased in the hippocampus of rats with high levels of innate anxiety

A genomic region neighboring the alpha-synuclein gene, on rat chromosome 4, has been associated with anxiety- and alcohol-related behaviors in different rat strains. In this study, we have investigated potential molecular and physiological links between alpha-synuclein and the behavioral differences observed between Lewis (LEW) and Spontaneously Hypertensive (SHR) inbred rats, a genetic model of anxiety. As expected, LEW rats appeared more fearful than SHR rats in three anxiety models: open field, elevated plus maze and light/dark box. Moreover, LEW rats displayed a higher preference for alcohol and consumed higher quantities of alcohol than SHR rats. alpha-Synuclein mRNA and protein concentrations were higher in the hippocampus, but not the hypothalamus of LEW rats. This result inversely correlated with differences in dopamine turnover in the hippocampus of LEW and SHR rats, supporting the hypothesis that alpha-synuclein is important in the downregulation of dopamine neurotransmission. A novel single nucleotide polymorphism was identified in the 30-untranslated region (3`-UTR) of the alpha-synuclein cDNA between these two rat strains. Plasmid constructs based on the LEW 3`-UTR sequence displayed increased expression of a reporter gene in transiently transfected PC12 cells, in accordance with in-vivo findings, suggesting that this nucleotide exchange might participate in the differential expression of alpha-synuclein between LEW and SHR rats. These results are consistent with a novel role for alpha-synuclein in modulating rat anxiety- like behaviors, possibly through dopaminergic mechanisms. Since the behavioral and genetic differences between these two strains are the product of independent evolutionary histories, the possibility that polymorphisms in the alpha-synuclein gene may be associated with vulnerability to anxiety- related disorders in humans requires further investigation. Molecular Psychiatry (2009) 14, 894-905; doi: 10.1038/mp.2008.43; published online 22 April 2008

Cloning and characterization of a zebrafish homologue of human AQP1: a bifunctional water and gas channel

Chen LM, Zhao J, Musa-Aziz R, Pelletier MF, Drummond IA, Boron WF. Cloning and characterization of a zebrafish homologue of human AQP1: a bifunctional water and gas channel. Am J Physiol Regul Integr Comp Physiol 299: R1163-R1174, 2010. First published August 25, 2010; doi:10.1152/ajpregu.00319.2010.-The mammalian aquaporins AQP1, AQP4, and AQP5 have been shown to function not only as water channels but also as gas channels. Zebrafish have two genes encoding an AQP1 homologue, aqp1a and aqp1b. In the present study, we cloned the cDNA that encodes the zebrafish protein Aqp1a from the 72-h postfertilization (hpf) embryo of Danio rerio, as well as from the swim bladder of the adult. The deduced amino-acid sequence of aqp1a consists of 260 amino acids and is 59% identical to human AQP1. By analyzing the genomic DNA sequence, we identified four exons in the aqp1a gene. By in situ hybridization, aqp1a is expressed transiently in the developing vasculature and in erythrocytes from 16 to 48 h of development. Later, at 72 hpf, aqp1a is expressed in dermal ionocytes and in the swim bladder. Western blot analysis of adult tissues reveals that Aqp1a is most highly expressed in the eye and swim bladder. Xenopus oocytes expressing aqp1a have a channel-dependent (*) osmotic water permeability (P(f)*) that is indistinguishable from that of human AQP1. On the basis of the magnitude of the transient change in surface pH (Delta pHS) that were recorded as the oocytes were exposed to either CO(2) or NH(3), we conclude that zebrafish Aqp1a is permeable to both CO(2) and NH(3). The ratio (Delta pHS*)CO2/P(f)* is about half that of human AQP1, and the ratio (Delta pHS*)NH3/P(f)* is about one-quarter that of human AQP1. Thus, compared with human AQP1, zebrafish Aqp1a has about twice the selectivity for CO(2) over NH(3).