The very early evolution of biological complexity.

All extant life is descended from a common ancestor, which, despite being very ancient, appears to have been a complex cellular organism. A new study by Moody et al. shows that this ancestor was not only a complex cell, but also lived within a microbial ecology likely inhabited by other complex cells.

Electron transport chains as a window into the earliest stages of evolution.

The origin and early evolution of life is generally studied under two different paradigms: bottom up and top down. Prebiotic chemistry and early Earth geochemistry allow researchers to explore possible origin of life scenarios. But for these "bottom-up" approaches, even successful experiments only amount to a proof of principle. On the other hand, "top-down" research on early evolutionary history is able to provide a historical account about ancient organisms, but is unable to investigate stages that occurred during and just after the origin of life. Here, we consider ancient electron transport chains (ETCs) as a potential bridge between early evolutionary history and a protocellular stage that preceded it. Current phylogenetic evidence suggests that ancestors of several extant ETC components were present at least as late as the last universal common ancestor of life. In addition, recent experiments have shown that some aspects of modern ETCs can be replicated by minerals, protocells, or organic cofactors in the absence of biological proteins. Here, we discuss the diversity of ETCs and other forms of chemiosmotic energy conservation, describe current work on the early evolution of membrane bioenergetics, and advocate for several lines of research to enhance this understanding by pairing top-down and bottom-up approaches.

Encyclopedia of Family A DNA Polymerases Localized in Organelles: Evolutionary Contribution of Bacteria Including the Proto-Mitochondrion.

DNA polymerases synthesize DNA from deoxyribonucleotides in a semiconservative manner and serve as the core of DNA replication and repair machinery. In eukaryotic cells, there are 2 genome-containing organelles, mitochondria, and plastids, which were derived from an alphaproteobacterium and a cyanobacterium, respectively. Except for rare cases of genome-lacking mitochondria and plastids, both organelles must be served by nucleus-encoded DNA polymerases that localize and work in them to maintain their genomes. The evolution of organellar DNA polymerases has yet to be fully understood because of 2 unsettled issues. First, the diversity of organellar DNA polymerases has not been elucidated in the full spectrum of eukaryotes. Second, it is unclear when the DNA polymerases that were used originally in the endosymbiotic bacteria giving rise to mitochondria and plastids were discarded, as the organellar DNA polymerases known to date show no phylogenetic affinity to those of the extant alphaproteobacteria or cyanobacteria. In this study, we identified from diverse eukaryotes 134 family A DNA polymerase sequences, which were classified into 10 novel types, and explored their evolutionary origins. The subcellular localizations of selected DNA polymerases were further examined experimentally. The results presented here suggest that the diversity of organellar DNA polymerases has been shaped by multiple transfers of the PolI gene from phylogenetically broad bacteria, and their occurrence in eukaryotes was additionally impacted by secondary plastid endosymbioses. Finally, we propose that the last eukaryotic common ancestor may have possessed 2 mitochondrial DNA polymerases, POP, and a candidate of the direct descendant of the proto-mitochondrial DNA polymerase I, rdxPolA, identified in this study.

The evolutionary history of the HUP domain.

Among the enzyme lineages that undoubtedly emerged prior to the last universal common ancestor is the so-called HUP, which includes Class I aminoacyl tRNA synthetases (AARSs) as well as enzymes mediating NAD, FAD, and CoA biosynthesis. Here, we provide a detailed analysis of HUP evolution, from emergence to structural and functional diversification. The HUP is a nucleotide binding domain that uniquely catalyzes adenylation via the release of pyrophosphate. In contrast to other ancient nucleotide binding domains with the αßα sandwich architecture, such as P-loop NTPases, the HUP's most conserved feature is not phosphate binding, but rather ribose binding by backbone interactions to the tips of ß1 and/or ß4. Indeed, the HUP exhibits unusual evolutionary plasticity and, while ribose binding is conserved, the location and mode of binding to the base and phosphate moieties of the nucleotide, and to the substrate(s) reacting with it, have diverged with time, foremost along the emergence of the AARSs. The HUP also beautifully demonstrates how a well-packed scaffold combined with evolvable surface elements promotes evolutionary innovation. Finally, we offer a scenario for the emergence of the HUP from a seed ßαß fragment, and suggest that despite an identical architecture, the HUP and the Rossmann represent independent emergences.

Evolution of factors shaping the endoplasmic reticulum.

Endomembrane system compartments are significant elements in virtually all eukaryotic cells, supporting functions including protein synthesis, post-translational modifications and protein/lipid targeting. In terms of membrane area the endoplasmic reticulum (ER) is the largest intracellular organelle, but the origins of proteins defining the organelle and the nature of lineage-specific modifications remain poorly studied. To understand the evolution of factors mediating ER morphology and function we report a comparative genomics analysis of experimentally characterized ER-associated proteins involved in maintaining ER structure. We find that reticulons, REEPs, atlastins, Ufe1p, Use1p, Dsl1p, TBC1D20, Yip3p and VAPs are highly conserved, suggesting an origin at least as early as the last eukaryotic common ancestor (LECA), although many of these proteins possess additional non-ER functions in modern eukaryotes. Secondary losses are common in individual species and in certain lineages, for example lunapark is missing from the Stramenopiles and the Alveolata. Lineage-specific innovations include protrudin, Caspr1, Arl6IP1, p180, NogoR, kinectin and CLIMP-63, which are restricted to the Opisthokonta. Hence, much of the machinery required to build and maintain the ER predates the LECA, but alternative strategies for the maintenance and elaboration of ER shape and function are present in modern eukaryotes. Moreover, experimental investigations for ER maintenance factors in diverse eukaryotes are expected to uncover novel mechanisms.

Evolution of Cellular Organization Along the First Branches of the Tree of Life.

Current evidence suggests that some form of cellular organization arose well before the time of the last universal common ancestor (LUCA). Standard phylogenetic analyses have shown that several protein families associated with membrane translocation, membrane transport, and membrane bioenergetics were very likely present in the proteome of the LUCA. Despite these cellular systems emerging prior to the LUCA, extant archaea, bacteria, and eukaryotes have significant differences in cellular infrastructure and the molecular functions that support it, leading some researchers to argue that true cellularity did not evolve until after the LUCA. Here, we use recently reconstructed minimal proteomes of the LUCA as well as the last archaeal common ancestor (LACA) and the last bacterial common ancestor (LBCA) to characterize the evolution of cellular systems along the first branches of the tree of life. We find that a broad set of functions associated with cellular organization were already present by the time of the LUCA. The functional repertoires of the LACA and LBCA related to cellular organization nearly doubled along each branch following the divergence of the LUCA. These evolutionary trends created the foundation for similarities and differences in cellular organization between the taxonomic domains that are still observed today.

Leucoma salicis nucleopolyhedrovirus (LesaNPV) genome sequence shed new light on the origin of the Alphabaculovirus orpseudotsugatae species.

LesaNPV (Leucoma salicis nucleopolyhedrovirus) is an alphabaculovirus group Ib. Potentially, it can be an eco-friendly agent to control the white satin moth Leucoma salicis population. In this study, we have established the relationship between LesaNPV and other closely related alphabaculoviruses. Environmental samples of late instar of white satin moth collected in Poland infected with baculovirus have been homogenized, polyhedra were purified and subjected to scanning and transmission electron microscopy. Viral DNA was sequenced using the Illumina platform and the whole-genome sequence was established by de novo assembly of paired reads. Genome annotation and phylogenetic analyses were performed with the use of bioinformatics tools. The genome of LesaNPV is 132 549 bp long with 154 ORFs and 54.9% GC content. Whole-genome sequencing revealed deletion of dUTPase as well as ribonucleoside reductases small and large subunits region in LesaNPV genome compared to Dasychira pudibunda nucleopolyhedrovirus (DapuNPV) and Orgyia pseudotsugata multiple nucleopolyhedrovirus (OpMNPV) where this region is complete. Phylogenetic analysis of Baculoviridae family members showed that LesaNPV is less divergent from a common ancestor than closely related species DapuNPV and OpMNPV. This is interesting because their hosts do not occur in the same area. The baculoviruses described in this manuscript are probably isolates of one species and could be assigned to recently denominated species Alphabaculovirus orpseudotsugatae, historically originating from OpMNPV. This finding could have significant implications for the classification and understanding of the phylogeographical spread of baculoviruses.

Structure, Evolution, and Mitochondrial Genome Analysis of Mussel Species (Bivalvia, Mytilidae).

Based on the nucleotide sequences of the mitochondrial genome (mitogenome) of specimens taken from two mussel species (Arcuatula senhousia and Mytilus coruscus), an investigation was performed by means of the complex approaches of the genomics, molecular phylogenetics, and evolutionary genetics. The mitogenome structure of studied mussels, like in many other invertebrates, appears to be much more variable than in vertebrates and includes changing gene order, duplications, and deletions, which were most frequent for tRNA genes; the mussel species' mitogenomes also have variable sizes. The results demonstrate some of the very important properties of protein polypeptides, such as hydrophobicity and its determination by the purine and pyrimidine nucleotide ratio. This fact might indirectly indicate the necessity of purifying natural selection for the support of polypeptide functionality. However, in accordance with the widely accepted and logical concept of natural cutoff selection for organisms living in nature, which explains its action against deleterious nucleotide substitutions in the nonsynonymous codons (mutations) and its holding of the active (effective) macromolecules of the polypeptides in a population, we were unable to get unambiguous evidence in favor of this concept in the current paper. Here, the phylogeny and systematics of mussel species from one of the largest taxons of bivalve mollusks are studied, the family known as Mytilidae. The phylogeny for Mytilidae (order Mytilida), which currently has no consensus in terms of systematics, is reconstructed using a data matrix of 26-27 mitogenomes. Initially, a set of 100 sequences from GenBank were downloaded and checked for their gender: whether they were female (F) or male (M) in origin. Our analysis of the new data confirms the known drastic differences between the F/M mitogenome lines in mussels. Phylogenetic reconstructions of the F-lines were performed using the combined set of genetic markers, reconstructing only protein-coding genes (PCGs), only rRNA + tRNA genes, and all genes. Additionally, the analysis includes the usage of nucleotide sequences composed of other data matrices, such as 20-68 mitogenome sequences. The time of divergence from MRCA, estimated via BEAST2, for Mytilidae is close to 293 Mya, suggesting that they originate in the Silurian Period. From all these data, a consensus for the phylogeny of the subfamily of Mytilinae and its systematics is suggested. In particular, the long-debated argument on mussel systematics was resolved as to whether Mytilidae, and the subfamily of Mytilinae, are monophyletic. The topology signal, which was strongly resolved in this paper and in the literature, has refuted the theory regarding the monophyly of Mytilinae.

Neolithic expansion and the 17q21.31 inversion in Iberia: an evolutionary approach to H2 haplotype distribution in the Near East and Europe.

The chromosomal region 17q21.31 harbors a 900 kb inversion polymorphism named after the microtubule-associated protein tau (MAPT) gene. Since no recombination occurs, two haplotypes are recognized: a directly oriented variant (H1) and an inverted variant (H2). The H2 haplotype features a distribution pattern with high frequencies in the Near East and Europe, medium levels in South Asia and North Africa, and low levels elsewhere. Studies of this genomic region are relevant owing to its likely association with numerous neurodegenerative diseases. However, the causes underlying the geographic distribution of the haplotype frequencies remain a bone of contention among researchers. With this work, we have intended to outline a plausible hypothesis on the origin of the high European H2 frequencies. To that end, we have analyzed an extensive population database (including three new Iberian populations) to explore potential clinal variations of H2 frequencies. We found a sigmoidal frequency cline with an upward trend from South Asia to Europe. The maximum value was detected in the Basques from Gipuzkoa province (0.494) with the curve's inflection point in the Near East. From our results, we suggest that the most likely scenario for high H2 frequencies in Europe would be a founding event in the Near East during the late Paleolithic or early Neolithic. Subsequently, such H2 overrepresentation would have reached Europe with the arrival of the first Neolithic farmers. The current frequencies and geographic distribution of the 17q21.31 inversion suggest that the founding events mainly affected the H2D subhaplotype.

Evolution of vertebral numbers in primates, with a focus on hominoids and the last common ancestor of hominins and panins.

The primate vertebral column has been extensively studied, with a particular focus on hominoid primates and the last common ancestor of humans and chimpanzees. The number of vertebrae in hominoids-up to and including the last common ancestor of humans and chimpanzees-is subject to considerable debate. However, few formal ancestral state reconstructions exist, and none include a broad sample of primates or account for the correlated evolution of the vertebral column. Here, we conduct an ancestral state reconstruction using a model of evolution that accounts for both homeotic (changes of one type of vertebra to another) and meristic (addition or loss of a vertebra) changes. Our results suggest that ancestral primates were characterized by 29 precaudal vertebrae, with the most common formula being seven cervical, 13 thoracic, six lumbar, and three sacral vertebrae. Extant hominoids evolved tail loss and a reduced lumbar column via sacralization (homeotic transition at the last lumbar vertebra). Our results also indicate that the ancestral hylobatid had seven cervical, 13 thoracic, five lumbar, and four sacral vertebrae, and the ancestral hominid had seven cervical, 13 thoracic, four lumbar, and five sacral vertebrae. The last common ancestor of humans and chimpanzees likely either retained this ancestral hominid formula or was characterized by an additional sacral vertebra, possibly acquired through a homeotic shift at the sacrococcygeal border. Our results support the 'short-back' model of hominin vertebral evolution, which postulates that hominins evolved from an ancestor with an African ape-like numerical composition of the vertebral column.

Investigation of tRNA-based relatedness within the Planctomycetes-Verrucomicrobia-Chlamydiae (PVC) superphylum: a comparative analysis.

The PVC superphylum is a diverse group of prokaryotes that require stringent growth conditions. RNA is a fascinating molecule to find evolutionary relatedness according to the RNA World Hypothesis. We conducted tRNA gene analysis to find evolutionary relationships in the PVC phyla. The analysis of genomic data (P = 9, V = 4, C = 8) revealed that the number of tRNA genes varied from 28 to 90 in Planctomycetes and Chlamydia, respectively. Verrucomicrobia has whole genomes and the longest scaffold (3 + 1), with tRNA genes ranging from 49 to 53 in whole genomes and 4 in the longest scaffold. Most tRNAs in the E. coli genome clustered with homologs, but approximately 43% clustered with tRNAs encoding different amino acids. Planctomyces, Akkermansia, Isosphaera, and Chlamydia were similar to E. coli tRNAs. In a phylum, tRNAs coding for different amino acids clustered at a range of 8 to 10%. Further analysis of these tRNAs showed sequence similarity with Cyanobacteria, Proteobacteria, Viridiplantae, Ascomycota and Basidiomycota (Eukaryota). This indicates the possibility of horizontal gene transfer or, otherwise, a different origin of tRNA in PVC bacteria. Hence, this work proves its importance for determining evolutionary relatedness and potentially identifying bacteria using tRNA. Thus, the analysis of these tRNAs indicates that primitive RNA may have served as the genetic material of LUCA before being replaced by DNA. A quantitative analysis is required to test these possibilities that relate the evolutionary significance of tRNA to the origin of life.

Flagellar export apparatus and ATP synthetase: Homology evidenced by synteny predating the Last Universal Common Ancestor.

We report evidence further supporting homology between proteins in the F1 FO -ATP synthetase and the bacterial flagellar motor (BFM). BFM proteins FliH, FliI, and FliJ have been hypothesized to be homologous to FO -b + F1 -δ, F1 -α/ß, and F1 -γ, with similar structure and interactions. We conduct a further test by constructing a gene order dataset, examining the order of fliH, fliI, and fliJ genes across the phylogenetic breadth of flagellar and nonflagellar type 3 secretion systems, and comparing this to published surveys of gene order in the F1 FO -ATP synthetase, its N-ATPase relatives, and the bacterial/archaeal V- and A-type ATPases. Strikingly, the fliHIJ gene order was deeply conserved, with the few exceptions appearing derived, and exactly matching the widely conserved F-ATPase gene order atpFHAG, coding for subunits b-δ-α-γ. The V/A-type ATPases have a similar conserved gene order. Our results confirm homology between these systems, and suggest a rare case of synteny conserved over billions of years, predating the Last Universal Common Ancestor (LUCA).

A compendium of covariances and correlation coefficients of coalescent tree properties.

Gene genealogies are frequently studied by measuring properties such as their height (H), length (L), sum of external branches (E), sum of internal branches (I), and mean of their two basal branches (B), and the coalescence times that contribute to the other genealogical features (T). These tree properties and their relationships can provide insight into the effects of population-genetic processes on genealogies and genetic sequences. Here, under the coalescent model, we study the 15 correlations among pairs of features of genealogical trees: Hn, Ln, En, In, Bn, and Tk for a sample of size n, with 2≤k≤n. We report high correlations among Hn, Ln, In, and Bn, with all pairwise correlations of these quantities having values greater than or equal to 6[6ζ(3)+6-π2]/(π18+9π2-π4)≈0.84930 in the limit as n→∞, where ζ is the Riemann zeta function. Although En has expectation 2 for all n and Hn has expectation 2 in the n→∞ limit, their limiting correlation is 0. The results contribute toward understanding features of the shapes of coalescent trees.

The bounded coalescent model: Conditioning a genealogy on a minimum root date.

The coalescent model represents how individuals sampled from a population may have originated from a last common ancestor. The bounded coalescent model is obtained by conditioning the coalescent model such that the last common ancestor must have existed after a certain date. This conditioned model arises in a variety of applications, such as speciation, horizontal gene transfer or transmission analysis, and yet the bounded coalescent model has not been previously analysed in detail. Here we describe a new algorithm to simulate from this model directly, without resorting to rejection sampling. We show that this direct simulation algorithm is more computationally efficient than the rejection sampling approach. We also show how to calculate the probability of the last common ancestor occurring after a given date, which is required to compute the probability density of realisations under the bounded coalescent model. Our results are applicable in both the isochronous (when all samples have the same date) and heterochronous (where samples can have different dates) settings. We explore the effect of setting a bound on the date of the last common ancestor, and show that it affects a number of properties of the resulting phylogenies. All our methods are implemented in a new R package called BoundedCoalescent which is freely available online.

The Sexual Ancestor of all Eukaryotes: A Defense of the "Meiosis Toolkit": A Rigorous Survey Supports the Obligate Link between Meiosis Machinery and Sexual Recombination.

The distribution pattern of the meiotic machinery in known eukaryotes is most parsimoniously explained by the hypothesis that all eukaryotes are ancestrally sexual. However, this assumption is questioned by preliminary results, in culture conditions. These suggested that Acanthamoeba, an organism considered to be largely asexual, constitutively expresses meiosis genes nevertheless-at least in the lab. This apparent disconnect between the "meiosis toolkit" and sexual processes in Acanthamoeba led to the conclusion that the eukaryotic ancestor is asexual. In this review, the "meiosis toolkit" is rigorously defended, drawing on numerous research articles. Additionally, the claim of constitutive meiotic gene expression is probed in Acanthamoeba via the same transcriptomics data. The results show that the expression of the meiotic machinery is not constitutive in Acanthamoeba as claimed before. Furthermore, it is argued that this would have no implications for understanding the nature of the eukaryotic ancestor, regardless of the result.

Reconstructing the Last Common Ancestor: Epistemological and Empirical Challenges.

Reconstructing the genetic traits of the Last Common Ancestor (LCA) and the Tree of Life (TOL) are two examples of the reaches of contemporary molecular phylogenetics. Nevertheless, the whole enterprise has led to paradoxical results. The presence of Lateral Gene Transfer poses epistemic and empirical challenges to meet these goals; the discussion around this subject has been enriched by arguments from philosophers and historians of science. At the same time, a few but influential research groups have aimed to reconstruct the LCA with rich-in-detail hypotheses and high-resolution gene catalogs and metabolic traits. We argue that LGT poses insurmountable challenges for detailed and rich in details reconstructions and propose, instead, a middle-ground position with the reconstruction of a slim LCA based on traits under strong pressures of Negative Natural Selection, and for the need of consilience with evidence from organismal biology and geochemistry. We defend a cautionary perspective that goes beyond the statistical analysis of gene similarities and assumes the broader consequences of evolving empirical data and epistemic pluralism in the reconstruction of early life.

Homozygosity Haplotype and Whole-Exome Sequencing Analysis to Identify Potentially Functional Rare Variants Involved in Multiple Sclerosis among Sardinian Families.

Multiple Sclerosis (MS) is a complex multifactorial autoimmune disease, whose sex- and age-adjusted prevalence in Sardinia (Italy) is among the highest worldwide. To date, 233 loci were associated with MS and almost 20% of risk heritability is attributable to common genetic variants, but many low-frequency and rare variants remain to be discovered. Here, we aimed to contribute to the understanding of the genetic basis of MS by investigating potentially functional rare variants. To this end, we analyzed thirteen multiplex Sardinian families with Immunochip genotyping data. For five families, Whole Exome Sequencing (WES) data were also available. Firstly, we performed a non-parametric Homozygosity Haplotype analysis for identifying the Region from Common Ancestor (RCA). Then, on these potential disease-linked RCA, we searched for the presence of rare variants shared by the affected individuals by analyzing WES data. We found: (i) a variant (43181034 T > G) in the splicing region on exon 27 of CUL9; (ii) a variant (50245517 A > C) in the splicing region on exon 16 of ATP9A; (iii) a non-synonymous variant (43223539 A > C), on exon 9 of TTBK1; (iv) a non-synonymous variant (42976917 A > C) on exon 9 of PPP2R5D; and v) a variant (109859349-109859354) in 3'UTR of MYO16.

Evolutionary analysis of rabies virus using the partial Nucleoprotein and Glycoprotein gene in Mumbai region of India.

Nearly 1.7 million cases of dog bites are reported every year in India and many cases of animal rabies are left unattended and undiagnosed. Therefore, a mere diagnosis of rabies is not sufficient to understand the epidemiology and the spread of the rabies virus (RV) in animals. There is a paucity of information about the evolutionary dynamics of RV in dogs and its biodiversity patterns in India. In total, 50 dog-brain samples suspected of rabies were screened by the nucleoprotein- (N) and glycoprotein- (G) gene PCR. The N and G genes were subsequently sequenced to understand the molecular evolution in these genes. The phylogenetic analysis of the N gene revealed that six isolates in the Mumbai region belonged to a single Arctic lineage. Time-scaled phylogeny by Bayesian coalescent analysis of the partial N gene revealed that the time to the most recent common ancestor (TMRCA) for the sequences belonged to the cluster from 2006.68 with a highest posterior density of 95 % betweeen 2005-2008, which is assigned to Indian lineage I. Migration pattern revealed a strong Bayes factor between Mumbai to Delhi, Panji to Hyderabad, Delhi to Chennai, and Chennai to Chandigarh. Phylogenetic analysis of the G gene revealed that the RVs circulating in the Mumbai region are divided into three lineages. Time-scaled phylogeny by the Bayesian coalescent analysis method estimated that the TMRCA for sequences under study was from 1993 and Indian clusters was from 1962. In conclusion, the phylogenetic analysis of the N gene revealed that six isolates belonged to single Arctic lineages along with other Indian isolates and they were clustered into a single lineage but divided into three clades based on the G-gene sequences. The present study highlights and enhances the current molecular epidemiology and evolution of RV and revealed strong location bias and geographical clustering within Indian isolates on the basis of N and G genes.

Protein-free ribosomal RNA scaffolds can assemble poly-lysine oligos from charged tRNA fragments.

Ribosomal protein synthesis is a central process of the modern biological world. Because the ribosome contains proteins itself, it is very important to understand its precursor and evolution. Small ribozymes have demonstrated the principle of "RNA world" hypothesis, but protein free peptide ligase remains elusive. In this report, we have identified two fragments in the peptidyl transfer center that can synthesize a 9-mer poly-lysine in a solution contains Mg2+. This result is deduced from isotope-shifting in high resolution MS. To our best knowledge, this is the longest peptide oligo that can be synthesized by a pure ribozyme. Via single molecule FRET experiments, we have demonstrated the ligase mechanism was probably by substrate proximity via dimerization. We prospect that these RNA fragments can be useful to synthesize template free natural and non-natural peptides, to be model system for peptidyl transfer reaction mechanism and can shed light to the evolution of ribosome.

Ancestral Eukaryotes Reproduced Asexually, Facilitated by Polyploidy: A Hypothesis.

The notion that eukaryotes are ancestrally sexual has been gaining attention. This idea comes in part from the discovery of sets of "meiosis-specific genes" in the genomes of protists. The existence of these genes has persuaded many that these organisms may be engaging in sex, even though this has gone undetected. The involvement of sex in protists is supported by the view that asexual reproduction results in the accumulation of mutations that would inevitably result in the decline and extinction of such lineages. It is argued that this phenomenon can be obviated by polyploidy and that the "meiosis-specific genes" are used in other processes, including polyploidy control and homologous recombination, independent of meiosis. These phenomena account for the finding that these genes are expressed in cultures devoid of apparent cell fusion events. Hence, it is also proposed that asexual, and not sexual, reproduction is the ancestral condition.