Trait‐based and phylogenetic associations between parasites and their hosts: a case study with small mammals and fleas in the Palearctic

We investigated the associations between ecological (density, shelter structure), morphological (body mass, hair morphology) and physiological traits (basal metabolic rate) of small mammals and ecological (seasonality of reproduction, microhabitat preferences, abundance, host specificity) and morphological (presence and number of combs) traits of their flea parasites that shape host selection processes by fleas. We adapted the extended version of the three‐table ordination and linked species composition of flea assemblages of host species with traits and phylogenies of both hosts and fleas. Fleas with similar trait values, independent of phylogenetic affinities, were clustered on the same host species. Fleas possessing certain traits selected hosts possessing certain traits. Fleas belonging to the same phylogenetic lineage were found on the same host more often than expected by chance. Certain phylogenetic lineages of hosts harbored certain phylogenetic lineages of fleas. The process of host selection by fleas appeared to be determined by reciprocal relationships between host and flea traits, as well as between host and flea phylogenies. We concluded that the connection between host and flea phylogenies, coupled with the connection between host and flea traits, suggests that the species compositions of the host spectra of fleas were driven by the interaction between historical processes and traits.     

New insights into the phylogenetic relationships among the oceanic dolphins (Cetacea: Delphinidae)

The phylogenetic relationships of oceanic dolphins (family Delphinidae) remain unclear. Several works using mitochondrial and/or nuclear DNA on different genera and species have been published, though no consensus exists regarding even the subfamilies that conform the family. Here, a new phylogeny for the family Delphinidae, including 36 different complete mitochondrial genomes (plus two outgroups), was constructed under Bayesian and maximum likelihood approaches. Results indicate identical tree topology in both cases, with almost all nodes fully supported independently of the reconstruction approach. This topology is different from those previously published and proposes new phylogenetic relationships among subfamilies, genera and species of the family. These findings are critically important for the study of oceanic dolphin taxonomy, ecology, evolution and conservation, and highlight the importance of revisiting and resolving uncertain phylogenies.     

Assessing the performance of single-copy genes for recovering robust phylogenies

Phylogenies involving nonmodel species are based on a few genes, mostly chosen following historical or practical criteria. Because gene trees are sometimes incongruent with species trees, the resulting phylogenies may not accurately reflect the evolutionary relationships among species. The increase in availability of genome sequences now provides large numbers of genes that could be used for building phylogenies. However, for practical reasons only a few genes can be sequenced for a wide range of species. Here we asked whether we can identify a few genes, among the single-copy genes common to most fungal genomes, that are sufficient for recovering accurate and well-supported phylogenies. Fungi represent a model group for phylogenomics because many complete fungal genomes are available. An automated procedure was developed to extract single-copy orthologous genes from complete fungal genomes using a Markov Clustering Algorithm (Tribe-MCL). Using 21 complete, publicly available fungal genomes with reliable protein predictions, 246 single-copy orthologous gene clusters were identified. We inferred the maximum likelihood trees using the individual orthologous sequences and constructed a reference tree from concatenated protein alignments. The topologies of the individual gene trees were compared to that of the reference tree using three different methods. The performance of individual genes in recovering the reference tree was highly variable. Gene size and the number of variable sites were highly correlated and significantly affected the performance of the genes, but the average substitution rate did not. Two genes recovered exactly the same topology as the reference tree, and when concatenated provided high bootstrap values. The genes typically used for fungal phylogenies did not perform well, which suggests that current fungal phylogenies based on these genes may not accurately reflect the evolutionary relationships among species. Analyses on subsets of species showed that the phylogenetic performance did not seem to depend strongly on the sample. We expect that the best-performing genes identified here will be very useful for phylogenetic studies of fungi, at least at a large taxonomic scale. Furthermore, we compare the method developed here for finding genes for building robust phylogenies with previous ones and we advocate that our method could be applied to other groups of organisms when more complete genomes are available.     

Co‐occurrence and phylogenetic distance in communities of mammalian ectoparasites: limiting similarity versus environmental filtering

Similarity between species plays a key role in the processes governing community assembly. The co‐occurrence of highly similar species may be unlikely if their similar needs lead to intense competition (limiting similarity). On the other hand, persistence in a particular habitat may require certain traits, such that communities end up consisting of species sharing the same traits (environmental filtering). Relatively little information exists on the relative importance of these processes in structuring parasite communities. Assuming that phylogenetic relatedness reflects ecological similarity, we tested whether the co‐occurrence of pairs of flea species (Siphonaptera) on the same host individuals was explained by the phylogenetic distance between them, among 40 different samples of mammalian hosts (rodents and shrews) from different species, areas or seasons. Our results indicate that frequency of co‐occurrence between flea species increased with decreasing phylogenetic distance between them in 37 out of 40 community samples, with 14 of these correlations being statistically significant. A meta‐analysis across all samples confirmed the overall trend for closely related species to co‐occur more frequently on the same individual hosts than expected by chance, independently of the identity of the host species or of environmental conditions. These findings suggest that competition between closely related, and therefore presumably ecologically similar, species is not important in shaping flea communities. Instead, if only fleas with certain behavioural, ecological and physiological properties can encounter and exploit a given host, and if phylogenetic relationships determine trait similarity among flea species, then a process akin to environmental filtering, or host filtering, could favour the co‐occurrence of related species on the same host.     

Molecular Phylogeny and Systematics of Anoplocephaline Cestodes in Rodents and Lagomorphs

A molecular phylogenetic hypothesis is presented for the anoplocephaline cestodes of placental mammals based on sequence data from the mitochondrial cytochrome c oxidase I (COI) gene, the nuclear-encoded 28S rRNA gene and the internal transcribed spacer region I of rRNA (ITS1). The material consists of 35 species representing nine genera of cestodes, with emphasis on taxa parasitising rodents and lagomorphs in the Holarctic region. The resulting phylogenies show considerable disagreement with earlier systematic and phylogenetic hypotheses derived from morphology. Specifically, the results contradict the view of uterine morphology being the primary determinant of deeper phylogenetic splits within Anoplocephalinae. Also, the role of genital duplication as a means of generic divergence was not found to follow consistently the pattern suggested by earlier hypotheses. Colonisation of novel host lineages has evidently been the predominant mode of diversification in anoplocephaline cestodes of placental mammals; evidence for phyletic co-evolution was obscure. The phylogenies consistently distinguished a large monophyletic group including all species from arvicoline rodents (voles and lemmings), primarily representing the genera Anoplocephaloides Baer, 1923 and Paranoplocephala Lühe, 1910. Phylogenetic relationships within the “arvicoline clade” of cestodes were generally poorly resolved. Consistent support for nodes above and below the unresolved polytomy indicates a rapid radiation involving a nearly simultaneous diversification of many lineages, a scenario also proposed for the arvicoline hosts.     

High-throughput sequencing of six bamboo chloroplast genomes: phylogenetic implications for temperate woody bamboos (Poaceae: Bambusoideae)

Background

Bambusoideae is the only subfamily that contains woody members in the grass family, Poaceae. In phylogenetic analyses, Bambusoideae, Pooideae and Ehrhartoideae formed the BEP clade, yet the internal relationships of this clade are controversial. The distinctive life history (infrequent flowering and predominance of asexual reproduction) of woody bamboos makes them an interesting but taxonomically difficult group. Phylogenetic analyses based on large DNA fragments could only provide a moderate resolution of woody bamboo relationships, although a robust phylogenetic tree is needed to elucidate their evolutionary history. Phylogenomics is an alternative choice for resolving difficult phylogenies.

Methodology/Principal Findings

Here we present the complete nucleotide sequences of six woody bamboo chloroplast (cp) genomes using Illumina sequencing. These genomes are similar to those of other grasses and rather conservative in evolution. We constructed a phylogeny of Poaceae from 24 complete cp genomes including 21 grass species. Within the BEP clade, we found strong support for a sister relationship between Bambusoideae and Pooideae. In a substantial improvement over prior studies, all six nodes within Bambusoideae were supported with ≥0.95 posterior probability from Bayesian inference and 5/6 nodes resolved with 100% bootstrap support in maximum parsimony and maximum likelihood analyses. We found that repeats in the cp genome could provide phylogenetic information, while caution is needed when using indels in phylogenetic analyses based on few selected genes. We also identified relatively rapidly evolving cp genome regions that have the potential to be used for further phylogenetic study in Bambusoideae.

Conclusions/Significance

The cp genome of Bambusoideae evolved slowly, and phylogenomics based on whole cp genome could be used to resolve major relationships within the subfamily. The difficulty in resolving the diversification among three clades of temperate woody bamboos, even with complete cp genome sequences, suggests that these lineages may have diverged very rapidly.     

Multigene analysis of phylogenetic relationships and divergence times of primate sucking lice (Phthiraptera: Anoplura)

Cospeciation between hosts and parasites offers a unique opportunity to use information from parasites to infer events in host evolutionary history. Although lice (Insecta: Phthiraptera) are known to cospeciate with their hosts and have frequently served as important markers to infer host evolutionary history, most molecular studies are based on only one or two markers. Resulting phylogenies may, therefore, represent gene histories (rather than species histories), and analyses of multiple molecular markers are needed to increase confidence in the results of phylogenetic analyses. Herein, we phylogenetically examine nine molecular markers in primate sucking lice (Phthiraptera: Anoplura) and we use these markers to estimate divergence times among louse lineages. Individual and combined analyses of these nine markers are, for the most part, congruent, supporting relationships hypothesized in previous studies. Only one marker, the nuclear protein-coding gene Histone 3, has a significantly different tree topology compared to the other markers. The disparate evolutionary history of this marker, however, has no significant effect on topology or nodal support in the combined phylogenetic analyses. Therefore, phylogenetic results from the combined data set likely represent a solid hypothesis of species relationships. Additionally, we find that simultaneous use of multiple markers and calibration points provides the most reliable estimates of louse divergence times, in agreement with previous studies estimating divergences among species. Estimates of phylogenies and divergence times also allow us to verify the results of [Reed, D.L., Light, J.E., Allen, J.M., Kirchman, J.J., 2007. Pair of lice lost or parasites regained: the evolutionary history of anthropoid primate lice. BMC Biol. 5, 7.]; there was probable contact between gorilla and archaic hominids roughly 3 Ma resulting in a host switch of Pthirus lice from gorillas to archaic hominids. Thus, these results provide further evidence that data from cospeciating organisms can yield important information about the evolutionary history of their hosts.     

The role of phylogeny in desert rodent community assembly

Recent advances in molecular genetics and phylogenetic reconstruction have the potential to transform ecology by providing new insights into the historical evolution of ecological communities. This study by Stevens and collaborators complements decades of previous research on desert rodents, by combining data from a field study and a phylogenetic tree for Mojave Desert rodents to address patterns and processes of community assembly. The number of coexisting rodent species is positively correlated, and the average phylogenetic distance among these species is negatively correlated with perennial plant species richness. As rodent species diversity increases along a gradient of increasing environmental heterogeneity, communities are composed of increasingly related species: there is a consistent pattern of phylogenetic structure from over-dispersed through random to clumped. I discuss this pattern in the light of complementary results of previous studies. This paper is noteworthy for calling attention to still unanswered questions about how the historical events of speciation, colonization, extinction, and trait evolution and their relationship to past climates and vegetation have given rise to current patterns of community organization.     

Assessing phylogenetic resolution among mitochondrial, nuclear, and morphological datasets in Nothonotus darters (Teleostei: Percidae)

External morphological characters are the basis of our understanding of diversity and species relationships in many darter clades. The past decade has seen the publication of many studies utilizing mtDNA sequence data to investigate darter phylogenetics, but only recently have nuclear genes been used to investigate darter relationships. Despite a long tradition of use in darter systematics few studies have examined the phylogenetic utility of external morphological characters in estimating relationships among species in darter clades. We present DNA sequence data from the mitochondrial cytochrome b (cytb) gene, the nuclear encoded S7 intron 1, and discretely coded external morphological characters for all 20 species in the darter clade Nothonotus. Bayesian phylogenetic analyses result in phylogenies that are in broad agreement with previous studies. The cytb gene tree is well resolved, while the nuclear S7 gene tree lacks phylogenetic resolution, node support, and is characterized by a lack of reciprocal monophyly for many of the Nothonotus species. The phylogenies resulting from analysis of the morphological dataset lack resolution, but nodes present are found in the cytb and S7 gene trees. The highest resolution and node support is found in the Bayesian combined data phylogeny. Based on our results we propose continued exploration of the phylogenetic utility of external morphological characters in other darter clades. Given the extensive lack of reciprocal monophyly of species observed in the S7 gene tree we predict that nuclear gene sequences may have limited utility in intraspecific phylogeographic studies of Nothonotus darters.     

The Tree versus the Forest: The Fungal Tree of Life and the Topological Diversity within the Yeast Phylome

A recurrent topic in phylogenomics is the combination of various sequence alignments to reconstruct a tree that describes the evolutionary relationships within a group of species. However, such approach has been criticized for not being able to properly represent the topological diversity found among gene trees. To evaluate the representativeness of species trees based on concatenated alignments, we reconstruct several fungal species trees and compare them with the complete collection of phylogenies of genes encoded in the Saccharomyces cerevisiae genome. We found that, despite high levels of among-gene topological variation, the species trees do represent widely supported phylogenetic relationships. Most topological discrepancies between gene and species trees are concentrated in certain conflicting nodes. We propose to map such information on the species tree so that it accounts for the levels of congruence across the genome. We identified the lack of sufficient accuracy of current alignment and phylogenetic methods as an important source for the topological diversity encountered among gene trees. Finally, we discuss the implications of the high levels of topological variation for phylogeny-based orthology prediction strategies.     

An investigation of the co-evolutionary relationships between onchobothriid tapeworms and their elasmobranch hosts

There is general consensus that the living elasmobranchs comprise a monophyletic taxon. There is evidence that, among tetraphyllidean tapeworms, the approximately 201 hooked species (Onchobothriidae) may also comprise a monophyletic group. Determinations of host specificity are contingent upon correct specific identifications. Since 1960, over 200 new elasmobranch species and over 100 new onchobothriid species have been described. Some confidence can be placed in host and parasite identifications of recent studies, but specific identifications provided in older literature in many cases are suspect. There is some consensus among published works on the phylogenetic relationships among elasmobranchs. Phylogenetic relationships among onchobothriids remain largely unresolved. Elasmobranchs have been poorly sampled for onchobothriids; records exist for approximately 20% of the 911 species and approximately 44% of the 170 elasmobranch genera. Onchobothriids are remarkably host specific, exhibiting essentially oioxenous specificity for their definitive hosts. Multiple onchobothriid species commonly parasitise the same host species; in some cases these are congeners, in other cases these are members of two different onchobothriid genera. There is substantial incongruence between available host and parasite phylogenies. For example, Acanthobothrium is by far the most ubiquitous onchobothriid genus, parasitising almost all orders of elasmobranchs known to host onchobothriids, yet, there is no evidence of major clades of Acanthobothrium corresponding to postulated major subgroupings of elasmobranchs (e.g. Galea and Squalea or sharks and rays). Potamotrygonocestus appears to be among the most basal onchobothriid groups, yet it parasitises one of the most derived elasmobranch groups (the freshwater stingray genus Potamotrygon). It appears that congeners parasitising the same host species are not necessarily each other's closest relatives. At this point the preliminary and limited available data suggest that, at least in this system, strict host specificity is not necessarily indicative of strict co-evolution. This study was extremely limited by the lack of available robust phylogenies for onchobothriids and elasmobranchs.     

Phylogenetic analysis of RAD‐seq data: examining the influence of gene genealogy conflict on analysis of concatenated data

One of the major issues in phylogenetic analysis is that gene genealogies from different gene regions may not reflect the true species tree or history of speciation. This has led to considerable debate about whether concatenation of loci is the best approach for phylogenetic analysis. The application of Next‐generation sequencing techniques such as RAD‐seq generates thousands of relatively short sequence reads from across the genomes of the sampled taxa. These data sets are typically concatenated for phylogenetic analysis leading to data sets that contain millions of base pairs per taxon. The influence of gene region conflict among so many loci in determining the phylogenetic relationships among taxa is unclear. We simulated RAD‐seq data by sampling 100 and 500 base pairs from alignments of over 6000 coding regions that each produce one of three highly supported alternative phylogenies of seven species of Drosophila. We conducted phylogenetic analyses on different sets of these regions to vary the sampling of loci with alternative gene trees to examine the effect on detecting the species tree. Irrespective of sequence length sampled per region and which subset of regions was used, phylogenetic analyses of the concatenated data always recovered the species tree. The results suggest that concatenated alignments of Next‐generation data that consist of many short sequences are robust to gene tree/species tree conflict when the goal is to determine the phylogenetic relationships among taxa.     

Evidence for phylogenetic congruence among sulfur-oxidizing chemoautotrophic bacterial endosymbionts and their bivalve hosts

Sulfur-oxidizing chemoautotrophic (thioautotrophic) bacteria are now known to occur as endosymbionts in phylogenetically diverse bivalve hosts found in a wide variety of marine environments. The evolutionary origins of these symbioses, however, have remained obscure. Comparative 16S rRNA sequence analysis was used to investigate whether thioautotrophic endosymbionts are monophyletic or polyphyletic in origin and to assess whether phylogenetic relationships inferred among these symbionts reflect those inferred among their hosts. 16S rRNA gene sequences determined for endosymbionts from nine newly examined bivalve species from three families (Vesicomyidae, Lucinidae, and Solemyidae) were compared with previously published 16S rRNA sequences of thioautotrophic symbionts and free-living bacteria. Distance and parsimony methods were used to infer phylogenetic relationships among these bacteria. All newly examined symbionts fall within the gamma subdivision of the Proteobacteria, in clusters containing previously examined symbiotic thioautotrophs. The closest free-living relatives of these symbionts are bacteria of the genus Thiomicrospira. Symbionts of the bivalve superfamily Lucinacea and the family Vesicomyidae each form distinct monophyletic lineages which are strongly supported by bootstrap analysis, demonstrating that host phylogenies inferred from morphological and fossil evidence are congruent with phylogenies inferred for their respective symbionts by molecular sequence analysis. The observed congruence between host and symbiont phylogenies indicates shared evolutionary history of hosts and symbiont lineages and suggests an ancient origin for these symbioses. Correspondence to: D.L. Distel     

Genealogies: Pedigrees and Phylogenies are Reticulating Networks Not Just Divergent Trees

Pedigrees illustrate the genealogical relationships among individuals, and phylogenies do the same for groups of organisms (such as species, genera, etc.). Here, I provide a brief survey of current concepts and methods for calculating and displaying genealogical relationships. These relationships have long been recognized to be reticulating, rather than strictly divergent, and so both pedigrees and phylogenies are correctly treated as networks rather than trees. However, currently most pedigrees are instead presented as “family trees”, and most phylogenies are presented as phylogenetic trees. Nevertheless, the historical development of concepts shows that networks pre-dated trees in most fields of biology, including the study of pedigrees, biology theory, and biology practice, as well as in historical linguistics in the social sciences. Trees were actually introduced in order to provide a simpler conceptual model for historical relationships, since trees are a specific type of simple network. Computationally, trees and networks are a part of graph theory, consisting of nodes connected by edges. In this mathematical context they differ solely in the absence or presence of reticulation nodes, respectively. There are two types of graphs that can be called phylogenetic networks: (1) rooted evolutionary networks, and (2) unrooted affinity networks. There are quite a few computational methods for unrooted networks, which have two main roles in phylogenetics: (a) they act as a generic form of multivariate data display; and (b) they are used specifically to represent haplotype networks. Evolutionary networks are more difficult to infer and analyse, as there is no mathematical algorithm for reconstructing unique historical events. There is thus currently no coherent analytical framework for computing such networks.     

Gene trees, species trees, and morphology converge on a similar phylogeny of living gars (Actinopterygii: Holostei: Lepisosteidae), an ancient clade of ray-finned fishes

Extant gars represent the remaining members of a formerly diverse assemblage of ancient ray-finned fishes and have been the subject of multiple phylogenetic analyses using morphological data. Here, we present the first hypothesis of phylogenetic relationships among living gar species based on molecular data, through the examination of gene tree heterogeneity and coalescent species tree analyses of a portion of one mitochondrial (COI) and seven nuclear (ENC1, myh6, plagl2, S7 ribosomal protein intron 1, sreb2, tbr1, and zic1) genes. Individual gene trees displayed varying degrees of resolution with regards to species-level relationships, and the gene trees inferred from COI and the S7 intron were the only two that were completely resolved. Coalescent species tree analyses of nuclear genes resulted in a well-resolved and strongly supported phylogenetic tree of living gar species, for which Bayesian posterior node support was further improved by the inclusion of the mitochondrial gene. Species-level relationships among gars inferred from our molecular data set were highly congruent with previously published morphological phylogenies, with the exception of the placement of two species, Lepisosteus osseus and L. platostomus. Re-examination of the character coding used by previous authors provided partial resolution of this topological discordance, resulting in broad concordance in the phylogenies inferred from individual genes, the coalescent species tree analysis, and morphology. The completely resolved phylogeny inferred from the molecular data set with strong Bayesian posterior support at all nodes provided insights into the potential for introgressive hybridization and patterns of allopatric speciation in the evolutionary history of living gars, as well as a solid foundation for future examinations of functional diversification and evolutionary stasis in a "living fossil" lineage.     

Phylogenetic relationships among Phytophthora species inferred from sequence analysis of mitochondrially encoded cytochrome oxidase I and II genes

The phylogenetic relationships of 51 isolates representing 27 species of Phytophthora were assessed by sequence alignment of 568 bp of the mitochondrially encoded cytochrome oxidase II gene. A total of 1299 bp of the cytochrome oxidase I gene also were examined for a subset of 13 species. The cox II gene trees constructed by a heuristic search, based on maximum parsimony for a bootstrap 50% majority-rule consensus tree, revealed 18 species grouping into seven clades and nine species unaffiliated with a specific clade. The phylogenetic relationships among species observed on cox II gene trees did not exhibit consistent similarities in groupings for morphology, pathogenicity, host range or temperature optima. The topology of cox I gene trees, constructed by a heuristic search based on maximum parsimony for a bootstrap 50% majority-rule consensus tree for 13 species of Phytophthora, revealed 10 species grouping into three clades and three species unaffiliated with a specific clade. The groupings in general agreed with what was observed in the cox II tree. Species relationships observed for the cox II gene tree were in agreement with those based on ITS regions, with several notable exceptions. Some of these differences were noted in species in which the same isolates were used for both ITS and cox II analysis, suggesting either a differential rate of evolutionary divergence for these two regions or incorrect assumptions about alignment of ITS sequences. Analysis of combined data sets of ITS and cox II sequences generated a tree that did not differ substantially from analysis of ITS data alone, however, the results of a partition homogeneity test suggest that combining data sets may not be valid.     

Phylogenetic reconstruction using secondary structures of Internal Transcribed Spacer 2 (ITS2, rDNA): finding the molecular and morphological gap in Caribbean gorgonian corals

Background  

Most phylogenetic studies using current methods have focused on primary DNA sequence information. However, RNA secondary structures are particularly useful in systematics because they include characteristics, not found in the primary sequence, that give "morphological" information. Despite the number of recent molecular studies on octocorals, there is no consensus opinion about a region that carries enough phylogenetic resolution to solve intrageneric or close species relationships. Moreover, intrageneric morphological information by itself does not always produce accurate phylogenies; intra-species comparisons can reveal greater differences than intra-generic ones. The search for new phylogenetic approaches, such as by RNA secondary structure analysis, is therefore a priority in octocoral research.     

PHYLOGENETIC RELATIONSHIPS WITHIN THE COLONIAL VOLVOCALES (CHLOROPHYTA) INFERRED FROM CLADISTIC ANALYSIS BASED ON MORPHOLOGICAL DATA1

A cladistic analysis was used to deduce the phylogenetic relationships within the colonial Volvocales. Forty-one pairs of characters related to gross morphology and ultrastructure of vegetative colonies as well as asexual and sexual reproduction were analyzed based on parsimony, using the PAUP 3.0 computer program, for 25 species belonging to nine volvocacean and goniacean genera of the colonial Volvocales. Chlamydomonas reinhardtii Dangeard was the outgroup. The strict consensus tree indicated the presence of two monophyletic groups, one composed of all the volvocacean species analyzed in this study and the other containing the goniacean species except for the four-celled species Gonium sociale (Dujardin) Warming. In addition, these two groups constitute a large monophyletic group, to which G. sociale is a sister group. A new combination Tetrabaena socialis (Dujardin) Nozaki et Itoh and a new family Tetrabaenaceae Nozaki et Itoh are thus proposed for G. sociale. In addition, the analysis suggests that the volvocacean genera Eudorina and Pleodorina are paraphyletic groups, respectively, and that the monotypic genus Yamagishiella has no autapomorphic characters and represents primitive features of the anisogamous and oogamous genera of the Volvocaceae. Phylogenetic relationships within the Volvocaceae and the Goniaceae, as well as the various modes of sexual reproduction exhibited by these organisms, are discussed on the basis of the analysis.