Genetic signatures of strong recent positive selection at the lactase gene

In most human populations, the ability to digest lactose contained in milk usually disappears in childhood, but in European-derived populations, lactase activity frequently persists into adulthood (Scrimshaw and Murray 1988). It has been suggested (Cavalli-Sforza 1973; Hollox et al. 2001; Enattah et al. 2002; Poulter et al. 2003) that a selective advantage based on additional nutrition from dairy explains these genetically determined population differences (Simoons 1970; Kretchmer 1971; Scrimshaw and Murray 1988; Enattah et al. 2002), but formal population-genetics-based evidence of selection has not yet been provided. To assess the population-genetics evidence for selection, we typed 101 single-nucleotide polymorphisms covering 3.2 Mb around the lactase gene. In northern European-derived populations, two alleles that are tightly associated with lactase persistence (Enattah et al. 2002) uniquely mark a common (~77%) haplotype that extends largely undisrupted for >1 Mb. We provide two new lines of genetic evidence that this long, common haplotype arose rapidly due to recent selection: (1) by use of the traditional F(ST) measure and a novel test based on p(excess), we demonstrate large frequency differences among populations for the persistence-associated markers and for flanking markers throughout the haplotype, and (2) we show that the haplotype is unusually long, given its high frequency--a hallmark of recent selection. We estimate that strong selection occurred within the past 5,000-10,000 years, consistent with an advantage to lactase persistence in the setting of dairy farming; the signals of selection we observe are among the strongest yet seen for any gene in the genome.     

Mammalian sperm proteins are rapidly evolving: evidence of positive selection in functionally diverse genes

A growing number of genes involved in sex and reproduction have been demonstrated to be rapidly evolving. Here, we show that genes expressed solely in spermatozoa represent a highly diverged subset among mouse and human tissue-specific orthologs. The average rate of nonsynonymous substitutions per site (K(a)) is significantly higher in sperm proteins (mean K(a) = 0.18; N = 35) than in proteins expressed specifically in all other tissues (mean K(a) = 0.074; N = 473). No differences, however, are found in the synonymous substitution rate (K(s)) between tissues, suggesting that selective forces, and not mutation rate, explain the high rate of replacement substitutions in sperm proteins. Four out of 19 sperm-specific genes with characterized function demonstrated evidence of strong positive Darwinian selection, including a protein involved in gene regulation, Protamine-1 (PRM1), a protein involved in glycolysis, GAPDS, and two egg-binding proteins, Adam-2 precursor (ADAM2) and sperm-adhesion molecule-1 (SAM1). These results demonstrate the rapid evolution of sperm-specific genes and highlight the molecular action of sexual selection on a variety of characters involved in mammalian sperm function.     

Apoptotic HPV positive cancer cells exhibit transforming properties

Previous studies have shown that DNA can be transferred from dying engineered cells to neighboring cells through the phagocytosis of apoptotic bodies, which leads to cellular transformation. Here, we provide evidence of an uptake of apoptotic-derived cervical cancer cells by human mesenchymal cells. Interestingly, HeLa (HPV 18+) or Ca Ski (HPV16+) cells, harboring integrated high-risk HPV DNA but not C-33 A cells (HPV-), were able to transform the recipient cells. Human primary fibroblasts engulfed the apoptotic bodies effectively within 30 minutes after co-cultivation. This mechanism is active and involves the actin cytoskeleton. In situ hybridization of transformed fibroblasts revealed the presence of HPV DNA in the nucleus of a subset of phagocytosing cells. These cells expressed the HPV16/18 E6 gene, which contributes to the disruption of the p53/p21 pathway, and the cells exhibited a tumorigenic phenotype, including an increased proliferation rate, polyploidy and anchorage independence growth. Such horizontal transfer of viral oncogenes to surrounding cells that lack receptors for HPV could facilitate the persistence of the virus, the main risk factor for cervical cancer development. This process might contribute to HPV-associated disease progression in vivo.     

Identification of positive selection signatures in pigs by comparing linkage disequilibrium variances

Selection affects the patterns of linkage disequilibrium (LD) around the site of a beneficial allele with an increase in LD among the hitchhiking alleles. Comparing the differences in regional LD between pig populations could help to identify putative genomic regions with potential adaptations for economic traits. In this study, using Illumina Porcine SNP60K BeadChip genotyping data from 207 Chinese indigenous, 117 South American village and 408 Large White pigs, we estimated the variation of genome‐wide LD between populations using the varld program. The top 0.1% standardized VarLD scores were used as a criterion for all comparisons, and compared with LD blocks, a total of four selection signatures on Sus scrofa chromosome (SSC) 7, 9, 13 and 14 were identified in all populations. These signatures overlapped with quantitative trait loci for linoleic acid content, age at puberty, number of muscle fibers per unit area, hip structure and body weight traits in pigs. Among them, one of the signatures (56.5–56.6 Mb on SSC7) in Large White pigs harbored the ADAMTSL3 gene, which is known to affect body length. The findings of this study seem to point toward recent selection in different pig populations. Further investigations are encouraged to confirm the selection signatures detected by varld in the present study.     

Widespread positive selection in synonymous sites of mammalian genes

Evolution of protein sequences is largely governed by purifying selection, with a small fraction of proteins evolving under positive selection. The evolution at synonymous positions in protein-coding genes is not nearly as well understood, with the extent and types of selection remaining, largely, unclear. A statistical test to identify purifying and positive selection at synonymous sites in protein-coding genes was developed. The method compares the rate of evolution at synonymous sites (Ks) to that in intron sequences of the same gene after sampling the aligned intron sequences to mimic the statistical properties of coding sequences. We detected purifying selection at synonymous sites in approximately 28% of the 1,562 analyzed orthologous genes from mouse and rat, and positive selection in approximately 12% of the genes. Thus, the fraction of genes with readily detectable positive selection at synonymous sites is much greater than the fraction of genes with comparable positive selection at nonsynonymous sites, i.e., at the level of the protein sequence. Unlike other genes, the genes with positive selection at synonymous sites showed no correlation between Ks and the rate of evolution in nonsynonymous sites (Ka), indicating that evolution of synonymous sites under positive selection is decoupled from protein evolution. The genes with purifying selection at synonymous sites showed significant anticorrelation between Ks and expression level and breadth, indicating that highly expressed genes evolve slowly. The genes with positive selection at synonymous sites showed the opposite trend, i.e., highly expressed genes had, on average, higher Ks. For the genes with positive selection at synonymous sites, a significantly lower mRNA stability is predicted compared to the genes with negative selection. Thus, mRNA destabilization could be an important factor driving positive selection in nonsynonymous sites, probably, through regulation of expression at the level of mRNA degradation and, possibly, also translation rate. So, unexpectedly, we found that positive selection at synonymous sites of mammalian genes is substantially more common than positive selection at the level of protein sequences. Positive selection at synonymous sites might act through mRNA destabilization affecting mRNA levels and translation.     

Search for signatures in miRNAs associated with cancer

Since the first discovery in the early 1990''s, the predicted and validated population of microRNAs (miRNAs or miRs) has grownsignificantly. These small (~22 nucleotides long) regulators of gene expression have been implicated and associated with severalgenes in the cancer pathway as well. Globally, the identification and verification of microRNAs as biomarkers for cancer cell typeshas been the area of thrust for most miRNA biologists. However, there has been a noticeable vacuum when it comes to identifyinga common signature or trademark that could be used to demarcate a miR to be associated with the development or suppression ofcancer. To answer these queries, we report an in silico study involving the identification of global signatures in experimentallyvalidated microRNAs which have been associated with cancer. This study has thrown light on the presence of significant commonsignatures, viz., - sequential and hybridization, which may distinguish a miR to be associated with cancer. Based on our analysis,we suggest the utility of such signatures in the design and development of algorithms for prediction of miRs involved in the cancerpathway.     

Genomewide signatures of selection in Epichloë reveal candidate genes for host specialization

Host specialization is a key process in ecological divergence and speciation of plant‐associated fungi. The underlying determinants of host specialization are generally poorly understood, especially in endophytes, which constitute one of the most abundant components of the plant microbiome. We addressed the genetic basis of host specialization in two sympatric subspecies of grass‐endophytic fungi from the Epichloë typhina complex: subsp. typhina and clarkii. The life cycle of these fungi entails unrestricted dispersal of gametes and sexual reproduction before infection of a new host, implying that the host imposes a selective barrier on viability of the progeny. We aimed to detect genes under divergent selection between subspecies, experiencing restricted gene flow due to adaptation to different hosts. Using pooled whole‐genome sequencing data, we combined F_ST and D_XY population statistics in genome scans and detected 57 outlier genes showing strong differentiation between the two subspecies. Genomewide analyses of nucleotide diversity (π), Tajima's D and dN/dS ratios indicated that these genes have evolved under positive selection. Genes encoding secreted proteins were enriched among the genes showing evidence of positive selection, suggesting that molecular plant–fungus interactions are strong drivers of endophyte divergence. We focused on five genes encoding secreted proteins, which were further sequenced in 28 additional isolates collected across Europe to assess genetic variation in a larger sample size. Signature of positive selection in these isolates and putative identification of pathogenic function supports our findings that these genes represent strong candidates for host specialization determinants in Epichloë endophytes. Our results highlight the role of secreted proteins as key determinants of host specialization.     

Multiple signatures of positive selection downstream of notch on the X chromosome in Drosophila melanogaster

To identify genomic regions affected by the rapid fixation of beneficial mutations (selective sweeps), we performed a scan of microsatellite variability across the Notch locus region of Drosophila melanogaster. Nine microsatellites spanning 60 kb of the X chromosome were surveyed for variation in one African and three non-African populations of this species. The microsatellites identified an approximately 14-kb window for which we observed relatively low levels of variability and/or a skew in the frequency spectrum toward rare alleles, patterns predicted at regions linked to a selective sweep. DNA sequence polymorphism data were subsequently collected within this 14-kb region for three of the D. melanogaster populations. The sequence data strongly support the initial microsatellite findings; in the non-African populations there is evidence of a recent selective sweep downstream of the Notch locus near or within the open reading frames CG18508 and Fcp3C. In addition, we observe a significant McDonald-Kreitman test result suggesting too many amino acid fixations species wide, presumably due to positive selection, at the unannotated open reading frame CG18508. Thus, we observe within this small genomic region evidence for both recent (skew toward rare alleles in non-African populations) and recurring (amino acid evolution at CG18508) episodes of positive selection.     

A scan for signatures of positive selection in candidate loci for skin pigmentation in humans

Although the combination of pale skin and intense sun exposure results in an important health risk for the individual, it is less clear if at the population level this risk has possessed an evolutionary meaning. In this sense, a number of adaptive hypotheses have been put forward to explain the evolution of human skin pigmentation, such as photoprotection against sun-induced cancer, sexual selection, vitamin D synthesis or photoprotection of photolabile compounds, among others. It is expected that if skin pigmentation is adaptive, we might be able to see the signature of positive selection on some of the genes involved. In order to detect this signature, we analyze a battery of 81 candidate loci by means of phylogenetic and population genetic tests. Our results indicate that both light and dark skin may possess adaptive value. Of the main loci presenting this signature, TP53BP1 shows clear evidence of adaptive selection in Africans, whereas TYRP1 and SLC24A5 show evidence of adaptive selection in Caucasians. Although we cannot offer a mechanism that based on these genes explains the advantage of light skin, if TP53BP1, and perhaps RAD50, have truly conferred an adaptive value to the African population analyzed, photoprotection against sun-induced skin damage/cancer might be proposed as a mechanism that has driven the evolution of human skin pigmentation.     

Large-scale search for genes on which positive selection may operate

We conducted a systematic search for the candidate genes on which positive selection may operate, on the premise that for such genes the number of nonsynonymous substitution is expected to be larger than that of synonymous substitutions when the nucleotide sequences of genes under investigation are compared with each other. By obtaining 3,595 groups of homologous sequences from the DDBJ, EMBL, and GenBank DNA sequence databases, we found that 17 gene groups can be the candidates for the genes on which positive selection may operate. Thus, such genes are found to occupy only about 0.5% of the vast number of gene groups so far available. Interestingly enough 9 out of the 17 gene groups were the surface antigens of parasites or viruses.      

Identifying candidate positive selection genes in Korean imported pig breeds

Domestication and artificial selection have modified the genome landscape of the pig. The identification of selection signatures in the genome can help to elucidate the selection mechanisms and uncover the causal genes related to the phenotypic variations between domestic pig breeds. We scanned the genomes of Korean imported pig breeds against native breeds using Z-transformed Fst (ZF_ST) and Heterozygosity (ZH_p) statistics to search for the signatures of positive selection. We identified 411 (ZF_ST?=?175; ZH_p?=?256) putatively selected genes in commercial breeds. The gene regions identified were harboring those related to immunity, coat color, reproduction function and other traits. Several genes (e.g., PLSCR4, AGTR1 and CORIN) were related to reproduction traits such as fertility, ovulation rate, and uterine function. This study revealed genes which improve our understanding of the biological mechanisms of higher litter sizes, the phenotype of interest, in higher litter pig breeds.     

Persistence of costly novel genes in the absence of positive selection

Many genetic changes that ultimately lead to adaptive evolution come with a short‐term cost expressed in terms of reduced survival and reproduction. In the absence of genetic drift, it is unclear how such costly mutations may persist. Here we experimentally demonstrate that parasites can promote the persistence of costly genetic variants. We employed a genetically engineered strain (GMMO) of the bacterium Pseudomonas fluorescens as a model of the acquisition of a new gene either through a major mutation or through horizontal transfer, and examined its persistence in different evolving communities comprising an ancestral strain and a lytic bacteriophage. Whereas competition resulted in the elimination of the GMMO, inclusion of the phage promoted GMMO persistence. We provide evidence for why this effect is due to the differential susceptibility of GMMO and ancestral bacteria to phage.     

Maximum‐likelihood approaches reveal signatures of positive selection in BMP15 and GDF9 genes modulating ovarian function in mammalian female fertility

Bone morphogenetic proteins (BMPs) and the growth factors (GDFs) play an important role in ovarian folliculogenesis and essential regulator of processes of numerous granulosa cells. BMP15 gene variations linked to various ovarian phenotypic consequences subject to the species, from infertility to improved prolificacy in sheep, primary ovarian insufficiency in women or associated with minor subfertility in mouse. To study the evolving role of BMP15 and GDF9, a phylogenetic analysis was performed. To find out the candidate gene associated with prolificacy in mammals, the nucleotide sequence of BMP15 and GDF9 genes was recognized under positive selection in various mammalian species. Maximum‐likelihood approaches used on BMP15 and GDF9 genes exhibited a robust divergence and a prompted evolution as compared to other TGFβ family members. Furthermore, among 32 mammalian species, we identified positive selection signals in the hominidae clade resulting to 132D, 147E, 163Y, 191W, and 236P codon sites of BMP15 and 162F, 188K, 206R, 240A, 244L, 246H, 248S, 251D, 253L, 254F and other codon sites of GDF9. The positively selected amino acid sites such as Alanine, Lucien, Arginine, and lysine are important for signaling. In conclusion, this study evidences that GDF9 and BMP15 genes have rapid evolution than other TGFß family members and was subjected to positive selection in the mammalian clade. Selected sites under the positive selection are of remarkable significance for the particular functioning of the protein and consequently for female fertility.     

Identifying significant genes and functionally enriched pathways in familial hypercholesterolemia using integrated gene co-expression network analysis

Familial hypercholesterolemia (FH) is a monogenic lipid disorder which promotes atherosclerosis and cardiovascular diseases. Owing to the lack of sufficient published information, this study aims to identify the potential genetic biomarkers for FH by studying the global gene expression profile of blood cells. The microarray expression data of FH patients and controls was analyzed by different computational biology methods like differential expression analysis, protein network mapping, hub gene identification, functional enrichment of biological pathways, and immune cell restriction analysis. Our results showed the dysregulated expression of 115 genes connected to lipid homeostasis, immune responses, cell adhesion molecules, canonical Wnt signaling, mucin type O-glycan biosynthesis pathways in FH patients. The findings from expanded protein interaction network construction with known FH genes and subsequent Gene Ontology (GO) annotations have also supported the above findings, in addition to identifying the involvement of dysregulated thyroid hormone and ErbB signaling pathways in FH patients. The genes like CSNK1A1, JAK3, PLCG2, RALA, and ZEB2 were found to be enriched under all GO annotation categories. The subsequent phenotype ontology results have revealed JAK3I, PLCG2, and ZEB2 as key hub genes contributing to the inflammation underlying cardiovascular and immune response related phenotypes. Immune cell restriction findings show that above three genes are highly expressed by T-follicular helper CD4⁺ T cells, naïve B cells, and monocytes, respectively. These findings not only provide a theoretical basis to understand the role of immune dysregulations underlying the atherosclerosis among FH patients but may also pave the way to develop genomic medicine for cardiovascular diseases.