Search

High-quality, chromosome-scale genomes are essential for genomic analyses. Analyses, including 3D genomics, epigenetics, and comparative genomics rely on a high-quality genome assembly, which is often accomplished with the assistance of Hi-C data. Curation of genomes reveal that current Hi-C-assisted scaffolding algorithms either generate ordering and orientation errors or fail to assemble high-quality chromosome-level scaffolds. Here, we offer the software Puzzle Hi-C, which uses Hi-C reads to accurately assign contigs or scaffolds to chromosomes. Puzzle Hi-C uses the triangle region instead of the square region to count interactions in a Hi-C heatmap. This strategy dramatically diminishes scaffolding interference caused by long-range interactions. This software also introduces a dynamic, triangle window strategy during assembly. Initially small, the window expands with interactions to produce more effective clustering. Puzzle Hi-C outperforms available scaffolding tools.

Caenophidian snakes include the file snake genus Acrochordus and advanced colubroidean snakes that radiated mainly during the Neogene. Although caenophidian snakes are a well-supported clade, their inferred affinities, based either on molecular or morphological data, remain poorly known or controversial. Here, we provide an expanded molecular phylogenetic analysis of Caenophidia and use three non-parametric measures of support-Shimodaira-Hasegawa-Like test (SHL), Felsentein (FBP) and transfer (TBE) bootstrap measures-to evaluate the robustness of each clade in the molecular tree. That very different alternative support values are common suggests that results based on only one support value should be viewed with caution. Using a scheme to combine support values, we find 20.9% of the 1265 clades comprising the inferred caenophidian tree are unambiguously supported by both SHL and FBP values, while almost 37% are unsupported or ambiguously supported, revealing the substantial extent of phylogenetic problems within Caenophidia. Combined FBP/TBE support values show similar results, while SHL/TBE result in slightly higher combined values. We consider key morphological attributes of colubroidean cranial, vertebral and hemipenial anatomy and provide additional morphological evidence supporting the clades Colubroides, Colubriformes, and Endoglyptodonta. We review and revise the relevant caenophidian fossil record and provide a time-calibrated tree derived from our molecular data to discuss the main cladogenetic events that resulted in present-day patterns of caenophidian diversification. Our results suggest that all extant families of Colubroidea and Elapoidea composing the present-day endoglyptodont fauna originated rapidly within the early Oligocene-between approximately 33 and 28 Mya-following the major terrestrial faunal turnover known as the "Grande Coupure" and associated with the overall climate shift at the Eocene-Oligocene boundary. Our results further suggest that the caenophidian radiation originated within the Caenozoic, with the divergence between Colubroides and Acrochordidae occurring in the early Eocene, at ~ 56 Mya.

[This corrects the article DOI: 10.1371/journal.pone.0216148.].

Red Junglefowls (RJFs), the wild progenitor of modern day chickens (DCs), are believed to be in genetic endangerment due to introgression of domestic genes through opportunistic matings with domestic or feral chickens. Previous studies from India reported rare hybridization of RJFs in the wild. However, RJF population genetic structure, pattern of gene flow and their admixture with DC populations are poorly understood at the landscape level. We conducted this study with a large sample size, covering the predicted natural distribution range of RJFs in India. We documented strong evidence of directional gene flow from DCs to free-ranging wild RJFs, with the Northeastern RJF population exhibiting the most genetic variants in their nuclear and mitochondrial genomes, indicating it to be the ancestral population from which early radiation may have occurred. The results provide evidence that landscape features do not act as a barrier to gene flow and the distribution pattern could not be explored due to physical sharing or exchange of wild birds in the past when forests were continuous across RJF range in India.

Reverse transcription-quantitative polymerase chain reaction (RT-qPCR) is a powerful technique for examining gene expression changes during tumorigenesis. Target gene expression is generally normalized by a stably expressed endogenous reference gene; however, reference gene expression may differ among tissues under various circumstances. Because no valid reference genes have been documented for human breast cancer cell lines containing different cancer subtypes treated with transient transfection, we identified appropriate and reliable reference genes from thirteen candidates in a panel of 10 normal and cancerous human breast cell lines under experimental conditions with/without transfection treatments with two transfection reagents. Reference gene expression stability was calculated using four algorithms (geNorm, NormFinder, BestKeeper and comparative delta Ct), and the recommended comprehensive ranking was provided using geometric means of the ranking values using the RefFinder tool. GeNorm analysis revealed that two reference genes should be sufficient for all cases in this study. A stability analysis suggests that 18S rRNA-ACTB is the best reference gene combination across all cell lines; ACTB-GAPDH is best for basal breast cancer cell lines; and HSPCB-ACTB is best for ER+ breast cancer cells. After transfection, the stability ranking of the reference gene fluctuated, especially with Lipofectamine 2000 transfection reagent in two subtypes of basal and ER+ breast cell lines. Comparisons of relative target gene (HER2) expression revealed different expressional patterns depending on the reference genes used for normalization. We suggest that identifying the most stable and suitable reference genes is critical for studying specific cell lines under certain circumstances.

BackgroundThe -308 G/A polymorphism in the tumor necrosis factor (TNF) gene has been implicated in the risk of acne vulgaris, but the results are inconclusive. The present meta-analysis aimed to investigate the overall association between the -308 G/A polymorphism and acne vulgaris risk.MethodsWe searched in Pubmed, Embase, Web of Science and CNKI for studies evaluating the association between the -308 G/A gene polymorphism and acne vulgaris risk. Data were extracted and statistical analysis was performed using STATA 12.0 software.ResultsA total of five publications involving 1553 subjects (728 acne vulgaris cases and 825 controls) were included in this meta-analysis. Combined analysis revealed a significant association between this polymorphism and acne vulgaris risk under recessive model (OR = 2.73, 95% CI: 1.37-5.44, p = 0.004 for AA vs. AG + GG). Subgroup analysis by ethnicity showed that the acne vulgaris risk associated with the -308 G/A gene polymorphism was significantly elevated among Caucasians under recessive model (OR = 2.34, 95% CI: 1.13-4.86, p = 0.023).ConclusionThis meta-analysis suggests that the -308 G/A polymorphism in the TNF gene contributes to acne vulgaris risk, especially in Caucasian populations. Further studies among different ethnicity populations are needed to validate these findings.

Recent analyses of ancient Mesopotamian mitochondrial genomes have suggested a genetic link between the Indian subcontinent and Mesopotamian civilization. There is no consensus on the origin of the ancient Mesopotamians. They may be descendants of migrants, who founded regional Mesopotamian groups like that of Terqa or they may be merchants who were involved in trans Mesopotamia trade. To identify the Indian source population showing linkage to the ancient Mesopotamians, we screened a total of 15,751 mitochondrial DNAs (11,432 from the literature and 4,319 from this study) representing all major populations of India. Our results although suggest that south India (Tamil Nadu) and northeast India served as the source of the ancient Mesopotamian mtDNA gene pool, mtDNA of these ancient Mesopotamians probably contributed by Tamil merchants who were involved in the Indo-Roman trade.

Anthropogenic global climate changes are one of the greatest threats to biodiversity. Distribution modeling can predict the effects of climate changes and potentially their effects on genetic diversity. DNA barcoding quickly identifies patterns of genetic diversity. As a case study, we use DNA barcodes and distribution models to predict threats under climate changes in the frog Nanorana parkeri, which is endemic to the Qinghai-Tibetan Plateau. Barcoding identifies major lineages W and E. Lineage W has a single origin in a refugium and Lineage E derives from three refugia. All refugia locate in river valleys and each greatly contributes to the current level of intraspecific genetic diversity. Species distribution models suggest that global climate changes will greatly influence N. parkeri, especially in the level of genetic diversity, because two former refugia will fail to provide suitable habitat. Our pipeline provides a novel application of DNA barcoding and has important implications for the conservation of biodiversity in southern areas of the Qinghai-Tibetan Plateau.

The phylogenetic relationships of species in the Phasianidae, Order Galliformes, are the object of intensive study. However, convergent morphological evolution and rapid species radiation result in much ambiguity in the group. Further, matrilineal (mtDNA) genealogies conflict with trees based on nuclear DNA retrotransposable elements. Herein, we analyze 39 nearly complete mitochondrial genomes (three new) and up to seven nuclear DNA segments. We combine these multiple unlinked, more informative genetic markers to infer historical relationships of the major groups of phasianids. The nuclear DNA tree is largely congruent with the tree derived from mt genomes. However, branching orders of mt/nuclear trees largely conflict with those based on retrotransposons. For example, Gallus/Bambusicola/Francolinus forms the sister-group of Coturnix/Alectoris in the nuclear/mtDNA trees, yet the tree based on retrotransposable elements roots the former at the base of the tree and not with the latter. Further, while peafowls cluster with Gallus/Coturnix in the mt tree, they root at the base of the phasianids following Gallus in the tree based on retrotransposable elements. The conflicting branch orders in nuclear/mtDNA and retrotransposons-based trees in our study reveal the complex topology of the Phasianidae.

The Belgian Malinois (BM) is an excellent working dog that typically shows a circling behavior when placed in a confined space. Moreover, individuals showing moderate running in circles (one kind of obsessive compulsive behavior) in confined spaces typically show better work performance compared to those without the circling behavior or to those with a serious circling behavior (which can be defined as an obsessive compulsive disorder (OCD)). To determine whether the candidate gene CDH2, Cadherin 2, which is associated with OCD in the Doberman pinscher breed of dogs and in humans, was linked with this behavioral character in the BM, population genetic analyses were performed on a BM population and a natural population of the Chinese indigenous dog (CID). Many genetic signals of balancing selection were detected for one specific region of the CDH2 gene, which suggests that a genomic block, which is included in the CDH2 gene, experienced balancing selection in the BM, and that the CDH2 gene might be associated with the behavioral characteristics of the BM dog (a balance between circling behavior and work performance). Moreover one specific variant, G63913941A, which creates a predicted transcription factor-binding site, may be the key mutation in the CDH2 gene affecting the behavior of BMs by allowing the binding of a transcription factor and increasing CDH2 expression.

BackgroundPliocene uplifting of the Qinghai-Tibetan Plateau (QTP) and Quaternary glaciation may have impacted the Asian biota more than any other events. Little is documented with respect to how the geological and climatological events influenced speciation as well as spatial and genetic structuring, especially in vertebrate endotherms. Macaca mulatta is the most widely distributed non-human primate. It may be the most suitable model to test hypotheses regarding the genetic consequences of orogenesis on an endotherm.Methodology and principal findingsUsing a large dataset of maternally inherited mitochondrial DNA gene sequences and nuclear microsatellite DNA data, we discovered two maternal super-haplogroups exist, one in western China and the other in eastern China. M. mulatta formed around 2.31 Ma (1.51-3.15, 95%), and divergence of the two major matrilines was estimated at 1.15 Ma (0.78-1.55, 95%). The western super-haplogroup exhibits significant geographic structure. In contrast, the eastern super-haplogroup has far greater haplotypic variability with little structure based on analyses of six variable microsatellite loci using Structure and Geneland. Analysis using Migrate detected greater gene flow from WEST to EAST than vice versa. We did not detect signals of bottlenecking in most populations.ConclusionsAnalyses of the nuclear and mitochondrial datasets obtained large differences in genetic patterns for M. mulatta. The difference likely reflects inheritance mechanisms of the maternally inherited mtDNA genome versus nuclear biparentally inherited STRs and male-mediated gene flow. Dramatic environmental changes may be responsible for shaping the matrilineal history of macaques. The timing of events, the formation of M. mulatta, and the divergence of the super-haplogroups, corresponds to both the uplifting of the QTP and Quaternary climatic oscillations. Orogenesis likely drove divergence of western populations in China, and Pleistocene glaciations are likely responsible for genetic structuring in the eastern super-haplogroup via geographic isolation and secondary contact.

Goldfish, Carassius auratus, have experienced strong anthropogenic selection during their evolutionary history, generating a tremendous extent of morphological variation relative to that in native Carassius. To locate the geographic origin of goldfish, we analyzed nucleotide sequences from part of the control region (CR) and the entire cytochrome b (Cytb) mitochondrial DNA genes for 234 goldfish and a large series of native specimens. Four important morphological characteristics used in goldfish taxonomy-body shape, dorsal fin, eye shape, and tailfin-were selected for hypothesis-testing to identify those that better correspond to evolutionary history.Haplotypes of goldfish rooted in two sublineages (C5 and C6), which contained the haplotypes of native C. a. auratus from southern China. Values of F(ST) and N(m) revealed a close relationship between goldfish and native C. a. auratus from the lower Yangtze River. An extraordinary, stepwise loss of genetic diversity was detected from native fish to goldfish and from Grass-goldfish relative to other breeds. Significantly negative results for the tests of Tajima's D and Fu and Li's D* and F* were identified in goldfish, including the Grass breed. The results identified eye-shape as being the least informative character for grouping goldfish with respect to their evolutionary history. Fisher's exact test identified matrilineal constraints on domestication.Chinese goldfish have a matrilineal origin from native southern Chinese C. a. auratus, especially the lineages from the lower Yangtze River. Anthropogenic selection of the native Carassius eliminated aesthetically unappealing goldfish and this action appeared to be responsible for the stepwise decrease in genetic diversity of domesticated goldfish, a process similar to that reported for the domestication of pigs, rice, and maize. The three-breed taxonomy--Grass-goldfish, Egg-goldfish, and Wen-goldfish--better reflected the history of domestication.

Bats and cetaceans (i.e., whales, dolphins, porpoises) are two kinds of mammals with unique locomotive styles and occupy novel niches. Bats are the only mammals capable of sustained flight in the sky, while cetaceans have returned to the aquatic environment and are specialized for swimming. Associated with these novel adaptations to their environment, various development changes have occurred to their body plans and associated structures. Given the importance of Hox genes in many aspects of embryonic development, we conducted an analysis of the coding regions of all Hox gene family members from bats (represented by Pteropus vampyrus, Pteropus alecto, Myotis lucifugus and Myotis davidii) and cetaceans (represented by Tursiops truncatus) for adaptive evolution using the available draft genome sequences. Differences in the selective pressures acting on many Hox genes in bats and cetaceans were found compared to other mammals. Positive selection, however, was not found to act on any of the Hox genes in the common ancestor of bats and only upon Hoxb9 in cetaceans. PCR amplification data from additional bat and cetacean species, and application of the branch-site test 2, showed that the Hoxb2 gene within bats had significant evidence of positive selection. Thus, our study, with genomic and newly sequenced Hox genes, identifies two candidate Hox genes that may be closely linked with developmental changes in bats and cetaceans, such as those associated with the pancreatic, neuronal, thymus shape and forelimb. In addition, the difference in our results from the genome-wide scan and newly sequenced data reveals that great care must be taken in interpreting results from draft genome data from a limited number of species, and deep genetic sampling of a particular clade is a powerful tool for generating complementary data to address this limitation.

Coat color in dog breeds is an excellent character for revealing the power of artificial selection, as it is extremely diverse and likely the result of recent domestication. Coat color is generated by melanocytes, which synthesize pheomelanin (a red or yellow pigment) or eumelanin (a black or brown pigment) through the pigment type-switching pathway, and is regulated by three genes in dogs: MC1R (melanocortin receptor 1), CBD103 (β-defensin 103), and ASIP (agouti-signaling protein precursor). The genotypes of these three gene loci in dog breeds are associated with coat color pattern. Here, we resequenced these three gene loci in two Kunming dog populations and analyzed these sequences using population genetic approaches to identify evolutionary patterns that have occurred at these loci during the recent domestication and breeding of the Kunming dog. The analysis showed that MC1R undergoes balancing selection in both Kunming dog populations, and that the Fst value for MC1R indicates significant genetic differentiation across the two populations. In contrast, similar results were not observed for CBD103 or ASIP. These results suggest that high heterozygosity and allelic differences at the MC1R locus may explain both the mixed color coat, of yellow and black, and the difference in coat colors in both Kunming dog populations.

Odd-nosed monkeys represent one of the two major groups of Asian colobines. Our knowledge about this primate group is still limited as it is highlighted by the recent discovery of a new species in Northern Myanmar. Although a common origin of the group is now widely accepted, the phylogenetic relationships among its genera and species, and the biogeographic processes leading to their current distribution are largely unknown. To address these issues, we have analyzed complete mitochondrial genomes and 12 nuclear loci, including one X chromosomal, six Y chromosomal and five autosomal loci, from all ten odd-nosed monkey species. The gene tree topologies and divergence age estimates derived from different markers were highly similar, but differed in placing various species or haplogroups within the genera Rhinopithecus and Pygathrix. Based on our data, Rhinopithecus represent the most basal lineage, and Nasalis and Simias form closely related sister taxa, suggesting a Northern origin of odd-nosed monkeys and a later invasion into Indochina and Sundaland. According to our divergence age estimates, the lineages leading to the genera Rhinopithecus, Pygathrix and Nasalis+Simias originated in the late Miocene, while differentiation events within these genera and also the split between Nasalis and Simias occurred in the Pleistocene. Observed gene tree discordances between mitochondrial and nuclear datasets, and paraphylies in the mitochondrial dataset for some species of the genera Rhinopithecus and Pygathrix suggest secondary gene flow after the taxa initially diverged. Most likely such events were triggered by dramatic changes in geology and climate within the region. Overall, our study provides the most comprehensive view on odd-nosed monkey evolution and emphasizes that data from differentially inherited markers are crucial to better understand evolutionary relationships and to trace secondary gene flow.

The Cham people are the major Austronesian speakers of Mainland Southeast Asia (MSEA) and the reconstruction of the Cham population history can provide insights into their diffusion. In this study, we analyzed non-recombining region of the Y chromosome markers of 177 unrelated males from four populations in MSEA, including 59 Cham, 76 Kinh, 25 Lao, and 17 Thai individuals. Incorporating published data from mitochondrial DNA (mtDNA), our results indicated that, in general, the Chams are an indigenous Southeast Asian population. The origin of the Cham people involves the genetic admixture of the Austronesian immigrants from Island Southeast Asia (ISEA) with the local populations in MSEA. Discordance between the overall patterns of Y chromosome and mtDNA in the Chams is evidenced by the presence of some Y chromosome lineages that prevail in South Asians. Our results suggest that male-mediated dispersals via the spread of religions and business trade might play an important role in shaping the patrilineal gene pool of the Cham people.

BACKGROUND: While prions play a central role in the pathogenesis of transmissible spongiform encephalopathies, the biology of these proteins and the pathophysiology of these diseases remain largely unknown. Since no case of bovine spongiform encephalopathy (BSE) has ever been reported in buffalo despite their phylogenetic proximity to cattle, genetic differences may be driving the different susceptibilities of these two species to BSE. We thus hypothesized that differences in expression of the most recently identified member of the prion family or Shadoo (SPRN) gene may relate to these species-specific differences. PRINCIPAL FINDINGS: We first analyzed and compared the polymorphisms of the SPRN gene (~4.4 kb), including the putative promoter, coding and 3' regions, and further verified the entire ORF and putative promoter. This yielded a total of 117 fixed differences, remarkably: 1) a 12-bp insertion/deletion polymorphism in the hydrophobic domain of the cattle but not buffalo gene, introducing a four amino acid expansion/contraction in a series of 5 tandem Ala/Gly-containing repeats; 2) two fixed missense mutations (102Ser→Gly and 119Thr→Ala), and three missense mutations (92Pro>Thr/Met, 122Thr>Ile and 139Arg>Trp) in the coding region presenting different (P

For most mammals, running is their major locomotive style, however, cetaceans and bats are two mammalian groups that have independently developed new locomotive styles (swimming and flying) from their terrestrial ancestors. In this study, we used a genome-wide comparative analysis in an attempt to identify the selective imprint of the development of new locomotive styles by cetaceans and bats to adapt to their new ecological niches. We found that an elevated proportion of mitochondrion-associated genes show evidence of adaptive evolution in cetaceans and on the common ancestral lineage leading to bats, compared to other terrestrial mammals. This result is consistent with the fact that during the independent developments of swimming and flying in these two groups, the changes of energy metabolism ratios would be among the most important factors to overcome elevated energy demands. Furthermore, genes that show evidence of sequence convergence or parallel evolution in these two lineages were overrepresented in the categories of energy metabolism, muscle contraction, heart, and glucose metabolism, genes that perform functions which are essential for locomotion. In conclusion, our analyses showed that on the dolphin and bat lineages, genes associated with locomotion not only both show a greater propensity to adaptively evolve, but also show evidence of sequence convergence, which likely reflects a response to a common requirement during their development of these two drastic locomotive styles.

The molecular basis of the evolution of phenotypic characters is very complex and is poorly understood with few examples documenting the roles of multiple genes. Considering that a single gene cannot fully explain the convergence of phenotypic characters, we choose to study the convergent evolution of rod vision in two divergent bats from a network perspective. The Old World fruit bats (Pteropodidae) are non-echolocating and have binocular vision, whereas the sheath-tailed bats (Emballonuridae) are echolocating and have monocular vision; however, they both have relatively large eyes and rely more on rod vision to find food and navigate in the night. We found that the genes CRX, which plays an essential role in the differentiation of photoreceptor cells, SAG, which is involved in the desensitization of the photoactivated transduction cascade, and the photoreceptor gene RH, which is directly responsible for the perception of dim light, have undergone parallel sequence evolution in two divergent lineages of bats with larger eyes (Pteropodidae and Emballonuroidea). The multiple convergent events in the network of genes essential for rod vision is a rare phenomenon that illustrates the importance of investigating pathways and networks in the evolution of the molecular basis of phenotypic convergence.

GDF5 is a member of the bone morphogenetic protein (BMP) gene family, and plays an important role in the development of the skeletal system. Variants of the gene are associated with osteoarthritis and height in some human populations. Here, we resequenced the gene in individuals from four geographically separated human populations, and found that the evolution of the promoter region deviated from neutral expectations, with the sequence evolution driven by positive selection in the East Asian population, especially the haplotypes carrying the derived alleles of 5' UTR SNPs rs143384 and rs143383. The derived alleles of rs143384 and rs143383, which are associated with a risk of osteoarthritis and decreased height, have high frequencies in non-Africans and show strong extended haplotype homozygosity and high population differentiation in East Asian. It is concluded that positive selection has driven the rapid evolution of the two osteoarthritis osteoarthritis-risk and decreased height associated variants of the human GDF5 gene, and supports the suggestion that the reduction in body size during the terminal Pleistocene and Holocene period might have been an adaptive process influenced by genetic factors.

CLEC4M is a C-type lectin gene serving as cell adhesion receptor and pathogen recognition receptor. It recognizes several pathogens of important public health concern. In particular, a highly polymorphic variable number tandem repeat (VNTR) at the neck-region of CLEC4M had been associated with genetic predisposition to some infectious diseases. To gain insight into the origin and evolution of this VNTR in CLEC4M, we studied 21 Africans, 20 Middle Easterns, 35 Europeans, 38 Asians, 13 Oceania, and 18 Americans (a total of 290 chromosomes) from the (Human Genome Diversity Panel) HGDP-CEPH panel; these samples covered most of alleles of this VNTR locus present in human populations. We identified a limited number of haplotypes among the basic repeat subunits that is 69 base pairs in length. Only 8 haplotypes were found. Their sequence identities were determined in the 290 chromosomes. VNTR alleles of different repeat length (from 4 to 9 repeats) were analyzed for composition and orientation of these subunits. Our results showed that the subunit configuration of the same repeat number of VNTR locus from different populations were, in fact, virtually identical. It implies that most of the VNTR alleles existed before dispersion of modern humans outside Africa. Further analyses indicate that the present diversity profile of this locus in worldwide populations is generated from the effect of migration of different tribes and neutral evolution. Our findings do not support the hypothesis that the origin of the VNTR alleles were arisen by independent (separate) mutation events and caused by differential allele advantage and natural selection as suggested by previous report based on SNP data.

BACKGROUND: Erinaceidae is a family of small mammals that include the spiny hedgehogs (Erinaceinae) and the silky-furred moonrats and gymnures (Galericinae). These animals are widely distributed across Eurasia and Africa, from the tundra to the tropics and the deserts to damp forests. The importance of these animals lies in the fact that they are the oldest known living placental mammals, which are well represented in the fossil record, a rarity fact given their size and vulnerability to destruction during fossilization. Although the Family has been well studied, their phylogenetic relationships remain controversial. To test previous phylogenetic hypotheses, we combined molecular and morphological data sets, including representatives of all the genera. METHODOLOGY AND PRINCIPAL FINDINGS: We included in the analyses 3,218 bp mitochondrial genes, one hundred and thirty-five morphological characters, twenty-two extant erinaceid taxa, and five outgroup taxa. Phylogenetic relationships were reconstructed using both partitioned and combined data sets. As in previous analyses, our results strongly support the monophyly of both subfamilies (Galericinae and Erinaceinae), the Hylomys group (to include Neotetracus and Neohylomys), and a sister-relationship of Atelerix and Erinaceus. As well, we verified that the extremely long branch lengths within the Galericinae are consistent with their fossil records. Not surprisingly, we found significant incongruence between the phylogenetic signals of the genes and the morphological characters, specifically in the case of Hylomys parvus, Mesechinus, and relationships between Hemiechinus and Paraechinus. CONCLUSIONS: Although we discovered new clues to understanding the evolutionary relationships within the Erinaceidae, our results nonetheless, strongly suggest that more robust analyses employing more complete taxon sampling (to include fossils) and multiple unlinked genes would greatly enhance our understanding of the Erinaceidae. Until then, we have left the nomenclature of the taxa unchanged; hence it does not yet precisely reflect their phylogenetic relationships or the depth of their genetic diversity.

Extensive studies in different fields have been performed to reconstruct the prehistory of populations in the Japanese archipelago. Estimates the ancestral population dynamics based on Japanese molecular sequences can extend our understanding about the colonization of Japan and the ethnogenesis of modern Japanese.We applied Bayesian skyline plot (BSP) with a dataset based on 952 Japanese mitochondrial DNA (mtDNA) genomes to depict the female effective population size (N(ef)) through time for the total Japanese and each of the major mtDNA haplogroups in Japanese. Our results revealed a rapid N(ef) growth since ∼5 thousand years ago had left ∼72% Japanese mtDNA lineages with a salient signature. The BSP for the major mtDNA haplogroups indicated some different demographic history.The results suggested that the rapid population expansion acted as a major force in shaping current maternal pool of Japanese. It supported a model for population dynamics in Japan in which the prehistoric population growth initiated in the Middle Jomon Period experienced a smooth and swift transition from Jomon to Yayoi, and then continued through the Yayoi Period. The confounding demographic backgrounds of different mtDNA haplogroups could also have some implications for some related studies in future.

BACKGROUND: A major reduction of genetic diversity in mtDNA occurred during the domestication of East Asian pigs. However, the extent to which genetic diversity has been lost in the nuclear genome is uncertain. To reveal levels and patterns of nucleotide diversity and to elucidate the genetic relationships and demographic history of domestic pigs and their ancestors, wild boars, we investigated 14 nuclear markers (including 8 functional genes, 2 pseudogenes and 4 intergenic regions) from 11 different chromosomes in East Asia-wide samples and pooled them with previously obtained mtDNA data for a combined analysis. PRINCIPAL FINDINGS: The results indicated that domestic pigs and wild boars possess comparable levels of nucleotide diversity across the nuclear genome, which is inconsistent with patterns that have been found in mitochondrial genome. CONCLUSIONS: This incongruence between the mtDNA and nuclear genomes is suggestive of a large-scale backcross between male wild boars and female domestic pigs in East Asia. Our data reveal the impacts of founder effects and backcross on the pig genome and help us better understand the complex demographic histories of East Asian pigs, which will be useful for future work on artificial selection.

In the past decade, a high incidence of somatic mitochondrial DNA (mtDNA) mutations has been observed, mostly based on a fraction of the molecule, in various cancerous tissues; nevertheless, some of them were queried due to problems in data quality. Obviously, without a comprehensive understanding of mtDNA mutational profile in the cancerous tissue of a specific patient, it is unlikely to disclose the genuine relationship between somatic mtDNA mutations and tumorigenesis. To achieve this objective, the most straightforward way is to directly compare the whole mtDNA genome variation among three tissues (namely, cancerous tissue, para-cancerous tissue, and distant normal tissue) from the same patient. Considering the fact that most of the previous studies on the role of mtDNA in colorectal tumor focused merely on the D-loop or partial segment of the molecule, in the current study we have collected three tissues (cancerous, para-cancerous and normal tissues) respectively recruited from 20 patients with colorectal tumor and completely sequenced the mitochondrial genome of each tissue. Our results reveal a relatively lower incidence of somatic mutations in these patients; intriguingly, all somatic mutations are in heteroplasmic status. Surprisingly, the observed somatic mutations are not restricted to cancer tissues, for the para-cancer tissues and distant normal tissues also harbor somatic mtDNA mutations with a lower frequency than cancerous tissues but higher than that observed in the general population. Our results suggest that somatic mtDNA mutations in cancerous tissues could not be simply explained as a consequence of tumorigenesis; meanwhile, the somatic mtDNA mutations in normal tissues might reflect an altered physiological environment in cancer patients.

The leptin gene has received intensive attention and scientific investigation for its importance in energy homeostasis and reproductive regulation in mammals. Furthermore, study of the leptin gene is of crucial importance for public health, particularly for its role in obesity, as well as for other numerous physiological roles that it plays in mammals. In the present work, we report the identification of novel leptin genes in 4 species of Cetacea, and a comparison with 55 publicly available leptin sequences from mammalian genome assemblies and previous studies. Our study provides evidence for positive selection in the suborder Odontoceti (toothed whales) of the Cetacea and the family Phocidae (earless seals) of the Pinnipedia. We also detected positive selection in several leptin gene residues in these two lineages. To test whether leptin and its receptor evolved in a coordinated manner, we analyzed 24 leptin receptor gene (LPR) sequences from available mammalian genome assemblies and other published data. Unlike the case of leptin, our analyses did not find evidence of positive selection for LPR across the Cetacea and Pinnipedia lineages. In line with this, positively selected sites identified in the leptin genes of these two lineages were located outside of leptin receptor binding sites, which at least partially explains why co-evolution of leptin and its receptor was not observed in the present study. Our study provides interesting insights into current understanding of the evolution of mammalian leptin genes in response to selective pressures from life in an aquatic environment, and leads to a hypothesis that new tissue specificity or novel physiologic functions of leptin genes may have arisen in both odontocetes and phocids. Additional data from other species encompassing varying life histories and functional tests of the adaptive role of the amino acid changes identified in this study will help determine the factors that promote the adaptive evolution of the leptin genes in marine mammals.

Rhodopsin, encoded by the gene Rhodopsin (RH1), is extremely sensitive to light, and is responsible for dim-light vision. Bats are nocturnal mammals that inhabit poor light environments. Megabats (Old-World fruit bats) generally have well-developed eyes, while microbats (insectivorous bats) have developed echolocation and in general their eyes were degraded, however, dramatic differences in the eyes, and their reliance on vision, exist in this group. In this study, we examined the rod opsin gene (RH1), and compared its evolution to that of two cone opsin genes (SWS1 and M/LWS). While phylogenetic reconstruction with the cone opsin genes SWS1 and M/LWS generated a species tree in accord with expectations, the RH1 gene tree united Pteropodidae (Old-World fruit bats) and Yangochiroptera, with very high bootstrap values, suggesting the possibility of convergent evolution. The hypothesis of convergent evolution was further supported when nonsynonymous sites or amino acid sequences were used to construct phylogenies. Reconstructed RH1 sequences at internal nodes of the bat species phylogeny showed that: (1) Old-World fruit bats share an amino acid change (S270G) with the tomb bat; (2) Miniopterus share two amino acid changes (V104I, M183L) with Rhinolophoidea; (3) the amino acid replacement I123V occurred independently on four branches, and the replacements L99M, L266V and I286V occurred each on two branches. The multiple parallel amino acid replacements that occurred in the evolution of bat RH1 suggest the possibility of multiple convergences of their ecological specialization (i.e., various photic environments) during adaptation for the nocturnal lifestyle, and suggest that further attention is needed on the study of the ecology and behavior of bats.

BACKGROUND: Human skeletal system has evolved rapidly since the dispersal of modern humans from Africa, potentially driven by selection and adaptation. Osteogenin (BMP3) plays an important role in skeletal development and bone osteogenesis as an antagonist of the osteogenic bone morphogenetic proteins, and negatively regulates bone mineral density. METHODOLOGY/PRINCIPAL FINDINGS: Here, we resequenced the BMP3 gene from individuals in four geographically separated modern human populations. Features supportive of positive selection in the BMP3 gene were found including the presence of an excess of nonsynonymous mutations in modern humans, and a significantly lower genetic diversity that deviates from neutrality. The prevalent haplotypes of the first exon region in Europeans demonstrated features of long-range haplotype homogeneity. In contrast with findings in European, the derived allele SNP Arg192Gln shows higher extended haplotype homozygosity in East Asian. The worldwide allele frequency distribution of SNP shows not only a high-derived allele frequency in Asians, but also in Americans, which is suggestive of functional adaptation. CONCLUSIONS/SIGNIFICANCE: In conclusion, we provide evidence for recent positive selection operating upon a crucial gene in skeletal development, which may provide new insight into the evolution of the skeletal system and bone development.

Fig trees are pollinated by fig wasps, which also oviposit in female flowers. The wasp larvae gall and eat developing seeds. Although fig trees benefit from allowing wasps to oviposit, because the wasp offspring disperse pollen, figs must prevent wasps from ovipositing in all flowers, or seed production would cease, and the mutualism would go extinct. In Ficus racemosa, we find that syconia ('figs') that have few foundresses (ovipositing wasps) are underexploited in the summer (few seeds, few galls, many empty ovules) and are overexploited in the winter (few seeds, many galls, few empty ovules). Conversely, syconia with many foundresses produce intermediate numbers of galls and seeds, regardless of season. We use experiments to explain these patterns, and thus, to explain how this mutualism is maintained. In the hot summer, wasps suffer short lifespans and therefore fail to oviposit in many flowers. In contrast, cooler temperatures in the winter permit longer wasp lifespans, which in turn allows most flowers to be exploited by the wasps. However, even in winter, only in syconia that happen to have few foundresses are most flowers turned into galls. In syconia with higher numbers of foundresses, interference competition reduces foundress lifespans, which reduces the proportion of flowers that are galled. We further show that syconia encourage the entry of multiple foundresses by delaying ostiole closure. Taken together, these factors allow fig trees to reduce galling in the wasp-benign winter and boost galling (and pollination) in the wasp-stressing summer. Interference competition has been shown to reduce virulence in pathogenic bacteria. Our results show that interference also maintains cooperation in a classic, cooperative symbiosis, thus linking theories of virulence and mutualism. More generally, our results reveal how frequency-dependent population regulation can occur in the fig-wasp mutualism, and how a host species can 'set the rules of the game' to ensure mutualistic behavior in its symbionts.

The intron 5 of gene LMBR1 is the cis-acting regulatory module for the sonic hedgehog (SHH) gene. Mutation in this non-coding region is associated with preaxial polydactyly, and may play crucial roles in the evolution of limb and skeletal system.We sequenced a region of the LMBR1 gene intron 5 in East Asian human population, and found a significant deviation of Tajima's D statistics from neutrality taking human population growth into account. Data from HapMap also demonstrated extended linkage disequilibrium in the region in East Asian and European population, and significantly low degree of genetic differentiation among human populations.We proposed that the intron 5 of LMBR1 was presumably subject to balancing selection during the evolution of modern human.