Evolution of the rbcS gene family in Solanaceae: concerted evolution and gain and loss of introns, with a description of new statistical guidelines for determining the number of unique gene copies.

Thesis (Ph.D.)--University of Washington, 2014

Detecting regular sound changes in linguistics as events of concerted evolution

Background: Concerted evolution is normally used to describe parallel changes at different sites in a genome, but it is also observed in languages where a specific phoneme changes to the same other phoneme in many words in the lexicon—a phenomenon known as regular sound change. We develop a general statistical model that can detect concerted changes in aligned sequence data and apply it to study regular sound changes in the Turkic language family. Results: Linguistic evolution, unlike the genetic substitutional process, is dominated by events of concerted evolutionary change. Our model identified more than 70 historical events of regular sound change that occurred throughout the evolution of the Turkic language family, while simultaneously inferring a dated phylogenetic tree. Including regular sound changes yielded an approximately 4-fold improvement in the characterization of linguistic change over a simpler model of sporadic change, improved phylogenetic inference, and returned more reliable and plausible dates for events on the phylogenies. The historical timings of the concerted changes closely follow a Poisson process model, and the sound transition networks derived from our model mirror linguistic expectations. Conclusions: We demonstrate that a model with no prior knowledge of complex concerted or regular changes can nevertheless infer the historical timings and genealogical placements of events of concerted change from the signals left in contemporary data. Our model can be applied wherever discrete elements—such as genes, words, cultural trends, technologies, or morphological traits—can change in parallel within an organism or other evolving group.

Documentation of reticulate evolution in peonies (Paeonia) using internal transcribed spacer sequences of nuclear ribosomal DNA: implications for biogeography and concerted evolution.

The internal transcribed spacers (ITS) of nuclear ribosomal DNA of 33 species of genus Paeonia (Paeoniaceae) were sequenced. In section Paeonia, different patterns of nucleotide additivity were detected in 14 diploid and tetraploid species at sites that are variable in the other 12 species of the section, suggesting that reticulate evolution has occurred. Phylogenetic relationships of species that do not show additivity, and thus ostensibly were not derived through hybridization, were reconstructed by parsimony analysis. The taxa presumably derived through reticulate evolution were then added to the phylogenetic tree according to additivity from putative parents. The study provides an example of successfully using ITS sequences to reconstruct reticulate evolution in plants and further demonstrates that the sequence data could be highly informative and accurate for detecting hybridization. Maintenance of parental sequences in the species of hybrid origin is likely due to slowing of concerted evolution caused by the long generation time of peonies. The partial and uneven homogenization of parental sequences displayed in nine species of putative hybrid origin may have resulted from gradients of gene conversion. The documented hybridizations may have occurred since the Pleistocene glaciations. The species of hybrid origin and their putative parents are now distantly allopatric. Reconstruction of reticulate evolution with sequence data, therefore, provides gene records for distributional histories of some of the parental species.

The 5S rDNA family evolves through concerted and birth-and-death evolution in fish genomes: an example from freshwater stingrays

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

松科nrDNA 的进化研究Evolution of nuclear ribosomal DNA in Pinaceae

核核糖体DNA（nrDNA）已被作为一个重要的标记，用于推断很多分类等级上的系统发育关系。相对于在被子植物中的快速致同进化，nrDNA在裸子植物中的致同进化速率低，且ITS和5S-NTS区有着较大的长度变异，这种现象在松科植物中尤为明显。在本研究中，我们克隆并测定了银杉属的5S rDNA以及冷杉属、银杉属、雪松属、油杉属、长苞铁杉属、金钱松属与铁杉属的ITS序列。基于获得的新数据，再结合前人报导的其它属的数据，我们探讨了如下四个问题： （1）松科 nrDNA ITS1 亚重复单位的组成、分布及进化;（2）ITS1区的长度变异与亚重复单位数目的关系以及它们的系统学意义;（3）松科ITS1的二级结构特征;（4）银杉5S rDNA编码区及非转录间隔区的结构特征。主要研究结果如下： 1. ITS区的序列分析ITS区的克隆及序列分析发现：（1） 松科ITS1的长度变异范围为 944-3271 bp, 这是目前已报导的真核生物中属间ITS变异最大的类群之一;（2） 所有松科植物的ITS区域都包含亚重复单位，亚重复单位的数目从2到9，并且这些亚重复单位可分为两种类型，即不含保守核心序列（5’-GGCCACCCTAGTC ） 的长亚重复单位（LSR）和含上述保守核心序列的短亚重复单位（SSR）;（3） ITS1区的巨大长度变异主要归因于亚重复单位的数量变异; （4） ITS1区的GC含量与 它的序列长度和亚重复单位的数目有一定关系，并能够提供一些系统发育信息，特别是支持云杉属、松属和银杉属三者具有很近的亲缘关系。 2. ITS1亚重复单位的系统发育分析为了研究亚重复单位的进化关系，我们用最大似然法和最大简约法构建了松科ITS1亚重复单位的系统发育树。结果表明：（1）在ML和MP树中可发现有共同的五个分支; （2） 银杉比松科其它属拥有更多的SSR，且该属的所有9个SSR在系统树中构成一个单系支，表明它们是在银杉属内发生重复的;（3）一些SSR在属间和种间具有同源性，可为nrDNA ITS 的进化历史以及松科的系统发育 研究提供重要信息;（4）亚重复单位的多次重复以及伴随的重组可能是导致LSR 和SSR在松科不同属中分布式样不同的原因。 3. 松科ITS1的二级结构 用 Mfold 3.2 软件对松科所有11个属的ITS1区进行了二级结构预测，共获得了563个最低自由能折叠。结合以前关于松科二级结构的报导，我们分析的结果表明：（1） 松科ITS1的二级结构主要由几个延展的发夹结构组成;（2） 构象的复杂性与亚重复的数目呈正相关;（3）配对的亚重复单位通常在保守核心区（5’-GGCCACCCTAGTC ） 处有部分重叠，并且构成一个长茎，而其它的亚重复单位通常会自身折叠，且保守核心区的部分出现在发夹结构的环中。 4. 银杉5S rDNA 序列分析 我们对来自银杉不同群体的3个个体的5S rDNA进行了克隆，共获得 45 条序列，分析结果表明：（1） 绝大多数银杉5S rDNA编码区长度为120 bp, 以GGG 开头，以CTC结尾，编码区出现的碱基替代主要为转换;（2） 银杉与其它裸子植物相比，5S rDNA基因编码区具很高的相似性（90-99%）; （3）间隔区含有一个poly-C和一个poly-T结构、两个TC丰富区以及五个GC丰富区。根据长度和序列特征，银杉的5S rDNA间隔区可分为三种类型：Type A 长751-764 bp，Type B 长770-807 bp （含一个32 bp的插入），Type C 长581-594 bp; （5）长间隔区（Type A，Type B ）中含有两个148-175 bp的串联亚重复单位，该亚重复单位与短间隔区（Type C ）中的一段143 bp的序列具有较高的相似性（56.0-66.8%）。 5. 银杉5S rRNA的二级结构 Mfold 3.2 预测结果表明：（1）银杉5S rRNA二级结构包括5个双螺旋区（干区）（Ⅰ-Ⅴ）、2个发夹结构环区（C和D）、3个中间环区（B1、B2 和 E）和1个铰链区（A）, 铰链区为三个双螺旋的结合处;（2） 二级结构中的环区通常比双螺旋区更加保守;（3）在5个双螺旋中，I 和 IV 区有较高的碱基替代率。

Comparative analysis of tandem repeats from hundreds of species reveals unique insights into centromere evolution

Background: Centromeres are essential for chromosome segregation, yet their DNA sequences evolve rapidly. In most animals and plants that have been studied, centromeres contain megabase-scale arrays of tandem repeats. Despite their importance, very little is known about the degree to which centromere tandem repeats share common properties between different species across different phyla. We used bioinformatic methods to identify high-copy tandem repeats from 282 species using publicly available genomic sequence and our own data.Results: Our methods are compatible with all current sequencing technologies. Long Pacific Biosciences sequence reads allowed us to find tandem repeat monomers up to 1,419 bp. We assumed that the most abundant tandem repeat is the centromere DNA, which was true for most species whose centromeres have been previously characterized, suggesting this is a general property of genomes. High-copy centromere tandem repeats were found in almost all animal and plant genomes, but repeat monomers were highly variable in sequence composition and length. Furthermore, phylogenetic analysis of sequence homology showed little evidence of sequence conservation beyond approximately 50 million years of divergence. We find that despite an overall lack of sequence conservation, centromere tandem repeats from diverse species showed similar modes of evolution.Conclusions: While centromere position in most eukaryotes is epigenetically determined, our results indicate that tandem repeats are highly prevalent at centromeres of both animal and plant genomes. This suggests a functional role for such repeats, perhaps in promoting concerted evolution of centromere DNA across chromosomes.

Purifying selection and birth-and-death evolution in the ubiquitin gene family

Ubiquitin is a highly conserved protein that is encoded by a multigene family. It is generally believed that this gene family is subject to concerted evolution, which homogenizes the member genes of the family. However, protein homogeneity can be attained also by strong purifying selection. We therefore studied the proportion (pS) of synonymous nucleotide differences between members of the ubiquitin gene family from 28 species of fungi, plants, and animals. The results have shown that pS is generally very high and is often close to the saturation level, although the protein sequence is virtually identical for all ubiquitins from fungi, plants, and animals. A small proportion of species showed a low level of pS values, but these values appeared to be caused by recent gene duplication. It was also found that the number of repeat copies of the gene family varies considerably with species, and some species harbor pseudogenes. These observations suggest that the members of this gene family evolve almost independently by silent nucleotide substitution and are subjected to birth-and-death evolution at the DNA level.

Evolution by the birth-and-death process in multigene families of the vertebrate immune system

Concerted evolution is often invoked to explain the diversity and evolution of the multigene families of major histocompatibility complex (MHC) genes and immunoglobulin (Ig) genes. However, this hypothesis has been controversial because the member genes of these families from the same species are not necessarily more closely related to one another than to the genes from different species. To resolve this controversy, we conducted phylogenetic analyses of several multigene families of the MHC and Ig systems. The results show that the evolutionary pattern of these families is quite different from that of concerted evolution but is in agreement with the birth-and-death model of evolution in which new genes are created by repeated gene duplication and some duplicate genes are maintained in the genome for a long time but others are deleted or become nonfunctional by deleterious mutations. We found little evidence that interlocus gene conversion plays an important role in the evolution of MHC and Ig multigene families.

Evolution of duplicate control regions in the mitochondrial genomes of metazoa: A case study with Australasian Ixodes ticks

To investigate the evolution pattern and phylogenetic utility of duplicate control regions (CRs) in mitochondrial (mt) genomes, we sequenced the entire mt genomes of three Ixodes species and part of the mt genomes of another I I species. All the species from the Australasian lineage have duplicate CRs, whereas the other species have one CR. Sequence analyses indicate that the two CRs of the Australasian Ixodes ticks have evolved in concert in each species. In addition to the Australasian Ixodes ticks, species from seven other lineages of metazoa also have mt genomes with duplicate CRs. Accumulated mtDNA sequence data from these metazoans and two recent experiments on replication of mt genomes in human cell lines with duplicate CRs allowed us to re-examine four intriguing questions about the presence of duplicate CRs in the mt genomes of metazoa: (1) Why do some mt genomes, but not others, have duplicate CRs? (2) How did mt genomes with duplicate CRs evolve? (3) How could the nucleotide sequences of duplicate CRs remain identical or very similar over evolutionary time? (4) Are duplicate CRs phylogenetic markers? It appears that mt genomes with duplicate CRs have a selective advantage in replication over mt genomes with one CR. Tandem duplication followed by deletion of genes is the most plausible mechanism for the generation of mt genomes with duplicate CRs. Once duplicate CRs occur in an mt genome, they tend to evolve in concert, probably by gene conversion. However, there are lineages where gene conversion may not always occur, and, thus, the two CRs may evolve independently in these lineages. Duplicate CRs have much potential as phylogenetic markers at low taxonomic levels, such as within genera, within families, or among families, but not at high taxonomic levels, such as among orders.

松属植物rRNA基因的变异模式及其进化生物学意义

被子植物的rRNA基因已经得到深入研究。二倍体被子植物一般拥有1-4对18S-5.8S-26S rDNA位点和1-2对5S rDNA位点。作为特殊的多基因家族成员，rDNA会受均一化力 (homogenizing forces) 的作用，通过基因转换、不等交换等机制，形成基因的致同进化 (concerted evolution)。长期以来，我们一直认为动植物rDNA致同进化水平很高，各种拷贝的序列几乎完全一致，因此可以直接应用PCR测序的方法进行分子系统学研究。但是在裸子植物中由于研究资料的匮乏，使我们对裸子植物rDNA的变异模式了解甚少。松属植物作为裸子植物的最大类群，它的rDNA变异和进化有何特点、与被子植物是否相同，是这个重要类群的进化研究中目前尚未解决的问题。本文的研究内容从三个方面进行： (1)rDNA的染色体定位 目前，松属的18S-5.8S-26S rDNA的染色体定位研究只包括5种植物，其中的3种同时涉及到5S rDNA定位。这些研究结果表明，不同种存在相异的rDNA位点数目，甚至不同的个体的rDNA位点均有变化。其共同点是，18S-5.8S-26S rDNA位点数平均较被子植物多，5S rDNA除Pinus radiata外，在其它种里则与被子植物相似。这种现象是松属或裸子植物的共同特征，亦或是特例呢？有限的研究限制了对裸子植物rDNA的了解。本研究的目的之一就是研究松属植物rDNA的染色体空间分布特征，希望借此了解松属植物间的关系，比较裸子植物和被子植物rDNA在染色体组水平的差异。 (2)5S rDNA的分子进化 5S rDNA的序列水平的进化研究在松属中尚属空白。5S rDNA在染色体数目上没有显示裸子植物与被子植物的差异，是否意味着松属乃至裸子植物的5S rDNA也同被子植物一样——致同进化完全，序列高度一致呢？利用克隆测序方法对松属植物5S rDNA的研究无疑是有开创性的工作，可以探讨裸子植物的5S rDNA的进化机制和种间关系。 (3)杂种基因组研究 杂交物种的起源演化是当前生物学研究的热点，通过杂种基因组的研究，可以了解杂种的的基因组构成，组织方式和进化历史，探讨杂交事件对成种过程的影响及意义。这项研究涉及到高山松、云南松和油松。之所以采用这三种植物，因为等位酶、cpDNA和mtDNA证据证明高山松为油松和云南松的自然杂交种。但这些证据不足以反映杂种核基因组的重组特征和构成及其进化规律。我们利用rDNA-FISH、5S rDNA和基因组原位杂交分析三种松树间的基因组关系，为揭示高山松的进化机制和历史提供新的依据。 本项研究得到以下结果： 一. rDNA荧光原位杂交 (FISH) 通过对华山松和白皮松两种单维管束亚属植物及油松、云南松、高山松、马尾松和南亚松等五种双维管束亚属植物的18S rDNA与5S rDNA的荧光原位杂交，结果表明： ⑴ 裸子植物的18S rDNA位点数目明显多于二倍体被子植物。其中主要位点数目，油松有7对，高山松5对，云南松8对，马尾松10对，南亚松6对，白皮松3对，华山松10对，平均在7对;另外，部分松树还存在弱位点。无论强弱位点都有部分存在于染色体的着丝粒区，除了赤松 (Pinus densiflora)，在其它松科植物中并没有发现这种现象。究竟是基因转移的结果或该位点是18S rDNA的原始起源位置还有待确证。 ⑵ 5S rDNA位点相对变异较小，与被子植物相当。除了华山松5S rDNA有4对位点，马尾松只有1对位点外，其它松树的5S rDNA位点数目均为2对，并且在双维管束亚属植物中有一对属于弱位点。 ⑶ 两种rDNA存在不同连锁模式。双维管束亚属植物中，5S与18S rDNA连锁在同一染色体的同一臂或两条臂上。在同一染色体臂时，18S rDNA在臂的远端。单维管束亚属植物的5S与18S rDNA或连锁于同一染色体的同一臂上，或分别处于不同染色体。前一情况，5S rDNA位于臂的远端。据此可以说明两个亚属的rDNA结构在染色体组水平的很大分化。 ⑷ 松属植物的关系及高山松核型特征。由于5S与18S rDNA连锁关系的不同，可以将单维管束亚属和双维管束亚属分开。各亚属的不同物种可以依据杂交位点的多少、位置、信号强弱构成的核型图加以区分，并且构成一定的系统关系。杂交起源的高山松在染色体组上，表现出对油松和云南松两亲本不同染色体特征的分别继承与重组，并产生独有的特征。其II同源染色体之一18S rDNA位点的缺失，可能是染色体重组的痕迹。 二. 5S rDNA的序列变异与分子进化 利用分子克隆和DNA测序分析了油松、云南松、马尾松、白皮松和不同遗传背景的高山松居群的5S rRNA基因序列变异及基因进化规律，得到以下主要结果： ⑴ 5S rDNA的结构特征。双维管束亚属植物长度在658-728 bp，白皮松则为499-521 bp。长度差异体现在基因间隔区，而基因区极端保守，基本为120 bp。基因转录区内部存在着转录控制区，决定了5S rRNA的转录起始与转录效率。5S rRNA基因能够折叠成正常的二级结构，其中，相对于干区来说，环区要保守，但环E却表现出异乎寻常的变异，转换/颠换比值高达7.1，这种突变可能是假基因的产物。基因间隔区存在一定的保守单元，其中一些与转录的起始和终止调控相关，有些是裸子植物未知功能的特异保守区。 ⑵ 松属植物5S rDNA存在着基因组内与种间的异质性。基因组内的各个克隆中有超过80%的特异的，彼此不相同。整个5S rDNA分化距离为0.042 - 0.051，其中，间隔区的分化比基因区高，其速度约是基因区的3-7倍。比较种间5S rDNA序列发现：在122个克隆中，基因区只有50个特异的序列。基因组间的序列变异度与基因组内 (个体内) 没有明显差别。白皮松的间隔区与双维管束亚属松树的5S rDNA间隔区差异极大，几乎不能排序，而四种双维管束亚属植物的5S rDNA间隔区种间种内差异不大。 ⑶ 松属植物5S rDNA进化。PAUP分析建立的5S rRNA基因树显示，5S rRNA基因在基因组内是多系的 (polyphyletic)，表明成种事件以前，祖先种就已经存在序列的分化。观测到的5S rRNA基因序列变异状况，并非完全是致同进化或独立进化的单一因素造成的，而是二者的相互作用的结果。致同进化确实存在，只是速度较慢而已。 ⑷ 高山松5S rDNA 组成。高山松拥有最高的基因组内的序列多样性，高山松的5S rDNA拷贝既有亲本类型，又有重组类型，并且不同地理及遗传来源的高山松显示一定的分化趋势，有更多的拷贝来自母系亲本。 三. 基因组原位杂交 以油松和云南松总DNA作为探针，相互进行基因组原位杂交，结果显示云南松和油松的染色体组可以完全被对方探针标记，在现有基因组原位杂交的分辨率下不能将两个基因组区分开。说明云南松和油松基因组之间存在高比例的同源序列，两种松树的基因组组成十分相似。利用油松和云南松总DNA作为探针，对高山松的染色体组进行双探针基因组原位杂交。结果表明，高山松全部基因组都能与两亲本探针完全杂交，说明三者间有着异乎寻常的亲缘关系。但在PH失调影响下，高山松只有部分基因组被杂交，并且两种探针的杂交信号有轻微差异。这可能是高度重复序列优先杂交的结果。这些情况表明，高山松虽然在基因组构成上与两个亲本基本一致，但基因在染色体组的空间排布上是存在差异的，这一点可以从rDNA-FISH中证明。

Role of the triad N46, S106 and T107 and the surface charges in the determination of the acidic pH optimum of digestive lysozymes from Musca domestica

Structures of digestive lysozymes 1 and 2 from housefly (MdL1 and MdL2) show that S106-T107 delimit a polar pocket around E32 (catalytic acid/base) and N46 contributes to the positioning of 050 (catalytic nucleophile), whereas those residues are replaced by V109-A110 and D48 in the non-digestive lysozyme from hen egg-white (HEWL). Further analyses revealed that MdL1 and MdL2 surfaces are less positively charged than HEWL surface. To verify the relevance of these differences to the acidic pH optimum of digestive lysozymes it was determined that pKas of the catalytic residues of the triple mutant MdL2 (N46D-S106V-T107A) are similar to HEWL pKas and higher than those for MdL2. In agreement, triple mutant MdL2 and HEWL exhibits the same pH optimum upon methylumbelliferylchitotrioside. In addition to that, the introduction of six basic residues on MdL1 surface increased by 1 unit the pH optimum for the activity upon bacterial walls. Thus, the acidic pH optimum for MdL2 and MdL1 activities upon methylumbelliferylchitotrioside is determined by the presence of N46, S106 and T107 in the environment of their catalytic residues, which favors pKas reduction. Conversely, acidic pH optimum upon bacterial walls is determined by a low concentration of positive charges on the MdL2 and MdL1 surfaces. (C) 2010 Elsevier Inc. All rights reserved.

High similarity of U2 snDNA sequence between A and B chromosomes in the grasshopper Abracris flavolineata

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Multiple independent origins of mitochondrial control region duplications in the order Psittaciformes

mitochondrial genomes are generally thought to be under selection for compactness, due to their small size, consistent gene content, and a lack of introns or intergenic spacers. As more animal mitochondrial genomes are fully sequenced, rearrangements and partial duplications are being identified with increasing frequency, particularly in birds (Class Ayes). In this study, we investigate the evolutionary history of mitochondrial control region states within the avian order Psittaciformes (parrots and cockatoos). To this aim, we reconstructed a comprehensive multi-locus phylogeny of parrots, used PCR of three diagnostic fragments to classify the mitochondrial control region state as single or duplicated, and mapped these states onto the phylogeny. We further sequenced 44 selected species to validate these inferences of control region state. Ancestral state reconstruction using a range of weighting schemes identified six independent origins of mitochondrial control region duplications within Psittaciformes. Analysis of sequence data showed that varying levels of mitochondrial gene and tRNA homology and degradation were present within a given clade exhibiting duplications. Levels of divergence between control regions within an individual varied from 0-10.9% with the differences occurring mainly between 51 and 225 nucleotides 3' of the goose hairpin in domain I. Further investigations into the fates of duplicated mitochondrial genes, the potential costs and benefits of having a second control region, and the complex relationship between evolutionary rates, selection, and time since duplication are needed to fully explain these patterns in the mitochondrial genome. (C) 2012 Elsevier Inc. All rights reserved.

Transposable elements: structure and dynamic of the autonomous retrotransposon R2 in Arthropoda with non-canonical reproduction

Eukaryotic ribosomal DNA constitutes a multi gene family organized in a cluster called nucleolar organizer region (NOR); this region is composed usually by hundreds to thousands of tandemly repeated units. Ribosomal genes, being repeated sequences, evolve following the typical pattern of concerted evolution. The autonomous retroelement R2 inserts in the ribosomal gene 28S, leading to defective 28S rDNA genes. R2 element, being a retrotransposon, performs its activity in the genome multiplying its copy number through a “copy and paste” mechanism called target primed reverse transcription. It consists in the retrotranscription of the element’s mRNA into DNA, then the DNA is integrated in the target site. Since the retrotranscription can be interrupted, but the integration will be carried out anyway, truncated copies of the element will also be present in the genome. The study of these truncated variants is a tool to examine the activity of the element. R2 phylogeny appears, in general, not consistent with that of its hosts, except some cases (e.g. Drosophila spp. and Reticulitermes spp.); moreover R2 is absent in some species (Fugu rubripes, human, mouse, etc.), while other species have more R2 lineages in their genome (the turtle Mauremys reevesii, the Japanese beetle Popilia japonica, etc). R2 elements here presented are isolated in 4 species of notostracan branchiopods and in two species of stick insects, whose reproductive strategies range from strict gonochorism to unisexuality. From sequencing data emerges that in Triops cancriformis (Spanish gonochoric population), in Lepidurus arcticus (two putatively unisexual populations from Iceland) and in Bacillus rossius (gonochoric population from Capalbio) the R2 elements are complete and encode functional proteins, reflecting the general features of this family of transposable elements. On the other hand, R2 from Italian and Austrian populations of T. cancriformis (respectively unisexual and hermaphroditic), Lepidurus lubbocki (two elements within the same Italian population, gonochoric but with unfunctional males) and Bacillus grandii grandii (gonochoric population from Ponte Manghisi) have sequences that encode incomplete or non-functional proteins in which it is possible to recognize only part of the characteristic domains. In Lepidurus couesii (Italian gonochoric populations) different elements were found as in L. lubbocki, and the sequencing is still in progress. Two hypothesis are given to explain the inconsistency of R2/host phylogeny: vertical inheritance of the element followed by extinction/diversification or horizontal transmission. My data support previous study that state the vertical transmission as the most likely explanation; nevertheless horizontal transfer events can’t be excluded. I also studied the element’s activity in Spanish populations of T. cancriformis, in L. lubbocki, in L. arcticus and in gonochoric and parthenogenetic populations of B. rossius. In gonochoric populations of T. cancriformis and B. rossius I found that each individual has its own private set of truncated variants. The situation is the opposite for the remaining hermaphroditic/parthenogenetic species and populations, all individuals sharing – in the so far analyzed samples - the majority of variants. This situation is very interesting, because it isn’t concordant with the Muller’s ratchet theory that hypothesizes the parthenogenetic populations being either devoided of transposable elements or TEs overloaded. My data suggest a possible epigenetic mechanism that can block the retrotransposon activity, and in this way deleterious mutations don’t accumulate.

Contrasting histories of avian and mammalian Mhc genes revealed by class II B sequences from songbirds.

To explore the evolutionary dynamics of genes in the major histocompatibility complex (Mhc) in nonmammalian vertebrates, we have amplified complete sequences of the polymorphic second (beta1) and third (beta2) exons of class II beta chain genes of songbirds. The pattern of nucleotide substitution in the antigen-binding site of sequences cloned from three behaviorally and phylogenetically divergent songbirds [scrub jays Aphelocoma coerulescens), red-winged blackbirds (Agelaius phoeniceus), and house finches (Carpodacus mexicanus) reveals that class II B genes of songbirds are subject to the same types of diversifying forces as those observed at mammalian class II loci. By contrast, the tree of avian class II B genes reveals that orthologous relationships have not been retained as in placental mammals and that, unlike class II genes in mammals, genes in songbirds and chickens have had very recent common ancestors within their respective groups. Thus, whereas the selective forces diversifying class II B genes of birds are likely similar to those in mammals, their long-term evolutionary dynamics appear to be characterized by much higher rates of concerted evolution.