Ovule Transcriptome Analysis Discloses Deregulation of Genes and Pathways in Sexual and Apomictic Limonium Species (Plumbaginaceae)

The genus Limonium Mill. (sea lavenders) includes species with sexual and apomixis reproductive strategies, although the genes involved in these processes are unknown. To explore the mechanisms beyond these reproduction modes, transcriptome profiling of sexual, male sterile, and facultative apomictic species was carried out using ovules from different developmental stages. In total, 15,166 unigenes were found to be differentially expressed with apomictic vs. sexual reproduction, of which 4275 were uniquely annotated using an Arabidopsis thaliana database, with different regulations according to each stage and/or species compared. Gene ontology (GO) enrichment analysis indicated that genes related to tubulin, actin, the ubiquitin degradation process, reactive oxygen species scavenging, hormone signaling such as the ethylene signaling pathway and gibberellic acid-dependent signal, and transcription factors were found among differentially expressed genes (DEGs) between apomictic and sexual plants. We found that 24% of uniquely annotated DEGs were likely to be implicated in flower development, male sterility, pollen formation, pollen-stigma interactions, and pollen tube formation. The present study identifies candidate genes that are highly associated with distinct reproductive modes and sheds light on the molecular mechanisms of apomixis expression in Limonium sp.

Integrative analysis reveals the divergence and speciation between sister Sooty Copper butterflies Lycaena bleusei and L. tityrus

The comparison of closely related taxa is cornerstone in biology, as understanding mechanisms leading up to differentiation in relation to extant shared characters are powerful tools in interpreting the evolutionary process. Hotspots of biodiversity such as the west-Mediterranean, where many lineages meet are ideal grounds to study these processes. We set to explore the interesting example of Sooty Copper butterflies: widespread Eurasian Lycaena tityrus (Poda, 1761) comes into contact in Iberia with closely related and local endemic, L. bleusei (Oberthür, 1884), which hasn't always been considered a distinct species. An integrative analysis was designed, combining the use of extensive molecular data (five genes), geometric morphometrics analyses, verified and up-to-date distribution data, and environmental niche modelling, aimed at deciphering their true relationship, their placement within European Lycaena and trace their evolutionary history. We revealed several levels of differentiation: L. bleusei and L. tityrus appear to be reciprocally monophyletic independent gene-pools, distinct in all genes analysed, having mutually diverged 4.8 Ma ago. L. tityrus but not L. bleusei, further displays a genetic structure compatible with several glacial refugia, where populations assignable to infraspecific taxa surface. Conversely, L. bleusei shows a loss in mtDNA diversity in relation to nuDNA. Morphological analyses differentiate both species according to size and shape but also discriminate strong seasonal and sexual traits and a geographical phenotype segregation in L. tityrus. Finally, updated distribution and its modelling for current and glacial timeframes reveal both species respond differently to environmental variables, defining a mostly parapatric distribution and an overlapping belt where sympatry was recovered. During the last glacial maximum, a wider expansion in L. bleusei distribution explains current isolated populations. Our study highlights the importance of gathering several lines of evidence when deciphering the relationships between closely related populations in the fringe of cryptic species realm.

Overexpression of Water-Responsive Genes Promoted by Elevated CO2 Reduces ROS and Enhances Drought Tolerance in Coffea Species

Drought is a major constraint to plant growth and productivity worldwide and will aggravate as water availability becomes scarcer. Although elevated air [CO₂] might mitigate some of these effects in plants, the mechanisms underlying the involved responses are poorly understood in woody economically important crops such as Coffea. This study analyzed transcriptome changes in Coffea canephora cv. CL153 and C. arabica cv. Icatu exposed to moderate (MWD) or severe water deficits (SWD) and grown under ambient (aCO₂) or elevated (eCO₂) air [CO₂]. We found that changes in expression levels and regulatory pathways were barely affected by MWD, while the SWD condition led to a down-regulation of most differentially expressed genes (DEGs). eCO₂ attenuated the impacts of drought in the transcripts of both genotypes but mostly in Icatu, in agreement with physiological and metabolic studies. A predominance of protective and reactive oxygen species (ROS)-scavenging-related genes, directly or indirectly associated with ABA signaling pathways, was found in Coffea responses, including genes involved in water deprivation and desiccation, such as protein phosphatases in Icatu, and aspartic proteases and dehydrins in CL153, whose expression was validated by qRT-PCR. The existence of a complex post-transcriptional regulatory mechanism appears to occur in Coffea explaining some apparent discrepancies between transcriptomic, proteomic, and physiological data in these genotypes.

Testing drivers of acoustic divergence in cicadas (Cicadidae: Tettigettalna)

Divergence in acoustic signals may have a crucial role in the speciation process of animals that rely on sound for intra-specific recognition and mate attraction. The acoustic adaptation hypothesis (AAH) postulates that signals should diverge according to the physical properties of the signalling environment. To be efficient, signals should maximize transmission and decrease degradation. To test which drivers of divergence exert the most influence in a speciose group of insects, we used a phylogenetic approach to the evolution of acoustic signals in the cicada genus Tettigettalna, investigating the relationship between acoustic traits (and their mode of evolution) and body size, climate and micro-/macro-habitat usage. Different traits showed different evolutionary paths. While acoustic divergence was generally independent of phylogenetic history, some temporal variables' divergence was associated with genetic drift. We found support for ecological adaptation at the temporal but not the spectral level. Temporal patterns are correlated with micro- and macro-habitat usage and temperature stochasticity in ways that run against the AAH predictions, degrading signals more easily. These traits are likely to have evolved as an anti-predator strategy in conspicuous environments and low-density populations. Our results support a role of ecological selection, not excluding a likely role of sexual selection in the evolution of Tettigettalna calling songs, which should be further investigated in an integrative approach.

A first look at sea-lavenders genomics - can genome wide SNP information tip the scales of controversy in the Limonium vulgare species complex?

BACKGROUND

Sea-lavenders (Limonium Mill., Plumbaginaceae) are a cosmopolitan group of diploid and polyploid plants often adapted to extreme saline environments, with a mostly Tethyan distribution, occurring in the Mediterranean, Irano-Turanian, Euro-Siberian and in the New World. The halophylic Limonium vulgare polyploid complex in particular, presents a large distribution throughout extreme salt-marsh habitats and shows little morphological but high taximetric variation, frequently blurring species delimitation. In this work we pursue three main goals: assert whether SNP data from polyploid individuals has the resolution to distinguish the seven sampled species, to better understand how genetically structured Limonium vulgare is, and attempt to identify specific molecular mechanisms for the differentiation between L. maritimum and L. vulgare. For this purpose, 95 individuals were genotyped using Genotyping by Sequencing (GBS), which were assembled as two independent datasets using IPYRAD. All analyses performed downstream of assembly were fully automated. Phylogenetic inference, PCA, and admixture plots were used to infer answers to the study's main goals.

RESULTS

Close to 10,000 SNPs were obtained for each dataset. Phylogenetic analyses reveal that polyploid data can be used to infer species relationships. Population structure analyses suggest a genetically structured L. vulgare. A set of 34 SNPs were found to be fully segregated between L. vulgare and L. maritimum, two of which are potentially linked to proteins that might be involved in the speciation process.

CONCLUSION

Despite polyploid data analyses shortcomings, GBS generated SNPs have the resolution to discern all seven included species. Limonium vulgare revealed pronounced genetic structure along a geographical north-south cline. L. maritimum always appears as a distinct genetic entity. Segregated SNPs between L. vulgare and L. maritimum indicate salinity response and morphological trait control genes as potentially interesting to follow up for studying these species' divergence process.

Salt Stress Tolerance in Casuarina glauca: Insights from the Branchlets Transcriptome

Climate change and the accelerated rate of population growth are imposing a progressive degradation of natural ecosystems worldwide. In this context, the use of pioneer trees represents a powerful approach to reverse the situation. Among others, N₂-fixing actinorhizal trees constitute important elements of plant communities and have been successfully used in land reclamation at a global scale. In this study, we have analyzed the transcriptome of the photosynthetic organs of Casuarina glauca (branchlets) to unravel the molecular mechanisms underlying salt stress tolerance. For that, C. glauca plants supplied either with chemical nitrogen (KNO₃⁺) or nodulated by Frankia (NOD⁺) were exposed to a gradient of salt concentrations (200, 400, and 600 mM NaCl) and RNA-Seq was performed. An average of ca. 25 million clean reads was obtained for each group of plants, corresponding to 86,202 unigenes. The patterns of differentially expressed genes (DEGs) clearly separate two groups: (i) control- and 200 mM NaCl-treated plants, and (ii) 400 and 600 mM NaCl-treated plants. Additionally, although the number of total transcripts was relatively high in both plant groups, the percentage of significant DEGs was very low, ranging from 6 (200 mM NaCl/NOD⁺) to 314 (600 mM NaCl/KNO₃⁺), mostly involving down-regulation. The vast majority of up-regulated genes was related to regulatory processes, reinforcing the hypothesis that some ecotypes of C. glauca have a strong stress-responsive system with an extensive set of constitutive defense mechanisms, complemented by a tight mechanism of transcriptional and post-transcriptional regulation. The results suggest that the robustness of the stress response system in C. glauca is regulated by a limited number of genes that tightly regulate detoxification and protein/enzyme stability, highlighting the complexity of the molecular interactions leading to salinity tolerance in this species.

Worldwide Population Structure of the Coffee Rust Fungus Hemileia vastatrix Is Strongly Shaped by Local Adaptation and Breeding History

The devastating disease coffee leaf rust, caused by Hemileia vastatrix, has been a major constraint to worldwide coffee production. Recently, H. vastatrix populations were shown to be structured into three divergent genetic lineages with marked host specialization (C1, C2, and C3). However, there is yet no overall understanding of the population dynamics and adaptation of the most widespread and epidemiological relevant H. vastatrix group (C3). We used restriction site-associated DNA sequencing to generate 13,804 single nucleotide polymorphisms (SNPs) across a worldwide collection of 99 H. vastatrix isolates. Phylogenetic analyses uncovered a well-supported structuring within C3, with three main subgroups (SGs; SGI, SGII, and SGIII), which seem to reflect the historical distribution, breeding, and exchange of coffee cultivars. SGI shows a ladder-like diversification pattern and occurs across all four continents sampled, SGII is mainly restricted to Africa, and SGIII is observed only in Timor, revealing a higher genetic differentiation. Outlier and association tests globally identified 112 SNPs under putative positive selection, which impacted population structure. In particular, 29 overlapping SNPs per se seemed to have an extremely strong effect on H. vastatrix population divergence. We also found exclusive and fixed alleles associated with the SGs supporting local adaptation. Functional annotation revealed that transposable elements may play a role in host adaptation. Our study provides a higher-resolution perspective on the evolutionary history of H. vastatrix on cultivated coffee, showing its strong ability to adapt and the strength of the selective force imposed by coffee hosts, which should be taken into account when designing strategies for pathogen dissemination control and selective breeding.

Population structure in Quercus suber L. revealed by nuclear microsatellite markers

Quercus suber L. is a sclerophyllous tree species native to the western Mediterranean, a region that is considered highly vulnerable to increased temperatures and severe dry conditions due to environmental changes. Understanding the population structure and demographics of Q. suber is essential in order to anticipate whether populations at greater risk and the species as a whole have the genetic background and reproductive dynamics to enable rapid adaptation. The genetic diversity of Q. suber has been subject to different studies using both chloroplast and nuclear data, but population structure patterns remain unclear. Here, we perform genetic analyses on Q. suber using 13 nuclear microsatellite markers, and analysed 17 distinct locations across the entire range of the species. Structure analyses revealed that Q. suber may contain three major genetic clusters that likely result from isolation in refugia combined with posterior admixture and putative introgression from other Quercus species. Our results show a more complex structure scenario than previously inferred for Q. suber using nuclear markers and suggest that different southern populations contain high levels of genetic variation that may contribute to the resilience of Q. suber in a context of environmental change and adaptive pressure.

Global urban environmental change drives adaptation in white clover

Urbanization transforms environments in ways that alter biological evolution. We examined whether urban environmental change drives parallel evolution by sampling 110,019 white clover plants from 6169 populations in 160 cities globally. Plants were assayed for a Mendelian antiherbivore defense that also affects tolerance to abiotic stressors. Urban-rural gradients were associated with the evolution of clines in defense in 47% of cities throughout the world. Variation in the strength of clines was explained by environmental changes in drought stress and vegetation cover that varied among cities. Sequencing 2074 genomes from 26 cities revealed that the evolution of urban-rural clines was best explained by adaptive evolution, but the degree of parallel adaptation varied among cities. Our results demonstrate that urbanization leads to adaptation at a global scale.

Population structure, adaptation and divergence of the meadow spittlebug, Philaenus spumarius (Hemiptera, Aphrophoridae), revealed by genomic and morphological data

Understanding patterns of population differentiation and gene flow in insect vectors of plant diseases is crucial for the implementation of management programs of disease. We investigated morphological and genome-wide variation across the distribution range of the spittlebug Philaenus spumarius (Linnaeus, 1758) (Hemiptera, Auchenorrhyncha, Aphrophoridae), presently the most important vector of the plant pathogenic bacterium Xylella fastidiosa Wells et al., 1987 in Europe. We found genome-wide divergence between P. spumarius and a very closely related species, P. tesselatus Melichar, 1899, at RAD sequencing markers. The two species may be identified by the morphology of male genitalia but are not differentiated at mitochondrial COI, making DNA barcoding with this gene ineffective. This highlights the importance of using integrative approaches in taxonomy. We detected admixture between P. tesselatus from Morocco and P. spumarius from the Iberian Peninsula, suggesting gene-flow between them. Within P. spumarius, we found a pattern of isolation-by-distance in European populations, likely acting alongside other factors restricting gene flow. Varying levels of co-occurrence of different lineages, showing heterogeneous levels of admixture, suggest other isolation mechanisms. The transatlantic populations of North America and Azores were genetically closer to the British population analyzed here, suggesting an origin from North-Western Europe, as already detected with mitochondrial DNA. Nevertheless, these may have been produced through different colonization events. We detected SNPs with signatures of positive selection associated with environmental variables, especially related to extremes and range variation in temperature and precipitation. The population genomics approach provided new insights into the patterns of divergence, gene flow and adaptation in these spittlebugs and led to several hypotheses that require further local investigation.

Distribution of alien tetrapods in the Iberian Peninsula

We present a dataset that assembles occurrence records of alien tetrapods (amphibians, reptiles, birds and mammals) in the Iberian Peninsula, a coherent biogeographically unit where introductions of alien species have occurred for millennia. These data have important potential applications for ecological research and management, including the assessment of invasion risks, formulation of preventive and management plans, and research at the biological community level on alien species. This dataset summarizes inventories and data sources on the taxonomy and distribution of alien tetrapods in the Iberia Peninsula, comprising known locations from published literature, expert knowledge and citizen science platforms. An expert-based assessment process allowed the identification of unreliable records (misclassification or natural dispersion from native range), and the classification of species according to their status of reproduction in the wild. Distributional data was harmonized into a common area unit, the 10 × 10 km Universal Transverse Mercator (UTM) system (n = 6,152 cells). The year of observation and/or year of publication were also assigned to the records. In total, we assembled 35,940 unique distribution records (UTM × species × Year) for 253 species (6 amphibians, 16 reptiles, 218 birds and 13 mammals), spanning between 1912 and 2020. The species with highest number of distribution records were the Mediterranean painted frog Discoglossus pictus (n = 59 UTM), the pond slider Trachemys scripta (n = 471), the common waxbill Estrilda astrild (n = 1,275) and the house mouse Mus musculus (n = 4,043), for amphibians, reptiles, birds and mammals, respectively. Most alien species recorded are native to Africa (33%), followed by South America (21%), Asia (19%), North America (12%) and Oceania (10%). Thirty-six species are classified by IUCN as threatened in their native range, namely 2 Critically Endangered (CR), 6 Endangered (EN), 8 Vulnerable (VU), and 20 species Near Threatened (NT). Species maps are provided in DataSet1, as well R code and GIS layers to update them as new records are obtained.

Transcriptomic Leaf Profiling Reveals Differential Responses of the Two Most Traded Coffee Species to Elevated [CO2]

As atmospheric [CO₂] continues to rise to unprecedented levels, understanding its impact on plants is imperative to improve crop performance and sustainability under future climate conditions. In this context, transcriptional changes promoted by elevated CO₂ (eCO₂) were studied in genotypes from the two major traded coffee species: the allopolyploid Coffea arabica (Icatu) and its diploid parent, C. canephora (CL153). While Icatu expressed more genes than CL153, a higher number of differentially expressed genes were found in CL153 as a response to eCO₂. Although many genes were found to be commonly expressed by the two genotypes under eCO₂, unique genes and pathways differed between them, with CL153 showing more enriched GO terms and metabolic pathways than Icatu. Divergent functional categories and significantly enriched pathways were found in these genotypes, which altogether supports contrasting responses to eCO₂. A considerable number of genes linked to coffee physiological and biochemical responses were found to be affected by eCO₂ with the significant upregulation of photosynthetic, antioxidant, and lipidic genes. This supports the absence of photosynthesis down-regulation and, therefore, the maintenance of increased photosynthetic potential promoted by eCO₂ in these coffee genotypes.

Ancient divergence, a crisis of salt and another of ice shaped the evolution of the west Mediterranean butterfly Euchloe tagis

The Mediterranean region is an extremely complex hotspot where, since the Miocene, extensive geological, habitat and climate changes have taken place, alternating between warm and cold periods. These phenomena have taken a toll on the genetic composition of species, and surviving lineages have often adapted locally and diverged to the point of (complete) speciation. To study these phenomena, in this study we used one of the most enigmatic butterflies, the Portuguese dappled white, Euchloe tagis, a west Mediterranean endemic with fragmented, morphologically differentiated populations whose status have long been disputed. Even its affiliations with other Anthocharidini are largely unresolved. We used mitochondrial and nuclear markers under a phylogenetic and phylogeographical framework to evaluate its placement among relatives and population differentiation, reconstructing its evolutionary history. We found that this species had a Miocene origin ~15 Mya and was nearest to Euchloe s.s. and Elphinstonia. Its populations showed high genetic diversity but all coalesced to 5.3 Mya, when European and all but one African population diverged. Our multiple findings concerning the evolution of E. tagis through a changing, narrow habitable area might provide a more general perspective on how species survive within this hotspot of paramount importance.

Highly regional population structure of Spondyliosoma cantharus depicted by nuclear and mitochondrial DNA data

Resolution of population structure represents an effective way to define biological stocks and inform efficient fisheries management. In the present study, the phylogeography of the protogynous sparid Spondyliosoma cantharus, in the East Atlantic and Mediterranean Sea, was investigated with nuclear (S7) and mitochondrial (cytochrome b) DNA markers. Significant divergence of four regional genetic groups was observed: North Eastern Atlantic, Mediterranean Sea, Western African Transition (Cape Verde) and Gulf of Guinea (Angola). The two southern populations (Cape Verde and Angola) each comprised reciprocally monophyletic mtDNA lineages, revealed low levels of diversity in Cape Verde and high diversity for Angola despite being represented by only 14 individuals. A complete divergence between North Atlantic and Mediterranean populations was depicted by the mitochondrial marker, but a highly shared nuclear haplotype revealed an incomplete lineage sorting between these regions. Bayesian skyline plots and associated statistics revealed different dynamics among the four regions. Cape Verde showed no expansion and the expansion time estimated for Angola was much older than for the other regions. Mediterranean region seems to have experienced an early population growth but has remained with a stable population size for the last 30000 years while the North Atlantic population has been steadily growing. The lack of genetic structuring within these regions should not be taken as evidence of demographic panmixia in light of potential resolution thresholds and previous evidence of intra-regional phenotypic heterogeneity.

Standardised inventories of spiders (Arachnida, Araneae) of Macaronesia II: The native forests and dry habitats of Madeira archipelago (Madeira and Porto Santo islands)

Background

Here we present the data obtained from the samples collected as part of a large research project (MACDIV) which aims at understanding the drivers of spider (Araneae) community assembly in Macaronesian islands. To obtain the data, we applied the sampling protocol COBRA (Conservation Oriented Biodiversity Rapid Assessment), in twelve 50 m x 50 m native forest plots and five dry habitat plots on the island of Madeiraand in 5 dry habitat plots on the island of Porto Santo. Through this publication, we contribute to the knowledge of the arachnofauna of the Madeiran archipelago.

New information

From the samples that we collected, we obtained a total of 14,902 specimens, of which 49% were adults (7,263). We identified these specimens to 87 species and 18 morphospecies (undescribed), belonging to 26 families. Species of the family Linyphiidae dominated the samples, with 24 (morpho)species. Out of the 105 recorded (morpho)species, 34 were endemic, 26 native non-endemic, 22 introduced and 23 species of unknown origin. We report seven new records of possibly recently introduced species in the Madeiran archipelago. We also present 21 new records for Madeira island and 32 for Porto Santo (33 for the whole archipelago).

Population genomics of Bombus terrestris reveals high but unstructured genetic diversity in a potential glacial refugium

Ongoing climate change is expected to cause an increase in temperature and a reduction of precipitation levels in the Mediterranean region, which might cause changes in many species distributions. These effects negatively influence species gene pools, decreasing genetic variability and adaptive potential. Here, we use mitochondrial DNA and RADseq to analyse population genetic structure and genetic diversity of the bumblebee species Bombus terrestris (subspecies Bombus terrestris lusitanicus), in the Iberian Peninsula. Although this subspecies shows a panmictic pattern of population structure across Iberia and beyond, we found differentiation between subspecies B. t. lusitanicus and B. t. africanus, probably caused by the existence of barriers to gene flow between Iberia and North Africa. Furthermore, the results revealed that the Iberian Peninsula harbours a large fraction of B. terrestris intraspecific genetic variation, with the highest number of mitochondrial haplotypes found when compared with any other region in Europe studied so far, suggesting a potential role for the Iberian Peninsula as a glacial refugium. Our findings strengthen the idea that Iberia is a very important source of diversity for the global genetic pool of this species, because rare alleles might play a role in population resilience against human- or climate-mediated changes.

Shortcomings of Phylogenetic Studies on Recent Radiated Insular Groups: A Meta-Analysis Using Cabo Verde Biodiversity

Over the previous decades, numerous studies focused on how oceanic islands have contributed to determine the phylogenetic relationships and times of origin and diversification of different endemic lineages. The Macaronesian Islands (i.e., Azores, Madeira, Selvagens, Canaries, and Cabo Verde), harbour biotas with exceptionally high levels of endemism. Within the region, the vascular plants and reptiles constitute two of the most important radiations. In this study we compare relevant published phylogenetic data and diversification rates retrieved within Cabo Verde endemic lineages and discuss the importance of choosing appropriate phylogeny-based methods to investigate diversification dynamics on islands. From this selective literature-based review, we summarize the software packages used in Macaronesian studies and discuss their adequacy considering the published data to obtain well-supported phylogenies in the target groups. We further debate the importance of Next Generation Sequencing (NGS), to investigate the evolutionary processes of diversification in the Macaronesian Islands. Analysis of genomic data provides phylogenetic resolution for rapidly evolving species radiations, suggesting a great potential to improve the phylogenetic signal and divergence time estimates in insular lineages. The most important Macaronesian reptile radiations provide good case-studies to compare classical phylogenetic methods with new tools, such as phylogenomics, revealing a high value for research on this hotspot area.

Genomic signatures of introgression between commercial and native bumblebees, Bombus terrestris, in western Iberian Peninsula-Implications for conservation and trade regulation

Human-mediated introductions of species may have profound impacts on native ecosystems. One potential impact with largely unforeseen consequences is the potential admixture of introduced with autochthonous species through hybridization. Throughout the world, bumblebees have been deliberately introduced for crop pollination with known negative impacts on native pollinators. Given the likely allochthonous origin of commercial bumblebees used in Portugal (subspecies Bombus terrestris terrestris and B. t. dalmatinus), our aim was to assess their putative introgression with the native Iberian subspecies B. terrestris lusitanicus. We analysed one mitochondrial gene, cytochrome c oxidase subunit I (COX1) and genomic data involving thousands of genome-wide restriction-site-associated DNA markers (RAD-seq). In the mitochondrial COX1 analyses, we detected one relatively common haplotype in commercial bumblebees, also present in wild samples collected nearby the greenhouses where the commercial hives are used. In the RAD-seq analysis, we found a clear genetic differentiation between native and commercial lineages. Furthermore, we detected candidate hybrids in the wild, as well as putatively escaped commercial bumblebees, some of which being potentially fertile males. Although we cannot assess directly the fitness effects of introgressed alleles, there is a risk of maladaptive allele introgression to the local bumblebee subspecies, which can negatively impact autochthon populations. One immediate recommendation to farmers is for the proper disposal of hive boxes, after their use in greenhouses, so as to minimize the risk of escapees contaminating native populations. On the other hand, the feasibility of using local subspecies B. t. lusitanicus, preferably with local production, should be evaluated.

New insights into adaptation and population structure of cork oak using genotyping by sequencing

Species respond to global climatic changes in a local context. Understanding this process, including its speed and intensity, is paramount due to the pace at which such changes are currently occurring. Tree species are particularly interesting to study in this regard due to their long generation times, sedentarism, and ecological and economic importance. Quercus suber L. is an evergreen forest tree species of the Fagaceae family with an essentially Western Mediterranean distribution. Despite frequent assessments of the species' evolutionary history, large-scale genetic studies have mostly relied on plastidial markers, whereas nuclear markers have been used on studies with locally focused sampling strategies. In this work, "Genotyping by sequencing" is used to derive 1,996 single nucleotide polymorphism markers to assess the species' evolutionary history from a nuclear DNA perspective, gain insights into how local adaptation is shaping the species' genetic background, and to forecast how Q. suber may respond to global climatic changes from a genetic perspective. Results reveal (a) an essentially unstructured species, where (b) a balance between gene flow and local adaptation keeps the species' gene pool somewhat homogeneous across its distribution, but still allowing (c) variation clines for the individuals to cope with local conditions. "Risk of Non-Adaptedness" (RONA) analyses suggest that for the considered variables and most sampled locations, (d) the cork oak should not require large shifts in allele frequencies to survive the predicted climatic changes. Future directions include integrating these results with ecological niche modeling perspectives, improving the RONA methodology, and expanding its use to other species. With the implementation presented in this work, the RONA can now also be easily assessed for other organisms.

Genetic and Genomic Tools to Asssist Sugar Beet Improvement: The Value of the Crop Wild Relatives

Sugar beet (Beta vulgaris L. ssp. vulgaris) is one of the most important European crops for both food and sugar production. Crop improvement has been developed to enhance productivity, sugar content or other breeder's desirable traits. The introgression of traits from Crop Wild Relatives (CWR) has been done essentially for lessening biotic stresses constraints, namely using Beta and Patellifolia species which exhibit disease resistance characteristics. Several studies have addressed crop-to-wild gene flow, yet, for breeding programs genetic variability associated with agronomically important traits remains unexplored regarding abiotic factors. To accomplish such association from phenotype-to-genotype, screening for wild relatives occurring in habitats where selective pressures are in play (i.e., populations in salt marshes for salinity tolerance; populations subjected to pathogen attacks and likely evolved resistance to pathogens) are the most appropriate streamline to identify causal genetic information. By selecting sugar beet CWR species based on genomic tools, rather than random variations, is a promising but still seldom explored route toward the development of improved crops. In this perspective, a viable streamline for sugar beet improvement is proposed through the use of different genomic tools by recurring to sugar beet CWRs and focusing on agronomic traits associated with abiotic stress tolerance. Overall, identification of genomic and epigenomic landscapes associated to adaptive ecotypes, along with the cytogenetic and habitat characterization of sugar beet CWR, will enable to identify potential hotspots for agrobiodiversity of sugar beet crop improvement toward abiotic stress tolerance.

Morphology, songs and genetics identify two new cicada species from Morocco: <i>Tettigettalna</i> <i>afroamissa</i> sp. nov. and <i>Berberigetta</i> <i>dimelodica</i> gen. nov. &amp; sp. nov. (Hemiptera: Cicadettini)

Morocco has been the subject of very few expeditions on the last century with the objective of studying small cicadas. In the summer of 2014 an expedition was carried out to Morocco to update our knowledge with acoustic recordings and genetic data of these poorly known species. We describe here two new small-sized cicadas that could not be directly assigned to any species of North African cicadas: Tettigettalna afroamissa sp. nov. and Berberigetta dimelodica gen. nov. & sp. nov. In respect to T. afroamissa it is the first species of the genus to be found outside Europe and we frame this taxon within the evolutionary history of the genus. Acoustic analysis of this species allows us to confidently separate T. afroamissa from its congeners. With B. dimelodica, a small species showing a remarkable calling song characterized by an abrupt frequency modulation, a new genus had to be erected. Bayesian inference and maximum likelihood phylogenetic analyses with DNA-barcode sequences of Cytochrome C Oxidase 1 support the monophyly of both species, their distinctness and revealed genetic structure within B. dimelodica. Alongside the descriptions we also provide GPS coordinates of collection points, distributions and habitat preferences.

Deep analysis of wild Vitis flower transcriptome reveals unexplored genome regions associated with sex specification

RNA-seq of Vitis during early stages of bud development, in male, female and hermaphrodite flowers, identified new loci outside of annotated gene models, suggesting their involvement in sex establishment. The molecular mechanisms responsible for flower sex specification remain unclear for most plant species. In the case of V. vinifera ssp. vinifera, it is not fully understood what determines hermaphroditism in the domesticated subspecies and male or female flowers in wild dioecious relatives (Vitis vinifera ssp. sylvestris). Here, we describe a de novo assembly of the transcriptome of three flower developmental stages from the three Vitis vinifera flower types. The validation of de novo assembly showed a correlation of 0.825. The main goals of this work were the identification of V. v. sylvestris exclusive transcripts and the characterization of differential gene expression during flower development. RNA from several flower developmental stages was used previously to generate Illumina sequence reads. Through a sequential de novo assembly strategy one comprehensive transcriptome comprising 95,516 non-redundant transcripts was assembled. From this dataset 81,064 transcripts were annotated to V. v. vinifera reference transcriptome and 11,084 were annotated against V. v. vinifera reference genome. Moreover, we found 3368 transcripts that could not be mapped to Vitis reference genome. From all the non-redundant transcripts that were assembled, bioinformatics analysis identified 133 specific of V. v. sylvestris and 516 transcripts differentially expressed among the three flower types. The detection of transcription from areas of the genome not currently annotated suggests active transcription of previously unannotated genomic loci during early stages of bud development.

Genetic and Morphological Variation of the Forkbeard, Phycis phycis (Pisces, Phycidae): Evidence of Panmixia and Recent Population Expansion along Its Distribution Area

The knowledge of population structure of a species is essential to effectively assess and manage fisheries. In the present study, genetics, by mitochondrial DNA cytochrome b sequence analysis, and body geometric morphometrics were used to evaluate the existence of distinct populations of the forkbeard (Phycis phycis), an important commercial species in several European countries, especially Portugal and Spain. For geometric morphometric analysis, specimens were collected in the Northeast Atlantic Ocean—Azores, Madeira and mainland Portugal, and for genetic analysis, these samples were complemented with samples collected in the Mediterranean Sea—Spain, Italy and Croatia, in order to cover the entire distribution area of the species. Body shape of the forkbeard from the Northeast Atlantic was found to be highly variable. This variation was probably associated with different environmental factors between the study areas. Despite morphological variation, a low genetic differentiation between samples from different areas was found, most likely due to gene flow that occurred in the past or with the demographic history of the species. Moreover, the presence of unique haplotypes in the Northeast Atlantic and in the Mediterranean suggests that recent gene flow between populations from these areas should be limited. Altogether, a high haplotype diversity, a low nucleotide diversity, a “star-like” network and the results of the mismatch distribution, indicate a possible signature of recent population expansions, which probably started during the end of the Last Glacial Maximum and led to the colonization of the Northeast Atlantic and the Mediterranean.

Revisiting Vitis vinifera Subtilase Gene Family: A Possible Role in Grapevine Resistance against Plasmopara viticola

Subtilisin-like proteases, also known as subtilases, are a very diverse family of serine peptidases present in many organisms. In grapevine, there are hints of the involvement of subtilases in defense mechanisms, but their role is not yet understood. The first characterization of the subtilase gene family was performed in 2014. However, simultaneously, the grapevine genome was re-annotated and several sequences were re-annotated or retrieved. We have performed a re-characterization of this family in grapevine and identified 82 genes coding for 97 putative proteins, as result of alternative splicing. All the subtilases identified present the characteristic S8 peptidase domain and the majority of them also have a pro-domain I9 inhibitor, a protease-associated (PA) domain, and a signal peptide for targeting to the secretory pathway. Phylogenetic studies revealed six subtilase groups denominated VvSBT1 to VvSBT6. As several evidences have highlighted the participation of plant subtilases in response to biotic stimulus, we have investigated subtilase participation in grapevine resistance to Plasmopara viticola, the causative agent of downy mildew. Fourteen grapevine subtilases presenting either high homology to P69C from tomato, SBT3.3 from Arabidopsis thaliana or located near the Resistance to P. viticola (RPV) locus were selected. Expression studies were conducted in the grapevine-P. viticola pathosystem with resistant and susceptible cultivars. Our results may indicate that some of grapevine subtilisins are potentially participating in the defense response against this biotrophic oomycete.

Assessing genotype-phenotype associations in three dorsal colour morphs in the meadow spittlebug Philaenus spumarius (L.) (Hemiptera: Aphrophoridae) using genomic and transcriptomic resources

Background

Colour polymorphisms are common among animal species. When combined with genetic and ecological data, these polymorphisms can be excellent systems in which to understand adaptation and the molecular changes underlying phenotypic evolution. The meadow spittlebug, Philaenus spumarius (L.) (Hemiptera, Aphrophoridae), a widespread insect species in the Holarctic region, exhibits a striking dorsal colour/pattern balanced polymorphism. Although experimental crosses have revealed the Mendelian inheritance of this trait, its genetic basis remains unknown. In this study we aimed to identify candidate genomic regions associated with the colour balanced polymorphism in this species.

Results

By using restriction site-associated DNA (RAD) sequencing we were able to obtain a set of 1,837 markers across 33 individuals to test for associations with three dorsal colour phenotypes (typicus, marginellus, and trilineatus). Single and multi-association analyses identified a total of 60 SNPs associated with dorsal colour morphs. The genome size of P. spumarius was estimated by flow cytometry, revealing a 5.3 Gb genome, amongst the largest found in insects. A partial genome assembly, representing 24% of the total size, and an 81.4 Mb transcriptome, were also obtained. From the SNPs found to be associated with colour, 35% aligned to the genome and 10% to the transcriptome. Our data suggested that major loci, consisting of multi-genomic regions, may be involved in dorsal colour variation among the three dorsal colour morphs analysed. However, no homology was found between the associated loci and candidate genes known to be responsible for coloration pattern in other insect species. The associated markers showed stronger differentiation of the trilineatus colour phenotype, which has been shown previously to be more differentiated in several life-history and physiological characteristics as well. It is possible that colour variation and these traits are linked in a complex genetic architecture.

Conclusions

The loci detected to have an association with colour and the genomic and transcriptomic resources developed here constitute a basis for further research on the genetic basis of colour pattern in the meadow spittlebug P. spumarius.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-016-0455-5) contains supplementary material, which is available to authorized users.

Phylogeny of the most species-rich freshwater bivalve family (Bivalvia: Unionida: Unionidae): Defining modern subfamilies and tribes

Freshwater mussels of the order Unionida are key elements of freshwater habitats and are responsible for important ecological functions and services. Unfortunately, these bivalves are among the most threatened freshwater taxa in the world. However, conservation planning and management are hindered by taxonomic problems and a lack of detailed ecological data. This highlights the urgent need for advances in the areas of systematics and evolutionary relationships within the Unionida. This study presents the most comprehensive phylogeny to date of the larger Unionida family, i.e., the Unionidae. The phylogeny is based on a combined dataset of 1032bp (COI+28S) of 70 species in 46 genera, with 7 of this genera being sequenced for the first time. The resulting phylogeny divided the Unionidae into 6 supported subfamilies and 18 tribes, three of which are here named for the first time (i.e., Chamberlainiini nomen novum, Cristariini nomen novum and Lanceolariini nomen novum). Molecular analyses were complemented by investigations of selected morphological, anatomical and behavioral characters used in traditional phylogenetic studies. No single morphological, anatomical or behavioral character was diagnostic at the subfamily level and few were useful at the tribe level. However, within subfamilies, many tribes can be recognized based on a subset of these characters. The geographical distribution of each of the subfamilies and tribes is also presented. The present study provides important advances in the systematics of these extraordinary taxa with implications for future ecological and conservation studies.

Evolutionary and Biogeographic Insights on the Macaronesian Beta-Patellifolia Species (Amaranthaceae) from a Time-Scaled Molecular Phylogeny

The Western Mediterranean Region and Macaronesian Islands are one of the top biodiversity hotspots of Europe, containing a significant native genetic diversity of global value among the Crop Wild Relatives (CWR). Sugar beet is the primary crop of the genus Beta (subfamily Betoideae, Amaranthaceae) and despite the great economic importance of this genus, and of the close relative Patellifolia species, a reconstruction of their evolutionary history is still lacking. We analyzed nrDNA (ITS) and cpDNA gene (matK, trnH-psbA, trnL intron, rbcL) sequences to: (i) investigate the phylogenetic relationships within the Betoideae subfamily, and (ii) elucidate the historical biogeography of wild beet species in the Western Mediterranean Region, including the Macaronesian Islands. The results support the Betoideae as a monophyletic group (excluding the Acroglochin genus) and provide a detailed inference of relationships within this subfamily, revealing: (i) a deep genetic differentiation between Beta and Patellifolia species, which may have occurred in Late Oligocene; and (ii) the occurrence of a West-East genetic divergence within Beta, indicating that the Mediterranean species probably differentiated by the end of the Miocene. This was interpreted as a signature of species radiation induced by dramatic habitat changes during the Messinian Salinity Crisis (MSC, 5.96-5.33 Mya). Moreover, colonization events during the Pleistocene also played a role in shaping the current diversity patterns among and within the Macaronesian Islands. The origin and number of these events could not be revealed due to insufficient phylogenetic resolution, suggesting that the diversification was quite recent in these archipelagos, and unravelling potential complex biogeographic patterns with hybridization and gene flow playing an important role. Finally, three evolutionary lineages were identified corresponding to major gene pools of sugar beet wild relatives, which provide useful information for establishing in situ and ex situ conservation priorities in the hotspot area of the Macaronesian Islands.

Comparative Validation of Conventional and RNA-Seq Data-Derived Reference Genes for qPCR Expression Studies of Colletotrichum kahawae

Colletotrichum kahawae is an emergent fungal pathogen causing severe epidemics of Coffee Berry Disease on Arabica coffee crops in Africa. Currently, the molecular mechanisms underlying the Coffea arabica—C. kahawae interaction are still poorly understood, as well as the differences in pathogen aggressiveness, which makes the development of functional studies for this pathosystem a crucial step. Quantitative real time PCR (qPCR) has been one of the most promising approaches to perform gene expression analyses. However, proper data normalization with suitable reference genes is an absolute requirement. In this study, a set of 8 candidate reference genes were selected based on two different approaches (literature and Illumina RNA-seq datasets) to assess the best normalization factor for qPCR expression analysis of C. kahawae samples. The gene expression stability of candidate reference genes was evaluated for four isolates of C. kahawae bearing different aggressiveness patterns (Ang29, Ang67, Zim12 and Que2), at different stages of fungal development and key time points of the plant-fungus interaction process. Gene expression stability was assessed using the pairwise method incorporated in geNorm and the model-based method used by NormFinder software. For C. arabica—C. kahawae interaction samples, the best normalization factor included the combination of PP1, Act and ck34620 genes, while for C. kahawae samples the combination of PP1, Act and ck20430 revealed to be the most appropriate choice. These results suggest that RNA-seq analyses can provide alternative sources of reference genes in addition to classical reference genes. The analysis of expression profiles of bifunctional catalase-peroxidase (cat2) and trihydroxynaphthalene reductase (thr1) genes further enabled the validation of the selected reference genes. This study provides, for the first time, the tools required to conduct accurate qPCR studies in C. kahawae considering its aggressiveness pattern, developmental stage and host interaction.

Molecular candidates for early-stage flower-to-fruit transition in stenospermocarpic table grape (Vitis vinifera L.) inflorescences ascribed by differential transcriptome and metabolome profiles

Flower-to-fruit transition depends of nutrient availability and regulation at the molecular level by sugar and hormone signalling crosstalk. However, in most species, the identities of fruit initiation regulators and their targets are largely unknown. To ascertain the main pathways involved in stenospermocarpic table grape fruit set, comprehensive transcriptional and metabolomic analyses were conducted specifically targeting the early phase of this developmental stage in 'Thompson Seedless'. The high-throughput analyses performed disclosed the involvement of 496 differentially expressed genes and 28 differently accumulated metabolites in the sampled inflorescences. Our data show broad transcriptome reprogramming of molecule transporters, globally down-regulating gene expression, and suggest that regulation of sugar- and hormone-mediated pathways determines the downstream activation of berry development. The most affected gene was the SWEET14 sugar transporter. Hormone-related transcription changes were observed associated with increased indole-3-acetic acid, stimulation of ethylene and gibberellin metabolisms and cytokinin degradation, and regulation of MADS-box and AP2-like ethylene-responsive transcription factor expression. Secondary metabolism, the most representative biological process at transcriptome level, was predominantly repressed. The results add to the knowledge of molecular events occurring in grapevine inflorescence fruit set and provide a list of candidates, paving the way for genetic manipulation aimed at model research and plant breeding.

Shared and divergent pathways for flower abscission are triggered by gibberellic acid and carbon starvation in seedless Vitis vinifera L

BACKGROUND

Abscission is a highly coordinated developmental process by which plants control vegetative and reproductive organs load. Aiming at get new insights on flower abscission regulation, changes in the global transcriptome, metabolome and physiology were analyzed in 'Thompson Seedless' grapevine (Vitis vinifera L.) inflorescences, using gibberellic acid (GAc) spraying and shading as abscission stimuli, applied at bloom.

RESULTS

Natural flower drop rates increased from 63.1% in non-treated vines to 83% and 99% in response to GAc and shade treatments, respectively. Both treatments had a broad effect on inflorescences metabolism. Specific impacts from shade included photosynthesis inhibition, associated nutritional stress, carbon/nitrogen imbalance and cell division repression, whereas GAc spraying induced energetic metabolism simultaneously with induction of nucleotide biosynthesis and carbon metabolism, therefore, disclosing alternative mechanisms to regulate abscission. Regarding secondary metabolism, changes in flavonoid metabolism were the most represented metabolic pathways in the samples collected following GAc treatment while phenylpropanoid and stilbenoid related pathways were predominantly affected in the inflorescences by the shade treatment. However, both GAc and shade treated inflorescences revealed also shared pathways, that involved the regulation of putrescine catabolism, the repression of gibberellin biosynthesis, the induction of auxin biosynthesis and the activation of ethylene signaling pathways and antioxidant mechanisms, although often the quantitative changes occurred on specific transcripts and metabolites of the pathways.

CONCLUSIONS

Globally, the results suggest that chemical and environmental cues induced contrasting effects on inflorescence metabolism, triggering flower abscission by different mechanisms and pinpointing the participation of novel abscission regulators. Grapevine showed to be considered a valid model to study molecular pathways of flower abscission competence acquisition, noticeably responding to independent stimuli.

4Pipe4--A 454 data analysis pipeline for SNP detection in datasets with no reference sequence or strain information

Background

Next-generation sequencing datasets are becoming more frequent, and their use in population studies is becoming widespread. For non-model species, without a reference genome, it is possible from a panel of individuals to identify a set of SNPs that can be used for further population genotyping. However the lack of a reference genome to which the sequenced data could be compared makes the finding of SNPs more troublesome. Additionally when the data sources (strains) are not identified (e.g. in datasets of pooled individuals), the problem of finding reliable variation in these datasets can become much more difficult due to the lack of specialized software for this specific task.

Results

Here we describe 4Pipe4, a 454 data analysis pipeline particularly focused on SNP detection when no reference or strain information is available. It uses a command line interface to automatically call other programs, parse their outputs and summarize the results. The variation detection routine is built-in in the program itself. Despite being optimized for SNP mining in 454 EST data, it is flexible enough to automate the analysis of genomic data or even data from other NGS technologies. 4Pipe4 will output several HTML formatted reports with metrics on many of the most common assembly values, as well as on all the variation found. There is also a module available for finding putative SSRs in the analysed datasets.

Conclusions

This program can be especially useful for researchers that have 454 datasets of a panel of pooled individuals and want to discover and characterize SNPs for subsequent individual genotyping with customized genotyping arrays. In comparison with other SNP detection approaches, 4Pipe4 showed the best validation ratio, retrieving a smaller number of SNPs but with a considerably lower false positive rate than other methods.

4Pipe4’s source code is available at https://github.com/StuntsPT/4Pipe4.

Phylogeography and modes of reproduction in diploid and tetraploid halophytes of Limonium species (Plumbaginaceae): evidence for a pattern of geographical parthenogenesis

BACKGROUND AND AIMS

The genus Limonium (Plumbaginaceae) has long been recognized to have sexual and apomictic (asexual seed formation) modes of reproduction. This study aimed to elucidate phylogeographical patterns and modes of reproduction in diploid and tetraploid Limonium species, namely three putative sexual diploid species with morphological affinities (L. nydeggeri, L. ovalifolium, L. lanceolatum) and three related, probably apomict tetraploid species (L. binervosum, L. dodartii, L. multiflorum).

METHODS

cpDNA diversity and differentiation between natural populations of the species were investigated using two chloroplast sequence regions (trnL intron and trnL-trnF intergenic spacer). Floral heteromorphies, ovule cytoembryological analyses and pollination and crossing tests were performed in representative species of each ploidy group, namely diploid L. ovalifolium and tetraploid L. multiflorum, using plants from greenhouse collections.

KEY RESULTS AND CONCLUSIONS

Genetic analyses showed that diploid species have a higher haplotype diversity and a higher number of unique (endemic) haplotypes than tetraploid species. Network analysis revealed correlations between cpDNA haplotype distribution and ploidy groups, species groups and geographical origin, and haplotype sharing within and among species with distinct ploidy levels. Reproductive biology analyses showed that diploid L. ovalifolium mainly forms meiotically reduced tetrasporic embryo sacs of Gagea ova, Adoxa and Drusa types. Limonium multiflorum, however, has only unreduced, diplosporic (apomictic) embryo sacs of Rudbeckia type, and autonomous apomictic development seems to occur. Taken together, the findings provide evidence of a pattern of 'geographical parthenogenesis' in which quaternary climatic oscillations appear to be involved in the geographical patterns of coastal diploid and tetraploid Limonium species.

Genetic Diversity and Physiological Performance of Portuguese Wild Beet (Beta vulgaris spp. maritima) from Three Contrasting Habitats

The establishment of stress resilient sugar beets (Beta vulgaris spp. vulgaris) is an important breeding goal since this cash crop is susceptible to drought and salinity. The genetic diversity in cultivated sugar beets is low and the beet wild relatives are useful genetic resources for tolerance traits. Three wild beet populations (Beta vulgaris spp. maritima) from contrasting environments, Vaiamonte (VMT, dry inland hill), Comporta (CMP, marsh) and Oeiras (OEI, coastland), and one commercial sugar beet (Isella variety, SB), are compared. At the genetic level, the use of six microsatellite allowed to detect a total of seventy six alleles. It was observed that CMP population has the highest value concerning the effective number of alleles and of expected heterozygosity. By contrast, sugar beet has the lowest values for all the parameters considered. Loci analysis with STRUCTURE allows defining three genetic clusters, the sea beet (OEI and CMP), the inland ruderal beet (VMT) and the sugar beet (SB). A screening test for progressive drought and salinity effects demonstrated that: all populations were able to recover from severe stress; drought impact was higher than that from salinity; the impact on biomass (total, shoot, root) was population specific. The distinct strategies were also visible at physiological level. We evaluated the physiological responses of the populations under drought and salt stress, namely at initial stress stages, late stress stages, and early stress recovery. Multivariate analysis showed that the physiological performance can be used to discriminate between genotypes, with a strong contribution of leaf temperature and leaf osmotic adjustment. However, the separation achieved and the groups formed are dependent on the stress type, stress intensity and duration. Each of the wild beet populations evaluated is very rich in genetic terms (allelic richness) and exhibited physiological plasticity, i.e., the capacity to physiologically adjust to changing environments. These characteristics emphasize the importance of the wild beet ecotypes for beet improvement programs. Two striking ecotypes are VMT, which is the best to cope with drought and salinity, and CMP which has the highest root to shoot ratio. These genotypes can supply breeding programs with distinct goals.

Genomic Patterns of Positive Selection at the Origin of Rust Fungi

Understanding the origin and evolution of pathogenicity and biotrophic life-style of rust fungi has remained a conundrum for decades. Research on the molecular mechanisms responsible for rust fungi evolution has been hampered by their biotrophic life-style until the sequencing of some rust fungi genomes. With the availability of multiple whole genomes and EST data for this group, it is now possible to employ genome-wide surveys and investigate how natural selection shaped their evolution. In this work, we employed a phylogenomics approach to search for positive selection and genes undergoing accelerated evolution at the origin of rust fungi on an assembly of single copy genes conserved across a broad range of basidiomycetes. Up to 985 genes were screened for positive selection on the phylogenetic branch leading to rusts, revealing a pervasive signal of positive selection throughout the data set with the proportion of positively selected genes ranging between 19.6–33.3%. Additionally, 30 genes were found to be under accelerated evolution at the origin of rust fungi, probably due to a mixture of positive selection and relaxation of purifying selection. Functional annotation of the positively selected genes revealed an enrichment in genes involved in the biosynthesis of secondary metabolites and several metabolism and transporter classes.

Tettigettalnajosei (Boulard, 1982) (Hemiptera: Cicadoidea): first record in Spain, with notes on the distribution, genetic variation and behaviour of the species

The small cicada Tettigettalnajosei (Boulard, 1982) was until recently only known from southern Portugal and was considered endemic to this country. Fieldwork in 2013 led to the first record of the species in Spain, expanding its known eastern range to Andalusia. The northern limits remain poorly defined but it appears that the distribution of Tettigettalnajosei is restricted to the south Atlantic coastline in the Iberian Peninsula, with the highest densities found in Algarve. Some notes on behaviour and genetic variation of Tettigettalnajosei are also given.

New mitochondrial and nuclear evidences support recent demographic expansion and an atypical phylogeographic pattern in the spittlebug Philaenus spumarius (Hemiptera, Aphrophoridae)

Philaenus spumarius is a widespread insect species in the Holarctic region. Here, by focusing on the mtDNA gene COI but also using the COII and Cyt b genes and the nuclear gene EF-1α, we tried to explain how and when its current biogeographic pattern evolved by providing time estimates of the main demographic and evolutionary events and investigating its colonization patterns in and out of Eurasia. Evidence of recent divergence and expansion events at less than 0.5 Ma ago indicate that climate fluctuations in the Mid-Late Pleistocene were important in shaping the current phylogeographic pattern of the species. Data support a first split and differentiation of P. spumarius into two main mitochondrial lineages: the "western", in the Mediterranean region and the "eastern", in Anatolia/Caucasus. It also supports a following differentiation of the "western" lineage into two sub-lineages: the "western-Mediterranean", in Iberia and the "eastern-Mediterranean" in the Balkans. The recent pattern seems to result from postglacial range expansion from Iberia and Caucasus/Anatolia, thus not following one of the four common paradigms. Unexpected patterns of recent gene-flow events between Mediterranean peninsulas, a close relationship between Iberia and North Africa, as well as high levels of genetic diversity being maintained in northern Europe were found. The mitochondrial pattern does not exactly match to the nuclear pattern suggesting that the current biogeographic pattern of P. spumarius may be the result of both secondary admixture and incomplete lineage sorting. The hypothesis of recent colonization of North America from both western and northern Europe is corroborated by our data and probably resulted from accidental human translocations. A probable British origin for the populations of the Azores and New Zealand was revealed, however, for the Azores the distribution of populations in high altitude native forests is somewhat puzzling and may imply a natural colonization of the archipelago.

Oak root response to ectomycorrhizal symbiosis establishment: RNA-Seq derived transcript identification and expression profiling

Ectomycorrhizal symbiosis is essential for the life and health of trees in temperate and boreal forests where it plays a major role in nutrient cycling and in functioning of the forest ecosystem. Trees with ectomycorrhizal root tips are more tolerant to environmental stresses, such as drought, and biotic stresses such as root pathogens. Detailed information on these molecular processes is essential for the understanding of symbiotic tissue development in order to optimize the benefits of this natural phenomenon. Next generation sequencing tools allow the analysis of non model ectomycorrhizal plant-fungal interactions that can contribute to find the "symbiosis toolkits" and better define the role of each partner in the mutualistic interaction. By using 454 pyrosequencing we compared ectomycorrhizal cork oak roots with non-symbiotic roots. From the two cDNA libraries sequenced, over 2 million reads were obtained that generated 19,552 cork oak root unique transcripts. A total of 2238 transcripts were found to be differentially expressed when ECM roots were compared with non-symbiotic roots. Identification of up- and down-regulated gens in ectomycorrhizal roots lead to a number of insights into the molecular mechanisms governing this important symbiosis. In cork oak roots, ectomycorrhizal colonization resulted in extensive cell wall remodelling, activation of the secretory pathway, alterations in flavonoid biosynthesis, and expression of genes involved in the recognition of fungal effectors. In addition, we identified genes with putative roles in symbiotic processes such as nutrient exchange with the fungal partner, lateral root formation or root hair decay. These findings provide a global overview of the transcriptome of an ectomycorrhizal host root, and constitute a foundation for future studies on the molecular events controlling this important symbiosis.

A comprehensive assessment of the transcriptome of cork oak (Quercus suber) through EST sequencing

Background

Cork oak (Quercus suber) is one of the rare trees with the ability to produce cork, a material widely used to make wine bottle stoppers, flooring and insulation materials, among many other uses. The molecular mechanisms of cork formation are still poorly understood, in great part due to the difficulty in studying a species with a long life-cycle and for which there is scarce molecular/genomic information. Cork oak forests are of great ecological importance and represent a major economic and social resource in Southern Europe and Northern Africa. However, global warming is threatening the cork oak forests by imposing thermal, hydric and many types of novel biotic stresses. Despite the economic and social value of the Q. suber species, few genomic resources have been developed, useful for biotechnological applications and improved forest management.

Results

We generated in excess of 7 million sequence reads, by pyrosequencing 21 normalized cDNA libraries derived from multiple Q. suber tissues and organs, developmental stages and physiological conditions. We deployed a stringent sequence processing and assembly pipeline that resulted in the identification of ~159,000 unigenes. These were annotated according to their similarity to known plant genes, to known Interpro domains, GO classes and E.C. numbers. The phylogenetic extent of this ESTs set was investigated, and we found that cork oak revealed a significant new gene space that is not covered by other model species or EST sequencing projects. The raw data, as well as the full annotated assembly, are now available to the community in a dedicated web portal at http://www.corkoakdb.org.

Conclusions

This genomic resource represents the first trancriptome study in a cork producing species. It can be explored to develop new tools and approaches to understand stress responses and developmental processes in forest trees, as well as the molecular cascades underlying cork differentiation and disease response.

Overview of the functional virulent genome of the coffee leaf rust pathogen Hemileia vastatrix with an emphasis on early stages of infection

Hemileia vastatrix is the causal agent of coffee leaf rust, the most important disease of coffee Arabica. In this work, a 454-pyrosequencing transcriptome analysis of H. vastatrix germinating urediniospores (gU) and appressoria (Ap) was performed and compared to previously published in planta haustoria-rich (H) data. A total of 9234 transcripts were identified and annotated. Ca. 50% of these transcripts showed no significant homology to international databases. Only 784 sequences were shared by the three conditions, and 75% were exclusive of either gU (2146), Ap (1479) or H (3270). Relative transcript abundance and RT-qPCR analyses for a selection of genes indicated a particularly active metabolism, translational activity and production of new structures in the appressoria and intense signaling, transport, secretory activity and cellular multiplication in the germinating urediniospores, suggesting the onset of a plant-fungus dialogue as early as at the germ tube stage. Gene expression related to the production of carbohydrate-active enzymes and accumulation of glycerol in germinating urediniospores and appressoria suggests that combined lytic and physical mechanisms are involved in appressoria-mediated penetration. Besides contributing to the characterization of molecular processes leading to appressoria-mediated infection by rust fungi, these results point toward the identification of new H. vastatrix candidate virulence factors, with 516 genes predicted to encode secreted proteins.

Conflicting patterns of DNA barcoding and taxonomy in the cicada genus Tettigettalna from Southern Europe (Hemiptera: Cicadidae)

DNA barcodes have great potential to assist in species identification, especially when high taxonomical expertise is required. We investigated the utility of the 5' mitochondrial cytochrome c oxidase I (COI) region to discriminate between 13 European cicada species. These included all nine species currently recognized under the genus Tettigettalna, from which seven are endemic to the southern Iberian Peninsula. These cicadas have species-specific male calling songs but are morphologically very similar. Mean COI divergence between congeners ranged from 0.4% to 10.6%, but this gene was proven insufficient to determine species limits within genus Tettigettalna because a barcoding gap was absent for several of its species, that is, the highest intraspecific distance exceeded the lowest interspecific distance. The genetic data conflicted with current taxonomic classification for T. argentata and T. mariae. Neighbour-joining and Bayesian analyses revealed that T. argentata is geographically structured (clades North and South) and might constitute a species complex together with T. aneabi and T. mariae. The latter diverges very little from the southern clade of T. argentata and shares with it its most common haplotype. T. mariae is often in sympatry with T. argentata but it remains unclear whether introgression or incomplete lineage sorting may be responsible for the sharing of haplotypes. T. helianthemi and T. defauti also show high intraspecific variation that might signal hidden cryptic diversity. These taxonomic conflicts must be re-evaluated with further studies using additional genes and extensive morphological and acoustic analyses.

Assemblathon 2: evaluating de novo methods of genome assembly in three vertebrate species

BACKGROUND

The process of generating raw genome sequence data continues to become cheaper, faster, and more accurate. However, assembly of such data into high-quality, finished genome sequences remains challenging. Many genome assembly tools are available, but they differ greatly in terms of their performance (speed, scalability, hardware requirements, acceptance of newer read technologies) and in their final output (composition of assembled sequence). More importantly, it remains largely unclear how to best assess the quality of assembled genome sequences. The Assemblathon competitions are intended to assess current state-of-the-art methods in genome assembly.

RESULTS

In Assemblathon 2, we provided a variety of sequence data to be assembled for three vertebrate species (a bird, a fish, and snake). This resulted in a total of 43 submitted assemblies from 21 participating teams. We evaluated these assemblies using a combination of optical map data, Fosmid sequences, and several statistical methods. From over 100 different metrics, we chose ten key measures by which to assess the overall quality of the assemblies.

CONCLUSIONS

Many current genome assemblers produced useful assemblies, containing a significant representation of their genes and overall genome structure. However, the high degree of variability between the entries suggests that there is still much room for improvement in the field of genome assembly and that approaches which work well in assembling the genome of one species may not necessarily work well for another.

Molecular evidence of polyandry in the citrus mealybug, Planococcus citri (Hemiptera: Pseudococcidae)

The occurrence of polyandry in Planococcuscitri, presumed by earlier observations of mating behavior, was confirmed using microsatellite genotyping of pools of over 400 eggs resulting from controlled crosses of one female with two males. The genetic contribution of both mated males was confirmed in 13 out of 43 crosses. In three crosses it was possible to determine that only the first male fertilized the eggs, which may be due to sperm competition or unviable sperm supply. The microsatellite analysis also allowed the confirmation of aspects of the chromosomal inheritance detected previously in cytogenetic studies in Planococcuscitri, namely that only one of the alleles is transmitted by the male, indicating that the males are functionally haploid, supporting the observation of Paternal Genome Elimination (PGE) in these insects.