Molecular variation of the shuttle craft and Lim3 genes, controlling the development of the nervous system, in a natural Drosophila melanogaster population

Comparative polymorphism of the first exon and first intron of the shuttle craft (stc) and Lim3 genes and their putative regulatory 5'-flanking sequences was analyzed using 20 sequenced natural alleles. A comparison of the stc and Lim3 genes showed that the extent of polymorphism was similar in their introns and corresponded to the variation level characteristic of Drosophila melanogaster, while the putative regulatory region and first intron of the stc gene proved to be more variable than the corresponding regions of the Lim3 gene. Since the genes under study occurred on the same chromosomes isolated from one population and were close together in a region having a high recombination rate, the difference in the extent of polymorphism between the regulatory and coding regions was explained by individual characteristics of each gene. The results made it possible to assume that the extent of polymorphism of the coding gene regions is maintained by balancing selection.     

Gene Promoter Evolution Targets the Center of the Human Protein Interaction Network

Assessing the contribution of promoters and coding sequences to gene evolution is an important step toward discovering the major genetic determinants of human evolution. Many specific examples have revealed the evolutionary importance of cis-regulatory regions. However, the relative contribution of regulatory and coding regions to the evolutionary process and whether systemic factors differentially influence their evolution remains unclear. To address these questions, we carried out an analysis at the genome scale to identify signatures of positive selection in human proximal promoters. Next, we examined whether genes with positively selected promoters (Prom⁺ genes) show systemic differences with respect to a set of genes with positively selected protein-coding regions (Cod⁺ genes). We found that the number of genes in each set was not significantly different (8.1% and 8.5%, respectively). Furthermore, a functional analysis showed that, in both cases, positive selection affects almost all biological processes and only a few genes of each group are located in enriched categories, indicating that promoters and coding regions are not evolutionarily specialized with respect to gene function. On the other hand, we show that the topology of the human protein network has a different influence on the molecular evolution of proximal promoters and coding regions. Notably, Prom⁺ genes have an unexpectedly high centrality when compared with a reference distribution (P = 0.008, for Eigenvalue centrality). Moreover, the frequency of Prom⁺ genes increases from the periphery to the center of the protein network (P = 0.02, for the logistic regression coefficient). This means that gene centrality does not constrain the evolution of proximal promoters, unlike the case with coding regions, and further indicates that the evolution of proximal promoters is more efficient in the center of the protein network than in the periphery. These results show that proximal promoters have had a systemic contribution to human evolution by increasing the participation of central genes in the evolutionary process.     

Constraints on Intron Evolution in the Gene Encoding the Myosin Alkali Light Chain in Drosophila

Interspecific comparisons of intron sequences reveal conserved blocks of invariant nucleotides and several other departures from the strictly neutral model of molecular evolution. To distinguish the past action of evolutionary forces in introns known to have regulatory information, we examined nucleotide sequence variation at 991 sites in a random sample of 16 Drosophila melanogaster alleles of the gene encoding the myosin alkali light chain (Mlc1). The Mlc1 gene of D. melanogaster encodes two MLC1 isoforms via developmentally regulated alternative pre-mRNA splicing. Analyses of these data reveal that introns 4 and 5, which flank the alternatively spliced exon 5, have reduced levels of both intraspecific polymorphism and interspecific divergence relative to intron 3. No polymorphism was observed in any of the exons examined in D. melanogaster. A genealogical analysis clearly demonstrates the occurrence of intragenic recombination in the ancestral history of Mlc1. Recombination events are estimated to be 13 times more likely than mutation events over the span of the sequenced region. Although there is little evidence for pairwise linkage disequilibrium in the Mlc1 region, higher order disequilibrium does seem to be present in the 5' half of the portion of the gene that was examined. Predictions of the folding free energy of the pre-mRNA reveal that sampled alleles have a significantly higher (less stable) free energy than do randomly permuted sequences. These results are consistent with the hypothesis that introns surrounding an alternatively spliced exon are subjected to additional constraints, perhaps due to specific aspects of secondary structure required for appropriate splicing of the pre-mRNA molecule.     

Nucleotide variation of the duplicated amylase genes in Drosophila kikkawai

We examined levels and patterns of the nucleotide polymorphism of the Amylase genes with a head-to-head duplication in Drosophila kikkawai. The levels of variation in D. kikkawai were comparable to those in Drosophila melanogaster. Tajima's test, Fu and Li's test, HKA test, and MK test did not show significant departure from neutrality. We found an excess of replacement changes in the within-locus class, representing polymorphism in one of the duplicated genes, compared with the between-locus class, representing polymorphism shared between the duplicated genes. Most replacement changes in the within-locus class were singletons. These results suggest that most replacement changes are deleterious. A contrasting evolutionary pattern, involving concerted evolution in the coding regions but differential evolution in the 5'-flanking regions, was observed. However, unlike the duplicated Amy genes of D. melanogaster, the coding regions of the duplicated genes in D. kikkawai tended to diverge. Using Ohta's model of the small multigene family, we found that recombination (interchromosomal equal crossing-over) rate was one order higher than gene conversion (unequal crossing-over) rate, resulting in a considerable but incomplete homogenization of the duplicated coding regions. Linkage disequilibria were found in the intron as well as within and around the regulatory cis-element sequences of one of the duplicated genes (Amy1). The possible causes of these linkage disequilibria were discussed.     

Gene duplication,allelic diversity,selection processes and adaptive value of MHC class II <Emphasis Type="Italic">DRB</Emphasis> genes of the bank vole, <Emphasis Type="Italic">Clethrionomys glareolus</Emphasis>

The generation and maintenance of allelic polymorphism in genes of the major histocompatibility complex (MHC) is a central issue in evolutionary genetics. Recently, the focus has changed from ex situ to in situ populations to understand the mechanisms that determine adaptive MHC polymorphism under natural selection. Birth-and-death evolution and gene conversion events are considered to generate sequence diversity in MHC genes, which subsequently is maintained by balancing selection through parasites. The ongoing arms race between the host and parasites leads to an adaptive selection pressure upon the MHC, evident in high rates of non-synonymous vs synonymous substitution rates. We characterised the MHC class II DRB exon 2 of free living bank voles, Clethrionomys glareolus by single-strand conformation polymorphism and direct sequencing. Unlike other arvicolid species, the DRB locus of the bank vole is at least quadruplicated. No evidence for gene conversion events in the Clgl-DRB sequences was observed. We found not only high allelic polymorphism with 26 alleles in 36 individuals but also high rates of silent polymorphism. Exceptional for MHC class II genes is a purifying selection pressure upon the majority of MHC-DRB sequences. Further, we analysed the association between certain DRB alleles and the parasite burden with gastrointestinal trichostrongyle nematodes Heligmosomum mixtum and Heligmosomoides glareoli and found significant quality differences between specific alleles with respect to infection intensity. Our findings suggest a snapshot in an evolutionary process of ongoing birth-and-death evolution. One allele cluster has lost its function and is already silenced, another is loosing its adaptive value in terms of gastrointestinal nematode resistance, while a third group of alleles indicates all signs of classical functional MHC alleles.     

Molecular variation of the shuttle craft and Lim3 genes, controlling the development of the nervous system, in a natural Drosophila melanogaster population

Comparative polymorphism of the first exon and first intron of the shuttle craft (stc) and Lim3 genes and their putative regulatory 5′-flanking sequences was analyzed using 20 sequenced natural alleles. A comparison of the stc and Lim3 genes showed that the extent of polymorphism was similar in their introns and corresponded to the variation level characteristic of Drosophila melanogaster, while the putative regulatory region and first intron of the stc gene proved to be more variable than the corresponding regions of the Lim3 gene. Since the genes under study occurred on the same chromosomes isolated from one population and were close together in a region having a high recombination rate, the difference in the extent of polymorphism between the regulatory and coding regions was explained by individual characteristics of each gene. The results made it possible to assume that the extent of polymorphism of the coding gene regions is maintained by balancing selection.     

Neutral mutation as the source of genetic variation in life history traits

The mechanism underlying the maintenance of adaptive genetic variation is a long-standing question in evolutionary genetics. There are two concepts (mutation-selection balance and balancing selection) which are based on the phenotypic differences between alleles. Mutation - selection balance and balancing selection cannot properly explain the process of gene substitution, i.e. the molecular evolution of quantitative trait loci affecting fitness. I assume that such loci have non-essential functions (small effects on fitness), and that they have the potential to evolve into new functions and acquire new adaptations. Here I show that a high amount of neutral polymorphism at these loci can exist in real populations. Consistent with this, I propose a hypothesis for the maintenance of genetic variation in life history traits which can be efficient for the fixation of alleles with very small selective advantage. The hypothesis is based on neutral polymorphism at quantitative trait loci and both neutral and adaptive gene substitutions. The model of neutral - adaptive conversion (NAC) assumes that neutral alleles are not neutral indefinitely, and that in specific and very rare situations phenotypic (relative fitness) differences between them can appear. In this paper I focus on NAC due to phenotypic plasticity of neutral alleles. The important evolutionary consequence of NAC could be the increased adaptive potential of a population. Loci responsible for adaptation should be fast evolving genes with minimally discernible phenotypic effects, and the recent discovery of genes with such characteristics implicates them as suitable candidates for loci involved in adaptation.     

Interpretation of Genomic Variants Using a Unified Biological Network Approach

The decreasing cost of sequencing is leading to a growing repertoire of personal genomes. However, we are lagging behind in understanding the functional consequences of the millions of variants obtained from sequencing. Global system-wide effects of variants in coding genes are particularly poorly understood. It is known that while variants in some genes can lead to diseases, complete disruption of other genes, called ‘loss-of-function tolerant’, is possible with no obvious effect. Here, we build a systems-based classifier to quantitatively estimate the global perturbation caused by deleterious mutations in each gene. We first survey the degree to which gene centrality in various individual networks and a unified ‘Multinet’ correlates with the tolerance to loss-of-function mutations and evolutionary conservation. We find that functionally significant and highly conserved genes tend to be more central in physical protein-protein and regulatory networks. However, this is not the case for metabolic pathways, where the highly central genes have more duplicated copies and are more tolerant to loss-of-function mutations. Integration of three-dimensional protein structures reveals that the correlation with centrality in the protein-protein interaction network is also seen in terms of the number of interaction interfaces used. Finally, combining all the network and evolutionary properties allows us to build a classifier distinguishing functionally essential and loss-of-function tolerant genes with higher accuracy (AUC = 0.91) than any individual property. Application of the classifier to the whole genome shows its strong potential for interpretation of variants involved in Mendelian diseases and in complex disorders probed by genome-wide association studies.     

Genomics of natural populations: Evolutionary forces that establish and maintain gene arrangements in Drosophila pseudoobscura

The evolution of complex traits in heterogeneous environments may shape the order of genes within chromosomes. Drosophila pseudoobscura has a rich gene arrangement polymorphism that allows one to test evolutionary genetic hypotheses about how chromosomal inversions are established in populations. D. pseudoobscura has >30 gene arrangements on a single chromosome that were generated through a series of overlapping inversion mutations with >10 inversions with appreciable frequencies and wide geographic distributions. This study analyses the genomic sequences of 54 strains of Drosophila pseudoobscura that carry one of six different chromosomal arrangements to test whether (i) genetic drift, (ii) hitchhiking with an adaptive allele, (iii) direct effects of inversions to create gene disruptions caused by breakpoints, or (iv) indirect effects of inversions in limiting the formation of recombinant gametes are responsible for the establishment of new gene arrangements. We found that the inversion events do not disrupt the structure of protein coding genes at the breakpoints. Population genetic analyses of 2,669 protein coding genes identified 277 outlier loci harbouring elevated frequencies of arrangement‐specific derived alleles. Significant linkage disequilibrium occurs among distant loci interspersed between regions with low levels of association indicating that distant allelic combinations are held together despite shared polymorphism among arrangements. Outlier genes showing evidence of genetic differentiation between arrangements are enriched for sensory perception and detoxification genes. The data presented here support the indirect effect of inversion hypothesis where chromosomal inversions are favoured because they maintain linked associations among multilocus allelic combinations among different arrangements.     

陕北白绒山羊DRB1基因外显子2遗传变异分析

本研究通过对123只陕北白绒山羊DRB1基因外显子2的遗传变异分析,旨在获得陕北白绒山羊DRB1基因的多态性及变异信息,为山羊抗病基因的挖掘研究提供基础资料。本研究共获得6条陕北白绒山羊DRB1基因外显子2序列,其中4条为首次发现。生物信息学分析表明DRB1位点具有较高的多态性,6条等位基因可能起源于2个祖先基因。在长期的进化过程中,DRB1位点受到了明显的选择压力作用,这种选择作用有助于陕北白绒山羊对当地气候的适应。蛋白质结构的预测证实了DRB1*1与其它等位基因间的差异性,说明核苷酸变异可能会引起蛋白质结构的改变,最终可能影响宿主对病原体的免疫应答。本次对陕北白绒山羊DRB1基因多态性的调查与分析有助于筛选疾病抗性和易感性MHC (Major histocompatibility complex)候选基因,进而可加速绒山羊抗病品系的改良与培育进程。     

Characterization and Evolution of MHC Class II B Genes in Ardeid Birds

Major histocompatibility complex (MHC) is a multi-gene family that is very suitable to investigate a wide range of open questions in evolutionary ecology. In this study, we characterized two expressed MHC class II B genes (DAB1 and DAB2) in the Grey Heron (Aves: Ardea cinerea). We further developed the primer pairs to amplify and sequence two MHC class II B loci in ten ardeid birds. Phylogenetic analysis revealed that different parts of the genes showed different evolutionary patterns. The exon 2 sequences tended to cluster two gene-specific lineages. In each lineage, exon 2 sequences from several species showed closer relationships than sequences within species, and two shared identical alleles were found between species (Egretta sacra and Nycticorax nycticorax; Egretta garzetta and Bubulcus ibis), supporting the hypothesis of trans-species polymorphism. In contrast, the species-specific intron 2 plus partial exon 3 tree suggested that DAB1 and DAB2 were subject to concerted evolution. GENECONV analyses showed the gene exchange played an important role in the ardeid MHC evolution.     

Evolution of proteins and gene expression levels are coupled in Drosophila and are independently associated with mRNA abundance, protein length, and number of protein-protein interactions

Organismic evolution requires that variation at distinct hierarchical levels and attributes be coherently integrated, often in the face of disparate environmental and genetic pressures. A central part of the evolutionary analysis of biological systems remains to decipher the causal connections between organism-wide (or genome-wide) attributes (e.g., mRNA abundance, protein length, codon bias, recombination rate, genomic position, mutation rate, etc) as well as their role-together with mutation, selection, and genetic drift-in shaping patterns of evolutionary variation in any of the attributes themselves. Here we combine genome-wide evolutionary analysis of protein and gene expression data to highlight fundamental relationships among genomic attributes and their associations with the evolution of both protein sequences and gene expression levels. Our results show that protein divergence is positively coupled with both gene expression polymorphism and divergence. We show moreover that although the number of protein-protein interactions in Drosophila is negatively associated with protein divergence as well as gene expression polymorphism and divergence, protein-protein interactions cannot account for the observed coupling between regulatory and structural evolution. Furthermore, we show that proteins with higher rates of amino acid substitutions tend to have larger sizes and tend to be expressed at lower mRNA abundances, whereas genes with higher levels of gene expression divergence and polymorphism tend to have shorter sizes and tend to be expressed at higher mRNA abundances. Finally, we show that protein length is negatively associated with both number of protein-protein interactions and mRNA abundance and that interacting proteins in Drosophila show similar amounts of divergence. We suggest that protein sequences and gene expression are subjected to similar evolutionary dynamics, possibly because of similarity in the fitness effect (i.e., strength of stabilizing selection) of disruptions in a gene's protein sequence or its mRNA expression. We conclude that, as more and better data accumulate, understanding the causal connections among biological traits and how they are integrated over time to constrain or promote structural and regulatory evolution may finally become possible.     

The utility of QTL-Linked markers to detect selective sweeps in natural populations--a case study of the EDA gene and a linked marker in threespine stickleback

Sequence polymorphisms in coding genes and variability in quantitative trait loci (QTL)-linked markers can be used to uncover the evolutionary mechanisms of traits involved in adaptive processes. We studied sequence variation in the EDA gene and allelic variation in 18 microsatellites - one of which (Gac4174) is linked with the EDA QTL - in low, partially and completely plated morphs from eight threespine stickleback European populations. The results agree with previous studies in that EDA polymorphism is closely related to plate number variation: EDA sequences grouped populations into low and completely plated morphs, whereas microsatellites failed to do so. Furthermore, partially plated fish were heterozygous with respect to the distinctive EDA alleles for completely and low plated morphs, indicating that completely plated morph alleles are not entirely dominant in controlling the expression of lateral plate number. An examination of population differentiation in plate number with quantitative genetic methods revealed that the degree of differentiation exceeded that expected from genetic drift alone (Q(ST) > F(ST)). Our results support the adaptive genetic differentiation of plate morphs and the view that distinctive EDA gene polymorphism occurs in similar sites across the distribution range of this species. Yet, allele frequency differentiation in the Gac4174 microsatellite locus, informative in experimental crosses for plate number variation, did not differ from that of neutral markers and, was therefore unable to detect the signature of natural selection responsible for population divergence.     

Evidence for the contribution of LTR retrotransposons to C. elegans gene evolution

LTR retrotransposons may be important contributors to host gene evolution because they contain regulatory and coding signals. In an effort to assess the possible contribution of LTR retrotransposons to C. elegans gene evolution, we searched upstream and downstream of LTR retrotransposon sequences for the presence of predicted genes. Sixty-three percent of LTR retrotransposon sequences (79/124) are located within 1 kb of a gene or within gene boundaries. Most gene-retrotransposon associations were located along the chromosome arms. Our results are consistent with the hypothesis that LTR retrotransposons have contributed to the structural and/or regulatory evolution of genes in C. elegans.     

Hierarchical evolution of animal body plans

An open question in animal evolution is why the phylum- and superphylum-level body plans have changed so little, while the class- and family-level body plans have changed so greatly since the early Cambrian. Davidson and Erwin (Davidson and Erwin, 2006; Erwin and Davidson, 2009) proposed that the hierarchical structure of gene regulatory networks leads to different observed evolutionary rates for terminal properties of the body plan versus major aspects of body plan morphology. Here, we calculated the speed of evolution of genes in these gene regulatory networks. We found that the genes which determine the phylum and superphylum characters evolve slowly, while those genes which determine the classes, families, and speciation evolve more rapidly. This result furnishes genetic support to the hypothesis that the hierarchical structure of developmental regulatory networks provides an organizing structure which guides the evolution of aspects of the body plan.