Application of SNPs for population genetics of nonmodel organisms:new opportunities and challenges

Recent improvements in the speed, cost and accuracy of next generation sequencing are revolutionizing the discovery of single nucleotide polymorphisms (SNPs). SNPs are increasingly being used as an addition to the molecular ecology toolkit in nonmodel organisms, but their efficient use remains challenging. Here, we discuss common issues when employing SNP markers, including the high numbers of markers typically employed, the effects of ascertainment bias and the inclusion of nonneutral loci in a marker panel. We provide a critique of considerations specifically associated with the application and population genetic analysis of SNPs in nonmodel taxa, focusing specifically on some of the most commonly applied methods.

Production of extended-spectrum β-lactamases and the potential indirect pathogenic role of Prevotella isolates from the cystic fibrosis respiratory microbiota

Extended-spectrum β-lactamase (ESBL) production and the prevalence of the β-lactamase-encoding gene blaTEM were determined in Prevotella isolates (n=50) cultured from the respiratory tract of adults and young people with cystic fibrosis (CF). Time-kill studies were used to investigate the concept of passive antibiotic resistance and to ascertain whether a β-lactamase-positive Prevotella isolate can protect a recognised CF pathogen from the action of ceftazidime in vitro. The results indicated that approximately three-quarters (38/50; 76%) of Prevotella isolates produced ESBLs. Isolates positive for ESBL production had higher minimum inhibitory concentrations (MICs) of β-lactam antibiotics compared with isolates negative for production of ESBLs (P<0.001). The blaTEM gene was detected more frequently in CF Prevotella isolates from paediatric patients compared with isolates from adults (P=0.002), with sequence analysis demonstrating that 21/22 (95%) partial blaTEM genes detected were identical to blaTEM-116. Furthermore, a β-lactamase-positive Prevotella isolate protected Pseudomonas aeruginosa from the antimicrobial effects of ceftazidime (P=0.03). Prevotella isolated from the CF respiratory microbiota produce ESBLs and may influence the pathogenesis of chronic lung infection via indirect methods, including shielding recognised pathogens from the action of ceftazidime.

Bursaphelenchus antoniae sp. n. (Nematoda: Parasitaphelenchidae) associated with Hylobius sp. from Pinus pinaster in Portugal

Bursaphelenchus antoniae sp. n. is described and illustrated. Dauer juveniles were isolated from the body of the large pine weevil, Hylobius sp., collected from maritime pine (Pinus pinaster) stumps, in Portugal. Bursaphelenchus antoniae sp. n. was reared and maintained in P. pinaster wood segments and on Petri dish cultures of the fungi Botrytis cinerea and Monilinia fructicola. The new species is characterised by a relatively small body length of ca 583 μm (females) and 578 μm (males), a lateral field with two incisures, presence of a small vulval flap and a conoid female tail with a rounded or pointed terminus. Males have stout spicules with a disc-like cucullus and seven caudal papillae arranged as a single midventral precloacal papilla, one precloacal pair and two postcloacal pairs. In the character of the lateral field, B. antoniae sp. n. comes close to B. abietinus, B. rainulfi and B. hylobianum, whilst spicule characters place it within the piniperdae-group sensu Ryss et al. Morphologically, B. antoniae sp. n. is closest to B. hylobianum; the spicules of these two species having flattened, wing-like, alae on the distal third of the lamina. Bursaphelenchus antoniae sp. n. is distinguished from B. hylobianum on the arrangement of the caudal papillae (two vs three pairs). ITS-RFLP profiles and the failure to hybridise support the separation of the two species. Phylogenetic analysis of the new species, based on the 18S rDNA sequence, supports the inclusion of this new species in the B. hylobianum-group sensu Braasch. Sequence analysis of the 28S rDNA D2/D3 domain did not place the new species in a definite group.

Macrofaunal assemblages from mud volcanoes in the Gulf of Cadiz

Desde a descoberta em 1999 do primeiro vulcão de lama no Golfo de Cádis, cerca de 40 locais, de profundidade variável entre os 200 e os 3900 m, com diferentes graus de emissão de hidrocarbonetos foram localizados e amostrados dentro do programa IOC-UNESCO “Training Through Research (TTR) “ e mais recentemente dentro do projecto europeu HERMES. Neste estudo investigamos as comunidades da macrofauna dos vulcões de lama do Golfo de Cádis utilizando uma diversidade de equipamento de amostragem quantitativo e não quantitativo. Mais de 14550 espécimes foram examinados e incluídos nos diferentes grupos taxonómicos, sendo fornecida uma lista taxonómica detalhada com o menor nível taxonómico possível. A biodiversidade, distribuição dos principais taxa, as espécies quimiossintéticas e a biodiversidade regional e substituição de espécies são apresentados e discutidos. Dentro da macrofauna, os bivalves (nomeadamente super-familia Thyasiroidea, espécies quimisimbióticos e comunidade de bivalves) e os ofiurideos são estudados em pormenor. Os Thyasiroidea colhidos nos vulcões de lama do Golfo de Cádis são revistos. Das sete espécies identificadas, apenas uma Thyasira vulcolutre. sp. nov se encontra associada a um ambiente quimiossintético. Esta espécie é restrita a locais activos, mas não se verificam padrões de distribuição para as outras espécies. Os bivalves quimiosimbióticos amostrados são revistos. Das 10 espécies fortemente associadas a ambientes quimiossintéticos duas Solemyidae, Petrasma elarraichensis sp. nov. e Acharax gadirae sp. nov., uma Lucinidae, Lucinoma asapheus sp. nov., e uma Vesicomyidae, Isorropodon megadesmus sp. nov. são descritas e comparadas com similares das respectivas famílias. As comunidades de bivalves foram analisadas em detalhe e do estudo de 759 espécimes (49 espécies em 21 familias) descreve-se a diversidade e padrões de distribuição. Os Ophiuroidea amostrados nos vulcões de lama e ambientes batiais adjacentes são revistos. Treze espécies são incluídas em 4 famílias, Ophiacanthidae, Ophiactidae, Amphiuridae e Ophiuridae e são identificadas, tendo sido descrita uma nova espécie Ophiopristis cadiza sp. nov. Rácios isotópicos (δ13C, δ15N, δ34S) foram determinados em várias espécies no intuito de investigar a ecologia trófica das comunidades bênticas dos vulcões do Golfo de Cádis. Os valores de δ13C para os bivalves Solemyidae, Lucinidae e Thyasiridae estão de acordo com os valores para outros bivalves conhecidos por possuírem simbiontes tiotróficos. Por outro lado os valores de δ13C e δ34S para Bathymodiolus mauritanicus sugerem a ocorrência de metanotrofia. A análise da fauna heterotrófica indica igualmente que as espécies habitantes da cratera dos vulcões de lama derivam a sua nutrição de fontes quimiossintéticas. A indicação pela análise isotópica que as bactérias autotróficas contribuem substancialmente para a nutrição dos bivalves hospedeiros, levou-nos a investigar os endossimbiontes e as suas relações filogenéticas relativamente a outros bivalves através da análise comparativa de análises de sequências de 16S ribossomal RNS. Análises moleculares PCR-DGGE (Denaturing Gradient Gel Electrophoresis) e clonagem de genes de bacterias 16S rRNA confirmaram a presença de simbiontes oxidantes de enxofre e colocam a possibilidade de uma simbiose dupla para o B. mauritanicus. A diversidade microbiana dentro dos Frenulata foi igualmente estudada recorrendo a métodos moleculares e revelou a não existência de padrão entre espécies, vulcões, profundidade e idade do animal sugerindo assim a não procura de simbiontes específicos.

Lichen biodiversity and biomonitoring of atmospheric pollution

Esta tese debruça-se sobre a biodiversidade de líquenes epífitos de pinhais dunares portugueses e sobre uso de líquenes como biomonitores de poluição atmosférica nesse habitat. A Mata Nacional das Dunas de Quiaios (Figueira da Foz) foi o ponto de partida dos estudos de biodiversidade efetuados nesta tese, mas alguns deles estenderam-se à maior parte da costa portuguesa. Como resultado, encontrou-se uma espécie nova para a ciência, Lecanora sorediomarginata Rodrigues, Terrón & Elix, epifítica sobre Pinus pinaster Aiton e P. pinea L, que se encontra distribuída na maior parte da costa. Esta espécie caracteriza-se morfologicamente por um talo crustáceo, de cor esbranquiçada a acinzentada ou esverdeada e que desenvolve sorálios a partir de pequenas verrugas marginais. Quimicamente caracteriza-se pela presença dos ácidos 3,5-dicloro-2'-O-metilnorestenospórico [maior], 3,5-dicloro-2 -O-metilanziaico [menor], 3,5-dicloro-2 -O-metilnordivaricático [menor], 5-cloro-2'-Ometilanziaico [traço] e úsnico [traço]; atranorina [menor] e cloroatranorina [menor]. É quimicamente semelhante a L. lividocinerea Bagl., com a qual apresenta afinidades filogenéticas com base na análise da sequência ITS do rDNA, e a L. sulphurella Hepp. Adicionalmente, espécies Chrysothrix flavovirens Tønsberg e Ochrolechia arborea (Kreyer) Almb, também se encontraram epifíticas sobre P. pinaster e P. pinea em vários pinhais ao longo da costa, representando novos registos para a flora liquénica portuguesa, bem como a espécie Lepraria elobata Tønsberg encontrada epifítica sobre P. pinaster apenas nas Dunas de Quiaios. Além disso, as espécies Hypotrachyna lividescens (Kurok.) Hale e H. pseudosinuosa (Asahina) Hale encontraram-se epifíticas sobre P. pinaster e outros forófitos nas Dunas de Quiaios, constituindo novos registos para a flora liquénica da Península Ibérica. Estes resultados põe em evidência a importância dos pinhais dunares como habitat para líquenes epífitos. Num estudo conduzido entre janeiro e julho de 2008 num pinhal dunar (Mata do Urso, Figueira da Foz), em cuja bordadura existe uma fábrica de celulose de papel, usaram-se transplantes de líquenes da espécie Flavoparmelia caperata (L.) Hale para avaliar a acumulação de trinta e três elementos putativamente emitidos por fábricas de papel e pasta de papel. A cinética da fluorescência da clorofila a foi estudada nos líquenes transplantados, através da análise dos parâmetros Fv/Fm, F0, Fm, qP, NPQ, PSII, e Exc, de forma a avaliar os efeitos decorrentes da acumulação de elementos na vitalidade dos líquenes. Pretendeu-se avaliar se a acumulação de elementos e a cinética da fluorescência da clorofila a variavam significativamente com o local e o tempo de exposição, tendo em consideração os resultados obtidos de transplantes colocados num local de referência (Dunas de Quiaios) durante o mesmo período de tempo. (Continua no verso) resumo A maior parte dos elementos — Al, B, Ba, Ca, Co, Cr, Cu, Fe, Hg, Li, Mg, Mn, Mo, Na, Ni, P, S, Sb, Sc, Sr, Ti e V — ocorreu em concentrações significativamente mais elevadas nos transplantes expostos a 500 m da fábrica. Cerca de metade dos elementos estudados — B, Ba, Cr, Fe, Hg, Mg, Mn, Mo, Na, P, Pb, S, Sb e V — encontraram-se em concentrações significativamente mais elevadas nos transplantes expostos durante 180 dias. O solo foi identificado como uma fonte parcial da maior parte dos elementos. Os parâmetros Fv/Fm, Fm, PSII e Exc variaram significativamente com o local e/ou com o tempo de exposição. Observou-se um decréscimo significativo nos parâmetros Fv/Fm e Fm nos transplantes expostos a 500 e 1000 da fábrica, e também naqueles expostos durante 135 e 180 dias. Observou-se também um decréscimo significativo nos parâmetros PSII e Exc expostos durante 180 dias. Estes parâmetros correlacionaram-se de forma negativa e significativa com a acumulação de elementos: Fv/Fm: B, Ba, Co, Fe, Hg, Mg, Mn, Mo, N, P, S, Sb e Zn; Fm: Ba, Co, Hg, Mn, Mo, N, P, S, Sb e Zn; PSII: N e P; Exc: Mn, N, P e S. Estudos acerca da diversidade liquénica efetuados nos mesmos locais onde os transplantes foram colocados no local impactado, revelaram um menor valor de diversidade liquénica a 500 m da fábrica, que foi também o único local onde se encontraram espécies nitrófilas, o que se poderá dever à deposição de amónia e/ou poeiras. À semelhança de outros estudos, este trabalho confirma que os líquenes podem ser usados com sucesso em estudos de biomonitorização, mesmo em locais florestados. Além disso, traz também informações adicionais sobre como a acumulação de elementos pode influenciar a cinética da fluorescência da clorofila a em líquenes.

Molecular characterization of Citrus tristeza virus isolates from Cyprus on the basis of the coat protein gene

Within the context of a program in Cyprus for the control of Citrus tristeza virus (CTV), the coat protein (CP) genes of 12 local isolates of the virus that induced different symptoms on host trees, were compared to those of known isolates. The CP genes were reverse-transcribed (RT) and amplified by polymerase chain reaction (PCR) and the resulting amplicons were cloned and sequenced. Nucleotide sequence analysis revealed no signs of geographic speciation. All the sequences obtained clustered close to those of previously known isolates of worldwide origin that are in five distinct groups. The nucleotide diversity was high compared to that found using a worldwide database of CP gene sequences. These data support the existence of different CTV introductions into Cyprus or an introduction from a location in which CTV is relatively diverse. Some of the isolates induced stem pitting on branches of grapefruit and sweet orange. Such isolates have not been noted often in the Mediterranean basin. They were close in CP sequence to isolate B249 from Venezuela, which induces stem pitting, and are of particular concern for the whole region.

Lancefield group C and G streptococci in human infection : molecular typing, virulence and antimicrobial susceptibility

Tese de doutoramento, Ciências e Tecnologias da Saúde (Microbiologia), Universidade de Lisboa, Faculdade de Medicina, 2014

Characterizing the effects of ocean acidification in larval and juvenile Manila clam, Ruditapes philippinarum, using a transcriptomic approach

Ocean acidification as a result of anthropogenic carbon dioxide (CO2) emissions and global climate change poses a risk to the ecological landscape of intertidal and shallow subtidal communities. The organisms that inhabit these waters will have to cope with changing environmental conditions through the appropriate modulation of physiological processes. Calcifying organisms are particularly at risk, as increased atmospheric levels of CO2 in the atmosphere increase the partial pressure of CO2 (pCO2) in the oceans. Increased pCO2 reduces the saturation of carbonate minerals required to form calcified structures. Being able to cope with the increased energetic demand of maintaining these structures, in addition to other vital physiological processes, will be the key driver that determines which organisms will persist. Assessment of larval and juvenile Manila clam mortality and physiology in this study suggests that this species is capable of coping with elevated pCO2 conditions. The use of high throughput sequencing and RNA sequence analysis in larval clams revealed several physiological processes that play important roles in the Manila clam’s ability to tolerate elevated pCO2 conditions during this life stage. Exposure of juvenile Manila clams, acclimated to elevated pCO2 conditions, to a thermal stress revealed that this species might also be capable of coping with multiple stressors associated with global climate change. Manila clams could therefore represent a model for studying physiological mechanisms associated with successful acclimation of populations to ocean acidification.

Klebsiella pneumoniae subsp. pneumoniae–bacteriophage combination from the caecal effluent of a healthy woman

A recent study characterizing bacteriophage populations within human caecal effluent demonstrated the presence of numerous Podoviridae, Siphoviridae and Myoviridae within this material (Hoyles et al., 2014, Res Microbiol 165, 803–812). Further to this work, anaerobic bacteria were isolated on fastidious anaerobe agar from the caecal effluent of a healthy 31-year-old woman. Ten colonies were selected at random, streaked to purity and screened against the remaining caecal effluent (filter-sterilized, 0.45 μm pore size) in an attempt to isolate lytic bacteriophages. Bacteriophages within the effluent [2×105 ± 2.65×103 (n=3) pfu/ml] were active against five of the isolates, all identified by 16S rRNA gene sequence analysis as Klebsiella pneumoniae. One of the five isolates, L4-FAA5, was characterized further and found to be K. pneumoniae subsp. pneumoniae capsule type K2 rmpA+, and was used to propagate a bacteriophage (which we named KLPN1) to purity. Bacteriophage KLPN1 was a member of the Siphoviridae with a rosette-like tail tip and exhibited depolymerase activity, demonstrated by the formation of plaque-surrounding haloes that increased in size over the course of incubation. When screened against a panel of 21 clinical strains representing unknown K. pneumoniae subsp. pneumoniae capsule types and types K1, K2, K5, K20, K54 and K57, KLPN1 infected only K2 strains, but did not exhibit depolymerase activity against these. Whole-genome sequence analysis of KLPN1 showed the bacteriophage to have a genome of 49,037 bp (50.53 GC mol%) comprising 73 predicted ORFs, of which 22 encoded genes associated with structure, host recognition, packaging, DNA replication and cell lysis. The host recognition-associated gene was a potential depolymerase. This is the first report of the isolation of a bacterium–bacteriophage combination from the human caecum, and only the third member of the Siphoviridae known to infect K. pneumoniae subsp. pneumoniae.

Kroppenstedtia pulmonis sp. nov. and Kroppenstedtia sanguinis sp. nov., isolated from human patients

Three human clinical strains (W9323T , X0209T and X0394) isolated from lung biopsy, blood and cerebral spinal fluid, respectively, were characterized using a polyphasic taxonomic approach. Comparative analysis of the 16S rRNA gene sequences showed the three strains belonged to two novel branches within the genus Kroppenstedtia : 16S rRNA gene sequence analysis of W9323T showed closest sequence similarity to Kroppenstedtia eburnea JFMB- ATET (95.3 %), Kroppenstedtia guangzhouensis GD02T (94.7 %) and strain X0209T (94.6 %); sequence analysis of strain X0209T showed closest sequence similarity to K . eburnea JFMB- ATET (96.4 %) and K. guangzhouensis GD02T (96.0 %). Strains X0209T and X0394 were 99.9 % similar to each other by 16S rRNA gene sequence analysis. The DNA- DNA relatedness was 94.6 %, confirming that X0209T and X0394 belong to the same species. Chemotaxonomic data for strains W9323T and X0209T were consistent with those described for the genus Kroppenstedtia : whole- cell peptidoglycan contained LL- diaminopimelic acid; the major cellular fatty acids were iso- C15 and anteiso- C15 ; and the major menaquinone was MK- 7. Different endospore morphology, carbon utilization profiles, and whole cell wall sugar patterns of strains W9323T and X0209T supported by phylogenetic analysis enabled us to conclude that the strains represent two new species within the genus Kroppenstedtia , for which the names Kroppenstedtia pulmonis sp. nov. (type strain W9323T = DSM 45752T = CCUG 68107T) and Kroppenstedtia sanguinis sp. nov. (type strain X0209T = DSM 45749T = CCUG 38657T) are proposed.

Taxonomia bacteriana - aspectos atuais e perspectivas.

Métodos genotípicos. Determinação da razão de bases do DNA (% em moles de G + C); Hibridação DNA-DNA (HDD); Estudos filogenéticos com base no gene ribossomal 16S; Métodos de caracterização baseados nos polimorfismos de DNA (tipagem molecular); Métodos fenotípicos; Delineamento de uma espécie; Técnicas atuais na taxonomia bacteriana: utilização da genômica na taxonomia bacteriana; Utilização da metodologia de "Multilocus Sequence Analysis" - MLSA; Utilização da espectrometria de massa na identificação e classificação bacteriana.

Nuclear factor I-C links platelet-derived growth factor and transforming growth factor beta1 signaling to skin wound healing progression.

Transforming growth factor beta (TGF-beta) and platelet-derived growth factor A (PDGFAlpha) play a central role in tissue morphogenesis and repair, but their interplay remain poorly understood. The nuclear factor I C (NFI-C) transcription factor has been implicated in TGF-beta signaling, extracellular matrix deposition, and skin appendage pathologies, but a potential role in skin morphogenesis or healing had not been assessed. To evaluate this possibility, we performed a global gene expression analysis in NFI-C(-/-) and wild-type embryonic primary murine fibroblasts. This indicated that NFI-C acts mostly to repress gene expression in response to TGF-beta1. Misregulated genes were prominently overrepresented by regulators of connective tissue inflammation and repair. In vivo skin healing revealed a faster inflammatory stage and wound closure in NFI-C(-/-) mice. Expression of PDGFA and PDGF-receptor alpha were increased in wounds of NFI-C(-/-) mice, explaining the early recruitment of macrophages and fibroblasts. Differentiation of fibroblasts to contractile myofibroblasts was also elevated, providing a rationale for faster wound closure. Taken together with the role of TGF-beta in myofibroblast differentiation, our results imply a central role of NFI-C in the interplay of the two signaling pathways and in regulation of the progression of tissue regeneration.

Side effects of genome structural changes.

The first extensive catalog of structural human variation was recently released. It showed that large stretches of genomic DNA that vary considerably in copy number were extremely abundant. Thus it is conceivable that they play a major role in functional variation. Consistently, genomic insertions and deletions were shown to contribute to phenotypic differences by modifying not only the expression levels of genes within the aneuploid segments but also of normal copy-number neighboring genes. In this report, we review the possible mechanisms behind this latter effect.

Modeling sequencing errors by combining Hidden Markov models.

Among the largest resources for biological sequence data is the large amount of expressed sequence tags (ESTs) available in public and proprietary databases. ESTs provide information on transcripts but for technical reasons they often contain sequencing errors. Therefore, when analyzing EST sequences computationally, such errors must be taken into account. Earlier attempts to model error prone coding regions have shown good performance in detecting and predicting these while correcting sequencing errors using codon usage frequencies. In the research presented here, we improve the detection of translation start and stop sites by integrating a more complex mRNA model with codon usage bias based error correction into one hidden Markov model (HMM), thus generalizing this error correction approach to more complex HMMs. We show that our method maintains the performance in detecting coding sequences.

Results of a collaborative study of the EDNAP group regarding mitochondrial DNA heteroplasmy and segregation in hair shafts.

A collaborative exercise was carried out by the European DNA Profiling Group (EDNAP) in order to evaluate the distribution of mitochondrial DNA (mtDNA) heteroplasmy amongst the hairs of an individual who displays point heteroplasmy in blood and buccal cells. A second aim of the exercise was to study reproducibility of mtDNA sequencing of hairs between laboratories using differing chemistries, further to the first mtDNA reproducibility study carried out by the EDNAP group. Laboratories were asked to type 2 sections from each of 10 hairs, such that each hair was typed by at least two laboratories. Ten laboratories participated in the study, and a total of 55 hairs were typed. The results showed that the C/T point heteroplasmy observed in blood and buccal cells at position 16234 segregated differentially between hairs, such that some hairs showed only C, others only T and the remainder, C/T heteroplasmy at varying ratios. Additionally, differential segregation of heteroplasmic variants was confirmed in independent extracts at positions 16093 and the poly(C) tract at 302-309, whilst a complete A-G transition was confirmed at position 16129 in one hair. Heteroplasmy was observed at position 16195 on both strands of a single extract from one hair segment, but was not observed in the extracts from any other segment of the same hair. Similarly, heteroplasmy at position 16304 was observed on both strands of a single extract from one hair. Additional variants at positions 73, 249 and the HVII poly(C) region were reported by one laboratory; as these were not confirmed in independent extracts, the possibility of contamination cannot be excluded. Additionally, the electrophoresis and detection equipment used by this laboratory was different to those of the other laboratories, and the discrepancies at position 249 and the HVII poly(C) region appear to be due to reading errors that may be associated with this technology. The results, and their implications for forensic mtDNA typing, are discussed in the light of the biology of hair formation.