Complexity of avian evolution revealed by family-level genomes.

Despite tremendous efforts in the past decades, relationships among main avian lineages remain heavily debated without a clear resolution. Discrepancies have been attributed to diversity of species sampled, phylogenetic method and the choice of genomic regions1-3. Here we address these issues by analysing the genomes of 363 bird species4 (218 taxonomic families, 92% of total). Using intergenic regions and coalescent methods, we present a well-supported tree but also a marked degree of discordance. The tree confirms that Neoaves experienced rapid radiation at or near the Cretaceous-Palaeogene boundary. Sufficient loci rather than extensive taxon sampling were more effective in resolving difficult nodes. Remaining recalcitrant nodes involve species that are a challenge to model due to either extreme DNA composition, variable substitution rates, incomplete lineage sorting or complex evolutionary events such as ancient hybridization. Assessment of the effects of different genomic partitions showed high heterogeneity across the genome. We discovered sharp increases in effective population size, substitution rates and relative brain size following the Cretaceous-Palaeogene extinction event, supporting the hypothesis that emerging ecological opportunities catalysed the diversification of modern birds. The resulting phylogenetic estimate offers fresh insights into the rapid radiation of modern birds and provides a taxon-rich backbone tree for future comparative studies.

A region of suppressed recombination misleads neoavian phylogenomics.

Genomes are typically mosaics of regions with different evolutionary histories. When speciation events are closely spaced in time, recombination makes the regions sharing the same history small, and the evolutionary history changes rapidly as we move along the genome. When examining rapid radiations such as the early diversification of Neoaves 66 Mya, typically no consistent history is observed across segments exceeding kilobases of the genome. Here, we report an exception. We found that a 21-Mb region in avian genomes, mapped to chicken chromosome 4, shows an extremely strong and discordance-free signal for a history different from that of the inferred species tree. Such a strong discordance-free signal, indicative of suppressed recombination across many millions of base pairs, is not observed elsewhere in the genome for any deep avian relationships. Although long regions with suppressed recombination have been documented in recently diverged species, our results pertain to relationships dating circa 65 Mya. We provide evidence that this strong signal may be due to an ancient rearrangement that blocked recombination and remained polymorphic for several million years prior to fixation. We show that the presence of this region has misled previous phylogenomic efforts with lower taxon sampling, showing the interplay between taxon and locus sampling. We predict that similar ancient rearrangements may confound phylogenetic analyses in other clades, pointing to a need for new analytical models that incorporate the possibility of such events.

Exploring Conflicts in Whole Genome Phylogenetics: A Case Study Within Manakins (Aves: Pipridae).

Some phylogenetic problems remain unresolved even when large amounts of sequence data are analyzed and methods that accommodate processes such as incomplete lineage sorting are employed. In addition to investigating biological sources of phylogenetic incongruence, it is also important to reduce noise in the phylogenomic dataset by using appropriate filtering approach that addresses gene tree estimation errors. We present the results of a case study in manakins, focusing on the very difficult clade comprising the genera Antilophia and Chiroxiphia. Previous studies suggest that Antilophia is nested within Chiroxiphia, though relationships among Antilophia+Chiroxiphia species have been highly unstable. We extracted more than 11,000 loci (ultra-conserved elements and introns) from whole genomes and conducted analyses using concatenation and multispecies coalescent methods. Topologies resulting from analyses using all loci differed depending on the data type and analytical method, with 2 clades (Antilophia+Chiroxiphia and Manacus+Pipra+Machaeopterus) in the manakin tree showing incongruent results. We hypothesized that gene trees that conflicted with a long coalescent branch (e.g., the branch uniting Antilophia+Chiroxiphia) might be enriched for cases of gene tree estimation error, so we conducted analyses that either constrained those gene trees to include monophyly of Antilophia+Chiroxiphia or excluded these loci. While constraining trees reduced some incongruence, excluding the trees led to completely congruent species trees, regardless of the data type or model of sequence evolution used. We found that a suite of gene metrics (most importantly the number of informative sites and likelihood of intralocus recombination) collectively explained the loci that resulted in non-monophyly of Antilophia+Chiroxiphia. We also found evidence for introgression that may have contributed to the discordant topologies we observe in Antilophia+Chiroxiphia and led to deviations from expectations given the multispecies coalescent model. Our study highlights the importance of identifying factors that can obscure phylogenetic signal when dealing with recalcitrant phylogenetic problems, such as gene tree estimation error, incomplete lineage sorting, and reticulation events. [Birds; c-gene; data type; gene estimation error; model fit; multispecies coalescent; phylogenomics; reticulation].

Variables associated with nonresponders to high-frequency (10 kHz) spinal cord stimulation.

INTRODUCTION: The use of spinal cord stimulation (SCS) therapy to treat chronic pain continues to rise. Optimal patient selection remains one of the most important factors for SCS success. However, despite increased utilization and the existence of general indications, predicting which patients will benefit from neuromodulation remains one of the main challenges for this therapy. Therefore, this study aims to identify the variables that may correlate with nonresponders to high-frequency (10 kHz) SCS to distinguish the subset of patients less likely to benefit from this intervention. MATERIALS AND METHODS: This was a retrospective single-center observational study of patients who underwent 10 kHz SCS implant. Patients were divided into nonresponders and responders groups. Demographic data and clinical outcomes were collected at baseline and statistical analysis was performed for all continuous and categorical variables between the two groups to calculate statistically significant differences. RESULTS: The study population comprised of 237 patients, of which 67.51% were responders and 32.49% were nonresponders. There was a statistically significant difference of high levels of kinesiophobia, high self-perceived disability, greater pain intensity, and clinically relevant pain catastrophizing at baseline in the nonresponders compared to the responders. A few variables deemed potentially relevant, such as age, gender, history of spinal surgery, diabetes, alcohol use, tobacco use, psychiatric illness, and opioid utilization at baseline were not statistically significant. CONCLUSION: Our study is the first in the neuromodulation literature to raise awareness to the association of high levels of kinesiophobia preoperatively in nonresponders to 10 kHz SCS therapy. We also found statistically significant differences with greater pain intensity, higher self-perceived disability, and clinically relevant pain catastrophizing at baseline in the nonresponders relative to responders. It may be appropriate to screen for these factors preoperatively to identify patients who are less likely to respond to SCS. If these modifiable risk factors are present, it might be prudent to consider a pre-rehabilitation program with pain neuroscience education to address these factors prior to SCS therapy, to enhance successful outcomes in neuromodulation.

Systems biology as a framework to understand the physiological and endocrine bases of behavior and its evolution-From concepts to a case study in birds.

Organismal behavior, with its tremendous complexity and diversity, is generated by numerous physiological systems acting in coordination. Understanding how these systems evolve to support differences in behavior within and among species is a longstanding goal in biology that has captured the imagination of researchers who work on a multitude of taxa, including humans. Of particular importance are the physiological determinants of behavioral evolution, which are sometimes overlooked because we lack a robust conceptual framework to study mechanisms underlying adaptation and diversification of behavior. Here, we discuss a framework for such an analysis that applies a "systems view" to our understanding of behavioral control. This approach involves linking separate models that consider behavior and physiology as their own networks into a singular vertically integrated behavioral control system. In doing so, hormones commonly stand out as the links, or edges, among nodes within this system. To ground our discussion, we focus on studies of manakins (Pipridae), a family of Neotropical birds. These species have numerous physiological and endocrine specializations that support their elaborate reproductive displays. As a result, manakins provide a useful example to help imagine and visualize the way systems concepts can inform our appreciation of behavioral evolution. In particular, manakins help clarify how connectedness among physiological systems-which is maintained through endocrine signaling-potentiate and/or constrain the evolution of complex behavior to yield behavioral differences across taxa. Ultimately, we hope this review will continue to stimulate thought, discussion, and the emergence of research focused on integrated phenotypes in behavioral ecology and endocrinology.

A Prospective, Randomized Single-Blind Crossover Study Comparing High-Frequency 10,000 Hz and Burst Spinal Cord Stimulation.

OBJECTIVES: Although both high-frequency and burst spinal cord stimulation (SCS) have shown improved efficacy and patient satisfaction compared with conventional tonic stimulation, there are limited data directly comparing the two. This study aimed to compare both high-frequency 10,000 Hz and burst SCS in the same patients in terms of pain relief and satisfaction in those with axial back pain with or without leg pain. MATERIALS AND METHODS: This prospective, single-blind, randomized controlled trial was conducted at an outpatient pain clinic within an academic medical center. Participants were randomly allocated to one of two groups in which they trialed either burst or high-frequency 10,000 Hz SCS over five days, followed by a 24- to 48-hour washout period with no stimulation, and the alternative therapy over the remaining four days. Visual analog scale (VAS) scores were collected immediately before and after both therapy trials. Secondary end points included percentage change in VAS score and patient preference. RESULTS: Of 25 participants, those receiving burst followed by high-frequency SCS (n = 11) had a mean VAS difference of 4.73 after the first trial period and 2.86 after the second. Of those receiving high-frequency followed by burst SCS (n = 14), mean VAS difference after the first trial period was 4.00 and 1.93 after the second trial period. Four participants were withdrawn owing to lead migration. Both therapies showed statistically significant differences in pre- minus post-VAS scores and percent relief. There were no significant differences in carryover or treatment effects between the two groups. There was a statistically significant association between trial sequence and stimulator type implanted because the first stimulator trialed was more likely to be chosen. CONCLUSIONS: There were no observed differences in VAS pain score decrease when comparing burst and high-frequency 10,000 Hz SCS programming therapies. Patient preference followed an order effect, favoring the first programming therapy in the trial sequence.

Treatment of Refractory Low Back Pain Using Passive Recharge Burst in Patients Without Options for Corrective Surgery: Findings and Results From the DISTINCT Study, a Prospective Randomized Multicenter Controlled Trial.

OBJECTIVE: Spinal cord stimulation (SCS) is effective for relieving chronic intractable pain conditions. The Dorsal spInal cord STImulatioN vs mediCal management for the Treatment of low back pain study evaluates the effectiveness of SCS compared with conventional medical management (CMM) in the treatment of chronic low back pain in patients who had not undergone and were not candidates for lumbar spine surgery. METHODS AND MATERIALS: Patients were randomized to passive recharge burst therapy (n = 162) or CMM (n = 107). They reported severe pain and disability for more than a decade and had failed a multitude of therapies. Common diagnoses included degenerative disc disease, spondylosis, stenosis, and scoliosis-yet not to a degree amenable to surgery. The six-month primary end point compared responder rates, defined by a 50% reduction in pain. Hierarchical analyses of seven secondary end points were performed in the following order: composite responder rate (numerical rating scale [NRS] or Oswestry Disability Index [ODI]), NRS, ODI, Pain Catastrophizing Scale responder rate, Patient Global Impression of Change (PGIC) responder rate, and Patient-Reported Outcome Measure Information System-29 in pain interference and physical function. RESULTS: Intention-to-treat analysis showed a significant difference in pain responders on NRS between SCS (72.6%) and CMM (7.1%) arms (p < 0.0001). Of note, 85.2% of those who received six months of therapy responded on NRS compared with 6.2% of those with CMM (p < 0.0001). All secondary end points indicated the superiority of burst therapy over CMM. A composite measure on function or pain relief showed 91% of subjects with SCS improved, compared with 16% of subjects with CMM. A substantial improvement of 30 points was observed on ODI compared with a

Categorical edge-based analyses of phylogenomic data reveal conflicting signals for difficult relationships in the avian tree.

Phylogenetic analyses fail to yield a satisfactory resolution of some relationships in the tree of life even with genome-scale datasets, so the failure is unlikely to reflect limitations in the amount of data. Gene tree conflicts are particularly notable in studies focused on these contentious nodes, and taxon sampling, different analytical methods, and/or data type effects can further confound analyses. Although many efforts have been made to incorporate biological conflicts, few studies have curated individual genes for their efficiency in phylogenomic studies. Here, we conduct an edge-based analysis of Neoavian evolution, examining the phylogenetic efficacy of two recent phylogenomic bird datasets and three datatypes (ultraconserved elements [UCEs], introns, and coding regions). We assess the potential causes for biases in signal-resolution for three difficult nodes: the earliest divergence of Neoaves, the position of the enigmatic Hoatzin (Opisthocomus hoazin), and the position of owls (Strigiformes). We observed extensive conflict among genes for all data types and datasets even after meticulous curation. Edge-based analyses (EBA) increased congruence and provided information about the impact of data type, GC content variation (GCCV), and outlier genes on each of nodes we examined. First, outlier gene signals appeared to drive different patterns of support for the relationships among the earliest diverging Neoaves. Second, the placement of Hoatzin was highly variable, although our EBA did reveal a previously unappreciated data type effect with an impact on its position. It also revealed that the resolution with the most support here was Hoatzin + shorebirds. Finally, GCCV, rather than data type (i.e., coding vs non-coding) per se, was correlated with a signal that supports monophyly of owls + Accipitriformes (hawks, eagles, and vultures). Eliminating high GCCV loci increased the signal for owls + mousebirds. Categorical EBA was able to reveal the nature of each edge and provide a way to highlight especially problematic branches that warrant a further examination. The current study increases our understanding about the contentious parts of the avian tree, which show even greater conflicts than appreciated previously.

Historical specimens and the limits of subspecies phylogenomics in the New World quails (Odontophoridae).

As phylogenomics focuses on comprehensive taxon sampling at the species and population/subspecies levels, incorporating genomic data from historical specimens has become increasingly common. While historical samples can fill critical gaps in our understanding of the evolutionary history of diverse groups, they also introduce additional sources of phylogenomic uncertainty, making it difficult to discern novel evolutionary relationships from artifacts caused by sample quality issues. These problems highlight the need for improved strategies to disentangle artifactual patterns from true biological signal as historical specimens become more prevalent in phylogenomic datasets. Here, we tested the limits of historical specimen-driven phylogenomics to resolve subspecies-level relationships within a highly polytypic family, the New World quails (Odontophoridae), using thousands of ultraconserved elements (UCEs). We found that relationships at and above the species-level were well-resolved and highly supported across all analyses, with the exception of discordant relationships within the two most polytypic genera which included many historical specimens. We examined the causes of discordance and found that inferring phylogenies from subsets of taxa resolved the disagreements, suggesting that analyzing subclades can help remove artifactual causes of discordance in datasets that include historical samples. At the subspecies-level, we found well-resolved geographic structure within the two most polytypic genera, including the most polytypic species in this family, Northern Bobwhites (Colinus virginianus), demonstrating that variable sites within UCEs are capable of resolving phylogenetic structure below the species level. Our results highlight the importance of complete taxonomic sampling for resolving relationships among polytypic species, often through the inclusion of historical specimens, and we propose an integrative strategy for understanding and addressing the uncertainty that historical samples sometimes introduce to phylogenetic analyses.

Earth history and the passerine superradiation.

Avian diversification has been influenced by global climate change, plate tectonic movements, and mass extinction events. However, the impact of these factors on the diversification of the hyperdiverse perching birds (passerines) is unclear because family level relationships are unresolved and the timing of splitting events among lineages is uncertain. We analyzed DNA data from 4,060 nuclear loci and 137 passerine families using concatenation and coalescent approaches to infer a comprehensive phylogenetic hypothesis that clarifies relationships among all passerine families. Then, we calibrated this phylogeny using 13 fossils to examine the effects of different events in Earth history on the timing and rate of passerine diversification. Our analyses reconcile passerine diversification with the fossil and geological records; suggest that passerines originated on the Australian landmass ∼47 Ma; and show that subsequent dispersal and diversification of passerines was affected by a number of climatological and geological events, such as Oligocene glaciation and inundation of the New Zealand landmass. Although passerine diversification rates fluctuated throughout the Cenozoic, we find no link between the rate of passerine diversification and Cenozoic global temperature, and our analyses show that the increases in passerine diversification rate we observe are disconnected from the colonization of new continents. Taken together, these results suggest more complex mechanisms than temperature change or ecological opportunity have controlled macroscale patterns of passerine speciation.

Antitumor necrosis factor alpha and an epidural abscess during a spinal cord stimulator trial: A case report.

Spinal cord stimulation (SCS) is commonly utilized for treatment and management of chronic intractable low back and lower extremity pain. Although SCS is an overall low-risk procedure, there are potential life-threatening complications, including surgical site infections, such as an epidural abscess. Immunosuppression, a risk factor for epidural abscess, is becoming more common as an increasing number of patients are being treated with biologics for a multitude of autoimmune disorders. One class of commonly utilized biologics is antitumor necrosis factor (anti-TNF) alpha. Whereas these drugs can provide tremendous benefit for treatment and management of autoimmune disorders, there is no clear understanding of the degree to which these medications increase a patient's risk for surgical site infection, including those associated with SCS-related procedures. We present a case of an epidural abscess that developed immediately following an SCS trial in a patient with multiple undisclosed risk factors, including the use of an anti-TNF alpha agent to treat ankylosing spondylitis. For an epidural abscess, early diagnosis is key to preventing devastating complications and the need for surgical intervention. Immunosuppression can be the result of multiple issues including cancer, HIV, and biologic agents, such as anti-TNF alpha for the management of autoimmune diseases. There is limited evidence pertaining to the development of epidural abscesses in patients on anti-TNF alpha medications who undergo SCS. Studies focused on infections in patients undergoing SCS trials and permanent implants while on anti-TNF alpha agents could provide recommendations and guidance.

A phylogenomic supermatrix of Galliformes (Landfowl) reveals biased branch lengths.

Building taxon-rich phylogenies is foundational for macroevolutionary studies. One approach to improve taxon sampling beyond individual studies is to build supermatricies of publicly available data, incorporating taxa sampled across different studies and utilizing different loci. Most existing supermatrix studies have focused on loci commonly sequenced with Sanger technology ("legacy" markers, such as mitochondrial data and small numbers of nuclear loci). However, incorporating phylogenomic studies into supermatrices allows problem nodes to be targeted and resolved with considerable amounts of data, while improving taxon sampling with legacy data. Here we estimate phylogeny from a galliform supermatrix which includes well-known model and agricultural species such as the chicken and turkey. We assembled a supermatrix comprising 4500 ultra-conserved elements (UCEs) collected as part of recent phylogenomic studies in this group and legacy mitochondrial and nuclear (intron and exon) sequences. Our resulting phylogeny included 88% of extant species and recovered well-accepted relationships with strong support. However, branch lengths, which are particularly important in down-stream macroevolutionary studies, appeared vastly skewed. Taxa represented only by rapidly evolving mitochondrial data had high proportions of missing data and exhibited long terminal branches. Conversely, taxa sampled for slowly evolving UCEs with low proportions of missing data exhibited substantially shorter terminal branches. We explored several branch length re-estimation methods with particular attention to terminal branches and conclude that re-estimation using well-sampled mitochondrial sequences may be a pragmatic approach to obtain trees suitable for macroevolutionary analysis.

Phylogenomics of manakins (Aves: Pipridae) using alternative locus filtering strategies based on informativeness.

Target capture sequencing effectively generates molecular marker arrays useful for molecular systematics. These extensive data sets are advantageous where previous studies using a few loci have failed to resolve relationships confidently. Moreover, target capture is well-suited to fragmented source DNA, allowing data collection from species that lack fresh tissues. Herein we use target capture to generate data for a phylogeny of the avian family Pipridae (manakins), a group that has been the subject of many behavioral and ecological studies. Most manakin species feature lek mating systems, where males exhibit complex behavioral displays including mechanical and vocal sounds, coordinated movements of multiple males, and high speed movements. We analyzed thousands of ultraconserved element (UCE) loci along with a smaller number of coding exons and their flanking regions from all but one species of Pipridae. We examined three different methods of phylogenetic estimation (concatenation and two multispecies coalescent methods). Phylogenetic inferences using UCE data yielded strongly supported estimates of phylogeny regardless of analytical method. Exon probes had limited capability to capture sequence data and resulted in phylogeny estimates with reduced support and modest topological differences relative to the UCE trees, although these conflicts had limited support. Two genera were paraphyletic among all analyses and data sets, with Antilophia nested within Chiroxiphia and Tyranneutes nested within Neopelma. The Chiroxiphia-Antilophia clade was an exception to the generally high support we observed; the topology of this clade differed among analyses, even those based on UCE data. To further explore relationships within this group, we employed two filtering strategies to remove low-information loci. Those analyses resulted in distinct topologies, suggesting that the relationships we identified within Chiroxiphia-Antilophia should be interpreted with caution. Despite the existence of a few continuing uncertainties, our analyses resulted in a robust phylogenetic hypothesis of the family Pipridae that provides a comparative framework for future ecomorphological and behavioral studies.

Comparative Genomics Reveals a Burst of Homoplasy-Free Numt Insertions.

Mitochondrial DNA sequences are frequently transferred into the nuclear genome, giving rise to numts (nuclear mitochondrial DNA segments). In the absence of whole genomes, avian numts have been suggested to be rare and relatively short. We examined 64 bird genomes to test hypotheses regarding numt frequency, distribution among taxa, and likelihood of homoplasy. We discovered 100-fold variation in numt number across species. Two songbirds, Geospiza fortis (Darwin's finch) and Zonotrichia albicollis (white-throated sparrow) had the largest number of numts. Ancestral state reconstruction of 957 numt insertions in these two species and their close relatives indicated a remarkable acceleration of numt insertion in the ancestor of Geospiza and Zonotrichia followed by slower, continued accumulation in each lineage. These numts appear to result primarily from de novo insertion with the duplication of existing numts representing a secondary pathway. Insertion events were essentially homoplasy-free and numts appear to represent perfect rare genomic changes.

An evolutionary model motivated by physicochemical properties of amino acids reveals variation among proteins.

Motivation: The relative rates of amino acid interchanges over evolutionary time are likely to vary among proteins. Variation in those rates has the potential to reveal information about constraints on proteins. However, the most straightforward model that could be used to estimate relative rates of amino acid substitution is parameter-rich and it is therefore impractical to use for this purpose. Results: A six-parameter model of amino acid substitution that incorporates information about the physicochemical properties of amino acids was developed. It showed that amino acid side chain volume, polarity and aromaticity have major impacts on protein evolution. It also revealed variation among proteins in the relative importance of those properties. The same general approach can be used to improve the fit of empirical models such as the commonly used PAM and LG models. Availability and implementation: Perl code and test data are available from https://github.com/ebraun68/sixparam. Supplementary information: Supplementary data are available at Bioinformatics online.

The MAPK Kinase Kinase GmMEKK1 Regulates Cell Death and Defense Responses.

MAPK signaling pathways play critical roles in plant immunity. Here, we silenced multiple genes encoding MAPKs using virus-induced gene silencing mediated by Bean pod mottle virus to identify MAPK genes involved in soybean (Glycine max) immunity. Surprisingly, a strong hypersensitive response (HR) cell death was observed when soybean MAPK KINASE KINASE1 (GmMEKK1), a homolog of Arabidopsis (Arabidopsis thaliana) MEKK1, was silenced. The HR was accompanied by the overaccumulation of defense signaling molecules, salicylic acid (SA) and hydrogen peroxide. Genes involved in primary metabolism, translation/transcription, photosynthesis, and growth/development were down-regulated in GmMEKK1-silenced plants, while the expression of defense-related genes was activated. Accordingly, GmMEKK1-silenced plants were more resistant to downy mildew (Peronospora manshurica) and Soybean mosaic virus compared with control plants. Silencing GmMEKK1 reduced the activation of GmMPK6 but enhanced the activation of GmMPK3 in response to flg22 peptide. Unlike Arabidopsis MPK4, GmMPK4 was not activated by either flg22 or SA. Interestingly, transient overexpression of GmMEKK1 in Nicotiana benthamiana also induced HR. Our results indicate that GmMEKK1 plays both positive and negative roles in immunity and appears to differentially activate downstream MPKs by promoting GmMPK6 activation but suppressing GmMPK3 activation in response to flg22. The involvement of GmMPK4 kinase activity in cell death and in flg22- or SA-triggered defense responses in soybean requires further investigation.

What are the roles of taxon sampling and model fit in tests of cyto-nuclear discordance using avian mitogenomic data?

Conflicts between nuclear and mitochondrial phylogenies have led to uncertainty for some relationships within the tree of life. These conflicts have led some to question the value of mitochondrial DNA in phylogenetics now that genome-scale nuclear data can be readily obtained. However, since mitochondrial DNA is maternally inherited and does not recombine, its phylogeny should be closer to the species tree. Additionally, its rapid evolutionary rate may drive accumulation of mutations along short internodes where relevant information from nuclear loci may be limited. In this study, we examine the mitochondrial phylogeny of Cavitaves to elucidate its congruence with recently published nuclear phylogenies of this group of birds. Cavitaves includes the orders Trogoniformes (trogons), Bucerotiformes (hornbills), Coraciiformes (kingfishers and allies), and Piciformes (woodpeckers and allies). We hypothesized that sparse taxon sampling in previously published mitochondrial trees was responsible for apparent cyto-nuclear discordance. To test this hypothesis, we assembled 27 additional Cavitaves mitogenomes and estimated phylogenies using seven different taxon sampling schemes ranging from five to 42 ingroup species. We also tested the role that partitioning and model choice played in the observed discordance. Our analyses demonstrated that improved taxon sampling could resolve many of the disagreements. Similarly, partitioning was valuable in improving congruence with the topology from nuclear phylogenies, though the model used to generate the mitochondrial phylogenies had less influence. Overall, our results suggest that the mitochondrial tree is trustworthy when partitioning is used with suitable taxon sampling.

A Small Cysteine-Rich Protein from the Asian Soybean Rust Fungus, Phakopsora pachyrhizi, Suppresses Plant Immunity.

The Asian soybean rust fungus, Phakopsora pachyrhizi, is an obligate biotrophic pathogen causing severe soybean disease epidemics. Molecular mechanisms by which P. pachyrhizi and other rust fungi interact with their host plants are poorly understood. The genomes of all rust fungi encode many small, secreted cysteine-rich proteins (SSCRP). While these proteins are thought to function within the host, their roles are completely unknown. Here, we present the characterization of P. pachyrhizi effector candidate 23 (PpEC23), a SSCRP that we show to suppress plant immunity. Furthermore, we show that PpEC23 interacts with soybean transcription factor GmSPL12l and that soybean plants in which GmSPL12l is silenced have constitutively active immunity, thereby identifying GmSPL12l as a negative regulator of soybean defenses. Collectively, our data present evidence for a virulence function of a rust SSCRP and suggest that PpEC23 is able to suppress soybean immune responses and physically interact with soybean transcription factor GmSPL12l, a negative immune regulator.

A simple strategy for recovering ultraconserved elements, exons, and introns from low coverage shotgun sequencing of museum specimens: Placement of the partridge genus Tropicoperdix within the galliformes.

Next-generation DNA sequencing (NGS) offers a promising way to obtain massive numbers of orthologous loci to understand phylogenetic relationships among organisms. Of particular interest are old museum specimens and other samples with degraded DNA, where traditional sequencing methods have proven to be challenging. Low coverage shotgun sequencing and sequence capture are two widely used NGS approaches for degraded DNA. Sequence capture can yield sequence data for large numbers of orthologous loci, but it can only be used to sequence genomic regions near conserved sequences that can be used as probes. Low coverage shotgun sequencing has the potential to yield different data types throughout the genome. However, many studies using this method have often generated mitochondrial sequences, and few nuclear sequences, suggesting orthologous nuclear sequences are likely harder to recover. To determine the phylogenetic position of the galliform genus Tropicoperdix, whose phylogenetic position is currently uncertain, we explored two strategies to maximize data extraction from low coverage shotgun sequencing from approximately 100-year-old museum specimens from two species of Tropicoperdix. One approach, a simple read mapping strategy, outperformed the other (a reduced complexity assembly approach), and allowed us to obtain a large number of ultraconserved element (UCE) loci, relatively conserved exons, more variable introns, as well as mitochondrial genomes. Additionally, we demonstrated some simple approaches to explore possible artifacts that may result from the use of degraded DNA. Our data placed Tropicoperdix within a clade that includes many taxa characterized with ornamental eyespots (peafowl, argus pheasants, and peacock pheasants), and established relationships among species within the genus. Therefore, our study demonstrated that low coverage shotgun sequencing can easily be leveraged to yield substantial amounts and varying types of data, which opens the door for many research questions that might require information from different data types from museum specimens.

Why Do Phylogenomic Data Sets Yield Conflicting Trees? Data Type Influences the Avian Tree of Life more than Taxon Sampling.

Phylogenomics, the use of large-scale data matrices in phylogenetic analyses, has been viewed as the ultimate solution to the problem of resolving difficult nodes in the tree of life. However, it has become clear that analyses of these large genomic data sets can also result in conflicting estimates of phylogeny. Here, we use the early divergences in Neoaves, the largest clade of extant birds, as a "model system" to understand the basis for incongruence among phylogenomic trees. We were motivated by the observation that trees from two recent avian phylogenomic studies exhibit conflicts. Those studies used different strategies: 1) collecting many characters [$\sim$ 42 mega base pairs (Mbp) of sequence data] from 48 birds, sometimes including only one taxon for each major clade; and 2) collecting fewer characters ($\sim$ 0.4 Mbp) from 198 birds, selected to subdivide long branches. However, the studies also used different data types: the taxon-poor data matrix comprised 68% non-coding sequences whereas coding exons dominated the taxon-rich data matrix. This difference raises the question of whether the primary reason for incongruence is the number of sites, the number of taxa, or the data type. To test among these alternative hypotheses we assembled a novel, large-scale data matrix comprising 90% non-coding sequences from 235 bird species. Although increased taxon sampling appeared to have a positive impact on phylogenetic analyses the most important variable was data type. Indeed, by analyzing different subsets of the taxa in our data matrix we found that increased taxon sampling actually resulted in increased congruence with the tree from the previous taxon-poor study (which had a majority of non-coding data) instead of the taxon-rich study (which largely used coding data). We suggest that the observed differences in the estimates of topology for these studies reflect data-type effects due to violations of the models used in phylogenetic analyses, some of which may be difficult to detect. If incongruence among trees estimated using phylogenomic methods largely reflects problems with model fit developing more "biologically-realistic" models is likely to be critical for efforts to reconstruct the tree of life. [Birds; coding exons; GTR model; model fit; Neoaves; non-coding DNA; phylogenomics; taxon sampling.].