Morphobiometric and molecular characterization of Bursaphelenchus Fuchs, 1937 (Nematoda: Aphelenchoididae) species associated with Pinus pinaster Ainton in Portugal

Abstract - The genus Bursaphelenchus comprises almost 100 species mainly from the northern hemisphere, with conifers as the most important hosts. Among the various nematode species, the pine wood nematode (PWN), Bursaphelenchus xylophilus, is the casual agent of pine wilt disease (PWD), and the most important forest pest for pines worldwide, classified as an A1 quarantine organism within the European Union. In 1999 this nematode was detected for the first time in Portugal and Europa associated with maritime pine, Pinus pinaster. Following detection, a national program denominated "Programa Nacional de Luta contra o Nemátodo da Madeira do Pinheiro" (PROLUNP) was created to, among other objectives, determine the distribution of the PWN and its associated vector(s) and host(s), and therefore intensive surveys covering the entire country were conducted with thousands of wood samples and suspected insects being analyzed. This thesis presents the listing, distribution, frequency and the insects associated with Bursaphelenchus species found associated with maritime pine in Portugal, identifying and characterizing the various species by morphological, biometrical and molecular biology (ITS-RFLP and rDNA sequencing analysis) techniques. To achieve the objectives, a total of 4813 maritime pine wood samples and 3294 insects from 22 species and six families were individually analyzed. A total of nine Bursaphelenchus species were found, namely: B. antoniae, B. hellenicus, B. leoni, B. mucronatus, B. pinasteri, B. sexdentati, B. teratospicularis, B. tusciae and B. xylophilus, all of them (with the exception of B. xylophilus) being new records for Portugal. Some of the species appear to have a widespread distribution, such as B. leoni, B. teratospicularis and B. tusciae while others were very rarely found and apparently have a localized distribution range within the country, namely B. antoniae and B. mucronatus. The majority of the species is characteristic of the Mediterranean region and can also be found in countries such as Spain, Italy and Greece, reflecting the affinity of our fauna with those locations. The association of B. hellenicus and B. tusciae with maritime pine is here reported for the first time. Six of the Bursaphelenchus species were also found associated with insects, mainly from the family Scolytidae (Coleoptera). Some of these interactions were described for the first time, namely: B. hellenicus with both Ips sexdentatus and Hylurgus ligniperda, B. sexdentati with both H. ligniperda and Orthotomicus erosus and B. tusciae with H. ligniperda. The exclusive association of B. xylophilus with the cerambycid Monochamus galloprovincialis was also confirmed. The nematode's dauer juveniles were usually found in low numbers in the insect vectors (ca 10-100 per insect), although for B. xylophilus a few thousand specimens per insect were sometimes found. The location of the dauer juveniles differed according to the species, although they were more common under the elytra and wings of the adult insects. A species new to science was detected and formally described as B. antoniae, associated with Hylobius sp. (Coleoptera; Curculionidae) beetles. Morphologically, this new species is very similar to B. hylobianum, although it's distinct ITS-RFLP molecular pattern (with only the enzyme Haelll producing comparable restriction bands) and the failure of hybridization supported the two species as distinct entities. Additional phylogenetic analysis of the 18S rDNA sequence further supported the taxonomical proximity of B. antoniae with B. hylobianum. Concerning the PWN, detailed studies on the development and morphology of B. xylophilus were conducted, and comparative measurements of field-collected and laboratory-maintained populations demonstrated that nematodes from the second group displayed larger size in all morphometric parameters, which could derive from more adequate conditions of nourishment and/or temperature. Taxonomical studies on the development stages of B. xylophilus confirmed the existence of four propagative juvenile stages (J1,J2,J3 and J4), an adult stage with both sexes and two dispersal stages (jIII e jIV), with the measurements of the gonad length allowing the separation of the propagative stages. It is hoped that the acquired knowledge will be useful on future surveys of nematodes of the Bursaphelenchus genus collected from either wood material or insect vectors, and facilitate the correct distinction and identification of the various species which are now known to occur. ### Resumo - 0 género Bursaphelenchus compreende quase 100 espécies, distribuídas sobretudo nos países do hemisfério norte do globo terrestre. Embora algumas espécies já tenham sido detectadas em plantas herbáceas, os hospedeiros vegetais mais comuns deste género são as coníferas, particularmente pinheiros. 0 nemátode da madeira do pinheiro (NMP), Bursaphelenchus xylophilus, é considerado a espécie mais importante deste género uma vez que é o agente causal da doença da murchidão dos pinheiros ("pine wilt disease"). Originário dos Estados Unidos, onde não causa grande impacte, o NMP foi introduzido em alguns países da Ásia (China, Japão, Coreia e Taiwan) e mais recentemente na Europa (Portugal). Nestas regiões é responsável pela destruição de milhares de hectares de coníferas, assumindo uma elevada importância económica. Em Portugal, depois da sua detecção em 1999, associado a Pinus pinaster, foi implementado um programa nacional "Programa Nacional de Luta contra o Nemátodo da Madeira do Pinheiro" (PROLUNP) que permitiu determinar a área afectada pela praga (a sul do rio Tejo, península de Setúbal) bem como definir e implementar estratégias de controlo e prevenção da disseminação do NMP a outras zonas de Portugal. Recentemente, em Junho de 2008, foi confirmada a presença de B. xylophilus em outras regiões de Portugal levando as autoridades oficiais a definir todo o território continental como zona afectada e de restrição. As prospecções intensivas realizadas nos últimos anos incluíram a recolha e análise de milhares de amostras de madeira de pinheiro bem como de insectos associados ao pinheiro bravo conduzindo à identificação de várias espécies de Bursaphelenchus. Assim, os estudos conduzidos neste trabalho tiveram como objectivos efectuar uma caracterização morfológica, biométrica e molecular das espécies associadas a P. pinaster em Portugal bem como a sua distribuição geográfica e abundância. Os estudos biométricos foram realizados com populações extraídas directamente do meio natural. Foi ainda realizada uma pesquisa que permitiu identificar os insectos a que estão associadas essas espécies, os seus possíveis vectores. Foram analisadas no total 4813 amostras de P. pinaster e 3294 insectos (22 espécies pertencentes a seis famílias diferentes). Foram identificadas um total de nove espécies: B. antoniae n. sp., B. hellenicus, B. leoni, B. mucronatus, B. pinasteri, B. sexdentati, B. teratospicularis, B. tusciae e B. xylophilus. Foram realizados estudos morfológicos e biométricos de todas as espécies com excepção de B. mucronatus; o reduzido número de exemplares encontrados em apenas uma amostra foram utilizados para efectuar o diagnóstico molecular desta espécie (ITS-RFLP). Apesar de ter havido, sempre que possível, a confirmação molecular, na maioria dos casos a caracterização morfológica e biométrica permitiu a correcta identificação das espécies. Contudo, foi imprescindível a análise molecular em algumas amostras, nomeadamente para a identificação de B. xylophilus e B. sexdentati; dada a grande semelhança entre B. xylophilus e B. mucronatus e tendo sido encontradas algumas populações de B. xylophilus que possuíam fêmeas com cauda mucronada, foi necessária a realização da confirmação molecular. Com excepção de B. xylophilus, todas as outras espécies foram reportadas pela primeira vez em Portugal. Juntamente com B. xylophilus, B. pinasteri foi a espécie encontrada nas amostras de madeira de pinheiro com maior frequência. Algumas destas espécies como B. leoni, B. teratospicularis e B. tusciae foram reportadas em diferentes localidades do norte, centro e sul de Portugal, apresentando uma vasta distribuição geográfica; este resultado está em consonância com a forte associação destas espécies a climas mediterrânicos tal como acontece em Espanha, França, Itália e Grécia. Em oposição, espécies como B. antoniae e B. mucronatus foram encontradas apenas numa ocasião na região centro (Leiria) e norte (Figueira da Foz) do país, respectivamente. Bursaphelenchus mucronatus é igualmente pouco frequente em Espanha onde ocorre sobretudo na região norte, na Galiza. Esta espécie preferirá climas mais frios, ocorrendo com uma maior frequência nas regiões de latitude norte; esta análise é corroborada pela presença constante em países como Alemanha, Finlândia, França, Noruega, Rússia e Suécia. A nível mundial são descritas neste trabalho pela primeira vez as associações das espécies B. hellenicus e B. tusciae ao hospedeiro vegetal P. pinaster.

NEW DATA ON THE CHRONOLOGY OF THE VALE DO FORNO SEDIMENTARY SEQUENCE (LOWER TAGUS RIVER TERRACE STAIRCASE) AND ITS RELEVANCE AS FLUVIAL ARCHIVE OF THE MIDDLE PLEISTOCENE IN WESTERN IBERIA

NEW DATA ON THE CHRONOLOGY OF THE VALE DO FORNO SEDIMENTARY SEQUENCE (LOWER TAGUS RIVER TERRACE STAIRCASE) AND ITS RELEVANCE AS FLUVIAL ARCHIVE OF THE MIDDLE PLEISTOCENE IN WESTERN IBERIA Pedro P. Cunha 1, António A. Martins 2, Jan-Pieter Buylaert 3,4, Andrew S. Murray 4, Luis Raposo 5, Paolo Mozzi 6, Martin Stokes 7 1 MARE - Marine and Environmental Sciences Centre, Department of Earth Sciences, University of Coimbra, Portugal: pcunha@dct.uc.pt 2 MARE - Marine and Environmental Sciences Centre, Dep. Geociências, University of Évora, Portugal; aam@uevora.pt 3 Centre for Nuclear Technologies, Technical University of Denmark, Risø Campus, Denmark; jabu@dtu.dk 4 Nordic Laboratory for Luminescence Dating, Aarhus University, Risø DTU, Denmark; anmu@dtu.dk 5 Museu Nacional de Arqueologia, Lisboa, Portugal; 3raposos@sapo.pt 6 Department of Geosciences, University of Padova, Italy; paolo.mozzi@unipd.it 7 School of Geography, Earth and Environmental Sciences, University of Plymouth, UK; m.stokes@plymouth.ac.uk The stratigraphic units that record the evolution of the Tagus River in Portugal (study area between Vila Velha de Ródão and Porto Alto villages; Fig. 1) have different sedimentary characteristics and lithic industries (Cunha et al., 2012): - a culminant sedimentary unit (the ancestral Tagus, before the drainage network entrenchment) – SLD13 (+142 to 262 m above river bed – a.r.b.; with probable age ca. 3,6 to 1,8 Ma), without artefacts; - T1 terrace (+84 to 180 m; ca. 1000? to 900 ka), without artefacts; - T2 terrace (+57 to 150 m; top deposits with a probable age ca. 600 ka), without artefacts; - T3 terrace (+43 to 113 m; ca. 460 to 360? ka), without artefacts; - T4 terrace (+26 to 55 m; ca. 335 a 155 ka), Lower Paleolithic (Acheulian) at basal and middle levels but early Middle Paleolithic at top levels; - T5 terrace (+5 to 34 m; 135 to 73 ka), Middle Paleolithic (Mousterian; Levallois technique); - T6 terrace (+3 to 14 m; 62 to 32 ka), late Middle Paleolithic (late Mousterian); - Carregueira Sands (aeolian sands) and colluvium (+3 a ca. 100 m; 32 to 12 ka), Upper Paleolithic to Epipaleolithic; - alluvial plain (+0 to 8 m; ca. 12 ka to present), Mesolithic and more recent industries. The differences in elevation (a.r.b.) of the several terrace staircases results from differential uplift due to active faults. Longitudinal correlation with the terrace levels indicates that a graded profile ca. 200 km long was achieved during terrace formation periods and a strong control by sea base level was determinant for terrace formation. The Neogene sedimentary units constituted the main source of sediments for the fluvial terraces (Fig. 2). Geomorphological mapping, coupled with lithostratigraphy, sedimentology and luminescence dating (quartz-OSL and K-feldspar post-IRIR290) were used in this study focused on the T4 terrace, which comprises a Lower Gravels (LG) unit and an Upper Sand (US) unit. The thick, coarse and dominantly massive gravels of the LG unit indicate deposition by a coarse bed-load braided river, with strong sediment supply, high gradient and fluvial competence, during conditions of rapidly rising sea level. Luminescence dating only provided minimum ages but it is probable that the LG unit corresponds to the earlier part of the MIS9 (ca. 335 to 325 ka), immediately postdating the incision promoted by the very low sea level (reaching ca. -140 m) during MIS10 (362 to 337 ka), a period of relatively cold climate conditions with weak vegetation cover on slopes and low sea level. Fig. 1. Main Portuguese reaches in which the Tagus River can be divided (Lower Tagus Basin): I – from the Spanish border to Arneiro (a general E–W trend, mainly consisting of polygonal segments); II – from Arneiro to Gavião (NE–SW); III – from Gavião to Arripiado (E–W); IV – from Arripiado to Vila Franca de Xira (NNE-SSW); V – from Vila Franca de Xira to the Atlantic shoreline. The faults considered to be the limit of the referred fluvial sectors are: F1 – Ponsul-Arneiro fault (WSW-ENE); F2 – Gavião fault (NW-SE); F3 – Ortiga fault (NW-SE); F4 – Vila Nova da Barquinha fault (W-E); F5 – Arripiado-Chamusca fault (NNE-SSW). 1 – estuary; 2 – terraces; 3 – faults; 4 – Tagus main channel. The main Iberian drainage basins are also represented (inset). The lower and middle parts of the US unit, comprising an alternation of clayish silts with paleosols and minor sands to the east (flood-plain deposits) and sand deposits to the west (channel belt), have a probable age of ca. 325 to 200 ka. This points to formation during MIS9 to MIS7, under conditions of high to medium sea levels and warm to mild conditions. The upper part of the US unit, dominated by sand facies and with OSL ages of ca. 200 to 154 ka, correlates with the early part of the MIS6. During this period, progradation resulted from climate deterioration and relative depletion of vegetation that promoted enhanced sediment production in the catchment, coupled with initiation of sea-level lowering that increased the longitudinal slope. The Vale do Forno and Vale da Atela archaeological sites (Alpiarça, central Portugal) document the earliest human occupation in the Lower Tagus River, well established in geomorphological and environmental terms, within the Middle Pleistocene. The Lower Palaeolithic sites were found on the T4 terrace (+26 m, a.r.b.). The oldest artefacts previously found in the LG unit, display crude bifacial forms that can be attributed to the Acheulian, with a probable age of ca. 335 to 325 ka. The T4 US unit has archaeological sites stratigraphically documenting successive phases of an evolved Acheulian, that probably date ca. 325 to 300 ka. Notably, these Lower Palaeolithic artisans were able to produce tools with different sophistication levels, simply by applying different strategies: more elaborated reduction sequences in case of bifaces and simple reduction sequences to obtain cleavers. Fig. 2. . Simplified geologic map of the Lower Tagus Cenozoic basin, adapted from the Carta Geológica de Portugal, 1/500000, 1992). The study area (comprising the Vale do Forno and Vale de Atela sites) is located on the more upstream sector of the Lower Tagus River reach IV, between Arripiado and Chamusca villages. 1 – alluvium (Holocene); 2 – terraces (Pleistocene); 3 – sands, silts and gravels (Paleogene to Pliocene); 4 – Sintra Massif (Cretaceous); 5 – limestones, marls, silts and sandstones (Mesozoic); 6 – quartzites (Ordovician); 7 – basement (Proterozoic to Palaeozoic); 8 – main fault. The main Portuguese reaches of the Tagus River are identified (I to V). The VF3 site (Milharós), containing a Final Acheulian industry, with fine and elaborated bifaces) found in a stratigraphic level located between the T4 terrace deposits and a colluvium associated with Late Pleistocene aeolian sands (32 to 12 ka), has an age younger than ca. 154 ka but much older than 32 ka. In the study area, the sedimentary units of the T4 terrace seem to record the river response to sea-level changes and climatically-driven fluctuations in sediment supply. REFERENCES Cunha P. P., Almeida N. A. C., Aubry T., Martins A. A., Murray A. S., Buylaert J.-P., Sohbati R., Raposo L., Rocha L., 2012, Records of human occupation from Pleistocene river terrace and aeolian sediments in the Arneiro depression (Lower Tejo River, central eastern Portugal). Geomorphology, vol. 165-166, pp. 78-90.

A new approach for computing a flood vulnerability index using cluster analysis

A Flood Vulnerability Index (FloodVI) was developed using Principal Component Analysis (PCA) and a new aggregation method based on Cluster Analysis (CA). PCA simplifies a large number of variables into a few uncorrelated factors representing the social, economic, physical and environmental dimensions of vulnerability. CA groups areas that have the same characteristics in terms of vulnerability into vulnerability classes. The grouping of the areas determines their classification contrary to other aggregation methods in which the areas' classification determines their grouping. While other aggregation methods distribute the areas into classes, in an artificial manner, by imposing a certain probability for an area to belong to a certain class, as determined by the assumption that the aggregation measure used is normally distributed, CA does not constrain the distribution of the areas by the classes. FloodVI was designed at the neighbourhood level and was applied to the Portuguese municipality of Vila Nova de Gaia where several flood events have taken place in the recent past. The FloodVI sensitivity was assessed using three different aggregation methods: the sum of component scores, the first component score and the weighted sum of component scores. The results highlight the sensitivity of the FloodVI to different aggregation methods. Both sum of component scores and weighted sum of component scores have shown similar results. The first component score aggregation method classifies almost all areas as having medium vulnerability and finally the results obtained using the CA show a distinct differentiation of the vulnerability where hot spots can be clearly identified. The information provided by records of previous flood events corroborate the results obtained with CA, because the inundated areas with greater damages are those that are identified as high and very high vulnerability areas by CA. This supports the fact that CA provides a reliable FloodVI.