Non-Mendelian inheritance of paralogs of 2 cytoskeletal genes in the ciliate Chilodonella uncinata

Recognition of the role of non-Mendelian inheritance is on the rise, particularly as epigenetic phenomena are shown to shape the transformation of genomes into phenotypes. Ciliates provide a model system in which to explore the role of epigenetics because ciliates have both a germ line (micronuclear) and somatic (macronuclear) genome within every cell. In the ciliate Chilodonella uncinata, the macronucleus is extensively fragmented such that many genes end up on their own chromosomes. Hence, it is possible to track the fate of unlinked genes within macronuclei of C. uncinata. Here we demonstrate that the pattern of inheritance in isolates of C. uncinata is complex and involves both Mendelian transmission between micronuclei and macronuclei and epigenetic phenomena. The macronuclei from 2 isolates of C. uncinata and their progeny share identical rDNA loci and 2 identical beta-tubulin paralogs, yet have different actin paralogs and some beta-tubulin paralogs that are not shared. We propose a model in which all the divergent paralogs are present in the ciliate micronuclei. Under this model, different paralogs are retained in developing macronuclei following conjugation. We further speculate that an epigenetic mechanism, such as RNA interference, is involved in selective retention of specific paralogs within lines. This system allows the exploration of epigenetic phenomena that shape somatic genomes and provides parallels to studies of the development of somatic nuclei within animals.     

Variation in modes and rates of evolution in nuclear and mitochondrial ribosomal DNA in the mushroom genus Amanita (Agaricales, Basidiomycota): phylogenetic implications

Modes and rates of molecular evolution, and congruence and combinability for phylogenetic reconstruction, of portions of the nuclear large ribosomal subunit (nLSU-rDNA) and mitochondrial small subunit (mtSSU-rDNA) genes were investigated in the mushroom genus Amanita. The AT content was higher in the mtSSU-rDNA than in the nLSU-rDNA. A transition bias in which AT substitutions were as frequent as transitions was present in the mtSSU-rDNA but not in the nLSU-rDNA. Among-sites rate variation in nucleotide substitutions at variable sites was present in the nLSU-rDNA but not in the mtSSU-rDNA. Likelihood ratio tests indicated very different models of evolution for the two molecules. A molecular clock could be rejected for both data sets. Rates of molecular evolution in the two molecules were uncoupled: faster evolutionary rates in the mtSSU-rDNA and nLSU-rDNA were not observed for the same taxa. In separate phylogenetic analyses, the nLSU-rDNA data set had higher phylogenetic resolution. The partition homogeneity test and statistical bootstrap support for branches indicated absence of conflict in the phylogenetic signal in the two data sets; however, tree topologies produced from the separate data sets were not congruent. Heterogeneity in modes and rates of evolution in the two molecules pose difficulties for a combined analysis of the two data sets: the use of equally weighted parsimony is not fully satisfactory when rate heterogeneity is present, and it is impractical to determine a model for maximum-likelihood analysis that fits simultaneously two heterogeneous data sets. Overall topologies produced from either the separated or the combined analyses using various tree reconstruction methods were identical for nearly all statistically significant branches.     

核基因序列在昆虫分子系统学上的应用

核基因中含有更加丰富的生物学信息,运用核基因序列或将核基因序列与线粒体基因序列相结合研究昆虫的系统发育正成为分子系统学领域的一种发展趋势.核糖体基因中18S rDNA、28S rDNA、ITS已在昆虫分子系统学中得到了广泛的应用.与核糖体基因相比,虽然编码蛋白的核基因应用于昆虫分子系统学的种类不少,但大部分都是应用于双翅目和鳞翅目昆虫的分子系统学研究中,能够成功地普遍用于多个目昆虫的系统学研究的核基因并不多.本文简要介绍了应用于昆虫分子系统学的核中核糖体基因和编码蛋白的核基因,并分析了核基因序列在分子系统学应用上的局限性和应用前景.     

Micronuclear and macronuclear forms of beta-tubulin genes in the ciliate Chilodonella uncinata reveal insights into genome processing and protein evolution

Chilodonella uncinata, like all ciliates, contains two distinct nuclei in every cell: a germline micronucleus and a somatic macronucleus. During development of the macronucleus from a zygotic nucleus, the genome is processed in several ways, including elimination of internal sequences. In this study, we analyze micronuclear and macronuclear copies of beta-tubulin in C. uncinata and find at least four divergent paralogs of beta-tubulin in the macronucleus. We characterize the micronuclear version of one paralog and compare its internally eliminated sequences (IESs) with previously described IESs in this species. These comparisons reveal the presence of a conserved sequence motif within IESs. In addition, we compare the sequences of beta-tubulin from C. uncinata with other ciliates and to other alveolates in order to test the hypothesis that the mode of molecular evolution in ciliates obscures phylogenetic signal in protein-coding genes. We find that heterogeneous rates of substitution in beta-tubulin across ciliates result in unstable genealogies that are inconsistent with phylogenies based on small subunit rDNA genes and on ultrastructure. We discuss the implications of our findings for genome processing and protein evolution in ciliates.     

Phylogenetic relationships among Acanthamoeba spp. based on PCR-RFLP analyses of mitochondrial small subunit rRNA gene.

We investigated the value of mitochondrial small subunit rRNA gene (mt SSU rDNA) PCR-RFLP as a taxonomic tool for Acanthamoeba isolates with close inter-relationships. Twenty-five isolates representing 20 species were included in the analysis. As in nuclear 18S rDNA analysis, two type strains (A. astronyxis and A. tubiashi) of morphological group 1 diverged earliest from the other strains, but the divergence between them was less than in 18S riboprinting. Acanthamoeba griffini of morphological group 2 branched between pathogenic (A. culbertsoni A-1 and A. healyi OC-3A) and nonpathogenic (A. palestinensis Reich, A. pustulosa GE-3a, A. royreba Oak Ridge, and A lenticulata PD2S) strains of morphological group 3. Among the remaining isolates of morphological group 2, the Chang strain had the identical mitochondrial riboprints as the type strain of A. hatchetti. AA2 and AA1, the type strains of A. divionensis and A. paradivionensis, respectively, had the identical riboprints as A. quina Vil3 and A. castellanii Ma. Although the branching orders of A. castellanii Neff, A. polyphaga P23, A. triangularis SH621, and A. lugdunensis L3a were different from those in 18S riboprinting analysis, the results obtained from this study generally coincided well with those from 18S riboprinting. Mitochondrial riboprinting may have an advantage over nuclear 18S rDNA riboprinting because the mt SSU rDNAs do not seem to have introns that are found in the 18S genes of Acanthamoeba and that distort phylogenetic analyses.     

Molecular differentiation of the microgastrine species commonly found in paddy fields from Southeast Asia,with additional data on their phylogeny （Hymenoptera：Braconidae）

Partial DNA sequences of three genes, that is, mitochondrial large ribosomal subunit (16S), nuclear large ribosomal subunit (28S D2) and mitochondrial NADH1 dehydrogenase (NADH1) gene, were sequenced from different microgas trine species(Braconidae: Microgastrinae) collected fresh from paddy fields. The DNA sequences were used to determine the extent of sequence variation among species in order to evaluate the specific status of each species. Cladistic analysis was also used to infer a phylogenetic relationship among these species. The results showed that sequence divergence among species of the same genus Cotesia was much lower than those among different genera, such as Cotesia, Exoryza and Apanteles; the sequence similarity of 16S rDNA and NADH 1 genes between Cotesia sp. and C. chilonis was higher than that between C. sp. and C. ruficrus.Phylogenetic analyses suggested that four species of Cotesia were always grouped in the same clade regardless of using different analysis methods; Cotesia sp. and C. chilonis are more closely related to each other than to C. ruficrus, different from previous morphological results. Additionally, sequence analyses indicated that NADH1 gene has more parsimony informative characters than 28S rDNA D2 and 16S rDNA at the species-level analysis,indicating that NADH1 gene might be a useful marker for species-level analysis.     

Ribosomal DNA: molecular evolution and phylogenetic inference.

Ribosomal DNA (rDNA) sequences have been aligned and compared in a number of living organisms, and this approach has provided a wealth of information about phylogenetic relationships. Studies of rDNA sequences have been used to infer phylogenetic history across a very broad spectrum, from studies among the basal lineages of life to relationships among closely related species and populations. The reasons for the systematic versatility of rDNA include the numerous rates of evolution among different regions of rDNA (both among and within genes), the presence of many copies of most rDNA sequences per genome, and the pattern of concerted evolution that occurs among repeated copies. These features facilitate the analysis of rDNA by direct RNA sequencing, DNA sequencing (either by cloning or amplification), and restriction enzyme methodologies. Constraints imposed by secondary structure of rRNA and concerted evolution need to be considered in phylogenetic analyses, but these constraints do not appear to impede seriously the usefulness of rDNA. An analysis of aligned sequences of the four nuclear and two mitochondrial rRNA genes identified regions of these genes that are likely to be useful to address phylogenetic problems over a wide range of levels of divergence. In general, the small subunit nuclear sequences appear to be best for elucidating Precambrian divergences, the large subunit nuclear sequences for Paleozoic and Mesozoic divergences, and the organellar sequences of both subunits for Cenozoic divergences. Primer sequences were designed for use in amplifying the entire nuclear rDNA array in 15 sections by use of the polymerase chain reaction; these "universal" primers complement previously described primers for the mitochondrial rRNA genes. Pairs of primers can be selected in conjunction with the analysis of divergence of the rRNA genes to address systematic problems throughout the hierarchy of life.     

蛹虫草线粒体DNA与细胞核DNA进化关系的比较

【目的】分析蛹虫草是否存在核内线粒体DNA片段,比较蛹虫草线粒体DNA与细胞核DNA的碱基变异程度及所反映的菌株间的系统发育关系。【方法】通过本地BLAST或LAST对蛹虫草线粒体基因组和核基因组进行序列相似性搜索;从10个已知线粒体基因组的蛹虫草菌株中分别扩增7个细胞核蛋白编码基因片段,并与其在14个线粒体蛋白编码基因上的碱基变异情况进行比较。【结果】蛹虫草核基因组中存在5处较短的核内线粒体DNA片段,总长只有278bp。蛹虫草核DNA的变异频率整体上高于线粒体DNA。核DNA和线粒体DNA所反映的蛹虫草菌株间的系统发育关系存在显著差异。【结论】蛹虫草线粒体DNA与核DNA间不存在长片段的基因交流,二者变异频率不同,所反映的蛹虫草菌株间的系统发育关系也有差异。本研究增加了对蛹虫草线粒体与细胞核DNA进化关系的认识。     

Next generation sequencing from Hepatozoon canis (Apicomplexa: Coccidia: Adeleorina): Complete apicoplast genome and multiple mitochondrion-associated sequences

Extrachromosomal genomes of the adeleorinid parasite Hepatozoon canis infecting an Israeli dog were investigated using next-generation and standard sequencing technologies. A complete apicoplast genome and several mitochondrion-associated sequences were generated. The apicoplast genome (31,869?bp) possessed two copies of both large subunit (23S) and small subunit (16S) ribosomal RNA genes (rDNA) within an inverted repeat region, as well as 22 protein-coding sequences, 25 transfer RNA genes (tDNA) and seven open reading frames of unknown function. Although circular-mapping, the apicoplast genome was physically linear according to next-generation data. Unlike other apicoplast genomes, genes encoding ribosomal protein S19 and tDNAs for alanine, aspartic acid, histidine, threonine and valine were not identified. No complete mitochondrial genome was recovered using next-generation data or directed PCR amplifications. Eight mitochondrion-associated (215–3523?bp) contigs assembled from next-generation data encoded a complete cytochrome c oxidase subunit I coding sequence, a complete cytochrome c oxidase subunit III coding sequence, two complete cytochrome B coding sequences, a non-coding, pseudogene for cytochrome B and multiple fragmented mitochondrial rDNA genes (SSUA, SSUB, SSUD, LSUC, LSUG, RNA6, RNA10, RNA14, RNA18). The paucity of NGS reads generating each of the mitochondrion-like sequences suggested that a complete mitochondrial genome at typically high copy number was absent in H. canis. In contrast, the complete nuclear rDNA unit sequence of H. canis (18S rDNA to 28S rDNA, 6977?bp) had >1000-fold next-generation coverage. Multiple divergent (from 93.6% to 99.9% pairwise identities) nuclear 18S rDNA contigs were generated (three types with 10 subtypes total). To our knowledge this is the first apicoplast genome sequenced from any adeleorinid coccidium and the first mitochondrion-associated sequences from this serious pathogen of wild and domestic canids. These newly generated sequences may provide useful genetic loci for high-resolution species-level genotyping that is currently impossible using existing nuclear rDNA targets.     

Cryptic diversity within the choanoflagellate morphospecies complex Codosiga botrytis - Phylogeny and morphology of ancient and modern isolates

Choanoflagellates are closely related to metazoans and fungi according to recent phylogenetic studies; therefore the systematics of these organisms is of particular interest. The choanoflagellate morphospecies Codosiga botrytis is the first described choanoflagellate, and is one of the most frequently reported choanoflagellate species. In this study we present phylogenetic and morphological data on eight different strains of Codosiga botrytis. Among these there are five ancient strains; these cultures have been established from up to 43,000 years old cysts from Siberian permafrost. We found that based on the variable V4 region of the small subunit (SSU) of the rDNA, all the investigated freshwater isolates of Codosiga botrytis, together with Sphaeroeca volvox, form a cluster at the base of all other choanoflagellate species. Moreover, the morphospecies described classically as Codosiga botrytis contains at least four different genotypes separated by considerably high genetic distance. All these 'cryptic species' have identical general morphology and cell structure. Strains have a similar life cycle with several different life forms and large morphological plasticity. For the first time we were able to establish cultures from cryo-conserved cysts of choanoflagellates. The ancient strains did not differ significantly in partial SSU rDNA from the modern ones. Besides, no biogeographically pattern could be established. This fact and the low genetic distances of some strains from remote locations support the distribution of dormant stages via air.     

Wood ingestion by passalid beetles in the presence of xylose-fermenting gut yeasts

During a survey of insect gut micro-organisms, we consistently isolated Pichia stipitis-like yeasts (Fungi: Ascomycota, Saccharomycetes) from the wood-ingesting beetles, Odontotaenius disjunctus and Verres sternbergianus (Coleoptera: Passalidae). The yeasts were isolated from passalid beetles over a wide area, including the eastern and midwestern USA and Panama. Phylogenetic analyses of the nuclear encoded small and large subunit rRNA gene (rDNA) sequences distinguished a well-supported clade consisting of the passalid yeasts and Pichia stipitis, P. segobiensis, Candida shehatae and C. ergatensis. Members of this clade have the ability to ferment and assimilate xylose or to hydrolyse xylan, major components of the polysaccharide, hemicellulose. Sexual reproduction was present in the passalid isolates but was rare among the gut yeasts of other beetles to which they were compared. Minor genetic and phenotypic variation among some of the passalid yeasts was detected using markers from the internal transcribed spacer region of the rDNA repeat unit, morphology, and in vitro metabolic tests. The consistent association of xylose-fermenting yeasts of almost identical genotypes with passalid beetles across a broad geographical distribution, suggests a significant symbiotic association.     

How do insect nuclear and mitochondrial gene substitution patterns differ? Insights from Bayesian analyses of combined datasets

We analyzed 12 combined mitochondrial and nuclear gene datasets in seven orders of insects using both equal weights parsimony (to evaluate phylogenetic utility) and Bayesian methods (to investigate substitution patterns). For the Bayesian analyses we used relatively complex models (e.g., general time reversible models with rate variation) that allowed us to quantitatively compare relative rates among genes and codon positions, patterns of rate variation among genes, and substitution patterns within genes. Our analyses indicate that nuclear and mitochondrial genes differ in a number of important ways, some of which are correlated with phylogenetic utility. First and most obviously, nuclear genes generally evolve more slowly than mitochondrial genes (except in one case), making them better markers for deep divergences. Second, nuclear genes showed universally high values of CI and (generally) contribute more to overall tree resolution than mitochondrial genes (as measured by partitioned Bremer support). Third, nuclear genes show more homogeneous patterns of among-site rate variation (higher values of alpha than mitochondrial genes). Finally, nuclear genes show more symmetrical transformation rate matrices than mitochondrial genes. The combination of low values of alpha and highly asymmetrical transformation rate matrices may explain the overall poor performance of mitochondrial genes when compared to nuclear genes in the same analysis. Our analyses indicate that some parameters are highly correlated. For example, A/T bias was positively and significantly associated with relative rate and CI was positively and significantly associated with alpha (the shape of the gamma distribution). These results provide important insights into the substitution patterns that might characterized high quality genes for phylogenetic analysis: high values of alpha, unbiased base composition, and symmetrical transformation rate matrices. We argue that insect molecular systematists should increasingly focus on nuclear rather than mitochondrial gene datasets because nuclear genes do not suffer from the same substitutional biases that characterize mitochondrial genes.     

Genetic mosaic in a marine species flock

We used molecular approaches to study the status of speciation in coral reef fishes known as hamlets (Serranidae: Hypoplectrus). Several hamlet morphospecies coexist on Caribbean reefs, and mate assortatively with respect to their strikingly distinct colour patterns. We provide evidence that, genetically, the hamlets display characteristics common in species flocks on land and in freshwaters. Substitutions within two mitochondrial DNA (mtDNA) protein-coding genes place hamlets within a monophyletic group relative to members of two related genera (Serranus and Diplectrum), and establish that the hamlet radiation must have been very recent. mtDNA distances separating hamlet morphospecies were slight (0.6 +/- 0.04%), yielding a coalescent estimate for the age of the hamlet flock of approximately 430 000 years. Morphospecies did not sort into distinct mtDNA haplotype phylogroups, and alleles at five hypervariable microsatellite loci were shared broadly across species boundaries. None the less, molecular variation was not distributed at random. Analyses of mtDNA haplotype frequencies and nested clades in haplotype networks revealed significant genetic differences between geographical regions and among colour morphospecies. We also observed significant microsatellite differentiation between geographical regions and in Puerto Rico, among colour morphospecies; the latter providing evidence for reproductive isolation between colour morphospecies at this locale. In our Panama collection, however, colour morphospecies were mostly genetically indistinguishable. This mosaic pattern of DNA differentiation implies a complex interaction between population history, mating behaviour and geography and suggests that porous boundaries separate species in this flock of brilliantly coloured coral reef fishes.     

Variation of nuclear and mitochondrial DNAs in Korean and Chinese isolates of Clonorchis sinensis

We compared the DNA sequence difference of isolates of Clonorchis sinensis from one Korean (Kimhae) and two Chinese areas (Guangxi and Shenyang). The sequences of nuclear rDNA (18S, internal transcribed spacer 1 and 2: ITS1 and ITS2) and mitochondrial DNA (cytochrome c oxidase subunit 1: cox1) were compared. A very few intraspecific nucleotide substitution of the 18S, ITS1, ITS2 and cox1 was found among three isolates of C. sinensis and a few nucleotide insertion and deletion of ITS1 were detected. The 18S, ITS1, ITS2 and cox1 sequences were highly conserved among three isolates. These findings indicated that the Korean and two Chinese isolates are similar at the DNA sequence level.     

Mitochondrial genome sequence evolution in Chlamydomonas

The mitochondrial genomes of the Chlorophyta exhibit significant diversity with respect to gene content and genome compactness; however, quantitative data on the rates of nucleotide substitution in mitochondrial DNA, which might help explain the origin of this diversity, are lacking. To gain insight into the evolutionary forces responsible for mitochondrial genome diversification, we sequenced to near completion the mitochondrial genome of the chlorophyte Chlamydomonas incerta, estimated the evolutionary divergence between Chlamydomonas reinhardtii and C. incerta mitochondrial protein-coding genes and rRNA-coding regions, and compared the relative evolutionary rates in mitochondrial and nuclear genes. Synonymous and nonsynonymous substitution rates do not differ significantly between the mitochondrial and nuclear protein-coding genes. The mitochondrial rRNA-coding regions, however, are evolving much faster than their nuclear counterparts, and this difference might be explained by relaxed functional constraints on the mitochondrial translational apparatus due to the small number of proteins synthesized in Chlamydomonas mitochondria. Substitution rates at synonymous sites in a nonstandard mitochondrial gene (rtl) and at intronic and synonymous sites in nuclear genes expressed at low levels suggest that the mutation rate is similar in these two genetic compartments. Potential evolutionary forces shaping mitochondrial genome evolution in Chlamydomonas are discussed.     

Investigating phylogenetic relationships of sunfishes and black basses (Actinopterygii: Centrarchidae) using DNA sequences from mitochondrial and nuclear genes

The 32 species of the Centrarchidae are ecologically important components of the diverse fish communities that characterize North American freshwater ecosystems. In spite of a rich history of systematic investigations of centrarchid fishes there is extensive conflict among previous hypotheses that may be due to restricted taxon or character sampling. We present the first phylogenetic analysis of the Centrarchidae that combines DNA sequence data from both the mitochondrial and nuclear genomes and includes all described species. Gene sequence data were collected from a complete mtDNA protein coding gene (NADH subunit 2), a nuclear DNA intron (S7 ribosomal protein intron 1), and a portion of a nuclear DNA protein-coding region (Tmo-4C4). Phylogenetic trees generated from analysis of the three-gene dataset were used to test alternative hypotheses of centrarchid relationships that were gathered from the literature. Four major centrarchid lineages are present in trees generated in maximum parsimony (MP) and Bayesian maximum likelihood analyses (BML). These lineages are Acantharchus pomotis, Micropterus, Lepomis, and a clade containing Ambloplites, Archoplites, Centrarchus, Enneacanthus, and Pomoxis. Phylogenetic trees resulting from MP and BML analyses are highly consistent but differ with regard to the placement of A. pomotis. Significant phylogenetic incongruence between mtDNA and nuclear genes appears to result from different placement of Micropterus treculi, and is not characteristic of relationships in all other parts of the centrarchid phylogeny. Slightly more than half of the 27 previously proposed hypotheses of centrarchid relationships were rejected based on the Shomodaira-Hasegawa test.     

Molecular genetic variation within and among isolates of QPX (Thraustochytridae), a parasite of the hard clam Mercenaria mercenaria

The thraustochytrid known as QPX (Quahog Parasite Unknown) has sporadically caused disease in the hard clam Mercenaria mercenaria along the east coast of North America since the 1960s. We hypothesized that genetically distinct QPX strains might be responsible for outbreaks of QPX disease in different areas and tested this hypothesis by comparing several QPX isolates recovered from the recent outbreak in Raritan Bay, New York with QPX strains isolated from 2 outbreaks in Massachusetts, USA. There was no variation in small subunit rDNA (SSU rDNA), 5.8S rDNA, or 4 mitochondrial gene sequences. In contrast, both of the ribosomal ribonucleic acid (rRNA) operon intergenic spacers, internal transcribed spacers 1 and 2 (ITS1 and ITS2), revealed substantial sequence variation. However, strain-specific sequences were not detected because the ITS sequence variation within QPX isolates was comparable to the variation between isolates. ITS1 sequences recovered from an infected clam by amplification with a QPX ITS2-specific primer were identical to those recovered from the QPX isolates.