CR1 retroposons provide a new insight into the phylogeny of Phasianidae species (Aves: Galliformes)

Chicken repeat 1 (CR1) elements, a class of retroposons belonging to non-long-terminal repeats, have been recognized as powerful tools for phylogenetic studies. Here we examine the phylogenetic relationships of 11 Phasianidae species based on CR1 retroposons. Together with 19 loci reported previously, a total of 99 CR1 loci were identified from chicken genome and turkey BAC clone sequences. 75 insertion events were used to address the branching order of 11 species in Phasianidae. The topology of our tree suggests that: 1) Gallus gallus possessed a basal phylogenetic position within Phasianidae and was related to Bambusicola thoracica (BSP=100%); 2) After the split of G. gallus and B. thoracica, Arborophila rufipectus diverged from Phasianidae (BSP=100%). Nine unambiguous insertion events supported a phylogenetic position of A. rufipectus different to previous mitochondrial data suggesting a hybrid origin or an ancient introgression of A. rufipectus; and 3) 22 CR1 insertion events strongly supported the eight phasianids under investigation sharing a common ancestor. Our study has revisited the phylogenetic position of G. gallus and A. rufipectus and provided a new insight into the phylogeny of Phasianidae birds. It showed that a CR1-based methodology has a great potential to be informative within Phasianidae in resolving relationships of closely related species whose radiation and speciation have occurred very recently.     

山鹧鸪属鸟类线粒体基因组的比较及系统发育研究

海南山鹧鸪(Arborophila ardens)对生境选择比较严格,种群数量稀少,属于濒危物种。为进一步研究山鹧鸪属的进化和系统发育关系,文章利用Illumina Hiseq2000高通量测序技术获得了海南山鹧鸪线粒体全基因组序列,从比较基因组学角度分析了4种山鹧鸪鸟类的线粒体基因组特征,并探讨了山鹧鸪属鸟类的系统发育地位。研究结果表明：(1) 海南山鹧鸪线粒体基因组长度为16 730 bp,编码13个蛋白质编码基因、2个核糖体RNA基因、22个转运RNA基因以及1个控制区;(2) 山鹧鸪属物种受到了纯化选择的作用,且在进化过程中积累了更多的非同义替换;(3) 山鹧鸪属位于雉科鸟类系统树的基部位置,其中白眉山鹧鸪与红喉山鹧鸪互为姐妹群,海南山鹧鸪位于山鹧鸪属的基部位置,与其他3种山鹧鸪鸟类的亲缘关系较远。     

The first complete mitochondrial genome sequence of Nanorana parkeri and Nanorana ventripunctata (Amphibia: Anura: Dicroglossidae), with related phylogenetic analyses

Members of the Nanorana genus (family Dicroglossidae) are often referred to as excellent model species with which to study amphibian adaptations to extreme environments and also as excellent keystone taxa for providing insights into the evolution of the Dicroglossidae. However, a complete mitochondrial genome is currently only available for Nanorana pleskei. Thus, we analyzed the complete mitochondrial genomes of Nanorana parkeri and Nanorana ventripunctata to investigate their evolutionary relationships within Nanorana and their phylogenetic position in the family Dicroglossidae. Our results showed that the genomes of N. parkeri (17,837 bp) and N. ventripunctata (18,373 bp) encode 13 protein‐coding genes (PCGs), two ribosomal RNA genes, 23 transfer RNA (tRNA) genes, and a noncoding control region. Overall sequences and genome structure of the two species showed high degree of similarity with N. pleskei, although the motif structures and repeat sequences of the putative control region showed clear differences among these three Nanorana species. In addition, a tandem repeat of the tRNA‐Met gene was found located between the tRNA‐Gln and ND2 genes. On both the 5′ and 3′‐sides, the control region possessed distinct repeat regions; however, the CSB‐2 motif was not found in N. pleskei. Based on the nucleotide sequences of 13 PCGs, our phylogenetic analyses, using Bayesian inference and maximum‐likelihood methods, illustrate the taxonomic status of Nanorana with robust support showing that N. ventripunctata and N. pleskei are more closely related than they are to N. parkeri. In conclusion, our analyses provide a more robust and reliable perspective on the evolutionary history of Dicroglossidae than earlier analyses, which used only a single species (N. pleskei).     

Waves of genomic hitchhikers shed light on the evolution of gamebirds (Aves: Galliformes)

Background  

The phylogenetic tree of Galliformes (gamebirds, including megapodes, currassows, guinea fowl, New and Old World quails, chicken, pheasants, grouse, and turkeys) has been considerably remodeled over the last decades as new data and analytical methods became available. Analyzing presence/absence patterns of retroposed elements avoids the problems of homoplastic characters inherent in other methodologies. In gamebirds, chicken repeats 1 (CR1) are the most prevalent retroposed elements, but little is known about the activity of their various subtypes over time. Ascertaining the fixation patterns of CR1 elements would help unravel the phylogeny of gamebirds and other poorly resolved avian clades.     

Mitochondrial genomes advance phylogenetic hypotheses for Tylenchina (Nematoda: Chromadorea)

Within the nematode class Chromadorea, the suborder Tylenchina is an ecologically and morphologically diverse assemblage of nematodes that includes free‐living microbivores, fungivores and various types of plant parasites. A recent nematode classification system based largely on SSU rDNA phylogenetic trees classified suborder Tylenchina to include four infraorders: Panagrolaimomorpha, Cephalobomorpha, Tylenchomorpha and Drilonematomorpha, and phylogenetic relationships among species of these infraorders have not always been robustly supported. In this study, we determined the complete mitochondrial genome sequences of three Tylenchina species (Aphelenchus avenae [Aphelenchidae, Tylenchomorpha], Halicephalobus gingivalis, Panagrellus redivivus [Panagrolaimomorpha]) and the partial genome sequence of Acrobeles complexus (Cephalobomorpha) and used these sequences to infer phylogenetic relationships among representatives of the Tylenchina and other nematodes. Phylogenetic analysis of amino acid sequences for 12 protein‐coding genes of 100 nematode species supports monophyly of: Chromadorea, Spiruromorpha, Oxyuridomorpha, Ascarididae + Toxocaridae + Anisakidae, Meloidogynidae + Pratylenchidae + Heteroderidae and Aphelenchoidea. Bayesian and maximum‐likelihood analyses also show the nested position of Diplogasteromorpha within Rhabditomorpha, and Rhigonematomorpha within Ascaridomorpha. These analyses also show non‐monophyly of: clade III, Ancylostomatidae, Panagrolaimomorpha, Tylenchina and Tylenchomorpha. Reconstructed mitochondrial genome phylogeny also revealed that among two main Tylenchomorpha groups, the monophyletic group representing Aphelenchoidea species was sister to the large clade consisting of Ascaridomorpha, Diplogasteromorpha, Rhabditomorpha and Rhigonematomorpha and some Panagrolaimomorpha species, whereas Tylenchoidea species were sister to the most inclusive assemblage containing all infraordinal groups of Chromadorea, except for P. redivivus (Panagrolaimomorpha) and Acrobeles complexus (Cephalobomorpha). The monophyly of Aphelenchoidea (i.e. sister relationship between Aphelenchidae and Aphelenchoididae) recovered in this study indicates that similarity in certain aspects of pharyngeal structure between these two families appears best explained by common ancestry, rather than convergent evolution.     

Avian papillomaviruses: the parrot <Emphasis Type="Italic">Psittacus erithacus</Emphasis> papillomavirus (PePV) genome has a unique organization of the early protein region and is phylogenetically related to the chaffinch papillomavirus

Background  

An avian papillomavirus genome has been cloned from a cutaneous exophytic papilloma from an African grey parrot (Psittacus erithacus). The nucleotide sequence, genome organization, and phylogenetic position of the Psittacus erithacus papillomavirus (PePV) were determined. This PePV sequence represents the first complete avian papillomavirus genome defined.     

The first complete mitochondrial genome of the Indian Tent Turtle,Pangshura tentoria (Testudines: Geoemydidae): Characterization and comparative analysis

The characterization of a complete mitogenome is widely used in genomics studies for systematics and evolutionary research. However, the sequences and structural motifs contained within the mitogenome of Testudines taxa have rarely been examined. The present study decodes the first complete mitochondrial genome of the Indian Tent Turtle, Pangshura tentoria (16,657 bp) by using next‐generation sequencing. This denovo assembly encodes 37 genes: 13 protein‐coding genes (PCGs), 22 transfer RNA (tRNAs), two ribosomal RNA, and one control region (CR). Most of the genes were encoded on majority strand, except for one PCG (NADH dehydrogenase subunit 6) and eight tRNAs. Most of the PCGs were started with an ATG initiation codon, except for Cytochrome oxidase subunit 1 with “GTG” and NADH dehydrogenase subunit 5 with “ATA.” The termination codons, “TAA” and “AGA” were observed in two subunits of NADH dehydrogenase gene. The relative synonymous codon usage analysis revealed the maximum abundance of alanine, isoleucine, leucine, and threonine. The nonsynonymous/synonymous ratios were <1 in all PCGs, which indicates strong negative selection among all Geoemydid species. The study also found the typical cloverleaf secondary structure in most of the tRNA genes, except for serine with the lack of the conventional DHU arm. The comparative study of Geoemydid mitogenomes revealed the occurrence of tandem repeats was frequent in the 3′ end of CR. Further, two copies of a unique tandem repeat “TTCTCTTT” were identified in P. tentoria. The Bayesian and maximum‐likelihood phylogenetic trees using concatenation of 13 PCGs revealed the close relationships of P. tentoria with Batagur trivittata in the studied dataset. All the Geoemydid species showed distinct clustering with high bootstrap support congruent with previous evolutionary hypotheses. We suggest that the generations of more mitogenomes of Geoemydid species are required, to improve our understanding of their in‐depth phylogenetic and evolutionary relationships.     

Evolution,biogeography and systematics of the western Palaearctic Zamenis ratsnakes

The phylogenetic relationships between western Palaearctic Zamenis and Rhinechis ratsnakes have been troubled, with recent estimates based on the supermatrix approach questioning their monophyly and providing contradictory results. In this study, we generated a comprehensive molecular data set for Zamenis and closely related ratsnakes to assess their phylogenetic and systematic relationships and infer their spatial and temporal modes of diversification. We obtained a fully resolved and well‐supported phylogeny, which is consistent across markers, taxon‐sets and phylogenetic methods. The close phylogenetic relationship between Rhinechis and Zamenis is well‐established. However, the early branching pattern within this clade, and the position of R. scalaris and Z. hohenackeri, remains poorly supported. The Persian ratsnake Z. persicus is sister to the Mediterranean species Z. situla, Z. longissimus and Z. lineatus, of which Z. situla is sister to a clade containing the latter two species. These results are consistent with a recent phylogenomic study on ratsnakes based on hundreds of loci. Whereas, topological tests based on our data and evidence from such phylogenomic study strongly rejected previous phylogenetic estimates based on the supermatrix approach and demonstrate that these “mega‐phylogenies”, with hundreds of taxa and high levels of missing data, have recovered inconsistent relationships with spurious nodal support. Biogeographical and molecular dating analyses suggest an origin of the ancestor of Rhinechis and Zamenis in the Aegean region with early cladogenesis during the Late Miocene associated with the Aegean arch formation and support a scenario of east‐to‐west diversification. Finally, while we have little morphological and phylogenetic evidence for the distinctiveness between Rhinechis and Zamenis, a classification of them in a single genus, and the designation of Zamenis scalaris (Schinz, 1822), reflects better their evolutionary relationships.     

Internal Transcribed Spacer Sequence Analysis of Puccinia helianthi Schw. and Its Application in Detection of Sunflower Rust

Two basidiomycete‐specific primers ITS1‐F and ITS4‐B were used in identification of the genus Puccinia. The primers showed good specificity for the genus with an 816‐bp product that was amplified exclusively. Twenty sequences of internal transcribed spacer (ITS) regions of Puccinia helianthi isolates from China remain unchanged. The whole ITS length (including ITS1 sequence 194 bp, 5.8S rRNA gene 156 bp, ITS2 sequence 206 bp) was 556 bp. By comparing the aligned ITS sequences of several Puccinia isolates from China, Spain and the United States, ITS homogeneity among these sunflower rust isolates was >99%. Genetic homology and phylogeny of P. helianthi with other Puccinia spp. was investigated. Nineteen sequences of rDNA ITS1 and ITS2 were determined and used as phylogenetic markers. Phylogenetic analysis of ITS regions showed that Puccinia spp. of sunflower was clustered in one clade with P. komarovii and P. violae, divergent from Puccinia spp. of Chrysanthemum, P. tenaceti of tansy (Tanacetum vulgare) and Puccina spp. of big sagebrush (Artemisia tridentate) indicating sunflower rust had distant phylogenetic relationships with other Compositae rusts. With the specified primers SR‐1 and SR‐2, either from purified urediniospores or symptomless (but infected) sunflower leaves could be examined specifically. Therefore, results of this study help in detection and polygenetic study of rust fungi occurring on sunflower.     

Mitogenomic phylogeny of Nassarius (Gastropoda: Neogastropoda)

Nassariids (Family Nassariidae) are a group of marine snails that are distributed worldwide, with their maximum species diversity in tropical regions, particularly the Indo‐Pacific. However, the traditional taxonomy of Nassariidae defined by shell or radula characters is usually inconsistent with little phylogenetic signal. In the present study, the complete mitochondrial (mt) genomes of nine Nassarius species were sequenced and compared with other eight nassariid species previously reported. All nassariid mt genomes showed the same gene order as in most caenogastropods and shared a very similar pattern with respect to genome size, nucleotide composition and AT contents. A deletion of three nucleotides in nad6 gene was detected in Nassarius jacksonianus and Nassarius acuticostus, and this feature also provided implications for nassariid phylogeny. The genetic distance analysis and reconstructed phylogeny revealed a distant relationship between N. jacksonianus or N. acuticostus and other members in Nassarius. The mitogenomic phylogeny recovered the evolutionary relationships within Nassarius with high statistical support. In addition, a chronogram was reconstructed under an uncorrelated relaxed molecular clock, which dated the divergence among main lineages of Nassarius during ~31 MYA.     

Mitochondrial genome evolution in fire ants (Hymenoptera: Formicidae)

Background  

Complete mitochondrial genome sequences have become important tools for the study of genome architecture, phylogeny, and molecular evolution. Despite the rapid increase in available mitogenomes, the taxonomic sampling often poorly reflects phylogenetic diversity and is often also biased to represent deeper (family-level) evolutionary relationships.     

Phylogenetic positioning of the Antarctic alga Prasiola crispa (Trebouxiophyceae) using organellar genomes and their structural analysis

Antarctica is one of the most difficult habitats for sustaining life on earth; organisms that live there have developed different strategies for survival. Among these organisms is the green alga Prasiola crispa, belonging to the class Trebouxiophyceae. The literature on P. crispa taxonomy is scarce, and many gaps in the evolutionary relationship with its closest relatives remain. The goal of this study was to analyze the evolutionary relationships between P. crispa and other green algae using plastid and mitochondrial genomes. In addition, we analyzed the synteny conservation of these genomes of P. crispa with those of closely related species. Based on the plastid genome, P. crispa grouped with Prasiolopsis sp. SAG 84.81, another Trebouxiophyceaen species from the Prasiola clade. Based on the mitochondrial genome analysis, P. crispa grouped with other Trebouxiophyceaen species but had a basal position. The structure of the P. crispa chloroplast genome had low synteny with Prasiolopsis sp. SAG 84.81, despite some conserved gene blocks. The same was observed in the mitochondrial genome compared with Coccomyxa subellipsoidea C‐169. We were able to establish the phylogenetic position of P. crispa with other species of Trebouxiophyceae using its genomes. In addition, we described the plasticity of these genomes using a structural analysis. The plastid and mitochondrial genomes of P. crispa will be useful for further genetic studies, phylogenetic analysis and resource protection of P. crispa as well as for further phylogenetic analysis of Trebouxiophyceaen green algae.     

Phylogeny of Veneridae (Bivalvia) based on mitochondrial genomes

Veneridae is one of the most diverse families of bivalve molluscs. However, their phylogenetic relationships among subfamilies have been debated for years. To explore phylogenetic relationships of Veneridae, we sequenced 13 complete mitochondrial genome sequences from eight subfamilies and compared with available complete mitochondrial genome of other Veneridae taxa (18 previously reported sequences). Phylogenetic analyses using probabilistic methods recovered two highly supported clades. In addition, the protein‐coding gene order revealed a highly conserved pattern among the same subclade lineages. According to our molecular analyses, Tapetinae should be recognized as a valid subfamily, but the genera formed para‐polyphyletic clades. Chioninae was recovered not monophyletic that differs from a previously molecular phylogeny. Furthermore, the reconstructed chronogram calibrated with fossils recovered the Veneridae have originated during the early Permian (about 290 million years ago). Noticeably, programmed frameshift was found in the nad4 gene of Leukoma jedoensis, Anomalodiscus squamosus and Antigona lamellaris and cob gene of L. jedoensis. This is the first time that the presence of the programmed frameshift has been found in the protein‐coding genes of Heterodonta species. Our results improved the phylogenetic resolution within Veneridae, and a more taxonomic sampling analysis of the subfamily Chioninae is supposed to construct.     

Flip‐flop organization in the chloroplast genome of Capsosiphon fulvescens (Ulvophyceae,Chlorophyta)

To better understand organelle genome evolution of the ulvophycean green alga Capsosiphon fulvescens, we sequenced and characterized its complete chloroplast genome. The circular chloroplast genome was 111,561 bp in length with 31.3% GC content that contained 108 genes including 77 protein‐coding genes, two copies of rRNA operons, and 27 tRNAs. In this analysis, we found the two types of isoform, called heteroplasmy, were likely caused by a flip‐flop organization. The flip‐flop mechanism may have caused structural variation and gene conversion in the chloroplast genome of C. fulvescens. In a phylogenetic analysis based on all available ulvophycean chloroplast genome data, including a new C. fulvescens genome, we found three major conflicting signals for C. fulvescens and its sister taxon Pseudoneochloris marina within 70 individual genes: (i) monophyly with Ulotrichales, (ii) monophyly with Ulvales, and (iii) monophyly with the clade of Ulotrichales and Ulvales. Although the 70‐gene concatenated phylogeny supported monophyly with Ulvales for both species, these complex phylogenetic signals of individual genes need further investigations using a data‐rich approach (i.e., organelle genome data) from broader taxon sampling.     

The use of DNA barcodes to estimate phylogenetic diversity in forest communities of southern China

To elucidate potential ecological and evolutionary processes associated with the assembly of plant communities, there is now widespread use of estimates of phylogenetic diversity that are based on a variety of DNA barcode regions and phylogenetic construction methods. However, relatively few studies consider how estimates of phylogenetic diversity may be influenced by single DNA barcodes incorporated into a sequence matrix (conservative regions vs. hypervariable regions) and the use of a backbone family‐level phylogeny. Here, we use general linear mixed‐effects models to examine the influence of different combinations of core DNA barcodes (rbcL, matK, ITS, and ITS2) and phylogeny construction methods on a series of estimates of community phylogenetic diversity for two subtropical forest plots in Guangdong, southern China. We ask: (a) What are the relative influences of single DNA barcodes on estimates phylogenetic diversity metrics? and (b) What is the effect of using a backbone family‐level phylogeny to estimate topology‐based phylogenetic diversity metrics? The combination of more than one barcode (i.e., rbcL + matK + ITS) and the use of a backbone family‐level phylogeny provided the most parsimonious explanation of variation in estimates of phylogenetic diversity. The use of a backbone family‐level phylogeny showed a stronger effect on phylogenetic diversity metrics that are based on tree topology compared to those that are based on branch lengths. In addition, the variation in the estimates of phylogenetic diversity that was explained by the top‐rank models ranged from 0.1% to 31% and was dependent on the type of phylogenetic community structure metric. Our study underscores the importance of incorporating a multilocus DNA barcode and the use of a backbone family‐level phylogeny to infer phylogenetic diversity, where the type of DNA barcode employed and the phylogenetic construction method used can serve as a significant source of variation in estimates of phylogenetic community structure.     

Molecular phylogeny of the Cladophoraceae (Cladophorales,Ulvophyceae), with the resurrection of Acrocladus Nägeli and Willeella Børgesen,and the description of Lurbica gen. nov. and Pseudorhizoclonium gen. nov.

The taxonomy of the Cladophoraceae, a large family of filamentous green algae, has been problematic for a long time due to morphological simplicity, parallel evolution, phenotypic plasticity, and unknown distribution ranges. Partial large subunit (LSU) rDNA sequences were generated for 362 isolates, and the analyses of a concatenated dataset consisting of unique LSU and small subunit (SSU) rDNA sequences of 95 specimens greatly clarified the phylogeny of the Cladophoraceae. The phylogenetic reconstructions showed that the three currently accepted genera Chaetomorpha, Cladophora, and Rhizoclonium are polyphyletic. The backbone of the phylogeny is robust and the relationships of the main lineages were inferred with high support, only the phylogenetic position of both Chaetomorpha melagonium and Cladophora rupestris could not be inferred unambiguously. There have been at least three independent switches between branched and unbranched morphologies within the Cladophoraceae. Freshwater environments have been colonized twice independently, namely by the freshwater Cladophora species as well as by several lineages of the Rhizoclonium riparium clade. In an effort to establish monophyletic genera, the genera Acrocladus and Willeella are resurrected and two new genera are described: Pseudorhizoclonium and Lurbica.     

Chromosome‐level genome assembly of the predator Propylea japonica to understand its tolerance to insecticides and high temperatures

The ladybird beetle Propylea japonica is an important natural enemy in agro‐ecological systems. Studies on the strong tolerance of P. japonica to high temperatures and insecticides, and its population and phenotype diversity have recently increased. However, abundant genome resources for obtaining insights into stress‐resistance mechanisms and genetic intra‐species diversity for P. japonica are lacking. Here, we constructed the P. japonica genome maps using Pacific Bioscience (PacBio) and Illumina sequencing technologies. The genome size was 850.90 Mb with a contig N50 of 813.13 kb. The Hi‐C sequence data were used to upgrade draft genome assemblies; 4,777 contigs were assembled to 10 chromosomes; and the final draft genome assembly was 803.93 Mb with a contig N50 of 813.98 kb and a scaffold N50 of 100.34 Mb. Approximately 495.38 Mb of repeated sequences was annotated. The 18,018 protein‐coding genes were predicted, of which 95.78% were functionally annotated, and 1,407 genes were species‐specific. The phylogenetic analysis showed that P. japonica diverged from the ancestor of Anoplophora glabripennis and Tribolium castaneum ~ 236.21 million years ago. We detected that some important gene families involved in detoxification of pesticides and tolerance to heat stress were expanded in P. japonica, especially cytochrome P450 and Hsp70 genes. Overall, the high‐quality draft genome sequence of P. japonica will provide invaluable resource for understanding the molecular mechanisms of stress resistance and will facilitate the research on population genetics, evolution and phylogeny of Coccinellidae. This genome will also provide new avenues for conserving the diversity of predator insects.     

<Emphasis Type="Italic">Gallus gallus</Emphasis> aggrecan gene-based phylogenetic analysis of selected avian taxonomic groups

Mitochondrial DNA (mtDNA) sequences remain the most widely used for phylogenetic analysis in birds. A major limitation of mtDNA sequences, however, is that mitochondria genes are inherited as a single linkage group. Here we describe the use of a 540-bp DNA sequence corresponding to the G3 domain of Gallus gallus nuclear aggrecan gene (AGC1) for phylogenetic analysis of the main groups of Galliformes including Phasianidae, Numididae, and Odontophoridae. We also included species from Cracidae and Megapodiidae which are considered by some as Craciformes and others, including here as Galliformes. The uncorrected sequence divergence of the G3 fragments ranges from 1 among the grouses to 36% between some of the distant groups within Galliformes. These sequences contain 39–48% AT nucleotides and the ratios of transition versus transversion are above 1.5 in majority of the comparisons. Using G3 sequences from an Anseriform, Oxyura jamaicensis, as out-groups, phylogenetic trees were obtained using maximum parsimony and distance algorithms and bootstrap analyses. These trees were consistent with those described using Avian sarcoma and leucosis virus gag genes and those from amino acid sequences of hemoglobin and lysozyme c. Our data also support relationships among Galliformes which were defined using mtDNA sequences. In addition to the general support of the five main families of Galliformes, our data are also consistent with previous work that showed Francolinus africanus and Gallus gallus are in the same clade and that Tetraoninae is a well-supported monophyletic subfamily within Phasianidae. The results presented here suggest that the AGC1 sequences meet the criterion of novel nuclear DNA sequences that can be used to help resolve the relationships among Galliformes.