Certain observations concerning the effects of epistasis on complex traits and the evolution of genomes.

Thesis (Ph.D.)--University of Washington, 2016-05

Drosophila 3′ UTRs are more complex than protein-coding sequences

The 3′ UTRs of eukaryotic genes participate in a variety of post-transcriptional (and some transcriptional) regulatory interactions. Some of these interactions are well characterised, but an undetermined number remain to be discovered. While some regulatory sequences in 3′ UTRs may be conserved over long evolutionary time scales, others may have only ephemeral functional significance as regulatory profiles respond to changing selective pressures. Here we propose a sensitive segmentation methodology for investigating patterns of composition and conservation in 3′ UTRs based on comparison of closely related species. We describe encodings of pairwise and three-way alignments integrating information about conservation, GC content and transition/transversion ratios and apply the method to three closely related Drosophila species: D. melanogaster, D. simulans and D. yakuba. Incorporating multiple data types greatly increased the number of segment classes identified compared to similar methods based on conservation or GC content alone. We propose that the number of segments and number of types of segment identified by the method can be used as proxies for functional complexity. Our main finding is that the number of segments and segment classes identified in 3′ UTRs is greater than in the same length of protein-coding sequence, suggesting greater functional complexity in 3′ UTRs. There is thus a need for sustained and extensive efforts by bioinformaticians to delineate functional elements in this important genomic fraction. C code, data and results are available upon request.

Evolutionary Genomics of Prokaryotic Viruses

Evolutionary history of biological entities is recorded within their nucleic acid sequences and can (sometimes) be deciphered by thorough genomic analysis. In this study we sought to gain insights into the diversity and evolution of bacterial and archaeal viruses. Our primary interest was pointed towards those virus groups/families for which comprehensive genomic analysis was not previously possible due to the lack of sufficient amount of genomic data. During the course of this work twenty-five putative proviruses integrated into various prokaryotic genomes were identified, enabling us to undertake a comparative genomics approach. This analysis allowed us to test the previously formulated evolutionary hypotheses and also provided valuable information on the molecular mechanisms behind the genome evolution of the studied virus groups.

Characterization of the molecular components and function of BARE-1, Hin-Mu and Mu transposition machineries

Transposable elements, transposons, are discrete DNA segments that are able to move or copy themselves from one locus to another within or between their host genome(s) without a requirement for DNA homology. They are abundant residents in virtually all the genomes studied, for instance, the genomic portion of TEs is approximately 3% in Saccharomyces cerevisiae, 45% in humans, and apparently more than 70% in some plant genomes such as maize and barley. Transposons plays essential role in genome evolution, in lateral transfer of antibiotic resistance genes among bacteria and in life cycle of certain viruses such as HIV-1 and bacteriophage Mu. Despite the diversity of transposable elements they all use a fundamentally similar mechanism called transpositional DNA recombination (transposition) for the movement within and between the genomes of their host organisms. The DNA breakage and joining reactions that underlie their transposition are chemically similar in virtually all known transposition systems. The similarity of the reactions is also reflected in the structure and function of the catalyzing enzymes, transposases and integrases. The transposition reactions take place within the context of a transposition machinery, which can be particularly complex, as in the case of the VLP (virus like particle) machinery of retroelements, which in vivo contains RNA or cDNA and a number of element encoded structural and catalytic proteins. Yet, the minimal core machinery required for transposition comprises a multimer of transposase or integrase proteins and their binding sites at the element DNA ends only. Although the chemistry of DNA transposition is fairly well characterized, the components and function of the transposition machinery have been investigated in detail for only a small group of elements. This work focuses on the identification, characterization, and functional studies of the molecular components of the transposition machineries of BARE-1, Hin-Mu and Mu. For BARE-1 and Hin-Mu transpositional activity has not been shown previously, whereas bacteriophage Mu is a general model of transposition. For BARE-1, which is a retroelement of barley (Hordeum vulgare), the protein and DNA components of the functional VLP machinery were identified from cell extracts. In the case of Hin-Mu, which is a Mu-like prophage in Haemophilus influenzae Rd genome, the components of the core machinery (transposase and its binding sites) were characterized and their functionality was studied by using an in vitro methodology developed for Mu. The function of Mu core machinery was studied for its ability to use various DNA substrates: Hin-Mu end specific DNA substrates and Mu end specific hairpin substrates. The hairpin processing reaction by MuA was characterized in detail. New information was gained of all three machineries. The components or their activity required for functional BARE-1 VLP machinery and retrotransposon life cycle were present in vivo and VLP-like structures could be detected. The Hin-Mu core machinery components were identified and shown to be functional. The components of the Mu and Hin-Mu core machineries were partially interchangeable, reflecting both evolutionary conservation and flexibility within the core machineries. The Mu core machinery displayed surprising flexibility in substrate usage, as it was able to utilize Hin-Mu end specific DNA substrates and to process Mu end DNA hairpin substrates. This flexibility may be evolutionarily and mechanistically important.

Transposable elements in a marginal plant population

Background
How new forms arise in nature has engaged evolutionary biologists since Darwin's seminal treatise on the origin of species. Transposable elements (TEs) may be among the most important internal sources for intraspecific variability. Thus, we aimed to explore the temporal dynamics of several TEs in individual genotypes from a small, marginal population of Aegilops speltoides. A diploid cross-pollinated grass species, it is a wild relative of the various wheat species known for their large genome sizes contributed by an extraordinary number of TEs, particularly long terminal repeat (LTR) retrotransposons. The population is characterized by high heteromorphy and possesses a wide spectrum of chromosomal abnormalities including supernumerary chromosomes, heterozygosity for translocations, and variability in the chromosomal position or number of 45S and 5S ribosomal DNA (rDNA) sites. We propose that variability on the morphological and chromosomal levels may be linked to variability at the molecular level and particularly in TE proliferation.

Results
Significant temporal fluctuation in the copy number of TEs was detected when processes that take place in small, marginal populations were simulated. It is known that under critical external conditions, outcrossing plants very often transit to self-pollination. Thus, three morphologically different genotypes with chromosomal aberrations were taken from a wild population of Ae. speltoides, and the dynamics of the TE complex traced through three rounds of selfing. It was discovered that: (i) various families of TEs vary tremendously in copy number between individuals from the same population and the selfed progenies; (ii) the fluctuations in copy number are TE-family specific; (iii) there is a great difference in TE copy number expansion or contraction between gametophytes and sporophytes; and (iv) a small percentage of TEs that increase in copy number can actually insert at novel locations and could serve as a bona fide mutagen.

Conclusions
We hypothesize that TE dynamics could promote or intensify morphological and karyotypical changes, some of which may be potentially important for the process of microevolution, and allow species with plastic genomes to survive as new forms or even species in times of rapid climatic change.

高山松的二倍体杂种起源-来自三个基因组的证据

利用松科植物特殊的遗传体系（叶绿体基因组一父系遗传、线粒体基因组—母系遗传、核基因组一双亲遗传），我们对高山松及其两个亲本种进行了广泛的群体取样，通过线粒体基因nadl、叶绿体基因rbcL和trnL-F基因间区以及低拷贝核基因4CL的序列分析或PCR-RFLP分析，为高山松同倍体杂种起源假说提供了翔实的遗传学证据，同时在个体水平上探讨了高山松不同群体的遗传组成、群体遗传结构、基因交流方向、群体建立过程以及杂种基因组的进化。具体结果如下： 1．细胞质基因组分析 1)线粒体基因nudl分析 本研究对油松、高山松和云南松的19个群体、295个个体的线粒体基因nadl的一个内含子进行了序列分析或PCR-RFLP分析，共检测到3种线粒体DNA单倍型-A、B和C。油松所有的取样群体仅含单倍型A;除BX群体外，所有的云南松群体仅含单倍型B; 10个高山松群体中，5个群体固定单倍型A，4个群体固定单倍型B，1个群体(ZD)分布有A和B两种单倍型。2)叶绿体rbcL基因分析 对同一组群体的rbcL基因进行序列分析或PCR-RFLP分析，共检测到两个变异位点和三种叶绿体单倍型（TT、TC和GC）。TT和GC分别是油松和云南松种特异性叶绿体单倍型，而在高山松群体里则三种单倍型均有分布，而且TC单倍型广泛地分布在7个杂种群体中，该单倍型很可能来源于点突变或第三个已灭绝的亲本。rbcL基因检测到的高山松群体分化系数很高(Gst=0.533)。 3)叶绿体trn L-F区序列分析 叶绿体trnL-F分子标记检测到的不同单倍型的差异主要是由引物“e”下游120碱基处一个多聚T结构的长度变异所致（叶绿体SSR位点）。10个高山松群体中共检测到5种叶绿体单倍型，其中两种主要的单倍型(9T和11T)分别为油松和云南松的种特异性单倍型，其他单倍型均为非典型单倍型。群体遗传结构分析表明：杂种群体表现最高的遗传多样性，而且trnL-F分析得到的高山松群体的分化系数也很高( Gst=0.443)。 总之，对高山松、油松和云南松的同一组群体取样进行的细胞质基因组分析表明：高山松群体分布有油松和云南松种特异性的线粒体和叶绿体单倍型，该细胞质DNA单倍型的地理分布为假说“高山松为油松和云南松的的二倍体杂种”提供了翔实的遗传学证据。油松和云南松在不同的杂种群体中分别做父本和母本，即两亲本在杂交过程中发生了双向基因交流。群体遗传结构分析发现高山松群体表现最高的遗传多样性，而且群体间的分化系数很高。不同的杂种群体在遗传组成上的差异表明他们经历过不同的建立和进化历史。从线粒体和叶绿体单倍型的地理分布可以看出杂种群体的建立曾经历强烈的奠基者效应和回交。青藏高原的隆升对高山松的起源、杂种群体的适应辐射以及保持产生了重要的影响。川西南和滇西北作为青藏高原的东边边界，很可能是当初云南松和油松分布的重叠区及杂交地带，即高山松的起源地。 2．核基因4CL分析 对高山松、油松和云南松的19个群体、32个个体的低拷贝核基因4CL进行了克隆及序列分析，获得的78条序列可分为两种类型（类型A和类型B）。这两种类型明显的差别是类型A相对于类型B在内含子区有- 20bp的缺失。以华山松的3条序列为外类群，对得到的78条序列进行基因谱系分析，发现所有的序列分成明显的两支，分别对应于类型A和类型B，而且每一支均包含三个种的部分序列，表明4CL基因在这三个种分化之前就已发生重复。另一个明显的特点是某个种的一条序列与另一个种的序列比其与同种的其他序列关系更近，可能因基因交流（杂交和渐渗）、非共祖、致同进化和重组等进化事件所致。三种松树中共检测到4CL基因序列的两种类型和六个亚类型，高山松群体中没有发现杂种独特的类型或亚类型。高山松和云南松共享三种序列亚类型以及最多的序列多态性，表明这两个种之间曾存在广泛的基因交流。

基因组的演化与人类起源

人类起源是生命演化的奇迹,也是长期以来生命科学领域最受关注的重大课题之一.作为目前已知的唯一高级智慧生物,人类对自身的起源问题一直十分关心.了解自身是人类最基本的愿望.在远古时代,人们就开始思考我们究竟是什么?我们从哪里来?历史上有很多关于人类起源的传说,比如女娲造人、亚当夏娃的故事、海猿说、杂交说、外星人说等.随着科学的发展,特别是自达尔文发表划时代的进化论以来,"人类跟其他动物一样,是自然进化的产物,从一个远古的似猿祖先进化而来"的思想已经被大家普遍接受.

Comparison of Pax1/9 locus reveals 500-Myr-old syntenic block and evolutionary conserved noncoding regions

Identification of conserved genomic regions within and between different genomes is crucial when studying genome evolution. Here, we described regions of strong synteny conservation between vertebrate deuterostomes (tetrapods and teleosts) and invertebrat

Stepwise loss of motilin and its specific receptor genes in rodents

Specific interactions among biomolecules drive virtually all cellular functions and underlie phenotypic complexity and diversity. Biomolecules are not isolated particles, but are elements of integrated interaction networks, and play their roles through specific interactions. Simultaneous emergence or loss of multiple interacting partners is unlikely. If one of the interacting partners is lost, then what are the evolutionary consequences for the retained partner? Taking advantages of the availability of the large number of mammalian genome sequences and knowledge of phylogenetic relationships of the species, we examined the evolutionary fate of the motilin (MLN) hormone gene, after the pseudogenization of its specific receptor, MLN receptor (MLNR), on the rodent lineage. We speculate that the MLNR gene became a pseudogene before the divergence of the squirrel and other rodents about 75 mya. The evolutionary consequences for the MLN gene were diverse. While an intact open reading frame for the MLN gene, which appears functional, was preserved in the kangaroo rat, the MLN gene became inactivated independently on the lineages leading to the guinea pig and the common ancestor of the mouse and rat. Gain and loss of specific interactions among biomolecules through the birth and death of genes for biomolecules point to a general evolutionary dynamic: gene birth and death are widespread phenomena in genome evolution, at the genetic level; thus, once mutations arise, a stepwise process of elaboration and optimization ensues, which gradually integrates and orders mutations into a coherent pattern.

The origin of new genes: Glimpses from the young and old

Genome data have revealed great variation in the numbers of genes in different organisms, which indicates that there is a fundamental process of genome evolution: the origin of new genes. However, there has been little opportunity to explore how genes with new functions originate and evolve. The study of ancient genes has highlighted the antiquity and general importance of some mechanisms of gene origination, and recent observations of young genes at early stages in their evolution have unveiled unexpected molecular and evolutionary processes.

PCR survey of Hox genes in the goldfish Carassius auratus auratus

A tetraploidization event took place in the cyprinid lineage leading to goldfishes about 15 million years ago. A PCR survey for Hox genes in the goldfish Carassius auratus auratus (Actinopterygii: Cyprinidae) was performed to assess the consequences of this genome duplication. Not surprisingly, the genomic organization of the Hox gene clusters of goldfish is similar to that of the closely related zebrafish (Danio rerio). However, the goldfish exhibits a much larger number of recent pseudogenes, which are characterized by indels. These findings are consistent with the hypothesis that dosage effects cause selection pressure to rapidly silence crucial developmental regulators after a tetraploidization event.

C-values of seven marine mammal species determined by flow cytometry

C-values, which estimate genome size, have puzzled geneticists for years because they bear no relationship to organismal complexity. Though C-values have been estimated for thousands of species, considerably more data are required in order to better understanding genome evolution. This is particularly true for mammals, in which C-values are known for less than 8% of the total number of mammalian species. Among marine mammals, a C-value has been estimated only for the bottlenose dolphin (Tursiops truncatus). Thus examination of additional species of marine mammals is necessary for comparative purposes. It will enable a better understanding of marine mammal genome evolution, and it is also relevant to conservation, because larger genome size has been linked to increased likelihood of extinction in some plant and animal groups. Our study presents C-values of seven marine mammal species, including five cetacean species that are endangered to varying degrees. Similarly to the results for other groups, our results suggest that larger genome size in cetaceans is related to an increased likelihood of extinction.

Between two fern genomes.

Ferns are the only major lineage of vascular plants not represented by a sequenced nuclear genome. This lack of genome sequence information significantly impedes our ability to understand and reconstruct genome evolution not only in ferns, but across all land plants. Azolla and Ceratopteris are ideal and complementary candidates to be the first ferns to have their nuclear genomes sequenced. They differ dramatically in genome size, life history, and habit, and thus represent the immense diversity of extant ferns. Together, this pair of genomes will facilitate myriad large-scale comparative analyses across ferns and all land plants. Here we review the unique biological characteristics of ferns and describe a number of outstanding questions in plant biology that will benefit from the addition of ferns to the set of taxa with sequenced nuclear genomes. We explain why the fern clade is pivotal for understanding genome evolution across land plants, and we provide a rationale for how knowledge of fern genomes will enable progress in research beyond the ferns themselves.

Identification of dopamine receptors across the extant avian family tree and analysis with other clades uncovers a polyploid expansion among vertebrates.

Dopamine is an important central nervous system transmitter that functions through two classes of receptors (D1 and D2) to influence a diverse range of biological processes in vertebrates. With roles in regulating neural activity, behavior, and gene expression, there has been great interest in understanding the function and evolution dopamine and its receptors. In this study, we use a combination of sequence analyses, microsynteny analyses, and phylogenetic relationships to identify and characterize both the D1 (DRD1A, DRD1B, DRD1C, and DRD1E) and D2 (DRD2, DRD3, and DRD4) dopamine receptor gene families in 43 recently sequenced bird genomes representing the major ordinal lineages across the avian family tree. We show that the common ancestor of all birds possessed at least seven D1 and D2 receptors, followed by subsequent independent losses in some lineages of modern birds. Through comparisons with other vertebrate and invertebrate species we show that two of the D1 receptors, DRD1A and DRD1B, and two of the D2 receptors, DRD2 and DRD3, originated from a whole genome duplication event early in the vertebrate lineage, providing the first conclusive evidence of the origin of these highly conserved receptors. Our findings provide insight into the evolutionary development of an important modulatory component of the central nervous system in vertebrates, and will help further unravel the complex evolutionary and functional relationships among dopamine receptors.