Positive Selection and Enhancer Evolution Shaped Lifespan and Body Mass in Great Apes

Within primates, the great apes are outliers both in terms of body size and lifespan, since they include the largest and longest-lived species in the order. Yet, the molecular bases underlying such features are poorly understood. Here, we leveraged an integrated approach to investigate multiple sources of molecular variation across primates, focusing on over 10,000 genes, including approximately 1,500 previously associated with lifespan, and additional approximately 9,000 for which an association with longevity has never been suggested. We analyzed dN/dS rates, positive selection, gene expression (RNA-seq), and gene regulation (ChIP-seq). By analyzing the correlation between dN/dS, maximum lifespan, and body mass, we identified 276 genes whose rate of evolution positively correlates with maximum lifespan in primates. Further, we identified five genes, important for tumor suppression, adaptive immunity, metastasis, and inflammation, under positive selection exclusively in the great ape lineage. RNA-seq data, generated from the liver of six species representing all the primate lineages, revealed that 8% of approximately 1,500 genes previously associated with longevity are differentially expressed in apes relative to other primates. Importantly, by integrating RNA-seq with ChIP-seq for H3K27ac (which marks active enhancers), we show that the differentially expressed longevity genes are significantly more likely than expected to be located near a novel “ape-specific” enhancer. Moreover, these particular ape-specific enhancers are enriched for young transposable elements, and specifically SINE–Vntr–Alus. In summary, we demonstrate that multiple evolutionary forces have contributed to the evolution of lifespan and body size in primates.     

Comparative Genomics Elucidates the Origin of a Supergene Controlling Floral Heteromorphism

Supergenes are nonrecombining genomic regions ensuring the coinheritance of multiple, coadapted genes. Despite the importance of supergenes in adaptation, little is known on how they originate. A classic example of supergene is the S locus controlling heterostyly, a floral heteromorphism occurring in 28 angiosperm families. In Primula, heterostyly is characterized by the cooccurrence of two complementary, self-incompatible floral morphs and is controlled by five genes clustered in the hemizygous, ca. 300-kb S locus. Here, we present the first chromosome-scale genome assembly of any heterostylous species, that of Primula veris (cowslip). By leveraging the high contiguity of the P. veris assembly and comparative genomic analyses, we demonstrated that the S-locus evolved via multiple, asynchronous gene duplications and independent gene translocations. Furthermore, we discovered a new whole-genome duplication in Ericales that is specific to the Primula lineage. We also propose a mechanism for the origin of S-locus hemizygosity via nonhomologous recombination involving the newly discovered two pairs of CFB genes flanking the S locus. Finally, we detected only weak signatures of degeneration in the S locus, as predicted for hemizygous supergenes. The present study provides a useful resource for future research addressing key questions on the evolution of supergenes in general and the S locus in particular: How do supergenes arise? What is the role of genome architecture in the evolution of complex adaptations? Is the molecular architecture of heterostyly supergenes across angiosperms similar to that of Primula?     

新一代生物防治及其在天敌瓢虫中的研究进展

生物防治是一种利用自然天敌控制有害生物的手段.当前生物防治的发展面临着使用成本高和效果不稳定的两大难题.这要求我们更精准地选育和释放天敌.随着测序技术的进步,学界提出了新一代生物防治(Next Generation Biological Control),旨在利用遗传学和基因组学来认识天敌昆虫种下水平的多样性和变异规律...     

Comparative Genomics and Evolution of Avian Specialized Traits

Genomic data are important for understanding the origin and evolution of traits. Under the context of rapidly developing of sequencing technologies and more widely available genome sequences, researchers are able to study evolutionary mechanisms of traits via comparative genomic methods. Compared with other vertebrates, bird genomes are relatively small and exhibit conserved synteny with few repetitive elements, which makes them suitable for evolutionary studies. Increasing genomic progress has been reported on the evolution of powered flight, body size variation, beak morphology, plumage colouration, high-elevation colonization, migration, and vocalization. By summarizing previous studies, we demonstrate the genetic bases of trait evolution, highlighting the roles of small-scale sequence variation, genomic structural variation, and changes in gene interaction networks. We suggest that future studies should focus on improving the quality of reference genomes, exploring the evolution of regulatory elements and networks, and combining genomic data with morphological, ecological, behavioural, and developmental biology data.     

Specificity of mRNA Folding and Its Association with Evolutionarily Adaptive mRNA Secondary Structures

The secondary structure is a fundamental feature of both non-coding RNAs (ncRNAs) and messenger RNAs (mRNAs). However, our understanding of the secondary structures of mRNAs, especially those of the coding regions, remains elusive, likely due to translation and the lack of RNA-binding proteins that sustain the consensus structure like those binding to ncRNAs. Indeed, mRNAs have recently been found to adopt diverse alternative structures, but the overall functional significance remains untested. We hereby approach this problem by estimating the folding specificity, i.e., the probability that a fragment of an mRNA folds back to the same partner once refolded. We show that the folding specificity of mRNAs is lower than that of ncRNAs and exhibits moderate evolutionary conservation. Notably, we find that specific rather than alternative folding is likely evolutionarily adaptive since specific folding is frequently associated with functionally important genes or sites within a gene. Additional analysis in combination with ribosome density suggests the ability to modulate ribosome movement as one potential functional advantage provided by specific folding. Our findings reveal a novel facet of the RNA structurome with important functional and evolutionary implications and indicate a potential method for distinguishing the mRNA secondary structures maintained by natural selection from molecular noise.     

An automated annotation tool for genomic DNA sequences using GeneScan and BLAST

Genomic sequence data are often available well before the annotated sequence is published. We present a method for analysis of genomic DNA to identify coding sequences using the GeneScan algorithm and characterize these resultant sequences by BLAST. The routines are used to develop a system for automated annotation of genome DNA sequences.     

Introduction into Plant Genomics

The success in complete sequencing of small genomes and development of new technologies that markedly speed up the cloning and sequencing processes open the way to intense development of plant genomics and complete sequencing of DNA of some species. It is assumed that success in plant genomics will result in revolutionary changes in biotechnology and plant breeding. However, the enormous size of genomes (tens of billions of base pairs), their extraordinary abundance of repetitive sequences, and allopolyploidy (the presence in a nucleus of several related but not identical genomes) force us to think that only few basic plant species will undergo complete sequencing, whereas genome investigations in other species will follow the principles of comparative genomics. By the present time, sequencing of the Arabidopsis genome (125 Mbp) is completed and that of the rice genome (about 430 Mbp) is close to its end. Studying the genomes of other plants, including economically valuable ones, already began on the basis of these works. The peculiarities of plant genomes make extraordinarily important our detailed knowledge on plant chromosomes which, in its turn, calls for expansion of research in this direction and development of new chromosome technologies, including the DNA-sparing methods of high-resolution banding.     

Successful lateral transfer requires codon usage compatibility between foreign genes and recipient genomes

We present evidence supporting the notion that codon usage (CU) compatibility between foreign genes and recipient genomes is an important prerequisite to assess the selective advantage of imported functions, and therefore to increase the fixation probability of horizontal gene transfer (HGT) events. This contrasts with the current tendency in research to predict recent HGTs in prokaryotes by assuming that acquired genes generally display poor CU. By looking at the CU level (poor, typical, or rich) exhibited by putative xenologs still resembling their original CU, we found that most alien genes predominantly present typical CU immediately upon introgression, thereby suggesting that the role of CU amelioration in HGT has been overemphasized. In our strategy, we first scanned a representative set of 103 complete prokaryotic genomes for all pairs of candidate xenologs (exported/imported genes) displaying similar CU. We applied additional filtering criteria, including phylogenetic validations, to enhance the reliability of our predictions. Our approach makes no assumptions about the CU of foreign genes being typical or atypical within the recipient genome, thus providing a novel unbiased framework to study the evolutionary dynamics of HGT.     

Automatic target selection for structural genomics on eukaryotes

A central goal of structural genomics is to experimentally determine representative structures for all protein families. At least 14 structural genomics pilot projects are currently investigating the feasibility of high-throughput structure determination; the National Institutes of Health funded nine of these in the United States. Initiatives differ in the particular subset of "all families" on which they focus. At the NorthEast Structural Genomics consortium (NESG), we target eukaryotic protein domain families. The automatic target selection procedure has three aims: 1) identify all protein domain families from currently five entirely sequenced eukaryotic target organisms based on their sequence homology, 2) discard those families that can be modeled on the basis of structural information already present in the PDB, and 3) target representatives of the remaining families for structure determination. To guarantee that all members of one family share a common foldlike region, we had to begin by dissecting proteins into structural domain-like regions before clustering. Our hierarchical approach, CHOP, utilizing homology to PrISM, Pfam-A, and SWISS-PROT chopped the 103,796 eukaryotic proteins/ORFs into 247,222 fragments. Of these fragments, 122,999 appeared suitable targets that were grouped into >27,000 singletons and >18,000 multifragment clusters. Thus, our results suggested that it might be necessary to determine >40,000 structures to minimally cover the subset of five eukaryotic proteomes.