ALBERT

Description: Due to their functional independence, proteins that comprise standalone metabolic units, which we name single-protein metabolic modules, may be particularly prone to gene duplication (GD) and horizontal gene transfer (HGT). Flavohemoglobins (flavoHbs) are prime examples of single-protein metabolic modules, detoxifying nitric oxide (NO), a ubiquitous toxin whose antimicrobial properties many life forms exploit, to nitrate, a common source of nitrogen for organisms. FlavoHbs appear widespread in bacteria and have been identified in a handful of microbial eukaryotes, but how the distribution of this ecologically and biomedically important protein family evolved remains unknown. Reconstruction of the evolutionary history of 3,318 flavoHb protein sequences covering the family’s known diversity showed evidence of recurrent HGT at multiple evolutionary scales including intrabacterial HGT, as well as HGT from bacteria to eukaryotes. One of the most striking examples of HGT is the acquisition of a flavoHb by the dandruff- and eczema-causing fungus Malassezia from Corynebacterium Actinobacteria, a transfer that growth experiments show is capable of mediating NO resistance in fungi. Other flavoHbs arose via GD; for example, many filamentous fungi possess two flavoHbs that are differentially targeted to the cytosol and mitochondria, likely conferring protection against external and internal sources of NO, respectively. Because single-protein metabolic modules such as flavoHb function independently, readily undergo GD and HGT, and are frequently involved in organismal defense and competition, we suggest that they represent "plug-and-play" proteins for ecological arms races.

Description: The largest living carnivorous marsupial, the Tasmanian devil ( Sarcophilus harrisii ), is the sole survivor of a lineage originating about 12 Ma. We set out to investigate the spectrum of transposable elements found in the Tasmanian devil genome, the first high-coverage genome of an Australian marsupial. Marsupial genomes have been shown to have the highest amount of transposable elements among vertebrates. We analyzed the horizontally transmitted DNA transposons OC1 and hAT-1_MEu in the Tasmanian devil genome. OC1 is present in all carnivorous marsupials, while having a very limited distribution among the remaining Australian marsupial orders. In contrast, hAT-1_MEu is present in all Australian marsupial orders, and has so far only been identified in a few placental mammals. We screened 158 introns for phylogenetically informative retrotransposons in the order Dasyuromorphia, and found that the youngest SINE (Short INterspersed Element), WSINE1, is no longer active in the subfamily Dasyuridae. The lack of detectable WSINE1 activity in this group may be due to a retrotransposon inactivation event approximately 30 Ma. We found that the Tasmanian devil genome contains a relatively low number of continuous full-length LINE-1 (Long INterspersed Element 1, L1) retrotransposons compared with the opossum genome. Furthermore, all L1 elements in the Tasmanian devil appeared to be nonfunctional. Hidden Markov Model approaches suggested that other potential sources of functional reverse transcriptase are absent from the genome. We discuss the issues associated with assembling long, highly similar L1 copies from short read Illumina data and describe how assembly artifacts can potentially lead to erroneous conclusions.

Description: Annelid disparity has resulted in morphological-based classifications that disagree with phylogenies based on Sanger sequencing and phylogenomic analyses. However, the data used for the latter studies came from various sources and technologies, involved poorly occupied matrices and lacked key lineages. Here, we generated a new Illumina-based data set to address annelid relationships from a fresh perspective, independent from previously generated data and with nearly fully occupied matrices. Our sampling reflects the span of annelid diversity, including two symbiotic annelid groups (Myzostomida and Spinther ) and five meiofaunal groups once referred to as part of Archiannelida (three from Protodrilida, plus Dinophilus and Polygordius ). As well as the placement of these unusual annelids, we sought to address the overall phylogeny of Annelida, and provide a new perspective for naming of major clades. Our results largely corroborate the phylogenomic results of Weigert et al. (2014; Illuminating the base of the annelid tree using transcriptomics. Mol Biol Evol. 31:1391-1401), with " Magelona + Owenia " and Chaetopteridae forming a grade with respect to all other annelids. Echiura and Sipuncula are supported as being annelid groups, with Sipuncula closest to amphinomids as sister group to Sedentaria and Errantia. We recovered the three Protodrilida terminals as sister clade to Phyllodocida and Eunicida (=clade Aciculata). We therefore place Protodrilida as part of Errantia. Polygordius was found to be sister group to the scaleworm terminal and the possibility that it is a simplified scaleworm clade, as has been shown for the former family Pisionidae, is discussed. Our results were equivocal with respect to Dinophilus , Myzostomida, and Spinther possibly owing to confounding long-branch effects.

Description: Convergent evolution occurs when the same trait arises independently in multiple lineages. In most cases of phenotypic convergence such transitions are adaptive, so finding the underlying molecular causes of convergence can provide insight into the process of adaptation. Convergent evolution at the genomic level also lends itself to study by comparative methods, although molecular convergence can also occur by chance, adding noise to this process. Parker et al. studied convergence across the genomes of several mammals, including echolocating bats and dolphins (Parker J, Tsagkogeorga G, Cotton JA, Liu Y, Provero P, Stupka E, Rossiter SJ. 2013. Genome-wide signatures of convergent evolution in echolocating mammals. Nature 502:228–231). On the basis of a null distribution of site-specific likelihood support (SSLS) generated using simulated topologies, they concluded that there was evidence for genome-wide adaptive convergence between echolocating taxa. Here, we demonstrate that methods based on SSLS do not adequately measure convergence, and reiterate the use of an empirical null model that directly compares convergent substitutions between all pairs of species. We find that when the proper comparisons are made there is no surprising excess of convergence between echolocating mammals, even in sensory genes.

Description: Protein interaction networks play central roles in biological systems, from simple metabolic pathways through complex programs permitting the development of organisms. Multicellularity could only have arisen from a careful orchestration of cellular and molecular roles and responsibilities, all properly controlled and regulated. Disease reflects a breakdown of this organismal homeostasis. To better understand the evolution of interactions whose dysfunction may be contributing factors to disease, we derived the human protein coevolution network using our MatrixMatchMaker algorithm and using the Orthologous MAtrix project (OMA) database as a source for protein orthologs from 103 eukaryotic genomes. We annotated the coevolution network using protein–protein interaction data, many functional data sources, and we explored the evolutionary rates and dates of emergence of the proteins in our data set. Strikingly, clustering based only on the topology of the coevolution network partitions it into two subnetworks, one generally representing ancient eukaryotic functions and the other functions more recently acquired during animal evolution. That latter subnetwork is enriched for proteins with roles in cell–cell communication, the control of cell division, and related multicellular functions. Further annotation using data from genetic disease databases and cancer genome sequences strongly implicates these proteins in both ciliopathies and cancer. The enrichment for such disease markers in the animal network suggests a functional link between these coevolving proteins. Genetic validation corroborates the recruitment of ancient cilia in the evolution of multicellularity.

Description: Substitution rates vary between species, and many explanations regarding the causes of this variation have been proposed. Here we consider how new genomic data on the per-generation mutation rate impinge on proposed hypotheses for substitution rate variation in primates. We propose that the generation-time effect as it is usually understood cannot explain the observed rate variation, but instead that selection for decreased somatic mutation rates can. By considering the disparate causes underlying mutation rate changes in recent human history, we also show that the per-generation mutation rate is increasing even as the per-cell-division rate is decreasing.

Description: Current sequencing methods produce large amounts of data, but genome assemblies constructed from these data are often fragmented and incomplete. Incomplete and error-filled assemblies result in many annotation errors, especially in the number of genes present in a genome. This means that methods attempting to estimate rates of gene duplication and loss often will be misled by such errors and that rates of gene family evolution will be consistently overestimated. Here, we present a method that takes these errors into account, allowing one to accurately infer rates of gene gain and loss among genomes even with low assembly and annotation quality. The method is implemented in the newest version of the software package CAFE, along with several other novel features. We demonstrate the accuracy of the method with extensive simulations and reanalyze several previously published data sets. Our results show that errors in genome annotation do lead to higher inferred rates of gene gain and loss but that CAFE 3 sufficiently accounts for these errors to provide accurate estimates of important evolutionary parameters.