Endless forms most frustrating: disentangling species boundaries in the Ramalina decipiens group (Lecanoromycetes, Ascomycota), with the description of six new species and a key to the group.

Oceanic islands have been recognized as natural laboratories in which to study a great variety of evolutionary processes. One such process is evolutionary radiations, the diversification of a single ancestor into a number of species that inhabit different environments and differ in the traits that allow them to exploit those environments. The factors that drive evolutionary radiations have been studied for decades in charismatic organisms such as birds or lizards, but are lacking in lichen-forming fungi, despite recent reports of some lineages showing diversification patterns congruent with radiation. Here we propose the Ramalina decipiens group as a model system in which to carry out such studies. This group is currently thought to be comprised of five saxicolous species, all of them endemic to the Macaronesian region (the Azores, Madeira, Selvagens, Canary and Cape Verde islands). Three species are single-island endemics (a rare geographic distribution pattern in lichens), whereas two are widespread and show extreme morphological variation. The latter are suspected to harbor unrecognized species-level lineages. In order to use the Ramalina decipiens group as a model system it is necessary to resolve the group's phylogeny and to clarify its species boundaries. In this study we attempt to do so following an integrative taxonomy approach. We constructed a phylogenetic tree based on six molecular markers, four of which are newly developed and generated competing species hypotheses based on molecular (species discovery strategies based on both single locus and multilocus datasets) and phenotypic data (unsupervised clustering algorithms based on morphology, secondary chemistry and geographic origin). We found that taxonomic diversity in the Ramalina decipiens group has been highly underestimated in previous studies. In consequence, we describe six new species, most of them single-island endemics and provide a key to the group. Phylogenetic relationships among species have been reconstructed with almost full support which, coupled with the endemic character of the group, makes it an excellent system for the study of island radiations in lichen-forming fungi. Citation: Blázquez M, Pérez-Vargas I, Garrido-Benavent I, et al. 2024. Endless forms most frustrating: disentangling species boundaries in the Ramalina decipiens group (Lecanoromycetes, Ascomycota), with the description of six new species and a key to the group. Persoonia 52: 44-93. https://doi.org/10.3767/persoonia.2024.52.03 .

Comparative genomic analysis of chemosensory-related gene families in gastropods.

Chemoreception is critical for the survival and reproduction of animals. Except for a reduced group of insects and chelicerates, the molecular identity of chemosensory proteins is poorly understood in invertebrates. Gastropoda is the extant mollusk class with the greatest species richness, including marine, freshwater, and terrestrial lineages, and likely, highly diverse chemoreception systems. Here, we performed a comprehensive comparative genome analysis taking advantage of the chromosome-level information of two Gastropoda species, one of which belongs to a lineage that underwent a whole genome duplication event. We identified thousands of previously uncharacterized chemosensory-related genes, the majority of them encoding G protein-coupled receptors (GPCR), mostly organized into clusters distributed across all chromosomes. We also detected gene families encoding degenerin epithelial sodium channels (DEG-ENaC), ionotropic receptors (IR), sensory neuron membrane proteins (SNMP), Niemann-Pick type C2 (NPC2) proteins, and lipocalins, although with a lower number of members. Our phylogenetic analysis of the GPCR gene family across protostomes revealed: (i) remarkable gene family expansions in Gastropoda; (ii) clades including members from all protostomes; and (iii) species-specific clades with a substantial number of receptors. For the first time, we provide new and valuable knowledge into the evolution of the chemosensory gene families in invertebrates other than arthropods.

Shallow Whole-Genome Sequencing of Aedes japonicus and Aedes koreicus from Italy and an Updated Picture of Their Evolution Based on Mitogenomics and Barcoding.

Aedes japonicus and Aedes koreicus are two invasive mosquitoes native to East Asia that are quickly establishing in temperate regions of Europe. Both species are vectors of arboviruses, but we currently lack a clear understanding of their evolution. Here, we present new short-read, shallow genome sequencing of A. japonicus and A. koreicus individuals from northern Italy, which we used for downstream phylogenetic and barcode analyses. We explored associated microbial DNA and found high occurrences of Delftia bacteria in both samples, but neither Asaia nor Wolbachia. We then assembled complete mitogenomes and used these data to infer divergence times estimating the split of A. japonicus from A. koreicus in the Oligocene, which was more recent than that previously reported using mitochondrial markers. We recover a younger age for most other nodes within Aedini and other Culicidae. COI barcoding and phylogenetic analyses indicate that A. japonicus yaeyamensis, A. japonicus amamiensis, and the two A. koreicus sampled from Europe should be considered as separate species within a monophyletic species complex. Our studies further clarify the evolution of A. japonicus and A. koreicus, and indicate the need to obtain whole-genome data from putative species in order to disentangle their complex patterns of evolution.

Bowel Preparation for Colonoscopy Changes Serum Composition as Detected by Thermal Liquid Biopsy and Fluorescence Spectroscopy.

(1) Background: About 50% of prescribed colonoscopies report no pathological findings. A secondary screening test after fecal immunochemical test positivity (FIT+) would be required. Considering thermal liquid biopsy (TLB) as a potential secondary test, the aim of this work was to study possible interferences of colonoscopy bowel preparation on TLB outcome on a retrospective study; (2) Methods: Three groups were studied: 1/514 FIT(+) patients enrolled in a colorectal screening program (CN and CP with normal and pathological colonoscopy, respectively), with blood samples obtained just before colonoscopy and after bowel preparation; 2/55 patients from the CN group with blood sample redrawn after only standard 8-10 h fasting and no bowel preparation (CNR); and 3/55 blood donors from the biobank considered as a healthy control group; (3) Results: The results showed that from the 514 patients undergoing colonoscopy, 247 had CN and 267 had CP. TLB parameters in these two groups were similar but different from those of the blood donors. The resampled patients (with normal colonoscopy and no bowel preparation) had similar TLB parameters to those of the blood donors. TLB parameters together with fluorescence spectra and other serum indicators (albumin and C-reactive protein) confirmed the statistically significant differences between normal colonoscopy patients with and without bowel preparation; (4) Conclusions: Bowel preparation seemed to alter serum protein levels and altered TLB parameters (different from a healthy subject). The diagnostic capability of other liquid-biopsy-based methods might also be compromised. Blood extraction after bowel preparation for colonoscopy should be avoided.

The chromosome-scale assembly of the Canary Islands endemic spider Dysdera silvatica (Arachnida, Araneae) sheds light on the origin and genome structure of chemoreceptor gene families in chelicerates.

Here, we present the chromosome-level genome assembly of Dysdera silvatica Schmidt, 1981, a nocturnal ground-dwelling spider endemic from the Canary Islands. The genus Dysdera has undergone a remarkable diversification in this archipelago mostly associated with shifts in the level of trophic specialization, becoming an excellent model to study the genomic drivers of adaptive radiations. The new assembly (1.37 Gb; scaffold N50 of 174.2 Mb), was performed using the chromosome conformation capture scaffolding technique, represents a continuity improvement of more than 4500 times with respect to the previous version. The seven largest scaffolds or pseudochromosomes, which cover 87% of the total assembly size, probably correspond with the seven chromosomes of the karyotype of this species, including a characteristic large X chromosome. To illustrate the value of this new resource we performed a comprehensive analysis of the two major arthropod chemoreceptor gene families (i.e., gustatory and ionotropic receptors). We identified 545 chemoreceptor sequences distributed across all pseudochromosomes, with a notable underrepresentation in the X chromosome. At least 54% of them localize in 83 genomic clusters with a significantly lower evolutionary distances between them than the average of the family, suggesting a recent origin of many of them. This chromosome-level assembly is the first high-quality genome representative of the Synspermiata clade, and just the third among spiders, representing a new valuable resource to gain insights into the structure and organization of chelicerate genomes, including the role that structural variants, repetitive elements and large gene families played in the extraordinary biology of spiders.

Genetic data from the extinct giant rat from Tenerife (Canary Islands) points to a recent divergence from mainland relatives.

Evolution of vertebrate endemics in oceanic islands follows a predictable pattern, known as the island rule, according to which gigantism arises in originally small-sized species and dwarfism in large ones. Species of extinct insular giant rodents are known from all over the world. In the Canary Islands, two examples of giant rats, †Canariomys bravoi and †Canariomys tamarani, endemic to Tenerife and Gran Canaria, respectively, disappeared soon after human settlement. The highly derived morphological features of these insular endemic rodents hamper the reconstruction of their evolutionary histories. We have retrieved partial nuclear and mitochondrial data from †C. bravoi and used this information to explore its evolutionary affinities. The resulting dated phylogeny confidently places †C. bravoi within the African grass rat clade (Arvicanthis niloticus). The estimated divergence time, 650 000 years ago (95% higher posterior densities: 373 000-944 000), points toward an island colonization during the Günz-Mindel interglacial stage. †Canariomys bravoi ancestors would have reached the island via passive rafting and then underwent a yearly increase of mean body mass calculated between 0.0015 g and 0.0023 g; this corresponds to fast evolutionary rates (in darwins (d), ranging from 7.09 d to 2.78 d) that are well above those observed for non-insular mammals.

Fluorescence Liquid Biopsy for Cancer Detection Is Improved by Using Cationic Dendronized Hyperbranched Polymer.

(1) Background: Biophysical techniques applied to serum samples characterization could promote the development of new diagnostic tools. Fluorescence spectroscopy has been previously applied to biological samples from cancer patients and differences from healthy individuals were observed. Dendronized hyperbranched polymers (DHP) based on bis(hydroxymethyl)propionic acid (bis-MPA) were developed in our group and their potential biomedical applications explored. (2) Methods: A total of 94 serum samples from diagnosed cancer patients and healthy individuals were studied (20 pancreatic ductal adenocarcinoma, 25 blood donor, 24 ovarian cancer, and 25 benign ovarian cyst samples). (3) Results: Fluorescence spectra of serum samples (fluorescence liquid biopsy, FLB) in the presence and the absence of DHP-bMPA were recorded and two parameters from the signal curves obtained. A secondary parameter, the fluorescence spectrum score (FSscore), was calculated, and the diagnostic model assessed. For pancreatic ductal adenocarcinoma (PDAC) and ovarian cancer, the classification performance was improved when including DHP-bMPA, achieving high values of statistical sensitivity and specificity (over 85% for both pathologies). (4) Conclusions: We have applied FLB as a quick, simple, and minimally invasive promising technique in cancer diagnosis. The classification performance of the diagnostic method was further improved by using DHP-bMPA, which interacted differentially with serum samples from healthy and diseased subjects. These preliminary results set the basis for a larger study and move FLB closer to its clinical application, providing useful information for the oncologist during patient diagnosis.

Thermal Liquid Biopsy (TLB) of Blood Plasma as a Potential Tool to Help in the Early Diagnosis of Multiple Sclerosis.

BACKGROUND: Multiple sclerosis (MS) is frequently characterized by a variety of clinical signs, often exhibiting little specificity. The diagnosis requires a combination of medical observations and instrumental tests, and any support for its objective assessment is helpful. OBJECTIVE: Herein, we describe the application of thermal liquid biopsy (TLB) of blood plasma samples, a methodology for predicting the occurrence of MS with a noninvasive, quick blood test. METHODS: TLB allows one to define an index (TLB score), which provides information about overall real-time alterations in plasma proteome that may be indicative of MS. RESULTS: This pilot study, based on 85 subjects (45 MS patients and 40 controls), showed good performance indexes (sensitivity and specificity both around 70%). The diagnostic methods better discriminate between early stage and low-burden MS patients, and it is not influenced by gender, age, or assumption of therapeutic drugs. TLB is more accurate for patients having low disability level (≤ 3.0, measured by the expanded disability status scale, EDSS) and a relapsing-remitting diagnosis. CONCLUSION: Our results suggest that TLB can be applied to MS, especially in an initial phase of the disease when diagnosis is difficult and yet more important (in such cases, accuracy of prediction is close to 80%), as well as in personalized patient periodic monitoring. The next step will be determining its utility in differentiating between MS and other disorders, in particular in inflammatory diseases.

Evolutionary History of Major Chemosensory Gene Families across Panarthropoda.

Chemosensory perception is a fundamental biological process of particular relevance in basic and applied arthropod research. However, apart from insects, there is little knowledge of specific molecules involved in this system, which is restricted to a few taxa with uneven phylogenetic sampling across lineages. From an evolutionary perspective, onychophorans (velvet worms) and tardigrades (water bears) are of special interest since they represent the closest living relatives of arthropods, altogether comprising the Panarthropoda. To get insights into the evolutionary origin and diversification of the chemosensory gene repertoire in panarthropods, we sequenced the antenna- and head-specific transcriptomes of the velvet worm Euperipatoides rowelli and analyzed members of all major chemosensory families in representative genomes of onychophorans, tardigrades, and arthropods. Our results suggest that the NPC2 gene family was the only family encoding soluble proteins in the panarthropod ancestor and that onychophorans might have lost many arthropod-like chemoreceptors, including the highly conserved IR25a receptor of protostomes. On the other hand, the eutardigrade genomes lack genes encoding the DEG-ENaC and CD36-sensory neuron membrane proteins, the chemosensory members of which have been retained in arthropods; these losses might be related to lineage-specific adaptive strategies of tardigrades to survive extreme environmental conditions. Although the results of this study need to be further substantiated by an increased taxon sampling, our findings shed light on the diversification of chemosensory gene families in Panarthropoda and contribute to a better understanding of the evolution of animal chemical senses.

Genome mining and sequence analysis of chemosensory soluble proteins in arthropods.

Identifying protein-coding genes from genome and transcriptome data is the first and one of the most important steps towards their comprehensive study. This chapter introduces both general procedures for sequence mining, and specific approaches for recognizing characteristic motives and chemical properties in soluble proteins potentially involved in arthropod chemical communication. We describe (i) the workflow to identify members of the OBP (Odorant-Binding Proteins) and CSP (Chemosensory Proteins) families in genomic and transcriptomic sequences using our recently developed bioinformatic solution, BITACORA, and (ii) the main further steps to visualize and to accurately annotate these genes in the Apollo genome browser. The success of further biochemical, functional and evolutionary analyses largely depends on the quality of these initial steps.

Feature selection based multivariate time series forecasting: An application to antibiotic resistance outbreaks prediction.

Antimicrobial resistance has become one of the most important health problems and global action plans have been proposed globally. Prevention plays a key role in these actions plan and, in this context, we propose the use of Artificial Intelligence, specifically Time Series Forecasting techniques, for predicting future outbreaks of Methicillin-resistant Staphylococcus aureus (MRSA). Infection incidence forecasting is approached as a Feature Selection based Time Series Forecasting problem using multivariate time series composed of incidence of Staphylococcus aureus Methicillin-sensible and MRSA infections, influenza incidence and total days of therapy of both of Levofloxacin and Oseltamivir antimicrobials. Data were collected from the University Hospital of Getafe (Spain) from January 2009 to January 2018, using months as time granularity. The main contributions of the work are the following: the applications of wrapper feature selection methods where the search strategy is based on multi-objective evolutionary algorithms (MOEA) along with evaluators based on the most powerful state-of-the-art regression algorithms. The performance of the feature selection methods has been measured using the root mean square error (RMSE) and mean absolute error (MAE) performance metrics. A novel multi-criteria decision-making process is proposed in order to select the most satisfactory forecasting model, using the metrics previously mentioned, as well as the slopes of model prediction lines in the 1, 2 and 3 steps-ahead predictions. The multi-criteria decision-making process is applied to the best models resulting from a ranking of databases and regression algorithms obtained through multiple statistical tests. Finally, to the best of our knowledge, this is the first time that a feature selection based multivariate time series methodology is proposed for antibiotic resistance forecasting. Final results show that the best model according to the proposed multi-criteria decision making process provides a RMSE = (0.1349, 0.1304, 0.1325) and a MAE = (0.1003, 0.096, 0.0987) for 1, 2, and 3 steps-ahead predictions.

bitacora: A comprehensive tool for the identification and annotation of gene families in genome assemblies.

Gene annotation is a critical bottleneck in genomic research, especially for the comprehensive study of very large gene families in the genomes of nonmodel organisms. Despite the recent progress in automatic methods, state-of-the-art tools used for this task often produce inaccurate annotations, such as fused, chimeric, partial or even completely absent gene models for many family copies, errors that require considerable extra efforts to be corrected. Here we present bitacora, a bioinformatics solution that integrates popular sequence similarity-based search tools and Perl scripts to facilitate both the curation of these inaccurate annotations and the identification of previously undetected gene family copies directly in genomic DNA sequences. We tested the performance of bitacora in annotating the members of two chemosensory gene families with different repertoire size in seven available genome sequences, and compared its performance with that of augustus-ppx, a tool also designed to improve automatic annotations using a sequence similarity-based approach. Despite the relatively high fragmentation of some of these drafts, bitacora was able to improve the annotation of many members of these families and detected thousands of new chemoreceptors encoded in genome sequences. The program creates general feature format (GFF) files, with both curated and newly identified gene models, and FASTA files with the predicted proteins. These outputs can be easily integrated in genomic annotation editors, greatly facilitating subsequent manual annotation and downstream evolutionary analyses.

Genomic adaptations to aquatic and aerial life in mayflies and the origin of insect wings.

The evolution of winged insects revolutionized terrestrial ecosystems and led to the largest animal radiation on Earth. However, we still have an incomplete picture of the genomic changes that underlay this diversification. Mayflies, as one of the sister groups of all other winged insects, are key to understanding this radiation. Here, we describe the genome of the mayfly Cloeon dipterum and its gene expression throughout its aquatic and aerial life cycle and specific organs. We discover an expansion of odorant-binding-protein genes, some expressed specifically in breathing gills of aquatic nymphs, suggesting a novel sensory role for this organ. In contrast, flying adults use an enlarged opsin set in a sexually dimorphic manner, with some expressed only in males. Finally, we identify a set of wing-associated genes deeply conserved in the pterygote insects and find transcriptomic similarities between gills and wings, suggesting a common genetic program. Globally, this comprehensive genomic and transcriptomic study uncovers the genetic basis of key evolutionary adaptations in mayflies and winged insects.

Understanding the Early Evolutionary Stages of a Tandem Drosophilamelanogaster-Specific Gene Family: A Structural and Functional Population Study.

Gene families underlie genetic innovation and phenotypic diversification. However, our understanding of the early genomic and functional evolution of tandemly arranged gene families remains incomplete as paralog sequence similarity hinders their accurate characterization. The Drosophila melanogaster-specific gene family Sdic is tandemly repeated and impacts sperm competition. We scrutinized Sdic in 20 geographically diverse populations using reference-quality genome assemblies, read-depth methodologies, and qPCR, finding that ∼90% of the individuals harbor 3-7 copies as well as evidence of population differentiation. In strains with reliable gene annotations, copy number variation (CNV) and differential transposable element insertions distinguish one structurally distinct version of the Sdic region per strain. All 31 annotated copies featured protein-coding potential and, based on the protein variant encoded, were categorized into 13 paratypes differing in their 3' ends, with 3-5 paratypes coexisting in any strain examined. Despite widespread gene conversion, the only copy present in all strains has functionally diverged at both coding and regulatory levels under positive selection. Contrary to artificial tandem duplications of the Sdic region that resulted in increased male expression, CNV in cosmopolitan strains did not correlate with expression levels, likely as a result of differential genome modifier composition. Duplicating the region did not enhance sperm competitiveness, suggesting a fitness cost at high expression levels or a plateau effect. Beyond facilitating a minimally optimal expression level, Sdic CNV acts as a catalyst of protein and regulatory diversity, showcasing a possible evolutionary path recently formed tandem multigene families can follow toward long-term consolidation in eukaryotic genomes.

Genomic Analysis of European Drosophila melanogaster Populations Reveals Longitudinal Structure, Continent-Wide Selection, and Previously Unknown DNA Viruses.

Genetic variation is the fuel of evolution, with standing genetic variation especially important for short-term evolution and local adaptation. To date, studies of spatiotemporal patterns of genetic variation in natural populations have been challenging, as comprehensive sampling is logistically difficult, and sequencing of entire populations costly. Here, we address these issues using a collaborative approach, sequencing 48 pooled population samples from 32 locations, and perform the first continent-wide genomic analysis of genetic variation in European Drosophila melanogaster. Our analyses uncover longitudinal population structure, provide evidence for continent-wide selective sweeps, identify candidate genes for local climate adaptation, and document clines in chromosomal inversion and transposable element frequencies. We also characterize variation among populations in the composition of the fly microbiome, and identify five new DNA viruses in our samples.

A methodology based on multiple criteria decision analysis for combining antibiotics in empirical therapy.

BACKGROUND: The current situation of critical progression in resistance to more effective antibiotics has forced the reuse of old highly toxic antibiotics and, for several reasons, the extension of the indications of combined antibiotic therapy as alternative options to broad spectrum empirical mono-therapy. A key aspect for selecting an appropriate and adequate antimicrobial therapy is that prescription must be based on local epidemiology and knowledge since many aspects, such as prevalence of microorganisms and effectiveness of antimicrobials, change from hospitals, or even areas and services within a single hospital. Therefore, the selection of combinations of antibiotics requires the application of a methodology that provides objectivity, completeness and reproducibility to the analysis of the detailed microbiological, epidemiological, pharmacological information on which to base a rational and reasoned choice. METHODS: We proposed a methodology for decision making that uses a multiple criteria decision analysis (MCDA) to support the clinician in the selection of an efficient combined empiric therapy. The MCDA includes a multi-objective constrained optimization model whose criteria are the maximum efficacy of therapy, maximum activity, the minimum activity overlapping, the minimum use of restricted antibiotics, the minimum toxicity of antibiotics and the activity against the most prevalent and virulent bacteria. The decision process can be defined in 4 steps: (1) selection of clinical situation of interest, (2) definition of local optimization criteria, (3) definition of constraints for reducing combinations, (4) manual sorting of solutions according to patient's clinical conditions, and (5) selection of a combination. EXPERIMENTS AND RESULTS: In order to show the application of the methodology to a clinical case, we carried out experiments with antibiotic susceptibility tests in blood samples taken during a five years period at a university hospital. The validation of the results consists of a manual review of the combinations and experiments carried out by an expert physician that has explained the most relevant solutions proposed according to current clinical knowledge and their use. CONCLUSION: We show that with the decision process proposed, the physician is able to select the best combined therapy according to different criteria such as maximum efficacy, activity and minimum toxicity. A method for the recommendation of combined antibiotic therapy developed on the basis of a multi-objective optimization model may assist the physicians in the search for alternatives to the use of broad-spectrum antibiotics or restricted antibiotics for empirical therapy. The decision proposed can be easily reproduced for any local epidemiology and any different clinical settings.

Thermal Liquid Biopsy (TLB) Focused on Benign and Premalignant Pancreatic Cyst Diagnosis.

BACKGROUND: Current efforts in the identification of new biomarkers are directed towards an accurate differentiation between benign and premalignant cysts. Thermal Liquid Biopsy (TLB) has been previously applied to inflammatory and tumor diseases and could offer an interesting point of view in this type of pathology. METHODS: In this work, twenty patients (12 males and 8 females, average ages 62) diagnosed with a pancreatic cyst benign (10) and premalignant (10) cyst lesions were recruited, and biological samples were obtained during the endoscopic ultrasonography procedure. RESULTS: Proteomic content of cyst liquid samples was studied and several common proteins in the different groups were identified. TLB cyst liquid profiles reflected protein content. Also, TLB serum score was able to discriminate between healthy and cysts patients (71% sensitivity and 98% specificity) and between benign and premalignant cysts (75% sensitivity and 67% specificity). CONCLUSIONS: TLB analysis of plasmatic serum sample, a quick, simple and non-invasive technique that can be easily implemented, reports valuable information on the observed pancreatic lesion. These preliminary results set the basis for a larger study to refine TLB serum score and move closer to the clinical application of TLB providing useful information to the gastroenterologist during patient diagnosis.

Cryptic species delineation in freshwater planarians of the genus Dugesia (Platyhelminthes, Tricladida): Extreme intraindividual genetic diversity, morphological stasis, and karyological variability.

The keystone of planarian taxonomy traditionally has been the anatomy of the copulatory apparatus. However, many planarian species comprise asexual fissiparous populations, with the fissiparous animals not developing a copulatory apparatus, thus precluding their morphological identification. Incorporation of molecular data into planarian systematics has been of great value, not only in the identification of fissiparous individuals but also as an additional source of information for determining species boundaries. Nevertheless, the discrepancy between morphological and molecular data has highlighted the need for extra sources of taxonomic information. Moreover, a recent study has pointed out that fissiparous reproduction may lead to high levels of intraindividual genetic diversity in planarians, which may mislead molecular analyses. In the present study we aim to test a new up-to-date integrative taxonomic procedure for planarians, including intraindividual genetic data and additional sources of taxonomic information, besides morphology and DNA, using Dugesia subtentaculata sensu lato as a model organism, a species with an intricate taxonomic history. First, we used three different methods for molecular species delimitation on single locus datasets, both with and without intraindividual information, for formulating Primary Species Hypotheses (PSHs). Subsequently, Secondary Species Hypotheses (SSHs) were formulated on the basis of three types of information: (1) a coalescent-based species delimitation method applied to multilocus data, (2) morphology of the copulatory apparatus, and (3) karyological metrics. This resulted in the delimitation of four morphologically cryptic species within the nominal species D. subtentaculata. Our results provide evidence that the analysis of intraindividual genetic data is essential for properly developing PSHs in planarians. Our study reveals also that karyological differentiation, rather than morphological differentiation, may play an important role in speciation processes in planarians, thus suggesting that the currently known diversity of the group could be highly underestimated.

The avocado genome informs deep angiosperm phylogeny, highlights introgressive hybridization, and reveals pathogen-influenced gene space adaptation.

The avocado, Persea americana, is a fruit crop of immense importance to Mexican agriculture with an increasing demand worldwide. Avocado lies in the anciently diverged magnoliid clade of angiosperms, which has a controversial phylogenetic position relative to eudicots and monocots. We sequenced the nuclear genomes of the Mexican avocado race, P. americana var. drymifolia, and the most commercially popular hybrid cultivar, Hass, and anchored the latter to chromosomes using a genetic map. Resequencing of Guatemalan and West Indian varieties revealed that ∼39% of the Hass genome represents Guatemalan source regions introgressed into a Mexican race background. Some introgressed blocks are extremely large, consistent with the recent origin of the cultivar. The avocado lineage experienced 2 lineage-specific polyploidy events during its evolutionary history. Although gene-tree/species-tree phylogenomic results are inconclusive, syntenic ortholog distances to other species place avocado as sister to the enormous monocot and eudicot lineages combined. Duplicate genes descending from polyploidy augmented the transcription factor diversity of avocado, while tandem duplicates enhanced the secondary metabolism of the species. Phenylpropanoid biosynthesis, known to be elicited by Colletotrichum (anthracnose) pathogen infection in avocado, is one enriched function among tandems. Furthermore, transcriptome data show that tandem duplicates are significantly up- and down-regulated in response to anthracnose infection, whereas polyploid duplicates are not, supporting the general view that collections of tandem duplicates contribute evolutionarily recent "tuning knobs" in the genome adaptive landscapes of given species.

The draft genome sequence of the spider Dysdera silvatica (Araneae, Dysderidae): A valuable resource for functional and evolutionary genomic studies in chelicerates.

BACKGROUND: We present the draft genome sequence of Dysdera silvatica, a nocturnal ground-dwelling spider from a genus that has undergone a remarkable adaptive radiation in the Canary Islands. RESULTS: The draft assembly was obtained using short (Illumina) and long (PaciBio and Nanopore) sequencing reads. Our de novo assembly (1.36 Gb), which represents 80% of the genome size estimated by flow cytometry (1.7 Gb), is constituted by a high fraction of interspersed repetitive elements (53.8%). The assembly completeness, using BUSCO and core eukaryotic genes, ranges from 90% to 96%. Functional annotations based on both ab initio and evidence-based information (including D. silvatica RNA sequencing) yielded a total of 48,619 protein-coding sequences, of which 36,398 (74.9%) have the molecular hallmark of known protein domains, or sequence similarity with Swiss-Prot sequences. The D. silvatica assembly is the first representative of the superfamily Dysderoidea, and just the second available genome of Synspermiata, one of the major evolutionary lineages of the "true spiders" (Araneomorphae). CONCLUSIONS: Dysderoids, which are known for their numerous instances of adaptation to underground environments, include some of the few examples of trophic specialization within spiders and are excellent models for the study of cryptic female choice. This resource will be therefore useful as a starting point to study fundamental evolutionary and functional questions, including the molecular bases of the adaptation to extreme environments and ecological shifts, as well of the origin and evolution of relevant spider traits, such as the venom and silk.