The discovery of Antarctic RNA viruses: a new game changer

Antarctic ecosystems are dominated by micro‐organisms, and viruses play particularly important roles in the food webs. Since the first report in 2009 (López‐Bueno et al. 2009 ), ‘omic’‐based studies have greatly enlightened our understanding of Antarctic aquatic microbial diversity and ecosystem function (Wilkins et al. 2013 ; Cavicchioli 2015 ). This has included the discovery of many new eukaryotic viruses (López‐Bueno et al. 2009 ), virophage predators of algal viruses (Yau et al. 2011 ), bacteria with resistance to phage (Lauro et al. 2011 ) and mechanisms of haloarchaeal evasion, defence and adaptation to viruses (Tschitschko et al. 2015 ). In this issue of Molecular Ecology, López‐Bueno et al. ( 2015 ) report the first discovery of RNA viruses from an Antarctic aquatic environment. High sequence coverage enabled genome variation to be assessed for four positive‐sense single‐stranded RNA viruses from the order Picornavirales. By examining the populations present in the water column and in the lake's catchment area, populations of ‘quasispecies’ were able to be linked to local environmental factors. In view of the importance of viruses in Antarctic ecosystems but lack of data describing them, this study represents a significant advance in the field.     

Using RAD‐seq to recognize sex‐specific markers and sex chromosome systems

Next‐generation sequencing methods have initiated a revolution in molecular ecology and evolution (Tautz et al. 2010 ). Among the most impressive of these sequencing innovations is restriction site‐associated DNA sequencing or RAD‐seq (Baird et al. 2008 ; Andrews et al. 2016 ). RAD‐seq uses the Illumina sequencing platform to sequence fragments of DNA cut by a specific restriction enzyme and can generate tens of thousands of molecular genetic markers for analysis. One of the many uses of RAD‐seq data has been to identify sex‐specific genetic markers, markers found in one sex but not the other (Baxter et al. 2011 ; Gamble & Zarkower 2014 ). Sex‐specific markers are a powerful tool for biologists. At their most basic, they can be used to identify the sex of an individual via PCR. This is useful in cases where a species lacks obvious sexual dimorphism at some or all life history stages. For example, such tests have been important for studying sex differences in life history (Sheldon 1998 ; Mossman & Waser 1999 ), the management and breeding of endangered species (Taberlet et al. 1993 ; Griffiths & Tiwari 1995 ; Robertson et al. 2006 ) and sexing embryonic material (Hacker et al. 1995 ; Smith et al. 1999 ). Furthermore, sex‐specific markers allow recognition of the sex chromosome system in cases where standard cytogenetic methods fail (Charlesworth & Mank 2010 ; Gamble & Zarkower 2014 ). Thus, species with male‐specific markers have male heterogamety (XY) while species with female‐specific markers have female heterogamety (ZW). In this issue, Fowler & Buonaccorsi ( 2016 ) illustrate the ease by which RAD‐seq data can generate sex‐specific genetic markers in rockfish (Sebastes). Moreover, by examining RAD‐seq data from two closely related rockfish species, Sebastes chrysomelas and Sebastes carnatus (Fig.  1 ), Fowler & Buonaccorsi ( 2016 ) uncover shared sex‐specific markers and a conserved sex chromosome system.     

Improving Bayesian isotope mixing models: a response to Jackson et al. (2009)

We recently described a Bayesian framework for stable isotope mixing models and provided a software tool, MixSIR, for conducting such analyses (Ecol. Lett., 2008; 11 :470). Jackson et al. (Ecol. Lett., 2009; 12:E1) criticized the performance of our software based on tests using simulated data. However, their simulation data were flawed, rendering claims of erroneous behaviour inaccurate. A re‐evaluation of the MixSIR source code did, however, uncover two minor coding errors, which we have fixed. When data are correctly simulated according to eqns  (1)–(4) in Jackson et al. (2009) , MixSIR consistently and accurately estimated the proportional contribution of prey to a predator diet, and was surprisingly robust to additional unquantified error. Jackson et al. (2009) also suggested we use a Dirichlet prior on the source proportion parameters, which we agree with. Finally, Jackson et al. (2009) propose adding additional error parameters to our mixing model framework. We caution that such increases in model complexity should be evaluated based on data support.     

Host plant richness explains diversity of ectomycorrhizal fungi: Response to the comment of Tedersoo et al. (2014)

Exploring the relationships between the biodiversity of groups of interacting organisms yields insight into ecosystem stability and function (Hooper et al. 2000 ; Wardle 2006 ). We demonstrated positive relationships between host plant richness and ectomycorrhizal (EM) fungal diversity both in a field study in subtropical China (Gutianshan) and in a meta‐analysis of temperate and tropical studies (Gao et al. 2013 ). However, based on re‐evaluation of our data sets, Tedersoo et al. ( 2014 ) argue that the observed positive correlation between EM fungal richness and EM plant richness at Gutianshan and also in our metastudies was based mainly from (i) a sampling design with inconsistent species pool and (ii) poor data compilation for the meta‐analysis. Accordingly, we checked our data sets and repeated the analysis performed by Tedersoo et al. ( 2014 ). In contrast to Tedersoo et al. ( 2014 ), our re‐analysis still confirms a positive effect of plant richness on EM fungal diversity in Gutianshan, temperate and tropical ecosystems, respectively.     

Demystifying the RAD fad

We are writing in response to the population and phylogenomics meeting review by Andrews & Luikart ( 2014 ) entitled ‘Recent novel approaches for population genomics data analysis’. Restriction‐site‐associated DNA (RAD) sequencing has become a powerful and useful approach in molecular ecology, with several different published methods now available to molecular ecologists, none of which can be considered the best option in all situations. A&L report that the original RAD protocol of Miller et al. ( 2007 ) and Baird et al. ( 2008 ) is superior to all other RAD variants because putative PCR duplicates can be identified (see Baxter et al. 2011 ), thereby reducing the impact of PCR artefacts on allele frequency estimates (Andrews & Luikart 2014 ). In response, we (i) challenge the assertion that the original RAD protocol minimizes the impact of PCR artefacts relative to that of other RAD protocols, (ii) present additional biases in RADseq that are at least as important as PCR artefacts in selecting a RAD protocol and (iii) highlight the strengths and weaknesses of four different approaches to RADseq which are a representative sample of all RAD variants: the original RAD protocol (mbRAD, Miller et al. 2007 ; Baird et al. 2008 ), double digest RAD (ddRAD, Peterson et al. 2012 ), ezRAD (Toonen et al. 2013 ) and 2bRAD (Wang et al. 2012 ). With an understanding of the strengths and weaknesses of different RAD protocols, researchers can make a more informed decision when selecting a RAD protocol.     

Massive parallel MHC genotyping: titanium that shines

The power of population genetic analyses is often limited by sample size resulting from constraints in financial resources and time to genotype large numbers of individuals. This particularly applies to nonmodel species where detailed genomic knowledge is lacking. Next‐generation sequencing technology using primers ‘tagged’ with an individual barcode of a few nucleotides offers the opportunity to genotype hundreds of individuals at several loci in parallel ( Binladen et al. 2007 ; Meyer et al. 2008 ). The large number of sequence reads can also be used to identify artefacts by frequency distribution thresholds intrinsically determined for each run and data set. In Babik et al. (2009 ), next‐generation deep sequencing was used to genotype several major histocompatibility complex (MHC) class IIB loci of the European bank vole ( Fig. 1 ). Their approach can be useful for many researchers working with complex multiallelic templates and large sample sizes.

Figure 1 Open in figure viewer PowerPoint Hypothetical example of parallel genotyping of two individuals using individually bar‐coded primers. Polymerase chain reactions (PCRs) are performed separately for each individual using a forward primer with a unique Tag‐sequence of four nucleotides. After sequencing of pooled PCR products, sequences can be sorted by their forward primer Tag (Tag‐sorting error rate was estimated < 0.1%). Rare sequences most likely represent artefacts and due to the large amount of sequences obtained (up to 10⁶) the artefact threshold can be determined intrinsically for each data set and was estimated to be around 3% in the case of bank vole MHC class IIB genes ( Babik et al. 2009 ). Photos by Gabriela Bydlon.     

Expanding the population genetic perspective of cnidarian‐Symbiodinium symbioses

The modern synthesis was a seminal period in the biological sciences, establishing many of the core principles of evolutionary biology that we know today. Significant catalysts were the contributions of R.A. Fisher, J.B.S. Haldane and Sewall Wright (and others) developing the theoretical underpinning of population genetics, thus demonstrating adaptive evolution resulted from the interplay of forces such as natural selection and mutation within groups of individuals occupying the same space and time (i.e. a population). Given its importance, it is surprising that detailed population genetic data remain lacking for numerous organisms vital to many ecosystems. For example, the coral reef ecosystem is well recognized for its high biodiversity and productivity, numerous ecological services and significant economic and societal values (Moberg & Folke 1999; Cinner 2014). Many coral reef invertebrates form symbiotic relationships with single‐celled dinoflagellates within the genus Symbiodinium Freudenthal (Taylor 1974), with hosts providing these (typically) intracellular symbionts with by‐products of metabolism and in turn receiving photosynthetically fixed carbon capable of meeting hosts’ respiratory demands (Falkowski et al. 1984; Muscatine et al. 1984). Unfortunately, the health and integrity of the coral reef ecosystem has been significantly and negatively impacted by onslaughts like anthropogenic eutrophication and disease in addition to global climate change, with increased incidences of ‘bleaching’ events (characterized as the loss of photosynthetic pigments from the algal cell or massive reduction of Symbiodinium density from hosts’ tissue) and host mortality leading to staggering declines in geographic coverage (Bruno & Selig 2007) that have raised questions on the viability of this ecosystem as we know it (Bellwood et al. 2004; Parmesan 2006). One avenue towards anticipating the future of the coral reef ecosystem is by developing a broader and deeper understanding of the current genotypic diversity encompassed within and between populations of their keystone species, the scleractinian corals and dinoflagellate symbionts, as they potentially possess functional variation (either singularly or in combination) that may come under selection due to the ongoing and rapid environmental changes they are experiencing. However, such studies, especially for members of the genus Symbiodinium, are sparse. In this issue, Baums et al. (2014) provide a significant contribution by documenting the range‐wide population genetics of Symbiodinium ‘fitti’ (Fig. 1 ) in the context of complementary data from its host, the endangered Caribbean elkhorn coral Acropora palmata (Fig. 1 ). Notable results of this study include a single S. ‘fitti’ genotype typically dominates an individual A. palmata colony both spatially and temporally, gene flow among coral host populations is a magnitude higher to that of its symbiont populations, and the partners possess disparate patterns of genetic differentiation across the Greater Caribbean. The implications of such findings are discussed herein.     

Will reduced host connectivity curb the spread of a devastating epidemic?

The white‐nose syndrome (WNS), caused by the fungal pathogen Pseudogymnoascus destructans, is threatening the cave‐dwelling bat fauna of North America by killing individuals by the thousands in hibernacula each winter since its appearance in New York State less than ten years ago. Epidemiological models predict that WNS will reach the western coast of the USA by 2035, potentially eliminating most populations of susceptible bat species in its path (Frick et al. 2015; O'Regan et al. 2015). These models were built and validated using distributional data from the early years of the epidemic, which spread throughout eastern North America following a route driven by cave density and winter severity (Maher et al. 2012). In this issue of Molecular Ecology, Wilder et al. (2015) refine these findings by showing that connectivity among host populations, as assessed by population genetic markers, is crucial in determining the spread of the pathogen. Because host connectivity is much reduced in the hitherto disease free western half of North America, Wilder et al. make the reassuring prediction that the disease will spread more slowly west of the Great Plains.     

Mock communities highlight the diversity of host‐associated eukaryotes

Host‐associated microbes are ubiquitous. Every multicellular eukaryote, and even many unicellular eukaryotes (protists), hosts a diverse community of microbes. High‐throughput sequencing (HTS) tools have illuminated the vast diversity of host‐associated microbes and shown that they have widespread influence on host biology, ecology and evolution (McFall‐Ngai et al. 2013 ). Bacteria receive most of the attention, but protists are also important components of microbial communities associated with humans (Parfrey et al. 2011 ) and other hosts. As HTS tools are increasingly used to study eukaryotes, the presence of numerous and diverse host‐associated eukaryotes is emerging as a common theme across ecosystems. Indeed, HTS studies demonstrate that host‐associated lineages account for between 2 and 12% of overall eukaryotic sequences detected in soil, marine and freshwater data sets, with much higher relative abundances observed in some samples (Ramirez et al. 2014 ; Simon et al. 2015 ; de Vargas et al. 2015 ). Previous studies in soil detected large numbers of predominantly parasitic lineages such as Apicomplexa, but did not delve into their origin [e.g. (Ramirez et al. 2014 )]. In this issue of Molecular Ecology, Geisen et al. ( 2015 ) use mock communities to show that many of the eukaryotic organisms detected by environmental sequencing in soils are potentially associated with animal hosts rather than free‐living. By isolating the host‐associated fraction of soil microbial communities, Geisen and colleagues help explain the surprisingly high diversity of parasitic eukaryotic lineages often detected in soil/terrestrial studies using high‐throughput sequencing (HTS) and reinforce the ubiquity of these host‐associated microbes. It is clear that we can no longer assume that organisms detected in bulk environmental sequencing are free‐living, but instead need to design studies that specifically enumerate the diversity and function of host‐associated eukaryotes. Doing so will allow the field to determine the role host‐associated eukaryotes play in soils and other environments and to evaluate hypotheses on assembly of host‐associated communities, disease ecology and more.     

Coalescing molecular evolution and DNA barcoding

The DNA barcoding concept (Woese et al. 1990 ; Hebert et al. 2003 ) has considerably boosted taxonomy research by facilitating the identification of specimens and discovery of new species. Used alone or in combination with DNA metabarcoding on environmental samples (Taberlet et al. 2012 ), the approach is becoming a standard for basic and applied research in ecology, evolution and conservation across taxa, communities and ecosystems (Scheffers et al. 2012 ; Kress et al. 2015 ). However, DNA barcoding suffers from several shortcomings that still remain overlooked, especially when it comes to species delineation (Collins & Cruickshank 2012 ). In this issue of Molecular Ecology, Barley & Thomson ( 2016 ) demonstrate that the choice of models of sequence evolution has substantial impacts on inferred genetic distances, with a propensity of the widely used Kimura 2‐parameter model to lead to underestimated species richness. While DNA barcoding has been and will continue to be a powerful tool for specimen identification and preliminary taxonomic sorting, this work calls for a systematic assessment of substitution models fit on barcoding data used for species delineation and reopens the debate on the limitation of this approach.     

Army ant invasions reveal phylogeographic processes across the Isthmus of Panama

Female army ants cannot fly, making them very poor dispersers across water barriers. This dependence on terrestrial corridors motivated the investigation by Winston et al. ( 2017 ), published in this issue of Molecular Ecology, into the role of Panamanian isthmus formation in the diversification of Eciton army ants. Complete closure of this isthmus occurred around three million years ago (3 Ma), but it has also been hypothesized that earlier, temporary land connections facilitated additional colonization events between South and Central America over the past 13 million years or more. The phylogenomic and population genomic analyses by Winston et al. ( 2017 ) uncovered multiple incursions of Eciton lineages into Central America between 4 and 7 Ma. Their study contributes to a growing body of evidence arguing that transitory land bridges predating 3 Ma supported substantial intercontinental biotic exchange.     

Process,correlation and parameter fitting in species distribution models: a response to Kriticos et al

In a recent article (Dormann et al., 2012, Journal of Biogeography, 39, 2119–2131), we compared different approaches to species distribution modelling and depicted modelling approaches along an axis from purely ‘correlative’ to ‘forward process‐based’ models. In their correspondence, Kriticos et al. (2013, Journal of Biogeography, doi: 10.1111/j.1365‐2699.2012.02791.x ) challenge this view, claiming that our continuum representation neglects differences among models and does not consider the ability of fitted process‐based models to combine the advantages of both process‐based and correlative modelling approaches. Here we clarify that the continuum view resulted from recognition of the manifold differences between models. We also reinforce the point that the current trend towards combining different modelling approaches may lead not only to the desired combination of the advantages but also to the accumulation of the disadvantages of those approaches. This point has not been made sufficiently clear previously.     

A harvest of weeds yields insight into a case of contemporary evolution

When Charles Darwin was exploring the idea of evolution via natural selection, he looked to domesticated species, with the opening chapter of The Origin of Species titled ‘Variation Under Domestication’ (Darwin 1859 ). Domesticated species such as crops are a great example of artificial selection, which Darwin realized was analogous to natural selection. But growing among those carefully selected crop varieties are the unwelcome and unwanted plants we call weeds. Despite the importance of weeds and long‐standing interest in their evolution (Baker 1974 ), we still know little about how agricultural weeds evolve, and we often fail to take evolution into account when attempting to manage them (Neve et al. 2009 ). Agricultural weeds are subjected to the unique conditions of farm fields, such as frequent soil disturbance and the addition of water and nutrients. They are also confronted with aggressive attempts at their removal via herbicides and mechanical means. As such, they are under intense demographic and selective pressure and can potentially rapidly evolve in response. In this issue of Molecular Ecology, Kuester and co‐authors make a rare attempt to understand contemporary evolution in an agricultural weed (Kuester et al. 2016 ). They do so using the powerful resurrection approach of comparing ancestors and descendants under common conditions (Franks et al. 2008 ). They sampled multiple populations of the weedy plant Ipomoea purpurea at two points in time. A comparison of these greenhouse‐grown ancestor and descendent populations showed that, over time, populations had lost significant levels of neutral genetic diversity, consistent with genetic bottlenecks. The authors also found a slight increase, on average, of resistance to the herbicide glyphosate, which is the active ingredient in Roundup^®. This work is one of a growing number of studies demonstrating rapid evolution in natural populations (Thompson 2013 ) and also reveals evidence of both selection and drift in populations of an agricultural weed.     

Novel aspects in the life cycle and biotrophic interactions in Pezizomycetes (Ascomycota,Fungi)

The ascomycete class Pezizomycetes (single order Pezizales) is known for its cup‐shaped fruit bodies and the evolution of edible truffles and morels, but little is known about the ontogeny and ecology of this large and ecologically diverse fungal group. In this issue of Molecular Ecology, Healy et al. ( 2013 ) make a great leap forward by describing and identifying asexual, anamorphic structures that produce mitotic spores in many ectomycorrhiza‐forming truffle and nontruffle species on soil surfaces worldwide (Fig.  1 ). Although such anamorphic forms have been reported sporadically from certain ectomycorrhizal and saprotrophic Pezizomycetes (e.g. Warcup 1990 ), Healy et al. ( 2013 ) demonstrate that these terricolous asexual forms are both taxonomically and geographically more widespread and, in fact, much more common than previously understood. We anticipate that deeper insight into other substrates, provided by molecular analyses of materials such as dead wood and seeds, is likely to reveal numerous anamorphs of saprotrophic and pathogenic Pezizomycetes as well (see Marek et al. 2009 ).     

Population genomic analyses reveal possible drivers of population divergence

Recent advances in sequencing technology and efficiency enable new and improved methods to investigate how populations diverge and species evolve. Fungi have relatively small and simple genomes and can often be cultured in the laboratory. Fungal populations can thus be sequenced for a relatively low cost, which makes them ideal for population genomic analyses. In several recent population genomic studies, wild populations of fungal model organisms and human pathogens have been analysed, for example Neurospora crassa (Ellison et al. 2011 ), Saccharomyces uvarum (Almeida et al. 2014 ), Coccidioides spp. (Neafsey et al. 2010 ) and Cryptococcus gatti (Engelthaler et al. 2014 ). In this issue of Molecular Ecology, Branco et al. ( 2015 ) apply population genomic tools to understand population divergence and adaptation in a symbiotic (mycorrhizal) fungus. This study exemplifies the possibilities of diving deeper into the genomic features involved in population divergence and speciation, also for nonmodel organisms, and how molecular and analytical tools will improve our understanding of the patterns and mechanisms that underlie adaptation to habitats, population divergence and dispersal limitation of fungi.     

Enalikter is not an annelid: homology,autapomorphies and the interpretation of problematic fossils

A megacheiran arthropod, Enalikter aphson, was recently described by Siveter et al. (2014) from the mid‐Silurian (late Wenlock) of Herefordshire. Previously, megacheirans had only been recognized from the Cambrian. Struck et al. (2015) considered the body plan of Enalikter to be incompatible with this affinity, arguing that many of the arthropod features were either not present or misinterpreted. Instead, they compared Enalikter to polychaete annelids, identifying characters from numerous polychaete lineages which they considered to be present in Enalikter. A reply to this critique by Siveter et al. (2015) reaffirmed arthropod affinities for Enalikter by presenting additional evidence for key arthropod features, such as arthropodized appendages. Here, we augment Siveter et al. by critically addressing the putative annelid characters of Enalikter presented by Struck et al. and additionally explore the morphological and phylogenetic implications of their hypothesis. We conclude that similarities between Enalikter and polychaetes are superficial and that character combinations proposed by Struck et al. are not present in any annelid, living or extinct. This taxon highlights the importance of using a phylogenetic framework for interpreting fossils that present unusual morphologies, such that proposed shared characters are hypotheses of homology rather than merely phenotypic similarities. Crucially, we argue that autapomorphic characters of subgroups of large taxa (like families or classes within phyla) should not be used to diagnose problematic fossils.     

Homomorphic plant sex chromosomes are coming of age

Sex chromosomes are a very peculiar part of the genome that have evolved independently in many groups of animals and plants (Bull 1983 ). Major research efforts have so far been focused on large heteromorphic sex chromosomes in a few animal and plant species (Chibalina & Filatov 2011 ; Zhou & Bachtrog 2012 ; Bellott et al. 2014 ; Hough et al. 2014 ; Zhou et al. 2014 ), while homomorphic (cytologically indistinguishable) sex chromosomes have largely been neglected. However, this situation is starting to change. In this issue, Geraldes et al. ( 2015 ) describe a small (~100 kb long) sex‐determining region on the homomorphic sex chromosomes of poplars (Populus trichocarpa and related species, Fig.  1 ). All species in Populus and its sister genus Salix are dioecious, suggesting that dioecy and the sex chromosomes, if any, should be relatively old. Contrary to this expectation, Geraldes et al. ( 2015 ) demonstrate that the sex‐determining region in poplars is of very recent origin and probably evolved within the genus Populus only a few million years ago.     

A tale of two lineages: unexpected,long‐term persistence of the amphibian‐killing fungus in Brazil

For the past 17 years, scientists have been compiling a list of amphibian species susceptible to infection by the amphibian‐killing chytrid fungus, Batrachochytrium dendrobatidis (Bd), all over the world, with >500 species infected on every continent except Antarctica (Olson et al. 2013 ). Where Bd has been found, the impacts on amphibians has been one of two types: either Bd arrives into a naïve amphibian population followed by a mass die‐off and population declines (e.g. Lips et al. 2006 ), or Bd is present at some moderate prevalence, usually infecting many species but at apparently nonlethal intensities for a long time. In this issue of Molecular Ecology, Rodriguez et al. ( 2014 ) discover that the Atlantic Coastal Forest of Brazil is home to two Bd lineages: the Global Pandemic Lineage (Bd‐GPL) – the strain responsible for mass die‐offs and population declines – and a lineage endemic to Brazil (Bd‐Bz). Even more surprising was that both lineages have been present in this area for the past 100 years, making these the oldest records of Bd infecting amphibians. The team also described a moderate but steady prevalence of ~20% across all sampled anuran families for over 100 years, indicating that Brazil has been in an enzootic disease state for over a century. Most amphibians were infected with Bd‐GPL, suggesting this lineage may be a better competitor than Bd‐Bz or may be replacing the Bd‐Bz lineage. Rodriguez et al. ( 2014 ) also detected likely hybridization of the two Bd lineages, as originally described by Schloegel et al. ( 2012 ).     

Morphological Characterizations of Four Species of Parallelostrombidium (Ciliophora,Oligotrichia), with a Note on the Phylogeny of the Genus

The morphology and phylogeny of four oligotrichid ciliates, Parallelostrombidium paraellipticum sp. n., P. dragescoi sp. n., P. jankowskii (Xu et al. 2009) comb. n., and P. kahli (Xu et al. 2009) comb. n., are described or redescribed based on live observation, protargol stained material, and SSU rRNA gene sequences. The new species P. paraellipticum sp. n. is characterized by its obovoidal cell shape, adoral zone composed of 17–21 collar, 9–11 buccal, and two thigmotactic membranelles, and extrusomes attached in one row along the girdle kinety. The new species P. dragescoi sp. n. is distinguished from its congeners by its obovoidal cell shape and a lack of thigmotactic membranelles. Based on ciliary patterns recognizable in the original slides, Omegastrombidium jankowskii Xu et al. 2009 and O. kahli Xu et al. 2009 should be transferred to the genus Parallelostrombidium Agatha 2004. Phylogenetic analyses based on SSU rRNA gene sequence data demonstrate that all four new sequences cluster with previously described congeners. The genus Parallelostrombidium is separated into two clusters, suggesting its non‐monophyly and probably corresponding to the two subgenera proposed by Agatha and Strüder‐Kypke (2014), as well as their morphological difference (cell dorsoventrally flattened vs. unflattened).