Phylogenetic utility and evolution of indels: a study in neognathous birds

Performance of intron 7 of the β-fibrinogen gene for phylogenetic analysis: An example using gladiator frogs, Boana Gray, 1825 (Anura, Hylidae, Cophomantinae)

Boana, the third largest genus of Hylinae, has cryptic morphological species. The potential applicability of b-fibrinogen intron 7 – FGBI7 is explored to propose a robust phylogeny of Boana. The phylogenetic potential of FGBI7 was evaluated using maximum parsimony, MrBayes, and maximum likelihood analysis. Comparison of polymorphic sites and topologies obtained with concatenated analysis of FGBI7 and other nuclear genes (CXCR4, CXCR4, RHO, SIAH1, TYR, and 28S) allowed evaluation of the phylogenetic signal of FGBI7. Mean evolutionary rates were calculated using the sequences of the mitochondrial genes ND1 and CYTB available for Boana in GenBank. Dating of Boana and some of its groups was performed using the RelTime method with secondary calibration. FGBI7 analysis revealed high values at informative sites for parsimony. The absolute values of the mean evolutionary rate were higher for mitochondrial genes than for FGBI7. Dating of congruent Boana groups for ND1, CYTB, and FGBI7 revealed closer values between mitochondrial genes and slightly different values from those of FGBI7. Divergence times of basal groups tended to be overestimated when mtDNA was used and were more accurate when nDNA was used. Although there is evidence of phylogenetic potential arising from concatenation of specific genes, FGBI7 provides well-resolved independent gene trees. These results lead to a paradigm for linking data in phylogenomics that focuses on the uniqueness of species histories and ignores the multiplicities of individual gene histories.

NGS-Indel Coder v2.0: A Streamlined Pipeline to Code Indel Characters in Phylogenomic Data

Hypothesized evolutionary insertions and deletions in nucleic acid sequences (indels) contain significant phylogenetic information and can be integrated in phylogenomic analyses. However, assemblies of short reads obtained from next-generation sequencing (NGS) technologies can contain errors that result in falsely inferred indels that need to be detected and omitted to avoid inclusion in phylogenetic analysis. Here, we detail the commands that comprise a new version of the NGS-Indel Coder pipeline, which was developed to validate indels using assembly read depth.

A 51 bp indel polymorphism within the PTH1R gene is significantly associated with chicken growth and carcass traits

Parathyroid hormone (PTH) is a crucial regulator of calcium homeostasis and bone remodeling, and the parathyroid hormone 1 receptor (PTH1R) belongs to a class II G-protein-coupled receptor. PTH activates PTH1R, which mediates catabolic and anabolic processes in the skeleton. However, the functional mechanism of PTH1R has not been thoroughly elucidated in organisms. This study identified a 51 bp indel mutation in the first intron of the PTH1R gene and elucidated the effect of this gene mutation on the growth and carcass traits in chickens. The results indicated that the 51 bp indel was significantly associated with subcutaneous fat thickness, abdominal fat weight, body weight and daily gain over 4-8 weeks. Furthermore, we found that PTH1R gene expression was highest in the kidney and liver tissues, and it showed a trend of decreasing in leg and breast muscle tissues at different embryonic stages. In addition, we examined the expression of the three genotypes of the PTH1R gene in the liver, breast muscle and abdominal fat and found that the II genotype was significantly higher than the DD and ID genotypes. In summary, these findings suggest that the PTH1R gene can serve as a potential molecular marker for chicken breeding.

Molecular characterization and an 80-bp indel polymorphism within the prolactin receptor (PRLR) gene and its associations with chicken growth and carcass traits

The prolactin receptor (PRLR), a type I cytokine receptor, must bind prolactin (PRL) to act on target cells to mediate various physiological functions, including reproduction and lactation. This study identified an 80-bp insertion/deletion (indel) polymorphism in the 3'-untranslated region (3'-UTR) of the chicken PRLR gene in 3736 individuals from 15 breeds and analyzed its associations with growth and carcass traits in an F₂ resource population. The results of the association analysis indicated that the 80-bp indel polymorphism was significantly (P < 0.05) or very significantly (P < 0.01) associated with multiple growth and carcass traits, such as body weight, leg weight, and shank length. In addition, we found that during the breeding process of commercial laying hens and commercial broilers, the 80-bp indel locus was artificially selected for the II genotype. Together, our findings reveal that this 80-bp indel polymorphism has potential as a new molecular marker for marker-assisted selection of chicken growth and carcass traits.

Phylogenetic Signal of Indels and the Neoavian Radiation

The early radiation of Neoaves has been hypothesized to be an intractable “hard polytomy”. We explore the fundamental properties of insertion/deletion alleles (indels), an under-utilized form of genomic data with the potential to help solve this. We scored >5 million indels from >7000 pan-genomic intronic and ultraconserved element (UCE) loci in 48 representatives of all neoavian orders. We found that intronic and UCE indels exhibited less homoplasy than nucleotide (nt) data. Gene trees estimated using indel data were less resolved than those estimated using nt data. Nevertheless, Accurate Species TRee Algorithm (ASTRAL) species trees estimated using indels were generally similar to nt-based ASTRAL trees, albeit with lower support. However, the power of indel gene trees became clear when we combined them with nt gene trees, including a striking result for UCEs. The individual UCE indel and nt ASTRAL trees were incongruent with each other and with the intron ASTRAL trees; however, the combined indel+nt ASTRAL tree was much more congruent with the intronic trees. Finally, combining indel and nt data for both introns and UCEs provided sufficient power to reduce the scope of the polytomy that was previously proposed for several supraordinal lineages of Neoaves.

NGS-Indel Coder: A pipeline to code indel characters in phylogenomic data with an example of its application in milkweeds (Asclepias)

Targeted genome sequencing approaches allow characterization of evolutionary relationships using a considerable number of nuclear genes and informative characters. However, most phylogenomic analyses only utilize single nucleotide polymorphisms (SNPs). Studies at the species level, especially in groups that have recently radiated, often recover low amounts of phylogenetically informative variation in coding regions, and require non-coding sequences, which are richer in indels, to resolve gene trees. Here, NGS-Indel Coder, a pipeline to detect and omit false positive indels inferred from assemblies of short read sequence data, was developed to resolve the relationships among and within major clades of the American milkweeds (Asclepias), which are the result of a rapid and recent evolutionary radiation, and whose phylogeny has been difficult to resolve. This pipeline was applied to a Hyb-Seq data set of 768 loci including targeted exons and flanking intron regions from 33 milkweed species. Robust species tree inference was improved by excluding small alignment partitions (<100 bp) that increased gene tree ambiguity and incongruence. To further investigate the robustness of indel coding, data sets that included small and large indels were explored, and species trees derived from concatenated loci versus coalescent methods based on gene trees were compared. The phylogeny of Asclepias obtained using nuclear data was well resolved, and phylogenetic information from indels improved resolution of specific nodes. The Temperate North American, Mexican Highland, and Incarnatae clades were well supported as monophyletic. Asclepias coulteri, which has been considered part of the Sonoran Desert clade based on plastome analyses, was placed as sister to all the other milkweed species studied here, rather than as a member of that clade. Two groups within the Temperate North American and Mexican clades were not resolved, and the inferred relationships strongly conflicted when comparing results based on data sets that did or did not include indel characters. This new pipeline represents a step forward in making maximal use of the information content in phylogenomic data sets.

Split-inducing indels in phylogenomic analysis

Background

Most phylogenetic studies using molecular data treat gaps in multiple sequence alignments as missing data or even completely exclude alignment columns that contain gaps.

Results

Here we show that gap patterns in large-scale, genome-wide alignments are themselves phylogenetically informative and can be used to infer reliable phylogenies provided the gap data are properly filtered to reduce noise introduced by the alignment method. We introduce here the notion of split-inducing indels (splids) that define an approximate bipartition of the taxon set. We show both in simulated data and in case studies on real-life data that splids can be efficiently extracted from phylogenomic data sets.

Conclusions

Suitably processed gap patterns extracted from genome-wide alignment provide a surprisingly clear phylogenetic signal and an allow the inference of accurate phylogenetic trees.

Electronic supplementary material

The online version of this article (10.1186/s13015-018-0130-7) contains supplementary material, which is available to authorized users.

Identification of the fungal endophyte of Ammophila breviligulata (American beachgrass) as Epichloë amarillans

The grass Ammophila breviligulata (American beachgrass) is known to host an endophyte of the genus Epichloë. Based on morphological characteristics it was originally identified as Acremonium typhinum var. ammophilae and is currently designated as Epichloë typhina var. ammophilae. However, the Epichloë species has not previously been identified based on DNA sequence data. Based on phylogenetic placement of beta-tubulin and translation elongation factor 1-alpha DNA sequences the endophyte is identified as a member of E. amarillans rather than E. typhina.

Never judge an iguana by its spines: Systematics of the Yucatan spiny tailed iguana, Ctenosaura defensor (Cope, 1866)

Spiny tailed iguanas are highly diverse clade of lizards in Mesoamerica, ranging from northern Mexico through Panama. Utilizing 2 regions of mitochondrial DNA (1948bp) and 4 nuclear loci (2232bp) we explored the relationships between these species and the phylogeographic history of the major clades. We discovered that the lineage endemic to the Yucatan Peninsula renders the genus Ctenosaura paraphyletic. To resolve this non-monophyly, we resurrect the taxon Cachryx Cope, 1866, and provide a new diagnosis for the genus. We also find that small body-size and highly spinose tails in the species previously referred to the subgenus Enyaliosaurus, have evolved independently 3 times. Cachryx were recovered as sister to the lineage of iguanines endemic to the Galapagos Islands, and we discuss biogeographic scenarios to explain this relationship as well as those among the primary clades of Ctenosaura in Mesoamerica.

Chicken rRNA Gene Cluster Structure

Ribosomal RNA (rRNA) genes, whose activity results in nucleolus formation, constitute an extremely important part of genome. Despite the extensive exploration into avian genomes, no complete description of avian rRNA gene primary structure has been offered so far. We publish a complete chicken rRNA gene cluster sequence here, including 5'ETS (1836 bp), 18S rRNA gene (1823 bp), ITS1 (2530 bp), 5.8S rRNA gene (157 bp), ITS2 (733 bp), 28S rRNA gene (4441 bp) and 3'ETS (343 bp). The rRNA gene cluster sequence of 11863 bp was assembled from raw reads and deposited to GenBank under KT445934 accession number. The assembly was validated through in situ fluorescent hybridization analysis on chicken metaphase chromosomes using computed and synthesized specific probes, as well as through the reference assembly against de novo assembled rRNA gene cluster sequence using sequenced fragments of BAC-clone containing chicken NOR (nucleolus organizer region). The results have confirmed the chicken rRNA gene cluster validity.

Parsimony and model-based analyses of indels in avian nuclear genes reveal congruent and incongruent phylogenetic signals

Insertion/deletion (indel) mutations, which are represented by gaps in multiple sequence alignments, have been used to examine phylogenetic hypotheses for some time. However, most analyses combine gap data with the nucleotide sequences in which they are embedded, probably because most phylogenetic datasets include few gap characters. Here, we report analyses of 12,030 gap characters from an alignment of avian nuclear genes using maximum parsimony (MP) and a simple maximum likelihood (ML) framework. Both trees were similar, and they exhibited almost all of the strongly supported relationships in the nucleotide tree, although neither gap tree supported many relationships that have proven difficult to recover in previous studies. Moreover, independent lines of evidence typically corroborated the nucleotide topology instead of the gap topology when they disagreed, although the number of conflicting nodes with high bootstrap support was limited. Filtering to remove short indels did not substantially reduce homoplasy or reduce conflict. Combined analyses of nucleotides and gaps resulted in the nucleotide topology, but with increased support, suggesting that gap data may prove most useful when analyzed in combination with nucleotide substitutions.

A phylogeny of birds based on over 1,500 loci collected by target enrichment and high-throughput sequencing

Evolutionary relationships among birds in Neoaves, the clade comprising the vast majority of avian diversity, have vexed systematists due to the ancient, rapid radiation of numerous lineages. We applied a new phylogenomic approach to resolve relationships in Neoaves using target enrichment (sequence capture) and high-throughput sequencing of ultraconserved elements (UCEs) in avian genomes. We collected sequence data from UCE loci for 32 members of Neoaves and one outgroup (chicken) and analyzed data sets that differed in their amount of missing data. An alignment of 1,541 loci that allowed missing data was 87% complete and resulted in a highly resolved phylogeny with broad agreement between the Bayesian and maximum-likelihood (ML) trees. Although results from the 100% complete matrix of 416 UCE loci were similar, the Bayesian and ML trees differed to a greater extent in this analysis, suggesting that increasing from 416 to 1,541 loci led to increased stability and resolution of the tree. Novel results of our study include surprisingly close relationships between phenotypically divergent bird families, such as tropicbirds (Phaethontidae) and the sunbittern (Eurypygidae) as well as between bustards (Otididae) and turacos (Musophagidae). This phylogeny bolsters support for monophyletic waterbird and landbird clades and also strongly supports controversial results from previous studies, including the sister relationship between passerines and parrots and the non-monophyly of raptorial birds in the hawk and falcon families. Although significant challenges remain to fully resolving some of the deep relationships in Neoaves, especially among lineages outside the waterbirds and landbirds, this study suggests that increased data will yield an increasingly resolved avian phylogeny.