Adaptive diversification of vomeronasal receptor 1 genes in rodents

The vomeronasal receptor 1 (V1R) are believed to be pheromone receptors in rodents. Here we used computational methods to identify 95 and 62 new putative V1R genes from the draft rat and mouse genome sequence, respectively. The rat V1R repertoire consists of 11 subfamilies, 10 of which are shared with the mouse, while rat appears to lack the H and I subfamilies found in mouse and possesses one unique subfamily (M). The estimations of the relative divergence times suggest that many subfamilies originated after the split of rodents and primates. The analysis also reveals that these clusters underwent an expansion very close to the split of mouse and rat. In addition, maximum likelihood analysis showed that the nonsynonymous and synonymous rate ratio for most of these clusters was much higher than one, suggesting the role of positive selection in the diversification of these duplicated V1R genes. Because V1R are thought to mediate the process of signal transduction in response to pheromone detection, we speculate that the V1R genes have evolved under positive Darwinian selection to maintain the ability to discriminate between large and complex pheromonal mixtures.

Molecular evolution of growth hormone gene family in old world monkeys and hominoids

Growth hormone is a classic molecule in the study of the molecular clock hypothesis as it exhibits a relatively constant rate of evolution in most mammalian orders except primates and artiodactyls, where dramatically enhanced rate of evolution (25-50-fold) has been reported. The rapid evolution of primate growth hormone occurred after the divergence of tarsiers and simians, but before the separation of old world monkeys (OWM) from new world monkeys (NWM). Interestingly, this event of rapid sequence evolution coincided with multiple duplications of the growth hormone gene, suggesting gene duplication as a possible cause of the accelerated sequence evolution. Here we determined 21 different GH-like sequences from four species of OWM and hominoids. Combining with published sequences from OWM and hominoids, our analysis demonstrates that multiple gene duplications and several gene conversion events both occurred in the evolutionary history of this gene family in OWM/hominoids. The episode of recent duplications of CSH-like genes in gibbon is accompanied with rapid sequence evolution likely resulting from relaxation of purifying selection. GHN genes in both hominoids and OWM are under strong purifying selection. In contrast, CSH genes in both lineages are probably not. GHV genes in OWM and hominoids evolved at different evolutionary rates and underwent different selective constraints. Our results disclosed the complex history of the primate growth hormone gene family and raised intriguing questions on the consequences of these evolutionary events. © 2005 Elsevier B.V. All rights reserved.

Phylogeny and biogeography of the Petaurista philippensis complex (Rodentia : Sciuridae), inter- and intraspecific relationships inferred from molecular and morphometric analysis

With modified DNA extraction and Purification protocols, the complete cytochrome b gene sequences (1140 bp) were determined from degraded museum specimens. Molecular analysis and morphological examination of cranial characteristics of the giant flying squirrels of Petaurista philippensis complex (P. grandis, P. hainana, and P. yunanensis) and other Petaurista species yielded new insights into long-standing controversies in the Petaurista systematics. Patterns of genetic variations and morphological differences observed in this study indicate that P. hainana, P. albiventer, and P. yunanensis can be recognized as distinct species, and P. grandis and P. petaurista are conspecific populations. Phylogenetic relationships reconstructed by using parsimony, likelihood, and Bayesian methods reveal that, with P. leucogenys as the basal branch, all Petaurista groups formed two distinct clades. Petaurista philippensis, P. hainana, P. yunanensis, and P. albiventer are clustered in the same clade, while P. grandis shows a close relationship to P. petaurista. Deduced divergence times based on Bayesian analysis and the transversional substitution at the third codon suggest that the retreating of glaciers and upheavals or movements of tectonic plates in the Pliocene-Pleistocene were the major factors responsible for the present geographical distributions of Petaurista groups. (c) 2005 Elsevier Inc. All rights reserved.

Positive selection drives population differentiation in the skeletal genes in modern humans

During the course of evolution, the human skeletal system has evolved rapidly leading to an incredible array of phenotypic diversity, including variations in height and bone mineral density. However, the genetic basis of this phenotypic diversity and the relatively rapid tempo of evolution have remained largely undocumented. Here, we discover that skeletal genes exhibit a significantly greater level of population differentiation among humans compared with other genes in the genome. The pattern is exceptionally evident at amino acid-altering sites within these genes. Divergence is greater between Africans and both Europeans and East Asians. In contrast, relatively weak differentiation is observed between Europeans and East Asians. SNPs with higher levels of differentiation have correspondingly higher derived allele frequencies in Europeans and East Asians. Thus, it appears that positive selection has operated on skeletal genes in the non-African populations and this may have been initiated with the human colonization of Eurasia. In conclusion, we provide genetic evidence supporting the rapid evolution of the human skeletal system and the associated diversity of phenotypes.

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.

A multi-locus phylogeny of Nectogalini shrews and influences of the paleoclimate on speciation and evolution

Nectogaline shrews are a major component of the small mammalian fauna of Europe and Asia, and are notable for their diverse ecology, including utilization of aquatic habitats. So far, molecular phylogenetic analyses including nectogaline species have been unable to infer a well-resolved, well-supported phylogeny, thus limiting the power of comparative evolutionary and ecological analyses of the group. Here, we employ Bayesian phylogenetic analyses of eight mitochondrial and three nuclear genes to infer the phylogenetic relationships of nectogaline shrews. We subsequently use this phylogeny to assess the genetic diversity within the genus Episoriculus, and determine whether adaptation to aquatic habitats evolved independently multiple times. Moreover, we both analyze the fossil record and employ Bayesian relaxed clock divergence dating analyses of DNA to assess the impact of historical global climate change on the biogeography of Nectogalini. We infer strong support for the polyphyly of the genus Episoriculus. We also find strong evidence that the ability to heavily utilize aquatic habitats evolved independently in both Neomys and Chimarrogale + Nectogale lineages. Our Bayesian molecular divergence analysis suggests that the early history of Nectogalini is characterized by a rapid radiation at the Miocene/Pliocene boundary, thus potentially explaining the lack of resolution at the base of the tree. Finally, we find evidence that nectogalines once inhabited northern latitudes, but the global cooling and desiccating events at the Miocene/Pliocene and Pliocene/Pleistocene boundaries and Pleistocene glaciation resulted in the migration of most Nectogalini lineages to their present day southern distribution. (C) 2010 Elsevier Inc. All rights reserved.

Rapid and parallel chromosomal number reductions in muntjac deer inferred from mitochondrial DNA phylogeny

Muntjac deer (Muntiacinae, Cervidae) are of great interest in evolutionary studies because of their dramatic chromosome variations and recent discoveries of several new species. In this paper, we analyze the evolution of karyotypes of muntjac deer in the context of a phylogeny which is based on 1,844-bp mitochondrial DNA sequences of seven generally recognized species in the muntjac subfamily. The phylogenetic results support the hypothesis that karyotypic evolution in muntjac deer has proceeded via reduction in diploid number. However, the reduction in number is not always linear, i.e., not strictly following the order: 46-->14/13-->8/9-->6/7. For example, Muntiacus muntjak (2n = 6/7) shares a common ancestor with Muntiacus feae (2n = 13/14), which indicates that its karyotype was derived in parallel with M. feae's from an ancestral karyotype of 2n greater than or equal to 13/14. The newly discovered giant muntjac (Muntiacus vuquangensis) may represent another pa;allel reduction lineage from the ancestral 2n = 46 karyotype. Our phylogenetic results indicate that the giant muntjac is relatively closer to Muntiacus reevesi than to other muntjacs and may be placed in the genus Muntiacus. Analyses of sequence divergence reveal that the rate of change in chromosome number in muntjac deer is one of the fastest in vertebrates. Within the muntjac subfamily, the fastest evolutionary rate is found in the Fea's lineage, in which two species with different karyotypes diverged in around 0.5 Myr.

Origin and evolution of new exons in the rodent zinc finger protein 39 gene

The origin of new structures and functions is an important process in evolution. In the past decades, we have obtained some preliminary knowledge of the origin and evolution of new genes. However, as the basic unit of genes, the origin and evolution of exons remain unclear. Because young exons retain the footprints of origination, they can be good materials for studying origin and evolution of new exons. In this paper, we report two young exons in a zinc finger protein gene of rodents. Since they are unique sequences in mouse and rat genome and no homologous sequences were found in the orthologous genes of human and pig, the young exons might originate after the divergence of primates and rodents through exonization of intronic sequences. Strong positive selection was detected in the new exons between mouse and rat, suggesting that these exons have undergone significant functional divergence after the separation of the two species. On the other hand, population genetics data of mouse demonstrate that the new exons have been subject to functional constraint, indicating an important function of the new exons in mouse. Functional analyses suggest that these new exons encode a nuclear localization signal peptide, which may mediate new ways of nuclear protein transport. To our knowledge, this is the first example of the origin and evolution of young exons.

The Genomes of Oryza sativa: A history of duplications

We report improved whole-genome shotgun sequences for the genomes of indica and japonica rice, both with multimegabase contiguity, or almost 1,000-fold improvement over the drafts of 2002. Tested against a nonredundant collection of 19,079 full-length cDNAs, 97.7% of the genes are aligned, without fragmentation, to the mapped superscaffolds of one or the other genome. We introduce a gene identification procedure for plants that does not rely on similarity to known genes to remove erroneous predictions resulting from transposable elements. Using the available EST data to adjust for residual errors in the predictions, the estimated gene count is at least 38,000 - 40,000. Only 2% - 3% of the genes are unique to any one subspecies, comparable to the amount of sequence that might still be missing. Despite this lack of variation in gene content, there is enormous variation in the intergenic regions. At least a quarter of the two sequences could not be aligned, and where they could be aligned, single nucleotide polymorphism ( SNP) rates varied from as little as 3.0 SNP/kb in the coding regions to 27.6 SNP/kb in the transposable elements. A more inclusive new approach for analyzing duplication history is introduced here. It reveals an ancient whole-genome duplication, a recent segmental duplication on Chromosomes 11 and 12, and massive ongoing individual gene duplications. We find 18 distinct pairs of duplicated segments that cover 65.7% of the genome; 17 of these pairs date back to a common time before the divergence of the grasses. More important, ongoing individual gene duplications provide a never-ending source of raw material for gene genesis and are major contributors to the differences between members of the grass family.

Duplication and independent selection of cell-wall invertase genes GIF1 and OsCIN1 during rice evolution and domestication

Background: Various evolutionary models have been proposed to interpret the fate of paralogous duplicates, which provides substrates on which evolution selection could act. In particular, domestication, as a special selection, has played important role in crop cultivation with divergence of many genes controlling important agronomic traits. Recent studies have indicated that a pair of duplicate genes was often sub-functionalized from their ancestral functions held by the parental genes. We previously demonstrated that the rice cell-wall invertase (CWI) gene GIF1 that plays an important role in the grain-filling process was most likely subjected to domestication selection in the promoter region. Here, we report that GIF1 and another CWI gene OsCIN1 constitute a pair of duplicate genes with differentiated expression and function through independent selection. Results: Through synteny analysis, we show that GIF1 and another cell-wall invertase gene OsCIN1 were paralogues derived from a segmental duplication originated during genome duplication of grasses. Results based on analyses of population genetics and gene phylogenetic tree of 25 cultivars and 25 wild rice sequences demonstrated that OsCIN1 was also artificially selected during rice domestication with a fixed mutation in the coding region, in contrast to GIF1 that was selected in the promoter region. GIF1 and OsCIN1 have evolved into different expression patterns and probable different kinetics parameters of enzymatic activity with the latter displaying less enzymatic activity. Overexpression of GIF1 and OsCIN1 also resulted in different phenotypes, suggesting that OsCIN1 might regulate other unrecognized biological process. Conclusion: How gene duplication and divergence contribute to genetic novelty and morphological adaptation has been an interesting issue to geneticists and biologists. Our discovery that the duplicated pair of GIF1 and OsCIN1 has experiencedsub-functionalization implies that selection could act independently on each duplicate towards different functional specificity, which provides a vivid example for evolution of genetic novelties in a model crop. Our results also further support the established hypothesis that gene duplication with sub-functionalization could be one solution for genetic adaptive conflict.

Recombination Yet Inefficient Selection along the Drosophila melanogaster Subgroup's Fourth Chromosome

A central goal of evolutionary genetics is an understanding of the forces responsible for the observed variation, both within and between species. Theoretical and empirical work have demonstrated that genetic recombination contributes to this variation by breaking down linkage between nucleotide sites, thus allowing them to behave independently and for selective forces to act efficiently on them. The Drosophila fourth chromosome, which is believed to experience no-or very low-rates of recombination has been an important model for investigating these effects. Despite previous efforts, central questions regarding the extent of recombination and the predominant modes of selection acting on it remain open. In order to more comprehensively test hypotheses regarding recombination and its potential influence on selection along the fourth chromosome, we have resequenced regions from most of its genes from Drosophila melanogaster, D. simulans, and D. yakuba. These data, along with available outgroup sequence, demonstrate that recombination is low but significantly greater than zero for the three species. Despite there being recombination, there is strong evidence that its frequency is low enough to have rendered selection relatively inefficient. The signatures of relaxed constraint can be detected at both the level of polymorphism and divergence.

The sequence and de novo assembly of the giant panda genome

Using next-generation sequencing technology alone, we have successfully generated and assembled a draft sequence of the giant panda genome. The assembled contigs (2.25 gigabases (Gb)) cover approximately 94% of the whole genome, and the remaining gaps (0.05 Gb) seem to contain carnivore-specific repeats and tandem repeats. Comparisons with the dog and human showed that the panda genome has a lower divergence rate. The assessment of panda genes potentially underlying some of its unique traits indicated that its bamboo diet might be more dependent on its gut microbiome than its own genetic composition. We also identified more than 2.7 million heterozygous single nucleotide polymorphisms in the diploid genome. Our data and analyses provide a foundation for promoting mammalian genetic research, and demonstrate the feasibility for using next-generation sequencing technologies for accurate, cost-effective and rapid de novo assembly of large eukaryotic genomes.

Natural selection on EPAS1 (HIF2 alpha) associated with low hemoglobin concentration in Tibetan highlanders

By impairing both function and survival, the severe reduction in oxygen availability associated with high-altitude environments is likely to act as an agent of natural selection. We used genomic and candidate gene approaches to search for evidence of such genetic selection. First, a genome-wide allelic differentiation scan (GWADS) comparing indigenous highlanders of the Tibetan Plateau (3,200 3,500 m) with closely related lowland Han revealed a genome-wide significant divergence across eight SNPs located near EPAS1. This gene encodes the transcription factor HIF2 alpha, which stimulates production of red blood cells and thus increases the concentration of hemoglobin in blood. Second, in a separate cohort of Tibetans residing at 4,200 m, we identified 31 EPAS1 SNPs in high linkage disequilibrium that correlated significantly with hemoglobin concentration. The sex-adjusted hemoglobin concentration was, on average, 0.8 g/dL lower in the major allele homozygotes compared with the heterozygotes. These findings were replicated in a third cohort of Tibetans residing at 4,300 m. The alleles associating with lower hemoglobin concentrations were correlated with the signal from the GWADS study and were observed at greatly elevated frequencies in the Tibetan cohorts compared with the Han. High hemoglobin concentrations are a cardinal feature of chronic mountain sickness offering one plausible mechanism for selection. Alternatively, as EPAS1 is pleiotropic in its effects, selection may have operated on some other aspect of the phenotype. Whichever of these explanations is correct, the evidence for genetic selection at the EPAS1 locus from the GWADS study is supported by the replicated studies associating function with the allelic variants.

On the use of combined compact scheme for numerical solution of the two-layer oceanic shallow water equations

Many types of oceanic physical phenomena have a wide range in both space and time. In general, simplified models, such as shallow water model, are used to describe these oceanic motions. The shallow water equations are widely applied in various oceanic and atmospheric extents. By using the two-layer shallow water equations, the stratification effects can be considered too. In this research, the sixth-order combined compact method is investigated and numerically implemented as a high-order method to solve the two-layer shallow water equations. The second-order centered, fourth-order compact and sixth-order super compact finite difference methods are also used to spatial differencing of the equations. The first part of the present work is devoted to accuracy assessment of the sixth-order super compact finite difference method (SCFDM) and the sixth-order combined compact finite difference method (CCFDM) for spatial differencing of the linearized two-layer shallow water equations on the Arakawa's A-E and Randall's Z numerical grids. Two general discrete dispersion relations on different numerical grids, for inertia-gravity and Rossby waves, are derived. These general relations can be used for evaluation of the performance of any desired numerical scheme. For both inertia-gravity and Rossby waves, minimum error generally occurs on Z grid using either the sixth-order SCFDM or CCFDM methods. For the Randall's Z grid, the sixth-order CCFDM exhibits a substantial improvement , for the frequency of the barotropic and baroclinic modes of the linear inertia-gravity waves of the two layer shallow water model, over the sixth-order SCFDM. For the Rossby waves, the sixth-order SCFDM shows improvement, for the barotropic and baroclinic modes, over the sixth-order CCFDM method except on Arakawa's C grid. In the second part of the present work, the sixth-order CCFDM method is used to solve the one-layer and two-layer shallow water equations in their nonlinear form. In one-layer model with periodic boundaries, the performance of the methods for mass conservation is compared. The results show high accuracy of the sixth-order CCFDM method to simulate a complex flow field. Furthermore, to evaluate the performance of the method in a non-periodic domain the sixth-order CCFDM is applied to spatial differencing of vorticity-divergence-mass representation of one-layer shallow water equations to solve a wind-driven current problem with no-slip boundary conditions. The results show good agreement with published works. Finally, the performance of different schemes for spatial differencing of two-layer shallow water equations on Z grid with periodic boundaries is investigated. Results illustrate the high accuracy of combined compact method.

Investigation of polymorphism by PCR-RFLP method in Cobia fish Rachycentron canadum in the Persian Gulf and Oman Sea

Cobia is a native fish species in Iranian waters in the Persian Gulf and Sea of Oman and has a good internal and foreign market. This fish is a fast growing species and for this reason Iranian Fisheries is considering to go for it culture practices. To go for any utilization such as fishing from wild stocks or culture activities, needs a better understanding of its peculiarities and genetic characteristics of its natural resources. Therefore, this project was discribed and conducted. In this investigation, cuts 2 or 3 cm of fin tissue of  specimen of Cobia obtained from Sistan and Bluchestan, Hormozgan, Bushehr and Khuzestan water provinces, were collected. DNA was extracted by Phenol-chlorophorm method and produced PCR product in length of 1060 and 1450 base pair of two mitochondrial genes COI and NADH2. Using 13 cutting enzymes (4 enzymes were subscriber for both of genes), 205 base pair (from 2510 base pair, equal with %3.8 from gene regains) were directly investigated. But binding patterns of enzymatic digestion of PCR products of both COI and ND genes from electrophoresis were monomorph in all samples and no polymorphism was observed. This may be attributed to the unsuitable choice of COI and ND2 genes for showing of intra specific divergence. But in general non-existence of genetic diversity or noticeable decrease of that among individuals has been reported in regions were fish migration exist and they can freely move between two regions. Therefore, non-observation of polymorphism in the study area might be the case and indicates represents the area. On the other hand, some scientists believe that the distributions of populations in different regions are greatly affected by environmental and physical and ecological factors. Althoug Cobia is a migratory fish, but with regard to the fact that the environmental conditions are different (specially temperature and salinity) between east and west of Persian Gulf and Oman sea, there is a possibility that different genetic groups of this species exist in the regions. Of course It is clear that using more samples and enzymes from other genetically regions could produce better results. Since none of the two investigated genes didn’t show genetic divergence or polymorphism amongst the individuals of one region or between different regions, therefore, statistic analysis for estimating of haplotype diversity or nucleotide diversity and drawing of relationship tree among individuals using available softwares was not possible.