Search

Background: The Mephitidae is a family of skunks and stink-badgers that includes 12 extant species in four genera, namely, Mydaus, Conepatus, Mephitis and Spilogale. Mydaus is the only genus within Mephitidae found outside the American continent, with its distribution limited to the islands of Borneo, Indonesia and Philippines. There are two extant species of Mydaus i.e., javanensis and marchei. Currently, complete mitogenomes are unavailable for either species. Here, we present the characterization of the first complete mitogenome for the Sunda stink-badger (Mydaus javanensis) from the island of Borneo. Methods: Muscle tissue was obtained and the DNA was sequenced using a combination of Illumina Barcode Tagged Sequence (BTSeq) and Sanger sequencing techniques. The genome was annotated with MITOS and manually checked for accuracy. A circular map of the mitogenome was constructed with Proksee. Relative synonymous codon usage (RSCU) and codon frequency were calculated using MEGA-X. The protein coding genes (PCGs) were aligned with reference sequences from GenBank and used for the construction of phylogenetic trees (maximum liklihood (ML) and Bayesian inference (BI)). Additionally, due to the lack of available complete genomes in public databases, we constructed another tree with the cyt b gene. Results: The complete circular mitogenome was 16,391 base pairs in length. It comprises the typical 13 protein-coding genes, 22 tRNAs, two ribosomal RNA genes, one control region (CR) and an L-strand replication origin (OL). The G+C content was 38.1% with a clear bias towards A and T nucleotides. Of the 13 PGCs, only ND6 was positioned in the reverse direction, along with five other tRNAs. Five PCGs had incomplete stop codons and rely on post-transcriptional polyadenylation (TAA) for termination. Based on the codon count, Leucine was the most common amino acid (589), followed by Threonine (332) and Isoleucine (325). The ML and BI phylogenetic trees, based on concatenated PCGs and the cyt b gene, respectively, correctly clustered the species with other members of the Mephitidae family but were unique enough to set it apart from Conepatus, Mephitis and Spilogale. The results confirm Mydaus as a member of the mephitids and the mitogenome will be useful for evolutionary analysis and conservation of the species. [ABSTRACT FROM AUTHOR]

Background The Mephitidae is a family of skunks and stink-badgers that includes 12 extant species in four genera, namely, Mydaus, Conepatus, Mephitis and Spilogale. Mydaus is the only genus within Mephitidae found outside the American continent, with its distribution limited to the islands of Borneo, Indonesia and Philippines. There are two extant species of Mydaus i.e., javanensis and marchei. Currently, complete mitogenomes are unavailable for either species. Here, we present the characterization of the first complete mitogenome for the Sunda stink-badger (Mydaus javanensis) from the island of Borneo. Methods Muscle tissue was obtained and the DNA was sequenced using a combination of Illumina Barcode Tagged Sequence (BTSeq) and Sanger sequencing techniques. The genome was annotated with MITOS and manually checked for accuracy. A circular map of the mitogenome was constructed with Proksee. Relative synonymous codon usage (RSCU) and codon frequency were calculated using MEGA-X. The protein coding genes (PCGs) were aligned with reference sequences from GenBank and used for the construction of phylogenetic trees (maximum liklihood (ML) and Bayesian inference (BI)). Additionally, due to the lack of available complete genomes in public databases, we constructed another tree with the cyt b gene. Results The complete circular mitogenome was 16,391 base pairs in length. It comprises the typical 13 protein-coding genes, 22 tRNAs, two ribosomal RNA genes, one control region (CR) and an L-strand replication origin (OL). The G+C content was 38.1% with a clear bias towards A and T nucleotides. Of the 13 PGCs, only ND6 was positioned in the reverse direction, along with five other tRNAs. Five PCGs had incomplete stop codons and rely on post-transcriptional polyadenylation (TAA) for termination. Based on the codon count, Leucine was the most common amino acid (589), followed by Threonine (332) and Isoleucine (325). The ML and BI phylogenetic trees, based on concatenated PCGs and the cyt b gene, respectively, correctly clustered the species with other members of the Mephitidae family but were unique enough to set it apart from Conepatus, Mephitis and Spilogale. The results confirm Mydaus as a member of the mephitids and the mitogenome will be useful for evolutionary analysis and conservation of the species.

Kerstersia gyiorum is a Gram-negative bacterium found in various animals, including humans, where it has been associated with various infections. Knowledge of the basic biology of K. gyiorum is essential to understand the evolutionary strategies of niche adaptation and how this organism contributes to infectious diseases; however, genomic data about K. gyiorum is very limited, especially from non-human hosts. In this work, we sequenced 12 K. gyiorum genomes isolated from healthy free-living brown-throated sloths (Bradypus variegatus) in the Parque Estadual das Fontes do Ipiranga (São Paulo, Brazil), and compared them with genomes from isolates of human origin, in order to gain insights into genomic diversity, phylogeny, and host specialization of this species. Phylogenetic analysis revealed that these K. gyiorum strains are structured according to host. Despite the fact that sloth isolates were sampled from a single geographic location, the intra-sloth K. gyiorum diversity was divided into three clusters, with differences of more than 1,000 single nucleotide polymorphisms between them, suggesting the circulation of various K. gyiorum lineages in sloths. Genes involved in mobilome and defense mechanisms against mobile genetic elements were the main source of gene content variation between isolates from different hosts. Sloth-specific K. gyiorum genome features include an IncN2 plasmid, a phage sequence, and a CRISPR-Cas system. The broad diversity of defense elements in K. gyiorum (14 systems) may prevent further mobile element flow and explain the low amount of mobile genetic elements in K. gyiorum genomes. Gene content variation may be important for the adaptation of K. gyiorum to different host niches. This study furthers our understanding of diversity, host adaptation, and evolution of K. gyiorum, by presenting and analyzing the first genomes of non-human isolates. [ABSTRACT FROM AUTHOR]

Kerstersia gyiorum is a Gram-negative bacterium found in various animals, including humans, where it has been associated with various infections. Knowledge of the basic biology of K. gyiorum is essential to understand the evolutionary strategies of niche adaptation and how this organism contributes to infectious diseases; however, genomic data about K. gyiorum is very limited, especially from non-human hosts. In this work, we sequenced 12 K. gyiorum genomes isolated from healthy free-living brown-throated sloths (Bradypus variegatus) in the Parque Estadual das Fontes do Ipiranga (São Paulo, Brazil), and compared them with genomes from isolates of human origin, in order to gain insights into genomic diversity, phylogeny, and host specialization of this species. Phylogenetic analysis revealed that these K. gyiorum strains are structured according to host. Despite the fact that sloth isolates were sampled from a single geographic location, the intra-sloth K. gyiorum diversity was divided into three clusters, with differences of more than 1,000 single nucleotide polymorphisms between them, suggesting the circulation of various K. gyiorum lineages in sloths. Genes involved in mobilome and defense mechanisms against mobile genetic elements were the main source of gene content variation between isolates from different hosts. Sloth-specific K. gyiorum genome features include an IncN2 plasmid, a phage sequence, and a CRISPR-Cas system. The broad diversity of defense elements in K. gyiorum (14 systems) may prevent further mobile element flow and explain the low amount of mobile genetic elements in K. gyiorum genomes. Gene content variation may be important for the adaptation of K. gyiorum to different host niches. This study furthers our understanding of diversity, host adaptation, and evolution of K. gyiorum, by presenting and analyzing the first genomes of non-human isolates.

Background Blackcurrant (Ribes nigrum), red currant (R. rubrum), white currant (R. rubrum), and gooseberry (R. uva-crispa) belong to Grossulariaceae and are popular small-berry crops worldwide. The lack of genomic data has severely limited their systematic classification and molecular breeding. Methods The complete chloroplast (cp) genomes of these four taxa were assembled for the first time using MGI-DNBSEQ reads, and their genome structures, repeat elements and protein-coding genes were annotated. By genomic comparison of the present four and previous released five Ribes cp genomes, the genomic variations were identified. By phylogenetic analysis based on maximum-likelihood and Bayesian methods, the phylogeny of Grossulariaceae and the infrageneric relationships of the Ribes were revealed. Results The four cp genomes have lengths ranging from 157,450 to 157,802 bp and 131 shared genes. A total of 3,322 SNPs and 485 Indels were identified from the nine released Ribes cp genomes. Red currant and white currant have 100% identical cp genomes partially supporting the hypothesis that white currant (R. rubrum) is a fruit color variant of red currant (R. rubrum). The most polymorphic genic and intergenic region is ycf1 and trnT-psbD, respectively. The phylogenetic analysis demonstrated the monophyly of Grossulariaceae in Saxifragales and the paraphyletic relationship between Saxifragaceae and Grossulariaceae. Notably, the Grossularia subgenus is well nested within the Ribes subgenus and shows a paraphyletic relationship with the co-ancestor of Calobotrya and Coreosma sections, which challenges the dichotomous subclassification of the Ribes genus based on morphology (subgenus Ribes and subgenus Grossularia). These data, results, and insights lay a foundation for the phylogenetic research and breeding of Ribes species.

Much genomic-scale, especially transcriptomic, data on spider phylogeny has accumulated in the last few years. These data have recently been used to investigate the diverse architectures and the origin of spider webs, concluding that the ancestral spider spun no foraging web, that spider webs evolved de novo 10-14 times, and that the orb web evolved at least three times. These findings in fact result from a particular phylogenetic character coding strategy, specifically coding the absence of webs as logically equivalent, and homologous to, 10 other observable (i.e., not absent) web architectures. "Absence" of webs should be regarded as inapplicable data. To be analyzed properly by character optimization algorithms, it should be coded as "?" because these codes-or their equivalent-are handled differently by such algorithms. Additional problems include critical misspellings of taxon names from one analysis to the next (misspellings cause some optimization algorithms to drop terminals, which affects taxon sampling and results), and mistakes in spider natural history. In sum, the method causes character optimization algorithms to produce counter-intuitive results, and does not distinguish absence from secondary loss. Proper treatment of missing entries and corrected data instead imply that foraging webs are primitive for spiders and that webs have been lost ∼5-7 times, not gained 10-14 times. The orb web, specifically, may be homologous (originated only once) although lost 2-6 times.

Psychotria viridis (Rubioideae: Rubiaceae), popularly known as chacrona, is commonly found as a shrub in the Amazon region and is well-known to produce psychoactive compounds, such as the N,N-dimethyltryptamine (DMT). Together with the liana Banisteropsis caapi, P. viridis is one of the main components of the Amerindian traditional, entheogenic beverage known as ayahuasca. In this work, we assembled and annotated the organellar genomes (ptDNA and mtDNA), presenting the first genomics resources for this species. The P. viridis ptDNA exhibits 154,106 bp, encoding all known ptDNA gene repertoire found in angiosperms. The Psychotria genus is a complex paraphyletic group, and according to phylogenomic analyses, P. viridis is nested in the Psychotrieae clade. Comparative ptDNA analyses indicate that most Rubiaceae plastomes present conserved ptDNA structures, often showing slight differences at the junction sites of the major four regions (LSC-IR-SSC). For the mitochondrion, assembly graph-based analysis supports a complex mtDNA organization, presenting at least two alternative and circular mitogenomes structures exhibiting two main repeats spanning 24 kb and 749 bp that may symmetrically isomerize the mitogenome into variable arrangements and isoforms. The circular mtDNA sequences (615,370 and 570,344 bp) encode almost all plant mitochondrial genes (except for the ccmC, rps7, rps10, rps14, rps19, rpl2 and rpl16 that appears as pseudogenes, and the absent genes sdh3, rps2, rsp4, rsp8, rps11, rpl6, and rpl10), showing slight variations related to exclusive regions, ptDNA integration, and relics of previous events of LTR-RT integration. The detection of two mitogenomes haplotypes is evidence of heteroplasmy as observed by the complex organization of the mitochondrial genome using graph-based analysis. Taken together, these results elicit the primary insights into the genome biology and evolutionary history of Psychotria viridis and may be used to aid strategies for conservation of this sacred, entheogenic species. [ABSTRACT FROM AUTHOR]

North American minnows of the Shiner Clade, within the family Leuciscidae, represent one of the most taxonomically complex clades of the order Cypriniformes due to the large number of taxa coupled with conserved morphologies. Species within this clade were moved between genera and subgenera until the community decided to lump many of the unclassified taxa with similar morphologies into one genus, Notropis, which has held up to 325 species. Despite phylogentic studies that began to re-elevate some genera merged into Notropis, such as Cyprinella, Luxilus, Lythrurus, and Pteronotropis, the large genus Notropis remained as a taxonomic repository for many shiners of uncertain placement. Recent molecular advances in sequencing technologies have provided the opportunity to re-examine the Shiner Clade using phylogenomic markers. Using a fish probe kit, we sequenced 90 specimens in 87 species representing 16 genera included in the Shiner Clade, with a resulting dataset of 1,004 loci and 286,455 base pairs. Despite the large dataset, only 32,349 bp (11.29%) were phylogenetically informative. In our maximum likelihood tree, 78% of nodes are 100% bootstrap supported demonstrating the utility of the phylogenomic markers at lower taxonomic levels. Unsurprisingly, species within Notropis as well as Hudsonius, Luxilus, and Alburnops are not resolved as monophyletic groups. Cyprinella is monophyletic if Cyprinella callistia is excluded, and Pteronotropis is monophyletic if it includes Hudsonius cummingsae. Taxonomic changes we propose are: restriction of species included in Alburnops and Notropis, elevation of the subgenus Hydrophlox, expansion of species included in Miniellus, movement of Hudsonius cummingsae to Pteronotropis, and resurrection of the genera Coccotis and Paranotropis. We additionally had two specimens of three species, Notropis atherinoides, Ericymba amplamala, and Pimephales vigilax and found signficant differences between the localities (1,086, 1,424, and 845 nucleotides respectively). [ABSTRACT FROM AUTHOR]

Psychotria viridis (Rubioideae: Rubiaceae), popularly known as chacrona, is commonly found as a shrub in the Amazon region and is well-known to produce psychoactive compounds, such as the N,N-dimethyltryptamine (DMT). Together with the liana Banisteropsis caapi, P. viridis is one of the main components of the Amerindian traditional, entheogenic beverage known as ayahuasca. In this work, we assembled and annotated the organellar genomes (ptDNA and mtDNA), presenting the first genomics resources for this species. The P. viridis ptDNA exhibits 154,106 bp, encoding all known ptDNA gene repertoire found in angiosperms. The Psychotria genus is a complex paraphyletic group, and according to phylogenomic analyses, P. viridis is nested in the Psychotrieae clade. Comparative ptDNA analyses indicate that most Rubiaceae plastomes present conserved ptDNA structures, often showing slight differences at the junction sites of the major four regions (LSC-IR-SSC). For the mitochondrion, assembly graph-based analysis supports a complex mtDNA organization, presenting at least two alternative and circular mitogenomes structures exhibiting two main repeats spanning 24 kb and 749 bp that may symmetrically isomerize the mitogenome into variable arrangements and isoforms. The circular mtDNA sequences (615,370 and 570,344 bp) encode almost all plant mitochondrial genes (except for the ccmC, rps7, rps10, rps14, rps19, rpl2 and rpl16 that appears as pseudogenes, and the absent genes sdh3, rps2, rsp4, rsp8, rps11, rpl6, and rpl10), showing slight variations related to exclusive regions, ptDNA integration, and relics of previous events of LTR-RT integration. The detection of two mitogenomes haplotypes is evidence of heteroplasmy as observed by the complex organization of the mitochondrial genome using graph-based analysis. Taken together, these results elicit the primary insights into the genome biology and evolutionary history of Psychotria viridis and may be used to aid strategies for conservation of this sacred, entheogenic species.

North American minnows of the Shiner Clade, within the family Leuciscidae, represent one of the most taxonomically complex clades of the order Cypriniformes due to the large number of taxa coupled with conserved morphologies. Species within this clade were moved between genera and subgenera until the community decided to lump many of the unclassified taxa with similar morphologies into one genus, Notropis, which has held up to 325 species. Despite phylogentic studies that began to re-elevate some genera merged into Notropis, such as Cyprinella, Luxilus, Lythrurus, and Pteronotropis, the large genus Notropis remained as a taxonomic repository for many shiners of uncertain placement. Recent molecular advances in sequencing technologies have provided the opportunity to re-examine the Shiner Clade using phylogenomic markers. Using a fish probe kit, we sequenced 90 specimens in 87 species representing 16 genera included in the Shiner Clade, with a resulting dataset of 1,004 loci and 286,455 base pairs. Despite the large dataset, only 32,349 bp (11.29%) were phylogenetically informative. In our maximum likelihood tree, 78% of nodes are 100% bootstrap supported demonstrating the utility of the phylogenomic markers at lower taxonomic levels. Unsurprisingly, species within Notropis as well as Hudsonius, Luxilus, and Alburnops are not resolved as monophyletic groups. Cyprinella is monophyletic if Cyprinella callistia is excluded, and Pteronotropis is monophyletic if it includes Hudsonius cummingsae. Taxonomic changes we propose are: restriction of species included in Alburnops and Notropis, elevation of the subgenus Hydrophlox, expansion of species included in Miniellus, movement of Hudsonius cummingsae to Pteronotropis, and resurrection of the genera Coccotis and Paranotropis. We additionally had two specimens of three species, Notropis atherinoides, Ericymba amplamala, and Pimephales vigilax and found signficant differences between the localities (1,086, 1,424, and 845 nucleotides respectively).

Baird's tapir, or the Central American Tapir Tapirus bairdii (family Tapiridae), is one of the largest mammals native to the forests and wetlands of southern North America and Central America, and is categorized as 'endangered' on the 2014 IUCN Red List of Threatened Species. This study reports, for the first time, the complete mitochondrial genome of T. bairdii and examines the phylogenetic position of T. bairdii amongst closely related species in the same family and order to which it belongs using mitochondrial protein-coding genes (PCG's). The circular, double-stranded, A-T rich mitochondrial genome of T. bairdii is 16,697 bp in length consisting of 13 proteincoding genes (PCG's), two ribosomal RNA genes (rrnS (12s ribosomal RNA and rrnL (16s ribosomal RNA)), and 22 transfer RNA (tRNA) genes. A 33 bp long region was identified to be the origin of replication for the light strand (OL), and a 1,247 bp long control region (CR) contains the origin of replication for the heavy strand (OH). A majority of the PCG's and tRNA genes are encoded on the positive, or heavy, strand. The gene order in T. baiirdi is identical to that of T. indicus and T. terrestris, the only two other species of extant tapirs with assembled mitochondrial genomes. An analysis of Ka/Ks ratios for all the PCG's show values<1, suggesting that all these PCGs experience strong purifying selection. A maximum-likelihood phylogenetic analysis supports the monophyly of the genus Tapirus and the order Perissodactyla. The complete annotation and analysis of the mitochondrial genome of T. bairdii will contribute to a better understanding of the population genomic diversity and structure of this species, and it will assist in the conservation and protection of its dwindling populations. [ABSTRACT FROM AUTHOR]

Baird’s tapir, or the Central American Tapir Tapirus bairdii (family Tapiridae), is one of the largest mammals native to the forests and wetlands of southern North America and Central America, and is categorized as ‘endangered’ on the 2014 IUCN Red List of Threatened Species. This study reports, for the first time, the complete mitochondrial genome of T. bairdii and examines the phylogenetic position of T. bairdii amongst closely related species in the same family and order to which it belongs using mitochondrial protein-coding genes (PCG’s). The circular, double-stranded, A-T rich mitochondrial genome of T. bairdii is 16,697 bp in length consisting of 13 protein-coding genes (PCG’s), two ribosomal RNA genes (rrnS (12s ribosomal RNA and rrnL (16s ribosomal RNA)), and 22 transfer RNA (tRNA) genes. A 33 bp long region was identified to be the origin of replication for the light strand (OL), and a 1,247 bp long control region (CR) contains the origin of replication for the heavy strand (OH). A majority of the PCG’s and tRNA genes are encoded on the positive, or heavy, strand. The gene order in T. baiirdi is identical to that of T. indicus and T. terrestris, the only two other species of extant tapirs with assembled mitochondrial genomes. An analysis of Ka/Ks ratios for all the PCG’s show values

Candidate phylum KSB1 is composed of uncultured bacteria and has been reported across various environments. However, the phylogeny and metabolic potential of KSB1 have not been studied comprehensively. In this study, phylogenomic analysis of KSB1 genomes from public databases and eleven metagenome-assembled genomes (MAGs) from marine and hydrothermal sediments revealed that those genomes were clustered into four clades. Isolation source and relative abundance of KSB1 genomes showed that clade I was particularly abundant in bioreactor sludge. Genes related to dissimilatory reduction of nitrate to ammonia (DNRA), the last step of denitrification converting nitrous oxide to nitrogen and assimilatory sulfur reduction were observed in the expanded genomes of clade I, which may due to horizontal gene transfer that frequently occurred in bioreactor. Annotation and metabolic reconstruction of clades II and IV showed flagellum assembly and chemotaxis genes in the genomes, which may indicate that exploration and sensing for nutrients and chemical gradients are critical for the two clades in deep-sea and hydrothermal sediment. Metabolic potentials of fatty acids and short-chain hydrocarbons utilization were predicted in clades I and IV of KSB1. Collectively, phylogenomic and metabolic analyses of KSB1 clades provide insight into their anaerobic heterotrophic lifestyle and differentiation in potential ecological roles.

The genus Escherichia has been extensively studied and it is known to encompass a range of commensal and pathogenic bacteria that primarily inhabit the gastrointestinal tracts of warm-blooded vertebrates. However, the presence of E. coli as a model organism and potential pathogen has diverted attention away from commensal strains and other species in the genus. To investigate the diversity of Escherichia in healthy chickens, we collected fecal samples from antibiotic-free Lohmann Brown layer hens and determined the genome sequences of 100 isolates, 81 of which were indistinguishable at the HC0 level of the Hierarchical Clustering of Core Genome Multi-Locus Sequence Typing scheme. Despite initial selection on CHROMagar Orientation medium, which is considered selective for E. coli, in silico phylotyping and core genome single nucleotide polymorphism analysis revealed the presence of at least one representative of all major clades of Escherichia, except for E. albertii, Shigella, and E. coli phylogroup B2 and cryptic clade I. The most frequent phylogenomic groups were E. coli phylogroups A and B1 and E. ruysiae (clades III and IV). We compiled a collection of reference strains isolated from avian sources (predominantly chicken), representing every Escherichia phylogroup and species, and used it to confirm the phylogeny and diversity of our isolates. Overall, the isolates carried low numbers of the virulence and antibiotic resistance genes typically seen in avian pathogenic E. coli. Notably, the clades not recovered are ones that have been most strongly associated with virulence by other studies. [ABSTRACT FROM AUTHOR]

The genus Escherichia has been extensively studied and it is known to encompass a range of commensal and pathogenic bacteria that primarily inhabit the gastrointestinal tracts of warm-blooded vertebrates. However, the presence of E. coli as a model organism and potential pathogen has diverted attention away from commensal strains and other species in the genus. To investigate the diversity of Escherichia in healthy chickens, we collected fecal samples from antibiotic-free Lohmann Brown layer hens and determined the genome sequences of 100 isolates, 81 of which were indistinguishable at the HC0 level of the Hierarchical Clustering of Core Genome Multi-Locus Sequence Typing scheme. Despite initial selection on CHROMagar Orientation medium, which is considered selective for E. coli, in silico phylotyping and core genome single nucleotide polymorphism analysis revealed the presence of at least one representative of all major clades of Escherichia, except for E. albertii, Shigella, and E. coli phylogroup B2 and cryptic clade I. The most frequent phylogenomic groups were E. coli phylogroups A and B1 and E. ruysiae (clades III and IV). We compiled a collection of reference strains isolated from avian sources (predominantly chicken), representing every Escherichia phylogroup and species, and used it to confirm the phylogeny and diversity of our isolates. Overall, the isolates carried low numbers of the virulence and antibiotic resistance genes typically seen in avian pathogenic E. coli. Notably, the clades not recovered are ones that have been most strongly associated with virulence by other studies.

North American watersheds contain a high diversity of freshwater mussels (Unionoida). During the long-lived, benthic phase of their life cycle, up to 40 species can co-occur in a single riffle and there is typically little evidence for major differences in their feeding ecology or microhabitat partitioning. In contrast, their brief parasitic larval phase involves the infection of a wide diversity of fish hosts and female mussels have evolved a spectrum of adaptations for infecting host fish with their offspring. Many species use a passive broadcast strategy: placing high numbers of larvae in the water column and relying on chance encounters with potential hosts. Many other species, including most members of the Lampsilini, have a proactive strategy that entails the use of prey-mimetic lures to change the behavior of the hosts, i.e., eliciting a feeding response through which they become infected. Two main lure types are collectively produced: mantle tissue lures (on the female’s body) and brood lures, containing infective larvae, that are released into the external environment. In this study, we used a phylogenomic approach (ddRAD-seq) to place the diversity of infection strategies used by 54 North American lampsiline mussels into an evolutionary context. Ancestral state reconstruction recovered evidence for the early evolution of mantle lures in this clade, with brood lures and broadcast infection strategies both being independently derived twice. The most common infection strategy, occurring in our largest ingroup clade, is a mixed one in which mimetic mantle lures are apparently the predominant infection mechanism, but gravid females also release simple, non-mimetic brood lures at the end of the season. This mixed infection strategy clade shows some evidence of an increase in diversification rate and most members use centrarchids (Micropterus & Lepomis spp.) as their predominant fish hosts. Broad linkage between infection strategies and predominant fish host genera is also seen in other lampsiline clades: worm-like mantle lures of Toxolasma spp. with sunfish (Lepomisspp.); insect larvae-like brood lures (Ptychobranchus spp.), or mantle lures (Medionidus spp., Obovaria spp.), or mantle lures combined with host capture (Epioblasma spp.) with a spectrum of darter (Etheostoma & Percina spp.) and sculpin (Cottus spp.) hosts, and tethered brood lures (Hamiota spp.) with bass (Micropterus spp.). Our phylogenetic results confirm that discrete lampsiline mussel clades exhibit considerable specialization in the primary fish host clades their larvae parasitize, and in the host infection strategies they employ to do so. They are also consistent with the hypothesis that larval resource partitioning of fish hosts is an important factor in maintaining species diversity in mussel assemblages. We conclude that, taking their larval ecology and host-infection mechanisms into account, lampsiline mussels may be legitimately viewed as an adaptive radiation. [ABSTRACT FROM AUTHOR]

Several main families of Ranunculales are rich in alkaloids and other medicinal compounds; many species of these families are used in traditional and folk medicine. Dichocarpum is a representative medicinal genus of Ranunculaceae, but the genetic basis of its metabolic phenotype has not been investigated, which hinders its sustainable conservation and utilization. We use the third-generation high-throughput sequencing and metabolomic techniques to decipher the full-length transcriptomes and metabolomes of five Dichocarpum species endemic in China, and 71,598 non-redundant full-length transcripts were obtained, many of which are involved in defense, stress response and immunity, especially those participating in the biosynthesis of specialized metabolites such as benzylisoquinoline alkaloids (BIAs). Twenty-seven orthologs extracted from trancriptome datasets were concatenated to reconstruct the phylogenetic tree, which was verified by the clustering analysis based on the metabolomic profile and agreed with the Pearson correlation between gene expression patterns of Dichocarpum species. The phylogenomic analysis of phytometabolite biosynthesis genes, e.g., (S)-norcoclaurine synthase, methyltransferases, cytochrome p450 monooxygenases, berberine bridge enzyme and (S)-tetrahydroprotoberberine oxidase, revealed the evolutionary trajectories leading to the chemodiversity, especially that of protoberberine type, aporphine type and bis-BIA abundant in Dichocarpum and related genera. The biosynthesis pathways of these BIAs are proposed based on full-length transcriptomes and metabolomes of Dichocarpum. Within Ranunculales, the gene duplications are common, and a unique whole genome duplication is possible in Dichocarpum. The extensive correlations between metabolite content and gene expression support the co-evolution of various genes essential for the production of different specialized metabolites. Our study provides insights into the transcriptomic and metabolomic landscapes of Dichocarpum, which will assist further studies on genomics and application of Ranunculales plants.

North American watersheds contain a high diversity of freshwater mussels (Unionoida). During the long-lived, benthic phase of their life cycle, up to 40 species can co-occur in a single riffle and there is typically little evidence for major differences in their feeding ecology or microhabitat partitioning. In contrast, their brief parasitic larval phase involves the infection of a wide diversity of fish hosts and female mussels have evolved a spectrum of adaptations for infecting host fish with their offspring. Many species use a passive broadcast strategy: placing high numbers of larvae in the water column and relying on chance encounters with potential hosts. Many other species, including most members of the Lampsilini, have a proactive strategy that entails the use of prey-mimetic lures to change the behavior of the hosts, i.e., eliciting a feeding response through which they become infected. Two main lure types are collectively produced: mantle tissue lures (on the female’s body) and brood lures, containing infective larvae, that are released into the external environment. In this study, we used a phylogenomic approach (ddRAD-seq) to place the diversity of infection strategies used by 54 North American lampsiline mussels into an evolutionary context. Ancestral state reconstruction recovered evidence for the early evolution of mantle lures in this clade, with brood lures and broadcast infection strategies both being independently derived twice. The most common infection strategy, occurring in our largest ingroup clade, is a mixed one in which mimetic mantle lures are apparently the predominant infection mechanism, but gravid females also release simple, non-mimetic brood lures at the end of the season. This mixed infection strategy clade shows some evidence of an increase in diversification rate and most members use centrarchids (Micropterus & Lepomis spp.) as their predominant fish hosts. Broad linkage between infection strategies and predominant fish host genera is also seen in other lampsiline clades: worm-like mantle lures of Toxolasma spp. with sunfish (Lepomis spp.); insect larvae-like brood lures (Ptychobranchus spp.), or mantle lures (Medionidus spp., Obovaria spp.), or mantle lures combined with host capture (Epioblasma spp.) with a spectrum of darter (Etheostoma & Percina spp.) and sculpin (Cottus spp.) hosts, and tethered brood lures (Hamiota spp.) with bass (Micropterus spp.). Our phylogenetic results confirm that discrete lampsiline mussel clades exhibit considerable specialization in the primary fish host clades their larvae parasitize, and in the host infection strategies they employ to do so. They are also consistent with the hypothesis that larval resource partitioning of fish hosts is an important factor in maintaining species diversity in mussel assemblages. We conclude that, taking their larval ecology and host-infection mechanisms into account, lampsiline mussels may be legitimately viewed as an adaptive radiation.

Higher-level relationships of the Hydrozoan subclass Hydroidolina, which encompasses the vast majority of medusozoan cnidarian species diversity, have been elusive to confidently infer. The most widely adopted phylogenetic framework for Hydroidolina based on ribosomal RNA data received low support for several higher level relationships. To address this issue, we developed a set of RNA baits to target more than a hundred loci from the genomes of a broad taxonomic sample of Hydroidolina for high-throughput sequencing. Using these data, we inferred the relationships of Hydroidolina using maximum likelihood and Bayesian approaches. Both inference methods yielded well-supported phylogenetic hypotheses that largely agree with each other. Using maximum likelihood and Baysian hypothesis testing frameworks, we found that several alternate topological hypotheses proposed previously may be rejected in light of the genomic data generated for this study. Both the maximum likelihood and Bayesian topologies inferred herein consistently score well across testing frameworks, suggesting that their consensus represents the most likely phylogenetic hypothesis of Hydroidolina. This phylogenetic framework places Aplanulata as sister lineage to the remainder of Hydroidolina. This is a strong deviation from previous phylogenetic analyses that placed Capitata or Siphonophorae as sister group to the remainder of Hydroidolina. Considering that Aplanulata represents a lineage comprised of species that for the most part possess a life cycle involving a solitary polyp and free-swimming medusa stage, the phylogenetic hypotheses presented herein have potentially large implications for clarifying the evolution of life cycles, coloniality, and the division of labor in Hydrozoa as taxon sampling for phylogenetic analyses becomes more complete.

Higher-level relationships of the Hydrozoan subclass Hydroidolina, which encompasses the vast majority of medusozoan cnidarian species diversity, have been elusive to confidently infer. The most widely adopted phylogenetic framework for Hydroidolina based on ribosomal RNA data received low support for several higher level relationships. To address this issue, we developed a set of RNA baits to target more than a hundred loci from the genomes of a broad taxonomic sample of Hydroidolina for high-throughput sequencing. Using these data, we inferred the relationships of Hydroidolina using maximum likelihood and Bayesian approaches. Both inference methods yielded wellsupported phylogenetic hypotheses that largely agree with each other. Using maximum likelihood and Baysian hypothesis testing frameworks, we found that several alternate topological hypotheses proposed previously may be rejected in light of the genomic data generated for this study. Both the maximum likelihood and Bayesian topologies inferred herein consistently score well across testing frameworks, suggesting that their consensus represents the most likely phylogenetic hypothesis of Hydroidolina. This phylogenetic framework places Aplanulata as sister lineage to the remainder of Hydroidolina. This is a strong deviation from previous phylogenetic analyses that placed Capitata or Siphonophorae as sister group to the remainder of Hydroidolina. Considering that Aplanulata represents a lineage comprised of species that for the most part possess a life cycle involving a solitary polyp and free-swimming medusa stage, the phylogenetic hypotheses presented herein have potentially large implications for clarifying the evolution of life cycles, coloniality, and the division of labor in Hydrozoa as taxon sampling for phylogenetic analyses becomes more complete. [ABSTRACT FROM AUTHOR]

Here, I review phylogenetic studies of the lizard family Pygopodidae, a group of 47 extant species that diversified in Australia and New Guinea. The goal of this study was to examine published phylogenetic and phylogenomic hypotheses on pygopodids to identify the strengths and weaknesses in our understanding of their phylogeny. Many parts of the pygopodid family tree are well established by multiple independent tree inferences including: (1) all multispecies genera (i.e., Aprasia, Delma, Lialis, Pletholax, and Pygopus) are monophyletic groups; (2) the root of the pygopodid tree is located along the branch leading to the Delma clade, thus showing that Delma is the sister group to all other pygopodid genera; (3) the Aprasia repens group, Delma tincta group, and several other groups of closely related species are demonstrated to be monophyletic entities; and (4) the monotypic Paradelma orientalis is the sister lineage to the Pygopus clade. Based on accumulated phylogenetic evidence, two taxonomic recommendations are given: Paradelma merits generic status rather than being subsumed into Pygopus as some earlier studies had suggested, and the monotypic Aclys concinna should be recognized as a member of Delma (following current practice) until future studies clarify its placement inside or outside the Delma clade. One chronic problem with phylogenetic studies of pygopodids, which has limited the explanatory power of many tree hypotheses, concerns the undersampling of known species. Although the continual addition of newly described species, especially over the past two decades, has been a major reason for these taxon sampling gaps, deficits in species sampling for ingroups and/or outgroups in several studies of pygopodid species complexes has confounded the testing of some ingroup monophyly hypotheses. Ancient hybridization between non-sister lineages may also be confounding attempts to recover the relationships among pygopodids using molecular data. Indeed, such a phenomenon can explain at least five cases of mito-nuclear discordance and conflicts among trees based on nuclear DNA datasets. Another problem has been the lack of consensus on the relationships among most pygopodid genera, an issue that may stem from rapid diversification of these lineages early in the group’s history. Despite current weaknesses in our understanding of pygopodid phylogeny, enough evidence exists to clarify many major and minor structural parts of their family tree. Accordingly, a composite tree for the Pygopodidae was able to be synthesized. This novel tree hypothesis contains all recognized pygopodid species and reveals that about half of the clades are corroborated by multiple independent tree hypotheses, while the remaining clades have less empirical support.

TQMD is a tool for high-performance computing clusters which downloads, stores and produces lists of dereplicated prokaryotic genomes. It has been developed to counter the ever-growing number of prokaryotic genomes and their uneven taxonomic distribution. It is based on word-based alignment-free methods (k-mers), an iterative single-linkage approach and a divide-and-conquer strategy to remain both efficient and scalable. We studied the performance of TQMD by verifying the influence of its parameters and heuristics on the clustering outcome. We further compared TQMD to two other dereplication tools (dRep and Assembly-Dereplicator). Our results showed that TQMD is primarily optimized to dereplicate at higher taxonomic levels (phylum/class), as opposed to the other dereplication tools, but also works at lower taxonomic levels (species/strain) like the other dereplication tools. TQMD is available from source and as a Singularity container at [https://bitbucket.org/phylogeno/tqmd ].

Phylogenomic studies can provide insights into speciation, adaptation, and extinction, while providing a roadmap for conservation. Hawaiian tree snails are a model system for an adaptive radiation facing an extinction crisis. In the last 5 years, nearly all populations of Hawaiian tree snails across the 30 remaining species in the subfamily Achatinellinae (Achatinellidae) have declined from hundreds or thousands in the wild down to undetectable levels. Nearly 100 species historically occurred across dramatic environmental gradients on five of the Hawaiian Islands, but habitat loss, overcollection, and predation by invasive species have decimated populations. As such, this system offers the opportunity to integrate efforts to conserve evolutionary potential into conservation planning for a rapidly declining subfamily. Here, we used genome-wide, restriction-site associated DNA sequencing (RADseq), along with mitochondrial genome reconstruction, to resolve evolutionary relationships to inform conservation efforts. Phylogenetic analysis of nearly 400k genome-wide SNPs from 59 populations and 25 species across six genera in the family Achatinellidae, was generally concordant with taxonomy, geography, and mtDNA with several notable exceptions; mtDNA was unable to resolve some deeper nodes (e.g., the monophyly of Achatinella), while SNP data did not resolve as many shallow nodes. Both phylogenetic and coalescent analysis revealed deep divergences between populations within Achatinella mustelina that were consistent with species-level differences. Given cryptic species-level divergence within populations that are geographically proximate, they are at higher risk of extirpation from invasive predators and climate change than previously assumed. This study clarifies evolutionary relationships within this model system for adaptive radiation, forming the basis for conservation strategies such as translocation, captive rearing, and hybridization trials to prevent the loss of capacity to adapt to rapidly changing environmental conditions.

Here, I review phylogenetic studies of the lizard family Pygopodidae, a group of 47 extant species that diversified in Australia and New Guinea. The goal of this study was to examine published phylogenetic and phylogenomic hypotheses on pygopodids to identify the strengths and weaknesses in our understanding of their phylogeny. Many parts of the pygopodid family tree are well established by multiple independent tree inferences including: (1) all multispecies genera (i.e., Aprasia, Delma, Lialis, Pletholax, and Pygopus) are monophyletic groups; (2) the root of the pygopodid tree is located along the branch leading to the Delma clade, thus showing that Delma is the sister group to all other pygopodid genera; (3) the Aprasia repens group, Delma tincta group, and several other groups of closely related species are demonstrated to be monophyletic entities; and (4) the monotypic Paradelma orientalis is the sister lineage to the Pygopus clade. Based on accumulated phylogenetic evidence, two taxonomic recommendations are given: Paradelma merits generic status rather than being subsumed into Pygopus as some earlier studies had suggested, and the monotypic Aclys concinna should be recognized as a member of Delma (following current practice) until future studies clarify its placement inside or outside the Delma clade. One chronic problem with phylogenetic studies of pygopodids, which has limited the explanatory power of many tree hypotheses, concerns the undersampling of known species. Although the continual addition of newly described species, especially over the past two decades, has been a major reason for these taxon sampling gaps, deficits in species sampling for ingroups and/or outgroups in several studies of pygopodid species complexes has confounded the testing of some ingroup monophyly hypotheses. Ancient hybridization between non-sister lineages may also be confounding attempts to recover the relationships among pygopodids using molecular data. Indeed, such a phenomenon can explain at least five cases of mito-nuclear discordance and conflicts among trees based on nuclear DNA datasets. Another problem has been the lack of consensus on the relationships among most pygopodid genera, an issue that may stem from rapid diversification of these lineages early in the group's history. Despite current weaknesses in our understanding of pygopodid phylogeny, enough evidence exists to clarify many major and minor structural parts of their family tree. Accordingly, a composite tree for the Pygopodidae was able to be synthesized. This novel tree hypothesis contains all recognized pygopodid species and reveals that about half of the clades are corroborated by multiple independent tree hypotheses, while the remaining clades have less empirical support. [ABSTRACT FROM AUTHOR]

TQMD is a tool for high-performance computing clusters which downloads, stores and produces lists of dereplicated prokaryotic genomes. It has been developed to counter the ever-growing number of prokaryotic genomes and their uneven taxonomic distribution. It is based on word-based alignment-free methods (k-mers), an iterative single-linkage approach and a divide-and-conquer strategy to remain both efficient and scalable. We studied the performance of TQMD by verifying the influence of its parameters and heuristics on the clustering outcome. We further compared TQMD to two other dereplication tools (dRep and Assembly-Dereplicator). Our results showed that TQMD is primarily optimized to dereplicate at higher taxonomic levels (phylum/class), as opposed to the other dereplication tools, but also works at lower taxonomic levels (species/strain) like the other dereplication tools. TQMD is available from source and as a Singularity container at [https://bitbucket.org/phylogeno/tqmd]. [ABSTRACT FROM AUTHOR]

Phylogenomic studies can provide insights into speciation, adaptation, and extinction, while providing a roadmap for conservation. Hawaiian tree snails are a model system for an adaptive radiation facing an extinction crisis. In the last 5 years, nearly all populations of Hawaiian tree snails across the 30 remaining species in the subfamily Achatinellinae (Achatinellidae) have declined from hundreds or thousands in the wild down to undetectable levels. Nearly 100 species historically occurred across dramatic environmental gradients on five of the Hawaiian Islands, but habitat loss, overcollection, and predation by invasive species have decimated populations. As such, this system offers the opportunity to integrate efforts to conserve evolutionary potential into conservation planning for a rapidly declining subfamily. Here, we used genome-wide, restriction-site associated DNA sequencing (RADseq), along with mitochondrial genome reconstruction, to resolve evolutionary relationships to inform conservation efforts. Phylogenetic analysis of nearly 400k genome-wide SNPs from 59 populations and 25 species across six genera in the family Achatinellidae, was generally concordant with taxonomy, geography, and mtDNA with several notable exceptions; mtDNA was unable to resolve some deeper nodes (e.g., the monophyly of Achatinella), while SNP data did not resolve as many shallow nodes. Both phylogenetic and coalescent analysis revealed deep divergences between populations within Achatinella mustelina that were consistent with species-level differences. Given cryptic species-level divergence within populations that are geographically proximate, they are at higher risk of extirpation from invasive predators and climate change than previously assumed. This study clarifies evolutionary relationships within this model system for adaptive radiation, forming the basis for conservation strategies such as translocation, captive rearing, and hybridization trials to prevent the loss of capacity to adapt to rapidly changing environmental conditions. [ABSTRACT FROM AUTHOR]

Background Computing genomic similarity between strains is a prerequisite for genome-based prokaryotic classification and identification. Genomic similarity was first computed as Average Nucleotide Identity (ANI) values based on the alignment of genomic fragments. Since this is computationally expensive, faster and computationally cheaper alignment-free methods have been developed to estimate ANI. However, these methods do not reach the level of accuracy of alignment-based methods. Methods Here we introduce LINflow, a computational pipeline that infers pairwise genomic similarity in a set of genomes. LINflow takes advantage of the speed of the alignment-free sourmash tool to identify the genome in a dataset that is most similar to a query genome and the precision of the alignment-based pyani software to precisely compute ANI between the query genome and the most similar genome identified by sourmash. This is repeated for each new genome that is added to a dataset. The sequentially computed ANI values are stored as Life Identification Numbers (LINs), which are then used to infer all other pairwise ANI values in the set. We tested LINflow on four sets, 484 genomes in total, and compared the needed time and the generated similarity matrices with other tools. Results LINflow is up to 150 times faster than pyani and pairwise ANI values generated by LINflow are highly correlated with those computed by pyani. However, because LINflow infers most pairwise ANI values instead of computing them directly, ANI values occasionally depart from the ANI values computed by pyani. In conclusion, LINflow is a fast and memory-efficient pipeline to infer similarity among a large set of prokaryotic genomes. Its ability to quickly add new genome sequences to an already computed similarity matrix makes LINflow particularly useful for projects when new genome sequences need to be regularly added to an existing dataset.

Despite many bioinformatic solutions for analyzing sequencing data, few options exist for targeted sequence retrieval from whole genomic sequencing (WGS) data with the ultimate goal of generating a phylogeny. Available tools especially struggle at deep phylogenetic levels and necessitate amino-acid space searches, which may increase rates of false positive results. Many tools are also difficult to install and may lack adequate user resources. Here, we describe a program that uses freely available similarity search tools to find homologs in assembled WGS data with unparalleled freedom to modify parameters. We evaluate its performance compared to other commonly used bioinformatics tools on two divergent insect species (>200 My) for which annotated genomes exist, and on one large set each of highly conserved and more variable loci. Our software is capable of retrieving orthologs from well-curated or unannotated, low or high depth shotgun, and target capture assemblies as well or better than other software as assessed by recovering the most genes with maximal coverage and with a low rate of false positives throughout all datasets. When assessing this combination of criteria, ALiBaSeq is frequently the best evaluated tool for gathering the most comprehensive and accurate phylogenetic alignments on all types of data tested. The software (implemented in Python), tutorials, and manual are freely available at https://github.com/AlexKnyshov/alibaseq.

Despite many bioinformatic solutions for analyzing sequencing data, few options exist for targeted sequence retrieval from whole genomic sequencing (WGS) data with the ultimate goal of generating a phylogeny. Available tools especially struggle at deep phylogenetic levels and necessitate amino-acid space searches, which may increase rates of false positive results. Many tools are also difficult to install and may lack adequate user resources. Here, we describe a program that uses freely available similarity search tools to find homologs in assembled WGS data with unparalleled freedom to modify parameters. We evaluate its performance compared to other commonly used bioinformatics tools on two divergent insect species (>200 My) for which annotated genomes exist, and on one large set each of highly conserved and more variable loci. Our software is capable of retrieving orthologs from well-curated or unannotated, low or high depth shotgun, and target capture assemblies as well or better than other software as assessed by recovering the most genes with maximal coverage and with a low rate of false positives throughout all datasets. When assessing this combination of criteria, ALiBaSeq is frequently the best evaluated tool for gathering the most comprehensive and accurate phylogenetic alignments on all types of data tested. The software (implemented in Python), tutorials, and manual are freely available at https://github.com/AlexKnyshov/alibaseq. [ABSTRACT FROM AUTHOR]

Maples (Acer) are among the most diverse and ecologically important tree genera of the north-temperate forests. They include species highly valued as ornamentals and as a source of timber and sugar products. Previous phylogenetic studies employing plastid markers have not provided sufficient resolution, particularly at deeper nodes, leaving the backbone of the maple plastid tree essentially unresolved. We provide the plastid genome sequences of 16 species of maples spanning the sectional diversity of the genus and explore the utility of these sequences as a source of information for genetic and phylogenetic studies in this group. We analyzed the distribution of different types of repeated sequences and the pattern of codon usage, and identified variable regions across the plastome. Maximum likelihood and Bayesian analyses using two partitioning strategies were performed with these and previously published sequences. The plastomes ranged in size from 155,212 to 157,023 bp and had structure and gene content except for Acer palmatum (sect. Palmata), which had longer inverted repeats and an additional copy of the rps19 gene. Two genes, rps2 and rpl22, were found to be truncated at different positions and might be non-functional in several species. Most dispersed repeats, SSRs, and overall variation were detected in the non-coding sequences of the LSC and SSC regions. Fifteen loci, most of which have not been used before in the genus, were identified as the most variable and potentially useful as molecular markers for barcoding and genetic studies. Both ML and Bayesian analyses produced similar results irrespective of the partitioning strategy used. The plastome-based tree largely supported the topology inferred in previous studies using cp markers while providing resolution to the backbone relationships but was highly incongruous with a recently published nuclear tree presenting an opportunity for further research to investigate the causes of discordance, and particularly the role of hybridization in the diversification of the genus. Plastome sequences are valuable tools to resolve deep-level relationships within Acer. The variable loci and SSRs identified in this study will facilitate the development of markers for ecological and evolutionary studies in the genus. This study underscores the potential of plastid genome sequences to improve our understanding of the evolution of maples.

Background With the rapid increase in availability of genomic resources offered by Next-Generation Sequencing (NGS) and the availability of free online genomic databases, efficient and standardized metadata curation approaches have become increasingly critical for the post-processing stages of biological data. Especially in organelle-based studies using circular chloroplast genome datasets, the assembly of the main structural regions in random order and orientation represents a major limitation in our ability to easily generate “ready-to-align” datasets for phylogenetic reconstruction, at both small and large taxonomic scales. In addition, current practices discard the most variable regions of the genomes to facilitate the alignment of the remaining coding regions. Nevertheless, no software is currently available to perform curation to such a degree, through simple detection, organization and positioning of the main plastome regions, making it a time-consuming and error-prone process. Here we introduce a fast and user friendly software ECuADOR, a Perl script specifically designed to automate the detection and reorganization of newly assembled plastomes obtained from any source available (NGS, sanger sequencing or assembler output). Methods ECuADOR uses a sliding-window approach to detect long repeated sequences in draft sequences, which then identifies the inverted repeat regions (IRs), even in case of artifactual breaks or sequencing errors and automates the rearrangement of the sequence to the widely used LSC–Irb–SSC–IRa order. This facilitates rapid post-editing steps such as creation of genome alignments, detection of variable regions, SNP detection and phylogenomic analyses. Results ECuADOR was successfully tested on plant families throughout the angiosperm phylogeny by curating 161 chloroplast datasets. ECuADOR first identified and reordered the central regions (LSC–Irb–SSC–IRa) for each dataset and then produced a new annotation for the chloroplast sequences. The process took less than 20 min with a maximum memory requirement of 150 MB and an accuracy of over 99%. Conclusions ECuADOR is the sole de novo one-step recognition and re-ordination tool that provides facilitation in the post-processing analysis of the extra nuclear genomes from NGS data. The program is available at https://github.com/BiodivGenomic/ECuADOR/.

Congruence among analyses of plant genomic data partitions (nuclear, chloroplast and mitochondrial) is a strong indicator of accuracy in plant molecular phylogenetics. Recent analyses of both nuclear and chloroplast genome data of land plants (embryophytes) have, controversially, been shown to support monophyly of both bryophytes (mosses, liverworts, and hornworts) and tracheophytes (lycopods, ferns, and seed plants), with mosses and liverworts forming the clade Setaphyta. However, relationships inferred from mitochondria are incongruent with these results, and typically indicate paraphyly of bryophytes with liverworts alone resolved as the earliest-branching land plant group. Here, we reconstruct the mitochondrial land plant phylogeny from a newly compiled data set. When among-lineage composition heterogeneity is accounted for in analyses of codon-degenerate nucleotide and amino acid data, the clade Setaphyta is recovered with high support, and hornworts are supported as the earliest-branching lineage of land plants. These new mitochondrial analyses demonstrate partial congruence with current hypotheses based on nuclear and chloroplast genome data, and provide further incentive for revision of how plants arose on land.

“Higher” termites have been able to colonize all tropical and subtropical regions because of their ability to digest lignocellulose with the aid of their prokaryotic gut microbiota. Over the last decade, numerous studies based on 16S rRNA gene amplicon libraries have largely described both the taxonomy and structure of the prokaryotic communities associated with termite guts. Host diet and microenvironmental conditions have emerged as the main factors structuring the microbial assemblages in the different gut compartments. Additionally, these molecular inventories have revealed the existence of termite-specific clusters that indicate coevolutionary processes in numerous prokaryotic lineages. However, for lack of representative isolates, the functional role of most lineages remains unclear. We reconstructed 589 metagenome-assembled genomes (MAGs) from the different gut compartments of eight higher termite species that encompass 17 prokaryotic phyla. By iteratively building genome trees for each clade, we significantly improved the initial automated assignment, frequently up to the genus level. We recovered MAGs from most of the termite-specific clusters in the radiation of, for example, Planctomycetes, Fibrobacteres, Bacteroidetes, Euryarchaeota, Bathyarchaeota, Spirochaetes, Saccharibacteria, and Firmicutes, which to date contained only few or no representative genomes. Moreover, the MAGs included abundant members of the termite gut microbiota. This dataset represents the largest genomic resource for arthropod-associated microorganisms available to date and contributes substantially to populating the tree of life. More importantly, it provides a backbone for studying the metabolic potential of the termite gut microbiota, including the key members involved in carbon and nitrogen biogeochemical cycles, and important clues that may help cultivating representatives of these understudied clades.

Maples (Acer) are among the most diverse and ecologically important tree genera of the north-temperate forests. They include species highly valued as ornamentals and as a source of timber and sugar products. Previous phylogenetic studies employing plastid markers have not provided sufficient resolution, particularly at deeper nodes, leaving the backbone of the maple plastid tree essentially unresolved. We provide the plastid genome sequences of 16 species of maples spanning the sectional diversity of the genus and explore the utility of these sequences as a source of information for genetic and phylogenetic studies in this group. We analyzed the distribution of different types of repeated sequences and the pattern of codon usage, and identified variable regions across the plastome. Maximum likelihood and Bayesian analyses using two partitioning strategies were performed with these and previously published sequences. The plastomes ranged in size from 155,212 to 157,023 bp and had structure and gene content except for Acer palmatum (sect. Palmata), which had longer inverted repeats and an additional copy of the rps19 gene. Two genes, rps2 and rpl22, were found to be truncated at different positions and might be non-functional in several species. Most dispersed repeats, SSRs, and overall variation were detected in the non-coding sequences of the LSC and SSC regions. Fifteen loci, most of which have not been used before in the genus, were identified as the most variable and potentially useful as molecular markers for barcoding and genetic studies. Both ML and Bayesian analyses produced similar results irrespective of the partitioning strategy used. The plastome-based tree largely supported the topology inferred in previous studies using cp markers while providing resolution to the backbone relationships but was highly incongruous with a recently published nuclear tree presenting an opportunity for further research to investigate the causes of discordance, and particularly the role of hybridization in the diversification of the genus. Plastome sequences are valuable tools to resolve deeplevel relationships within Acer. The variable loci and SSRs identified in this study will facilitate the development of markers for ecological and evolutionary studies in the genus. This study underscores the potential of plastid genome sequences to improve our understanding of the evolution of maples. [ABSTRACT FROM AUTHOR]

Background: With the rapid increase in availability of genomic resources offered by Next-Generation Sequencing (NGS) and the availability of free online genomic databases, efficient and standardized metadata curation approaches have become increasingly critical for the post-processing stages of biological data. Especially in organelle-based studies using circular chloroplast genome datasets, the assembly of the main structural regions in random order and orientation represents a major limitation in our ability to easily generate “ready-to-align” datasets for phylogenetic reconstruction, at both small and large taxonomic scales. In addition, current practices discard the most variable regions of the genomes to facilitate the alignment of the remaining coding regions. Nevertheless, no software is currently available to perform curation to such a degree, through simple detection, organization and positioning of the main plastome regions, making it a time-consuming and error-prone process. Here we introduce a fast and user friendly software ECuADOR, a Perl script specifically designed to automate the detection and reorganization of newly assembled plastomes obtained from any source available (NGS, sanger sequencing or assembler output). Methods: ECuADOR uses a sliding-window approach to detect long repeated sequences in draft sequences, which then identifies the inverted repeat regions (IRs), even in case of artifactual breaks or sequencing errors and automates the rearrangement of the sequence to the widely used LSC–Irb–SSC–IRa order. This facilitates rapid post-editing steps such as creation of genome alignments, detection of variable regions, SNP detection and phylogenomic analyses. Results: ECuADOR was successfully tested on plant families throughout the angiosperm phylogeny by curating 161 chloroplast datasets. ECuADOR first identified and reordered the central regions (LSC–Irb–SSC–IRa) for each dataset and then produced a new annotation for the chloroplast sequences. The process took less. than 20 min with a maximum memory requirement of 150 MB and an accuracy of over 99%. Conclusions: ECuADOR is the sole de novo one-step recognition and re-ordination tool that provides facilitation in the post-processing analysis of the extra nuclear genomes from NGS data. The program is available at https://github.com/ BiodivGenomic/ECuADOR/. [ABSTRACT FROM AUTHOR]

Congruence among analyses of plant genomic data partitions (nuclear, chloroplast and mitochondrial) is a strong indicator of accuracy in plant molecular phylogenetics. Recent analyses of both nuclear and chloroplast genome data of land plants (embryophytes) have, controversially, been shown to support monophyly of both bryophytes (mosses, liverworts, and hornworts) and tracheophytes (lycopods, ferns, and seed plants), with mosses and liverworts forming the clade Setaphyta. However, relationships inferred from mitochondria are incongruent with these results, and typically indicate paraphyly of bryophytes with liverworts alone resolved as the earliest-branching land plant group. Here, we reconstruct the mitochondrial land plant phylogeny from a newly compiled data set. When among-lineage composition heterogeneity is accounted for in analyses of codon-degenerate nucleotide and amino acid data, the clade Setaphyta is recovered with high support, and hornworts are supported as the earliest-branching lineage of land plants. These new mitochondrial analyses demonstrate partial congruence with current hypotheses based on nuclear and chloroplast genome data, and provide further incentive for revision of how plants arose on land. [ABSTRACT FROM AUTHOR]

"Higher" termites have been able to colonize all tropical and subtropical regions because of their ability to digest lignocellulose with the aid of their prokaryotic gut microbiota. Over the last decade, numerous studies based on 16S rRNA gene amplicon libraries have largely described both the taxonomy and structure of the prokaryotic communities associated with termite guts. Host diet and microenvironmental conditions have emerged as the main factors structuring the microbial assemblages in the different gut compartments. Additionally, these molecular inventories have revealed the existence of termite-specific clusters that indicate coevolutionary processes in numerous prokaryotic lineages. However, for lack of representative isolates, the functional role of most lineages remains unclear. We reconstructed 589 metagenome-assembled genomes (MAGs) from the different gut compartments of eight higher termite species that encompass 17 prokaryotic phyla. By iteratively building genome trees for each clade, we significantly improved the initial automated assignment, frequently up to the genus level. We recovered MAGs from most of the termite-specific clusters in the radiation of, for example, Planctomycetes, Fibrobacteres, Bacteroidetes, Euryarchaeota, Bathyarchaeota, Spirochaetes, Saccharibacteria, and Firmicutes, which to date contained only few or no representative genomes. Moreover, the MAGs included abundant members of the termite gut microbiota. This dataset represents the largest genomic resource for arthropod-associated microorganisms available to date and contributes substantially to populating the tree of life. More importantly, it provides a backbone for studying the metabolic potential of the termite gut microbiota, including the key members involved in carbon and nitrogen biogeochemical cycles, and important clues that may help cultivating representatives of these understudied clades. [ABSTRACT FROM AUTHOR]

Evolutionary relationships among plants have been inferred primarily using chloroplast data. To date, no study has comprehensively examined the plastome for gene tree conflict. Using a broad sampling of angiosperm plastomes, we characterize gene tree conflict among plastid genes at various time scales and explore correlates to conflict (e.g., evolutionary rate, gene length, molecule type). We uncover notable gene tree conflict against a backdrop of largely uninformative genes. We find alignment length and tree length are strong predictors of concordance, and that nucleotides outperform amino acids. Of the most commonly used markers, matK, greatly outperforms rbcL; however, the rarely used gene rpoC2 is the top-performing gene in every analysis. We find that rpoC2 reconstructs angiosperm phylogeny as well as the entire concatenated set of protein-coding chloroplast genes. Our results suggest that longer genes are superior for phylogeny reconstruction. The alleviation of some conflict through the use of nucleotides suggests that stochastic and systematic error is likely the root of most of the observed conflict, but further research on biological conflict within plastome is warranted given documented cases of heteroplasmic recombination. We suggest that researchers should filter genes for topological concordance when performing downstream comparative analyses on phylogenetic data, even when using chloroplast genomes.

A phylogenomic analysis of 42 complete plastid genomes (plastomes), including 16 that were newly sequenced, was conducted. Plastomes were sampled from 19 subtribes of Pooideae, to investigate relationships within and between Chloroplast Group 1 (Aveneae) and Group 2 (Poeae) species. Two data partitions: complete plastomes, and a combined plastome and rare genomic change (RGC) data matrix, were analyzed. Overall, 156 non-ambiguous RGC were identified, of which homology was inferred for 38 RGC. Among the 38 RGC identified, six were synapomorphic among the Group 1 subtribes: Aveninae, Agrostidinae, and Anthoxanthinae, (Phalaridinae + Torreyochloinae), and 27 were synapomorphic among the Group 2 subtribes: Loliinae, (Ammochloinae + Parapholiinae + Dactylidinae), Parapholiinae, Dactylidinae, Poinae, and Coleanthinae. Four RGC were determined to be homoplasious in Groups 1 and 2. Two other RGC originated through intrastrand deletion events. The remaining RGC events likely originated through recombination given their size and lack of sequence evidence for other types of mutations. This study also determined that relationships between taxa, even those only weakly supported in previous studies, could be inferred with strong support when utilizing complete plastomes.

Background Chlorophyceae is one of three most species-rich green algal classes and also the only class in core Chlorophyta whose monophyly remains uncontested as gene and taxon sampling improves. However, some key relationships within Chlorophyceae are less clear-cut and warrant further investigation. The present study combined genome-scale chloroplast data and rich sampling in an attempt to resolve the ordinal classification in Chlorophyceae. The traditional division into Sphaeropleales and Volvocales (SV), and a clade containing Oedogoniales, Chaetopeltidales, and Chaetophorales (OCC) was of particular interest with the addition of deeply branching members of these groups, as well as the placement of several incertae sedis taxa. Methods We sequenced 18 chloroplast genomes across Chlorophyceae to compile a data set of 58 protein-coding genes of a total of 68 chlorophycean taxa. We analyzed the concatenated nucleotide and amino acid datasets in the Bayesian and Maximum Likelihood frameworks, supplemented by analyses to examine potential discordant signal among genes. We also examined gene presence and absence data across Chlorophyceae. Results Concatenated analyses yielded at least two well-supported phylogenies: nucleotide data supported the traditional classification with the inclusion of the enigmatic Treubarinia into Sphaeropleales sensu lato. However, amino acid data yielded equally strong support for Sphaeropleaceae as sister to Volvocales, with the rest of the taxa traditionally classified in Sphaeropleales in a separate clade, and Treubarinia as sister to all of the above. Single-gene and other supplementary analyses indicated that the data have low phylogenetic signal at these critical nodes. Major clades were supported by genomic structural features such as gene losses and trans-spliced intron insertions in the plastome. Discussion While the sequence and gene order data support the deep split between the SV and OCC lineages, multiple phylogenetic hypotheses are possible for Sphaeropleales s.l. Given this uncertainty as well as the higher-taxonomic disorder seen in other algal groups, dwelling on well-defined, strongly supported Linnaean orders is not currently practical in Chlorophyceae and a less formal clade system may be more useful in the foreseeable future. For example, we identify two strongly and unequivocally supported clades: Treubarinia and Scenedesminia, as well as other smaller groups that could serve a practical purpose as named clades. This system does not preclude future establishment of new orders, or emendment of the current ordinal classification if new data support such conclusions.

The atypoid mygalomorphs include spiders from three described families that build a diverse array of entrance web constructs, including funnel-and-sheet webs, purse webs, trapdoors, turrets and silken collars. Molecular phylogenetic analyses have generally supported the monophyly of Atypoidea, but prior studies have not sampled all relevant taxa. Here we generated a dataset of ultraconserved element loci for all described atypoid genera, including taxa (Mecicobothrium and Hexurella) key to understanding familial monophyly, divergence times, and patterns of entrance web evolution. We show that the conserved regions of the arachnid UCE probe set target exons, such that it should be possible to combine UCE and transcriptome datasets in arachnids. We also show that different UCE probes sometimes target the same protein, and under the matching parameters used here show that UCE alignments sometimes include non-orthologs. Using multiple curated phylogenomic matrices we recover a monophyletic Atypoidea, and reveal that the family Mecicobothriidae comprises four separate and divergent lineages. Fossil-calibrated divergence time analyses suggest ancient Triassic (or older) origins for several relictual atypoid lineages, with late Cretaceous/early Tertiary divergences within some genera indicating a high potential for cryptic species diversity. The ancestral entrance web construct for atypoids, and all mygalomorphs, is reconstructed as a funnel-and-sheet web.

Much genomic-scale, especially transcriptomic, data on spider phylogeny has accumulated in the last few years. These data have recently been used to investigate the diverse architectures and the origin of spider webs, concluding that the ancestral spider spun no foraging web, that spider webs evolved de novo 10–14 times, and that the orb web evolved at least three times. These findings in fact result from a particular phylogenetic character coding strategy, specifically coding the absence of webs as logically equivalent, and homologous to, 10 other observable (i.e., not absent) web architectures. “Absence” of webs should be regarded as inapplicable data. To be analyzed properly by character optimization algorithms, it should be coded as “?” because these codes—or their equivalent—are handled differently by such algorithms. Additional problems include critical misspellings of taxon names from one analysis to the next (misspellings cause some optimization algorithms to drop terminals, which affects taxon sampling and results), and mistakes in spider natural history. In sum, the method causes character optimization algorithms to produce counter-intuitive results, and does not distinguish absence from secondary loss. Proper treatment of missing entries and corrected data instead imply that foraging webs are primitive for spiders and that webs have been lost ∼5–7 times, not gained 10–14 times. The orb web, specifically, may be homologous (originated only once) although lost 2–6 times.

With the increased availability of genome sequences for bacteria, it has become routine practice to construct genome-based phylogenies. These phylogenies have formed the basis for various taxonomic decisions, especially for resolving problematic relationships between taxa. Despite the popularity of concatenating shared genes to obtain well-supported phylogenies, various issues regarding this combined-evidence approach have been raised. These include the introduction of phylogenetic error into datasets, as well as incongruence due to organism-level evolutionary processes, particularly horizontal gene transfer and incomplete lineage sorting. Because of the huge effect that this could have on phylogenies, we evaluated the impact of phylogenetic conflict caused by organism-level evolutionary processes on the established species phylogeny for Pantoea, a member of the Enterobacterales. We explored the presence and distribution of phylogenetic conflict at the gene partition and nucleotide levels, by identifying putative inter-lineage recombination events that might have contributed to such conflict. Furthermore, we determined whether smaller, randomly constructed datasets had sufficient signal to reconstruct the current species tree hypothesis or if they would be overshadowed by phylogenetic incongruence. We found that no individual gene tree was fully congruent with the species phylogeny of Pantoea, although many of the expected nodes were supported by various individual genes across the genome. Evidence of recombination was found across all lineages within Pantoea, and provides support for organism-level evolutionary processes as a potential source of phylogenetic conflict. The phylogenetic signal from at least 70 random genes recovered robust, well-supported phylogenies for the backbone and most species relationships of Pantoea, and was unaffected by phylogenetic conflict within the dataset. Furthermore, despite providing limited resolution among taxa at the level of single gene trees, concatenated analyses of genes that were identified as having no signal resulted in a phylogeny that resembled the species phylogeny of Pantoea. This distribution of signal and noise across the genome presents the ideal situation for phylogenetic inference, as the topology from a ≥70-gene concatenated species phylogeny is not driven by single genes, and our data suggests that this finding may also hold true for smaller datasets. We thus argue that, by using a concatenation-based approach in phylogenomics, one can obtain robust phylogenies due to the synergistic effect of the combined signal obtained from multiple genes.

The reported Agrobacterium radiobacter DSM 30174T genome is highly fragmented, hindering robust comparative genomics and genome-based taxonomic analysis. We re-sequenced the Agrobacterium radiobacter type strain, generating a dramatically improved genome with high contiguity. In addition, we sequenced the genome of Agrobacterium tumefaciens B6T, enabling for the first time, a proper comparative genomics of these contentious Agrobacterium species. We provide concrete evidence that the previously reported Agrobacterium radiobacter type strain genome (Accession Number: ASXY01) is contaminated which explains its abnormally large genome size and fragmented assembly. We propose that Agrobacterium tumefaciens be reclassified as Agrobacterium radiobacter subsp. tumefaciens and that Agrobacterium radiobacter retains it species status with the proposed name of Agrobacterium radiobacter subsp. radiobacter. This proposal is based, first on the high pairwise genome-scale average nucleotide identity supporting the amalgamation of both Agrobacterium radiobacter and Agrobacterium tumefaciens into a single species. Second, maximum likelihood tree construction based on the concatenated alignment of shared genes (core genes) among related strains indicates that Agrobacterium radiobacter NCPPB3001 is sufficiently divergent from Agrobacterium tumefaciens to propose two independent sub-clades. Third, Agrobacterium tumefaciens demonstrates the genomic potential to synthesize the L configuration of fucose in its lipid polysaccharide, fostering its ability to colonize plant cells more effectively than Agrobacterium radiobacter.

The advent of Next Generation Sequencing has reduced sequencing costs and increased genomic projects from a huge amount of organismal taxa, generating an unprecedented amount of genomic datasets publicly available. Often, only a tiny fraction of outstanding relevance of the genomic data produced by researchers is used in their works. This fact allows the data generated to be recycled in further projects worldwide. The assembly of complete mitogenomes is frequently overlooked though it is useful to understand evolutionary relationships among taxa, especially those presenting poor mtDNA sampling at the level of genera and families. This is exactly the case for ants (Hymenoptera:Formicidae) and more specifically for the subfamily Pseudomyrmecinae, a group of arboreal ants with several cases of convergent coevolution without any complete mitochondrial sequence available. In this work, we assembled, annotated and performed comparative genomics analyses of 14 new complete mitochondria from Pseudomyrmecinae species relying solely on public datasets available from the Sequence Read Archive (SRA). We used all complete mitogenomes available for ants to study the gene order conservation and also to generate two phylogenetic trees using both (i) concatenated set of 13 mitochondrial genes and (ii) the whole mitochondrial sequences. Even though the tree topologies diverged subtly from each other (and from previous studies), our results confirm several known relationships and generate new evidences for sister clade classification inside Pseudomyrmecinae clade. We also performed a synteny analysis for Formicidae and identified possible sites in which nucleotidic insertions happened in mitogenomes of pseudomyrmecine ants. Using a data mining/bioinformatics approach, the current work increased the number of complete mitochondrial genomes available for ants from 15 to 29, demonstrating the unique potential of public databases for mitogenomics studies. The wide applications of mitogenomes in research and presence of mitochondrial data in different public dataset types makes the “no budget mitogenomics” approach ideal for comprehensive molecular studies, especially for subsampled taxa.

Evolutionary relationships among plants have been inferred primarily using chloroplast data. To date, no study has comprehensively examined the plastome for gene tree conflict. Using a broad sampling of angiosperm plastomes, we characterize gene tree conflict among plastid genes at various time scales and explore correlates to conflict (e.g., evolutionary rate, gene length, molecule type). We uncover notable gene tree conflict against a backdrop of largely uninformative genes. We find alignment length and tree length are strong predictors of concordance, and that nucleotides outperform amino acids. Of the most commonly used markers, matK, greatly outperforms rbcL; however, the rarely used gene rpoC2 is the top-performing gene in every analysis. We find that rpoC2 reconstructs angiosperm phylogeny as well as the entire concatenated set of protein-coding chloroplast genes. Our results suggest that longer genes are superior for phylogeny reconstruction. The alleviation of some conflict through the use of nucleotides suggests that stochastic and systematic error is likely the root of most of the observed conflict, but further research on biological conflict within plastome is warranted given documented cases of heteroplasmic recombination. We suggest that researchers should filter genes for topological concordance when performing downstream comparative analyses on phylogenetic data, even when using chloroplast genomes. [ABSTRACT FROM AUTHOR]

Arthropods often have obligate relationships with symbiotic microbes, and recent investigations have demonstrated that such host-microbe relationships could be exploited to suppress natural populations of vector carrying mosquitos. Strategies that target the interplay between agricultural pests and their symbionts could decrease the burden caused by agricultural pests; however, the lack of comprehensive genomic insights into naturally occurring microbial symbionts presents a significant bottleneck. Here we employed amplicon surveys, genome-resolved metagenomics, and scanning electron microscopy to investigate symbionts of the wheat stem sawfly (Cephus cinctus), a major pest that causes an estimated $350 million dollars or more in wheat yield losses in the northwestern United States annually. Through 16S rRNA gene sequencing of two major haplotypes and life stages of wheat stem sawfly, we show a novel Spiroplasma species is ever-present and predominant, with phylogenomic analyses placing it as a member of the ixodetis clade of mollicutes. Using state-of-the-art metagenomic assembly and binning strategies we were able to reconstruct a 714 Kb, 72.7%-complete Spiroplasma genome, which represents just the second draft genome from the ixodetis clade of mollicutes. Functional annotation of the Spiroplasma genome indicated carbohydrate-metabolism involved PTS-mediated import of glucose and fructose followed by glycolysis to lactate, acetate, and propionoate. The bacterium also encoded biosynthetic pathways for essential vitamins B2, B3, and B9. We identified putative Spiroplasma virulence genes: cardiolipin and chitinase. These results identify a previously undescribed symbiosis between wheat stem sawfly and a novel Spiroplasma sp., availing insight into their molecular relationship, and may yield new opportunities for microbially-mediated pest control strategies. [ABSTRACT FROM AUTHOR]

Background. Common phylogenomic approaches for recovering phylogenies are often time-consuming and require annotations for orthologous gene relationships that are not always available. In contrast, alignment-free phylogenomic approaches typically use structure and oligomer frequencies to calculate pairwise distances between species. We have developed an approach to quickly calculate distances between species based on codon aversion. Methods. Utilizing a novel alignment-free character state, we present CAM, an alignment-free approach to recover phylogenies by comparing differences in codon aversion motifs (i.e., the set of unused codons within each gene) across all genes within a species. Synonymous codon usage is non-random and differs between organisms, between genes, and even within a single gene, and many genes do not use all possible codons. We report a comprehensive analysis of codon aversion within 229,742,339 genes from 23,428 species across all kingdoms of life, and we provide an alignment-free framework for its use in a phylogenetic construct. For each species, we first construct a set of codon aversion motifs spanning all genes within that species. We define the pairwise distance between two species, A and B, as one minus the number of shared codon aversion motifs divided by the total codon aversion motifs of the species, A or B, containing the fewest motifs. This approach allows us to calculate pairwise distances even when substantial differences in the number of genes or a high rate of divergence between species exists. Finally, we use neighbor-joining to recover phylogenies. Results. Using the Open Tree of Life and NCBI Taxonomy Database as expected phylogenies, our approach compares well, recovering phylogenies that largely match expected trees and are comparable to trees recovered using maximum likelihood and other alignment-free approaches. Our technique is much faster than maximum likelihood and similar in accuracy to other alignment-free approaches. Therefore, we propose that codon aversion be considered a phylogenetically conserved character that may be used in future phylogenomic studies. [ABSTRACT FROM AUTHOR]

A phylogenomic analysis of 42 complete plastid genomes (plastomes), including 16 that were newly sequenced, was conducted. Plastomes were sampled from 19 subtribes of Pooideae, to investigate relationships within and between Chloroplast Group 1 (Aveneae) and Group 2 (Poeae) species. Two data partitions: complete plastomes, and a combined plastome and rare genomic change (RGC) data matrix, were analyzed. Overall, 156 non-ambiguous RGC were identified, of which homology was inferred for 38 RGC. Among the 38 RGC identified, six were synapomorphic among the Group 1 subtribes: Aveninae, Agrostidinae, and Anthoxanthinae, (Phalaridinae + Torreyochloinae), and 27 were synapomorphic among the Group 2 subtribes: Loliinae, (Ammochloinae + Parapholiinae + Dactylidinae), Parapholiinae, Dactylidinae, Poinae, and Coleanthinae. Four RGC were determined to be homoplasious in Groups 1 and 2. Two other RGC originated through intrastrand deletion events. The remaining RGC events likely originated through recombination given their size and lack of sequence evidence for other types of mutations. This study also determined that relationships between taxa, even those only weakly supported in previous studies, could be inferred with strong support when utilizing complete plastomes. [ABSTRACT FROM AUTHOR]

Background. Chlorophyceae is one of three most species-rich green algal classes and also the only class in core Chlorophyta whose monophyly remains uncontested as gene and taxon sampling improves. However, some key relationships within Chlorophyceae are less clear-cut and warrant further investigation. The present study combined genome-scale chloroplast data and rich sampling in an attempt to resolve the ordinal classification in Chlorophyceae. The traditional division into Sphaeropleales and Volvocales (SV), and a clade containing Oedogoniales, Chaetopeltidales, and Chaetophorales (OCC) was of particular interest with the addition of deeply branching members of these groups, as well as the placement of several incertae sedis taxa. Methods. We sequenced 18 chloroplast genomes across Chlorophyceae to compile a data set of 58 protein-coding genes of a total of 68 chlorophycean taxa. We analyzed the concatenated nucleotide and amino acid datasets in the Bayesian and Maximum Likelihood frameworks, supplemented by analyses to examine potential discordant signal among genes. We also examined gene presence and absence data across Chlorophyceae. Results. Concatenated analyses yielded at least two well-supported phylogenies: nucleotide data supported the traditional classification with the inclusion of the enigmatic Treubarinia into Sphaeroplealessensu lato. However, amino acid data yielded equally strong support for Sphaeropleaceae as sister to Volvocales, with the rest of the taxa traditionally classified in Sphaeropleales in a separate clade, and Treubarinia as sister to all of the above. Single-gene and other supplementary analyses indicated that the data have low phylogenetic signal at these critical nodes. Major clades were supported by genomic structural features such as gene losses and trans-spliced intron insertions in the plastome. Discussion. While the sequence and gene order data support the deep split between the SV and OCC lineages, multiple phylogenetic hypotheses are possible for Sphaeropleales s.l. Given this uncertainty as well as the higher-taxonomic disorder seen in other algal groups, dwelling on well-defined, strongly supported Linnaean orders is not currently practical in Chlorophyceae and a less formal clade system may be more useful in the foreseeable future. For example, we identify two strongly and unequivocally supported clades: Treubarinia and Scenedesminia, as well as other smaller groups that could serve a practical purpose as named clades. This system does not preclude future establishment of new orders, or emendment of the current ordinal classification if new data support such conclusions. [ABSTRACT FROM AUTHOR]