Unravelling phylogenetic relationships of the tribe Cereeae using target enrichment sequencing

BACKGROUND AND AIMS

Cactaceae are succulent plants, quasi-endemic to the American continent, and one of the most endangered plant groups in the world. Molecular phylogenies have been key to unravelling phylogenetic relationships among major cactus groups, previously hampered by high levels of morphological convergence. Phylogenetic studies using plastid markers have not provided adequate resolution for determining generic relationships within cactus groups. This is the case for the tribe Cereeae s.l., a highly diverse group from tropical America. Here we aimed to reconstruct a well-resolved phylogenetic tree of tribe Cereeae and update the circumscription of suprageneric and generic groups in this tribe.

METHODS

We integrated sequence data from public gene and genomic databases with new target sequences (generated using the customized Cactaceae591 probe set) across representatives of this tribe, with a denser taxon sampling of the subtribe Cereinae. We inferred concatenated and coalescent phylogenetic trees and compared the performance of both approaches.

KEY RESULTS

Six well-supported suprageneric clades were identified using different datasets. However, only genomic datasets, especially the Cactaceae591, were able to resolve the contentious relationships within the subtribe Cereinae.

CONCLUSIONS

We propose a new taxonomic classification within Cereeae based on well-resolved clades, including new subtribes (Aylosterinae subtr. nov., Uebelmanniinae subtr. nov. and Gymnocalyciinae subtr. nov.) and revised subtribes (Trichocereinae, Rebutiinae and Cereinae). We emphasize the importance of using genomic datasets allied with coalescent inference to investigate evolutionary patterns within the tribe Cereeae.

Revisiting phylogeny, systematics, and biogeography of a Pleistocene radiation

PREMISE

Pilosocereus (Cactaceae) is an important dry forest element in all subregions and transitional zones of the neotropics, with the highest diversity in eastern Brazil. The genus is subdivided into informal taxonomic groups; however, most of these are not supported by recent molecular phylogenetic inferences. This lack of confidence is probably due to the use of an insufficient number of loci and the complexity of cactus diversification. Here, we explored the species relationships in Pilosocereus in more detail, integrating multilocus phylogenetic approaches with the assessment of the ancestral range and the effect of geography on diversification shifts.

METHODS

We used 28 nuclear, plastid, and mitochondrial loci from 54 plant samples of 31 Pilosocereus species for phylogenetic analyses. We used concatenated and coalescent phylogenetic trees and Bayesian models to estimate the most likely ancestral range and diversification shifts.

RESULTS

All Pilosocereus species were clustered in the same branch, except P. bohlei. The phylogenetic relationships were more associated with the geographic distribution than taxonomic affinities among taxa. The genus began diversifying during the Plio-Pleistocene transition in the Caatinga domain and experienced an increased diversification rate during the Calabrian age.

CONCLUSIONS

We recovered a well-supported multispecies coalescent phylogeny. Our results refine the pattern of rapid diversification of Pilosocereus species across neotropical drylands during the Pleistocene and highlight the need for taxonomic rearrangements in the genus. We recovered a pulse of diversification during the Pleistocene that was likely driven by multiple dispersal and vicariance events within and among the Caatinga, Cerrado, and Atlantic Forest domains.

Coalescent-based species delimitation meets deep learning: Insights from a highly fragmented cactus system

Delimiting species boundaries is a major goal in evolutionary biology. An increasing volume of literature has focused on the challenges of investigating cryptic diversity within complex evolutionary scenarios of speciation, including gene flow and demographic fluctuations. New methods based on model selection, such as approximate Bayesian computation, approximate likelihoods, and machine learning are promising tools arising in this field. Here, we introduce a framework for species delimitation using the multispecies coalescent model coupled with a deep learning algorithm based on convolutional neural networks (CNNs). We compared this strategy with a similar ABC approach. We applied both methods to test species boundary hypotheses based on current and previous taxonomic delimitations as well as genetic data (sequences from 41 loci) in Pilosocereus aurisetus, a cactus species complex with a sky-island distribution and taxonomic uncertainty. To validate our method, we also applied the same strategy on data from widely accepted species from the genus Drosophila. The results show that our CNN approach has a high capacity to distinguish among the simulated species delimitation scenarios, with higher accuracy than ABC. For the cactus data set, a splitter hypothesis without gene flow showed the highest probability in both CNN and ABC approaches, a result agreeing with previous taxonomic classifications and in line with the sky-island distribution and low dispersal of P. aurisetus. Our results highlight the cryptic diversity within the P. aurisetus complex and show that CNNs are a promising approach for distinguishing complex evolutionary histories, even outperforming the accuracy of other model-based approaches such as ABC.

Comparative and predictive phylogeography in the South American diagonal of open formations: Unravelling the biological and environmental influences on multitaxon demography

Phylogeography investigates historical drivers of the geographical distribution of intraspecific lineages. Particular attention has been given to ecological, climatic and geological processes in the diversification of the Neotropical biota. Several species sampled across the South American diagonal of open formations (DOF), comprising the Caatinga, Cerrado and Chaco biomes, experienced range shifts coincident with Quaternary climatic changes. However, comparative studies across different spatial, temporal and biological scales on DOF species are still meagre. Here, we combine phylogeographical model selection and machine learning predictive frameworks to investigate the influence of Pleistocene climatic changes on several plant and animal species from the DOF. We assembled mitochondrial/chloroplastic DNA sequences in public repositories and inferred the demographic responses of 44 species, comprising 70 intraspecific lineages of plants, lizards, frogs, spiders and insects. We then built a random forest model using biotic and abiotic information to identify the best predictors of demographic responses in the Pleistocene. Finally, we assessed the temporal synchrony of species demographic responses with hierarchical approximate Bayesian computation. Biotic variables related to population connectivity, gene flow and habitat preferences largely predicted how species responded to Pleistocene climatic changes, and demographic changes were synchronous primarily during the Middle Pleistocene. Although 22 (~31%) lineages underwent demographic expansion, presumably associated with the spread of aridity during the glacial Pleistocene periods, our findings suggest that nine lineages (~13%) exhibited the opposite response due to taxon-specific attributes.

Maintaining genetic integrity with high promiscuity: Frequent hybridization with low introgression in multiple hybrid zones of Melocactus (Cactaceae)

Hybridization and introgression between species in contact/hybrid zones provide important insight into the genetic and ecological mechanisms of speciation. Cactaceae represents the most important radiation of true succulent angiosperms in the New World. This diversification continues to date, with species experiencing few intrinsic barriers to gene flow and the frequent occurrence of natural hybridization. Here, we used RAD-Seq single-nucleotide polymorphism (SNP) data to investigate the genetic architecture of hybridization in four hybrid zones hosting Melocactus concinnus and four congeneric species (M. ernestii, M. glaucescens, M. paucispinus, and M. zehntneri). Our results revealed that M. concinnus is highly promiscuous in sympatric areas and hybridizes with various species distributed in Morro do Chapéu (Diamantina Plateau, Bahia), eastern Brazil. However, the contemporary genomic introgression among the investigated species is very low (c. 2-5%), confirming that even in the face of hybridization, Melocactus species maintain their genetic integrity. The genomic cline approach showed a large fraction of loci deviating from a model of neutral introgression, where most of the loci are consistent with selection favoring parental genotypes. Our results suggest the occurrence of weak premating but strong postmating reproductive isolation in the analyzed cactus species. Furthermore, as most of the Melocactus species are restricted in distribution, hybridization might negatively affect their integrity if hybrids replace the parental species.

Model-based analysis supports interglacial refugia over long-dispersal events in the diversification of two South American cactus species

Pilosocereus machrisii and P. aurisetus are cactus species within the P. aurisetus complex, a group of eight cacti that are restricted to rocky habitats within the Neotropical savannas of eastern South America. Previous studies have suggested that diversification within this complex was driven by distributional fragmentation, isolation leading to allopatric differentiation, and secondary contact among divergent lineages. These events have been associated with Quaternary climatic cycles, leading to the hypothesis that the xerophytic vegetation patches which presently harbor these populations operate as refugia during the current interglacial. However, owing to limitations of the standard phylogeographic approaches used in these studies, this hypothesis was not explicitly tested. Here we use Approximate Bayesian Computation to refine the previous inferences and test the role of different events in the diversification of two species within P. aurisetus group. We used molecular data from chloroplast DNA and simple sequence repeats loci of P. machrisii and P. aurisetus, the two species with broadest distribution in the complex, in order to test if the diversification in each species was driven mostly by vicariance or by long-dispersal events. We found that both species were affected primarily by vicariance, with a refuge model as the most likely scenario for P. aurisetus and a soft vicariance scenario most probable for P. machrisii. These results emphasize the importance of distributional fragmentation in these species, and add support to the hypothesis of long-term isolation in interglacial refugia previously proposed for the P. aurisetus species complex diversification.

Anonymous nuclear markers data supporting species tree phylogeny and divergence time estimates in a cactus species complex in South America

Supportive data related to the article “Anonymous nuclear markers reveal taxonomic incongruence and long-term disjunction in a cactus species complex with continental-island distribution in South America” (Perez et al., 2016) [1]. Here, we present pyrosequencing results, primer sequences, a cpDNA phylogeny, and a species tree phylogeny.

Lineage-specific evolutionary rate in plants: Contributions of a screening for Cereus (Cactaceae)

PREMISE OF THE STUDY

Predictable chloroplast DNA (cpDNA) sequences have been listed for the shallowest taxonomic studies in plants. We investigated whether plastid regions that vary between closely allied species could be applied for intraspecific studies and compared the variation of these plastid segments with two nuclear regions.

METHODS

We screened 16 plastid and two nuclear intronic regions for species of the genus Cereus (Cactaceae) at three hierarchical levels (species from different clades, species of the same clade, and allopatric populations).

RESULTS

Ten plastid regions presented interspecific variation, and six of them showed variation at the intraspecific level. The two nuclear regions showed both inter- and intraspecific variation, and in general they showed higher levels of variability in almost all hierarchical levels than the plastid segments.

DISCUSSION

Our data suggest no correspondence between variation of plastid regions at the interspecific and intraspecific level, probably due to lineage-specific variation in cpDNA, which appears to have less effect in nuclear data. Despite the heterogeneity in evolutionary rates of cpDNA, we highlight three plastid segments that may be considered in initial screenings in plant phylogeographic studies.

Anonymous nuclear markers reveal taxonomic incongruence and long-term disjunction in a cactus species complex with continental-island distribution in South America

The Pilosocereus aurisetus complex consists of eight cactus species with a fragmented distribution associated to xeric enclaves within the Cerrado biome in eastern South America. The phylogeny of these species is incompletely resolved, and this instability complicates evolutionary analyses. Previous analyses based on both plastid and microsatellite markers suggested that this complex contained species with inherent phylogeographic structure, which was attributed to recent diversification and recurring range shifts. However, limitations of the molecular markers used in these analyses prevented some questions from being properly addressed. In order to better understand the relationship among these species and make a preliminary assessment of the genetic structure within them, we developed anonymous nuclear loci from pyrosequencing data of 40 individuals from four species in the P. aurisetus complex. The data obtained from these loci were used to identify genetic clusters within species, and to investigate the phylogenetic relationship among these inferred clusters using a species tree methodology. Coupled with a palaeodistributional modelling, our results reveal a deep phylogenetic and climatic disjunction between two geographic lineages. Our results highlight the importance of sampling more regions from the genome to gain better insights on the evolution of species with an intricate evolutionary history. The methodology used here provides a feasible approach to develop numerous genealogical molecular markers throughout the genome for non-model species. These data provide a more robust hypothesis for the relationship among the lineages of the P. aurisetus complex.

Using Next Generation RAD Sequencing to Isolate Multispecies Microsatellites for Pilosocereus (Cactaceae)

Microsatellite markers (also known as SSRs, Simple Sequence Repeats) are widely used in plant science and are among the most informative molecular markers for population genetic investigations, but the development of such markers presents substantial challenges. In this report, we discuss how next generation sequencing can replace the cloning, Sanger sequencing, identification of polymorphic loci, and testing cross-amplification that were previously required to develop microsatellites. We report the development of a large set of microsatellite markers for five species of the Neotropical cactus genus Pilosocereus using a restriction-site-associated DNA sequencing (RAD-seq) on a Roche 454 platform. We identified an average of 165 microsatellites per individual, with the absolute numbers across individuals proportional to the sequence reads obtained per individual. Frequency distribution of the repeat units was similar in the five species, with shorter motifs such as di- and trinucleotide being the most abundant repeats. In addition, we provide 72 microsatellites that could be potentially amplified in the sampled species and 22 polymorphic microsatellites validated in two populations of the species Pilosocereus machrisii. Although low coverage sequencing among individuals was observed for most of the loci, which we suggest to be more related to the nature of the microsatellite markers and the possible bias inserted by the restriction enzymes than to the genome size, our work demonstrates that an NGS approach is an efficient method to isolate multispecies microsatellites even in non-model organisms.

Polymorphic microsatellite markers for the rare and endangered cactus Uebelmannia pectinifera (Cactaceae) and its congeneric species

The cactus genus Uebelmannia includes 3 narrow endemic species associated with rocky savanna habitats in eastern South America. Because of their rarity and illegal over-collection, all of these species are endangered. Taxonomic uncertainties resulting from dramatic local variation in morphology within Uebelmannia species preclude effective conservation efforts, such as the reintroduction or translocation of plants, to restore declining populations. In this study, we developed and characterized 18 perfect, dinucleotide simple-sequence repeat markers for U. pectinifera, the most widely distributed species in the genus, and tested the cross-amplification of these markers in the remaining congeneric species and subspecies. All markers were polymorphic in a sample from 2 U. pectinifera populations. The effective number of alleles ranged from 1.6 to 8.7, with an average per population of 3.3 (SE ± 0.30) and 4.5 (SE ± 0.50). Expected heterozygosity ranged from 0.375 to 0.847 and 8-10 loci showed departures from Hardy- Weinberg equilibrium in the analyzed populations. Based on the observed polymorphism level of each marker, as well as the analysis of null allele presence and evidence of amplification of duplicate loci, a subset of 12 loci can be used as reliable markers to investigate the genetic structure, diversity, and species limits of the Uebelmannia genus.

Cross-species amplification of microsatellites reveals incongruence in the molecular variation and taxonomic limits of the Pilosocereus aurisetus group (Cactaceae)

The Pilosocereus aurisetus group contains eight cactus species restricted to xeric habitats in eastern and central Brazil that have an archipelago-like distribution. In this study, 5-11 microsatellite markers previously designed for Pilosocereus machrisii were evaluated for cross-amplification and polymorphisms in ten populations from six species of the P. aurisetus group. The genotypic information was subsequently used to investigate the genetic relationships between the individuals, populations, and species analyzed. Only the Pmac101 locus failed to amplify in all of the six analyzed species, resulting in an 88 % success rate. The number of alleles per polymorphic locus ranged from 2 to 12, and the most successfully amplified loci showed at least one population with a larger number of alleles than were reported in the source species. The population relationships revealed clear genetic clustering in a neighbor-joining tree that was partially incongruent with the taxonomic limits between the P. aurisetus and P. machrisii species, a fact which parallels the problematic taxonomy of the P. aurisetus group. A Bayesian clustering analysis of the individual genotypes confirmed the observed taxonomic incongruence. These microsatellite markers provide a valuable resource for facilitating large-scale genetic studies on population structures, systematics and evolutionary history in this group.

Isolation, characterization, and cross-species amplification of polymorphic microsatellite markers for Pilosocereus machrisii (Cactaceae)

PREMISE OF THE STUDY

Microsatellite primers were developed in Pilosocereus machrisii, a columnar cactus with a patchy distribution in eastern tropical South America, to assess its genetic diversity and population structure.

METHODS AND RESULTS

Eleven microsatellite markers were developed, of which one was monomorphic among 51 individuals from two populations. The number of alleles per polymorphic locus ranged from two to eight, and the total number of alleles was 57. From the 11 isolated loci, nine were successfully amplified in the other four Pilosocereus species.

CONCLUSIONS

The results showed that these markers will be useful for studies of genetic diversity, patterns of gene flow, and population genetic structure in P. machrisii, as well as across other congeneric species.

Pentachlorophenol (PCP) dechlorination in horizontal-flow anaerobic immobilized biomass (HAIB) reactors

This study verifies the potential applicability of horizontal-flow anaerobic immobilized biomass (HAIB) reactors to pentachlorophenol (PCP) dechlorination. Two bench-scale HAIB reactors (R1 and R2) were filled with cubic polyurethane foam matrices containing immobilized anaerobic sludge. The reactors were then continuously fed with synthetic wastewater consisting of PCP, glucose, acetic acid, and formic acid as co-substrates for PCP anaerobic degradation. Before being immobilized in polyurethane foam matrices, the biomass was exposed to wastewater containing PCP in reactors fed at a semi-continuous rate of 2.0 microg PCP g(-1) VS. The applied PCP loading rate was increased from 0.05 to 2.59 mg PCP l(-1)day(-1) for R1, and from 0.06 to 4.15 mg PCP l(-1)day(-1) for R2. The organic loading rates (OLR) were 1.1 and 1.7 kg COD m(-3)day(-1) at hydraulic retention times (HRT) of 24h for R1 and 18 h for R2. Under such conditions, chemical oxygen demand (COD) removal efficiencies of up to 98% were achieved in the HAIB reactors. Both reactors exhibited the ability to remove 97% of the loaded PCP. Dichlorophenol (DCP) was the primary chlorophenol detected in the effluent. The adsorption of PCP and metabolites formed during PCP degradation in the packed bed was negligible for PCP removal efficiency.

Monitoring the trihalomethanes present in water after treatment with chlorine under laboratory condition

In this work assays involving chlorinated water samples, which were previous spiked with humic substances or algae blue green and following the production of the THMs for 30 days is described. To implement the assays, five portions of 1,000 ml of water were stored in glass bottles. The water samples were treated with solutions containing 2, 3, 4 and 5 mg l(-1) chlorine. The samples aliquots (60 ml) were transferred into the glass vials, 10 ml were removed to have a headspace and 100 microl of the 10 mg l(-1) pentafluortoluene bromide solution was added to each vial. The extraction step was performed by adding 10 g of Na(2)SO(4) followed by 5 ml of n-pentane. The vials were stopped with a TFE-faced septum and sealed with aluminum caps. The generated THMs were determined by gas chromatography with electron capture detector using reference solutions with concentration ranging from 8 to 120 microg l(-1) THMs. Three assays were monitored during 30 days and chloroform was the predominant compound found in the water samples, while other species of THMs were not detected. The results showed that when the chlorine concentration was increased in water samples containing algae the concentration of THM varied randomly. Nevertheless, in water samples containing humic substances the increase of the THM concentration presented a relationship with the chlorine concentration. It was also observed that chloroform concentration increased with the elapsed time up to one and six days to water samples spiked with humic substances and algae blue green, respectively and decreased along 30 days. By other hand, assays performed using water samples containing decanted algae material showed that THM was not generated by the chlorine addition.

Population genetic structure of two columnar cacti with a patchy distribution in eastern Brazil

The genetic variability and population genetic structure of six populations of Praecereus euchlorus and Pilosocereus machrisii were investigated. The genetic variability in single populations of Pilosocereus vilaboensis, Pilosocereus aureispinus, and Facheiroa squamosa was also examined. All of these cacti species have a patchy geographic distribution in which they are restricted to small areas of xeric habitats in eastern Brazil. An analysis of genetic structure was used to gain insights into the historical mechanisms responsible for the patchy distribution of P. euchlorus and P. machrisii. High genetic variability was found at the populational level in all species (P=58.9-92.8%, A(p)=2.34-3.33, H(e)=0.266-0.401), and did not support our expectations of low variability based on the small population size. Substantial inbreeding was detected within populations (F(IS)=0.370-0.623). In agreement with their insular distribution patterns, P. euchlorus and P. machrisii had a high genetic differentiation (F(ST)=0.484 and F(ST)=0.281, respectively), with no evidence of isolation by distance. Accordingly, estimates of gene flow (N(m)) calculated from F(ST) and private alleles were below the level of N(m)=1 in P. machrisii and P. euchlorus. These results favored historical fragmentation as the mechanism responsible for the patchy distribution of these two species. The genetic distance between P. machrisii and P. vilaboensis was not compatible with their taxonomic distinction, indicating a possible local speciation event in this genus, or the occurrence of introgression events.

Wing shape heritability and morphological divergence of the sibling species Drosophila mercatorum and Drosophila paranaensis

The fruit-flies Drosophila paranaensis and Drosophila mercatorum pararepleta are sibling species belonging to the repleta group. Females of these two species are normally considered to be morphologically indistinguishable while males only differ consistently in the morphology of their genitalia. These species are sympatric throughout a large area of their geographic distribution. In this study, we investigated the degree of morphological divergence between D. paranaensis and D. mercatorum pararepleta based on morphometric analysis of their wings. The ellipse method was used to describe the placement of the longitudinal and transversal wing veins as well as the size of the wing and the shape of its outline. The heritability under laboratory and field conditions was also estimated from the parameters generated. Multivariate analysis showed that wing morphology possessed sufficient differences to discriminate between the two species with a successful classification rate of 95-98% for females and 82-87% for males. The results of the autoclassification were confirmed by a cross-validation test for females (92-96%). Most measurements possessed significant natural heritability (a mean of 0.48 for D. mercatorum and 0.88 for D. paranaensis), indicating that the variation observed was related to differences in genes acting additively. The principal difference between the two species was in the placement of the posterior transverse wing vein. However, the pattern of morphological variation in the wings of both species was similar, possibly because of shared restrictions in wing development pathways.

Formaldehyde degradation in an anaerobic packed-bed bioreactor

The development of appropriate technologies for the treatment of formaldehyde discharged into the environment is important to minimize its impact. Aerobic systems have been employed, although alternative anaerobic treatments have also been widely studied, mainly due to their low energy consumption and sludge production. However, toxic substances can lead to disturbances in anaerobic reactors. Some research has already been developed on formaldehyde anaerobic biological treatment, but no consensus has yet been reached about its behavior nor has the most efficient system been identified. Aiming at finding supporting evidence for this issue, therefore, this study investigated the degradation and toxicity of formaldehyde in a Horizontal-Flow Anaerobic Immobilized Sludge Reactor. Formaldehyde concentrations of 26.2-1158.6 mg HCHO/L were applied in the reactor, resulting in formaldehyde and chemical oxygen demand removal efficiencies of 99.7% and 92%, respectively. Volatile fatty acids with up to five carbons, found during the degradation of formaldehyde, are believed to indicate that the degradation followed routes unlike those suggested in the literature, which reports the formation of intermediates such as methanol and formic acid. The Monod kinetic model adhered to the experimental data well, with apparent kinetic parameters estimated as r(app)max) = 2.79 x 10(-3) mg HCHO/mg SSVh and K (app)(s) = 242.8 mg HCHO/L.

Relationships between necrotic cactus availability and population size in a cactophilic Drosophila (Diptera, Drosophilidae) located on a sandstone table hill in Brazil

Drosophila gouveai is a cactophilic species endemic to South America. In southeast Brazil it is found on summits of isolated hills, which apparently are current refugia resulting from climatic changes during the Quaternary Period. It breeds only in necrotic cactus cladodes of Pilosocereus machrisii. Temporal differences in necrotic cactus availability could have a great impact upon D. gouveai population size, and could thus influence its evolutionary history. We analyzed the relationship between necrotic cactus availability and population size of D. gouveai. The fluctuation in the population size, variation in necrotic cactus availability and exploitation of this resource by larvae were surveyed bimonthly for one year on a sandstone table hill in central-south Brazil. Temporal necrotic cactus availability did not vary significantly, though in June there was a moderate decrease Larval populations were highest in October and December. The D. gouveai population size was highest in February and remained relatively stable the rest of the year. The observed fluctuation in population size was not a function of temporal necrotic cactus availability in quantitative terms.

Hybrid reactor performance in pentachlorophenol (pcp) removal by anaerobic granules

The present research aimed at evaluating pentachlorophenol (PCP) degradation in a hybrid reactor supplied with a mixture of fatty acids (propionic, butyric, acetic and lactic) and methanol. The performance of the reactor is remarkably stable and efficient during PCP additions at range of 2.0 to 21.0 mg/L. The reduction of chemical oxygen demand (COD) was around 97% and methane was found to be 88% in the biogas production. The efficiency of volatile fatty acids breakdown was 93%, 64% and 74% respectively for butyric, propionic and acetic. PCP total removal of more than 99% was reached by granular sludge activities formed during 21 months of reactor operation. Methanogenic microorganisms predominance was noticed with 10(5) to 10(6) cells/mL during enumeration on methanol or lactate added to sulfate culture media. The removal rate was 1.07 mg PCP.g-1 VS.d-1 during the highest PCP concentration addition.