The role of positive selection in determining the molecular cause of species differences in disease

Wakid M, A.M. Alkhaldi A, Hussain S, Malik MA, Saqib M, Qamar W, Alvi MA. Genetic diversity and population structure of Echinococcus multilocularis: An in-silico global analysis

Objective

Alveolar echinococcosis is caused by Echinococcus multilocularis, a parasite of zoonotic significance with a wide range of intermediate and final hosts, and the parasite survives successfully in diversified conditions. Plentiful studies have been done to study the genetic structure of the population of the parasite and the level of intimate kinship using mitochondrial (mt) DNA. The present study was conducted to investigate the population structure, genetic variation, and phylogenetic relationship of various isolates of E. multiocularis submitted to GenBank worldwide. Sequences of mt genes (mt-cytochrome c oxidase (cox1), mt-NADH dehydrogenase (nad1)) of E. multilocularis were analyzed to achieve the set goals.

Materials and Methods

A total of 275 and 124 gene sequences of mt-cox1 and mt-nad1 belonging to E. multilocularis, respectively, were retrieved from the National Center for Biotechnology Information GenBank. The retrieved sequences were subjected to alignment with respective reference sequences using MEGA software. The PopArt software was used to establish median-joining networks, while DnaSp was used to calculate neutrality and diversity indices. MrBayes software was used to investigate the phylogenetic association between haplotypes based on Bayesian phylogeny.

Results

Approximately 13 and 20 distinctive haplotypes of nad1 and cox1 genes, respectively, were observed in the present study. In both of the mt genes, diversity indices indicated low haplotype (mt-cox1 = 0.140; mt-nad1 = 0.374) and nucleotide (mt-cox1 = 0.00111; mt-nad1 = 0.00287) diversities. The values of Tajima's D and Fu Fs for a population of both of the genes under study were found to be negative.

Conclusion

This study is a maiden attempt to provide insights into the population structure and genetic variation of E. multilocularis on a global scale. However, it is suggested that to better understand the population structure and genetic diversity of E. multilocularis, more geographical locations and amplifications of full-length gene sequences should be considered, which could be helpful in widening the insights into the genetic diversity of E. multilocularis.

Global scenario of genetic diversity in cox1 and nad1 genes of Moniezia expansa

Monieziasis is a parasite-borne production-limiting disease of livestock. Moniezia expansa is the most important species having cosmopolitan distribution. Despite of numerous prevalence reports, very little information is available about the evolutionary biology and population genetics of M. expansa. To close this research gap, this study was undertaken to recognize and inspect the genetic variation of M. expansa populations around the world using the cox1 and nad1 genes and deduce phylogenetic relationships with M. expansa populations. The cox1 and nad1 gene sequences were downloaded from the NCBI GenBank database. Followed by sequence alignment, median-joining networks were constructed using PopArt software. Diversity and neutrality indices were computed through DnaSp software while MEGA software was used to draw the maximum-likelihood phylogenetic tree. Thirty-two cox1 sequences, from five different countries, and 9 nad1 sequences from three different countries, were among the sequences used in this study. The cox1 and nad1 gene sequences had mutations in 97 and 36 different places, respectively. Twenty and 7 unique haplotypes were discovered for the cox1 and nad1 gene sequences, respectively. Comparable haplotype diversities were observed for both the genes under study (cox1 = 0.950; nad1 = 0.944). Negative Tajima's D and Fu Fs were found for the cox1 gene while these indices were positive for the nad1 gene. Phylogenetic analysis also showed the existence of unique haplotypes for both the cox1 and nad1 genes. The results of this study indicate that there is the existence of a huge genetic diversity in M. expansa isolates. For future studies, it is recommended that longer gene sequences should be used to describe genetic variation among M. expansa isolates as the length of the gene under study affects the genetic variation. Moreover, additional mitochondrial markers should also be investigated because the assertive strength of a group of gene targets is superior to defining genetic diversity.

Ribonucleotide reductase as a therapeutic target for drug repurposing as anthelmintics

Visceral cestodiases, like echinococcoses and cysticercoses, are zoonoses of worldwide distribution and are responsible for public health problems in many countries, especially in underdeveloped regions. Current treatments have low efficiency and there are few drugs currently in use for chemotherapy, making the development of new anthelmintics an urgent matter. The nucleotide salvage pathways are the only ones available for nucleotide synthesis in cestodes and other parasitic helminths, and, here, we used in silico approaches to assess the potential of the enzymes in these pathways as targets for drug repurposing as anthelminthics. First, a genomic survey allowed to identify a repertoire of 28 enzymes of the purine and pyrimidine salvage pathways from the cestode Echinococcus granulosus sensu stricto. Regarding purines, the parasite relies on salvaging free bases rather than salvaging nucleosides. Pyrimidines, on the other hand, can be salvaged from both bases and nucleosides. Druggability of the parasite enzymes was assessed, as well as the availability of commercial inhibitors for them. Druggable enzymes were then ranked according to their potential for drug repurposing and the 17 most promising enzymes were selected for evolutionary analyses. The constructed phylogenetic trees allowed to assess the degree of conservation among ortholog enzymes from parasitic helminths and their mammalian hosts. Positive selection is absent in all assessed flatworm enzymes. A potential target enzyme for drug repurposing, ribonucleotide reductase (RNR), was selected for further assessment. RNR 3D-modelling showed structural similarities between the E. granulosus and the human orthologs suggesting that inhibitors of the human RNR should be effective against the E. granulosus enzyme. In line with that, E. granulosus protoscolices treated in vitro with the inhibitor hydroxyurea had their viability and DNA synthesis reduced. These results are consistent with nucleotide synthesis inhibition and confirm the potential of a nucleotide salvage inhibitors for repurposing as an anthelmintic.

Genetic variation and population structure of Fasciola hepatica: an in silico analysis

Fasciola hepatica is a trematode leading to heavy economic setbacks to the livestock sector globally. The population's genetic information and intimate kinship level are frequently assessed using analysis of mitochondrial DNA. In this analysis, we retrieved cox1 (n = 247) and nad1 (n = 357) sequences of F. hepatica from the NCBI GenBank database and aligned the sequences with the respective reference sequences using MEGA software. The median joining network was drawn using PopArt software while neutrality and diversity indices were estimated with the help of DnaSp software. Neighbor-joining phylogenetic tree was constructed using the MEGA software package. A total of 46 and 98 distinctive haplotypes were observed for cox1 and nad1 genes, respectively. Diversity indices indicated high haplotype and nucleotide diversities in both genes. Positive Tajima's D and Fu's Fs values were found for the entire population of both the genes under study. The cox1 and nad1 gene segments in this study showed high Tajima's D values, suggesting a low likelihood of future population growth. The Tajima's D value of the nad1 gene sequence is lower (2.14910) than that of the cox1 gene sequence (3.40314), which suggests that the former is growing at a slower rate. However, the region-wise analysis revealed that both the cox1 and nad1 genes showed deviation from neutrality suggesting a recent population expansion as a result of an excess of low-frequency polymorphism. Furthermore, the overall host-wise analysis showed positive and significant Tajima's D values for the cox1 and nad1 gene sequences. To the best of our knowledge, this is the first attempt to provide insights into genetic variations and population structure of F. hepatica at a global scale using cox1 and nad1 genes. Our findings suggest the existence of specific variants of F. hepatica in different parts of the world and provide information on the molecular ecology of F. hepatica. The results of this study also mark a critical development in upcoming epidemiological investigations on F. hepatica and will also contribute to understanding the global molecular epidemiology and population structure of F. hepatica.

The adaptive evolution of cancer driver genes

BACKGROUND

Cancer is a life-threatening disease in humans; yet, cancer genes are frequently reported to be under positive selection. This suggests an evolutionary-genetic paradox in which cancer evolves as a secondary product of selection in human beings. However, systematic investigation of the evolution of cancer driver genes is sparse.

RESULTS

Using comparative genomics analysis, population genetics analysis and computational molecular evolutionary analysis, the evolution of 568 cancer driver genes of 66 cancer types were evaluated at two levels, selection on the early evolution of humans (long timescale selection in the human lineage during primate evolution, i.e., millions of years), and recent selection in modern human populations (~ 100,000 years). Results showed that eight cancer genes covering 11 cancer types were under positive selection in the human lineage (long timescale selection). And 35 cancer genes covering 47 cancer types were under positive selection in modern human populations (recent selection). Moreover, SNPs associated with thyroid cancer in three thyroid cancer driver genes (CUX1, HERC2 and RGPD3) were under positive selection in East Asian and European populations, consistent with the high incidence of thyroid cancer in these populations.

CONCLUSIONS

These findings suggest that cancer can be evolved, in part, as a by-product of adaptive changes in humans. Different SNPs at the same locus can be under different selection pressures in different populations, and thus should be under consideration during precision medicine, especially for targeted medicine in specific populations.

The influence of 17q21.31 and APOE genetic ancestry on neurodegenerative disease risk

Advances in genomic research over the last two decades have greatly enhanced our knowledge concerning the genetic landscape and pathophysiological processes involved in multiple neurodegenerative diseases. However, current insights arise almost exclusively from studies on individuals of European ancestry. Despite this, studies have revealed that genetic variation differentially impacts risk for, and clinical presentation of neurodegenerative disease in non-European populations, conveying the importance of ancestry in predicting disease risk and understanding the biological mechanisms contributing to neurodegeneration. We review the genetic influence of two important disease-associated loci, 17q21.31 (the "MAPT locus") and APOE, to neurodegenerative disease risk in non-European populations, touching on global population differences and evolutionary genetics by ancestry that may underlie some of these differences. We conclude there is a need to increase representation of non-European ancestry individuals in genome-wide association studies (GWAS) and biomarker analyses in order to help resolve existing disparities in understanding risk for, diagnosis of, and treatment for neurodegenerative diseases in diverse populations.

Identification of multiple TAR DNA binding protein retropseudogene lineages during the evolution of primates

The TAR DNA Binding Protein (TARDBP) gene has become relevant after the discovery of its several pathogenic mutations. The lack of evolutionary history is in contrast to the amount of studies found in the literature. This study investigated the evolutionary dynamics associated with the retrotransposition of the TARDBP gene in primates. We identified novel retropseudogenes that likely originated in the ancestors of anthropoids, catarrhines, and lemuriformes, i.e. the strepsirrhine clade that inhabit Madagascar. We also found species-specific retropseudogenes in the Philippine tarsier, Bolivian squirrel monkey, capuchin monkey and vervet. The identification of a retropseudocopy of the TARDBP gene overlapping a lncRNA that is potentially expressed opens a new avenue to investigate TARDBP gene regulation, especially in the context of TARDBP associated pathologies.

Systems identification and characterization of β-glucuronosyltransferase genes involved in arabinogalactan-protein biosynthesis in plant genomes

Utilizing plant biomass for bioethanol production requires an understanding of the molecular mechanisms involved in plant cell wall assembly. Arabinogalactan-proteins (AGPs) are glycoproteins that interact with other cell wall polymers to influence plant growth and developmental processes. Glucuronic acid, which is transferred to the AGP glycan by β-glucuronosyltransferases (GLCATs), is the only acidic sugar in AGPs with the ability to bind calcium. We carried out a comprehensive genome-wide analysis of a putative GLCAT gene family involved in AGP biosynthesis by examining its sequence diversity, genetic architecture, phylogenetic and motif characteristics, selection pressure and gene expression in plants. We report the identification of 161 putative GLCAT genes distributed across 14 plant genomes and a widely conserved GLCAT catalytic domain. We discovered a phylogenetic clade shared between bryophytes and higher land plants of monocot grass and dicot lineages and identified positively selected sites that do not result in functional divergence of GLCATs. RNA-seq and microarray data analyses of the putative GLCAT genes revealed gene expression signatures that likely influence the assembly of plant cell wall polymers which is critical to the overall growth and development of edible and bioenergy crops.

Cestode strobilation: prediction of developmental genes and pathways

BACKGROUND

Cestoda is a class of endoparasitic worms in the flatworm phylum (Platyhelminthes). During the course of their evolution cestodes have evolved some interesting aspects, such as their increased reproductive capacity. In this sense, they have serial repetition of their reproductive organs in the adult stage, which is often associated with external segmentation in a developmental process called strobilation. However, the molecular basis of strobilation is poorly understood. To assess this issue, an evolutionary comparative study among strobilated and non-strobilated flatworm species was conducted to identify genes and proteins related to the strobilation process.

RESULTS

We compared the genomic content of 10 parasitic platyhelminth species; five from cestode species, representing strobilated parasitic platyhelminths, and five from trematode species, representing non-strobilated parasitic platyhelminths. This dataset was used to identify 1813 genes with orthologues that are present in all cestode (strobilated) species, but absent from at least one trematode (non-strobilated) species. Development-related genes, along with genes of unknown function (UF), were then selected based on their transcriptional profiles, resulting in a total of 34 genes that were differentially expressed between the larval (pre-strobilation) and adult (strobilated) stages in at least one cestode species. These 34 genes were then assumed to be strobilation related; they included 12 encoding proteins of known function, with 6 related to the Wnt, TGF-β/BMP, or G-protein coupled receptor signaling pathways; and 22 encoding UF proteins. In order to assign function to at least some of the UF genes/proteins, a global gene co-expression analysis was performed for the cestode species Echinococcus multilocularis. This resulted in eight UF genes/proteins being predicted as related to developmental, reproductive, vesicle transport, or signaling processes.

CONCLUSIONS

Overall, the described in silico data provided evidence of the involvement of 34 genes/proteins and at least 3 developmental pathways in the cestode strobilation process. These results highlight on the molecular mechanisms and evolution of the cestode strobilation process, and point to several interesting proteins as potential developmental markers and/or targets for the development of novel antihelminthic drugs.

APOE Alleles and Diet in Brain Aging and Alzheimer's Disease

The APOE gene alleles modify human aging and the response to the diet at many levels with diverse pleotropic effects from gut to brain. To understand the interactions of APOE isoforms and diet, we analyze how cellular trafficking of apoE proteins affects energy metabolism, the immune system, and reproduction. The age-accelerating APOE4 allele alters the endosomal trafficking of cell surface receptors that mediate lipid and glucose metabolism. The APOE4 allele is the ancestral human allele, joined by APOE3 and then APOE2 in the human species. Under conditions of high infection, uncertain food, and shorter life expectancy, APOE4 may be adaptive for reducing mortality. As humans transitioned into modern less-infectious environments and longer life spans, APOE4 increased risks of aging-related diseases, particularly impacting arteries and the brain. The association of APOE4 with glucose dysregulation and body weight promotes many aging-associated diseases. Additionally, the APOE gene locus interacts with adjacent genes on chromosome 19 in haplotypes that modify neurodegeneration and metabolism, for which we anticipate complex gene-environment interactions. We summarize how diet and Alzheimer's disease (AD) risk are altered by APOE genotype in both animal and human studies and identify gaps. Much remains obscure in how APOE alleles modify nutritional factors in human aging. Identifying risk variant haplotypes in the APOE gene complex will clarify homeostatic adaptive responses to environmental conditions.

A logical relationship for schizophrenia, bipolar, and major depressive disorder. Part 1: Evidence from chromosome 1 high density association screen

Familial clustering of schizophrenia (SCZ), bipolar disorder (BPD), and major depressive disorder (MDD) was investigated systematically (Aukes et al., Genetics in Medicine, 2012, 14, 338-341) and any two or even three of these disorders could coexist in some families. Furthermore, evidence from symptomatology and psychopharmacology also imply the existence of intrinsic connections between these three major psychiatric disorders. A total of 71,445 SNPs on chromosome 1 were genotyped on 119 SCZ, 253 BPD (type-I), 177 MDD cases and 1000 controls and further validated in 986 SCZ patients in the population of Shandong province of China. Outstanding psychosis genes are systematically revealed( ATP1A4, ELTD1, FAM5C, HHAT, KIF26B, LMX1A, NEGR1, NFIA, NR5A2, NTNG1, PAPPA2, PDE4B, PEX14, RYR2, SYT6, TGFBR3, TTLL7, and USH2A). Unexpectedly, flanking genes for up to 97.09% of the associated SNPs were also replicated in an enlarged cohort of 986 SCZ patients. From the perspective of etiological rather than clinical psychiatry, bipolar, and major depressive disorder could be subtypes of schizophrenia. Meanwhile, the varied clinical feature and prognosis might be the result of interaction of genetics and epigenetics, for example, irreversible or reversible shut down, and over or insufficient expression of certain genes, which may gives other aspects of these severe mental disorders.

Antagonistic Pleiotropy in Human Disease

Between the 1930s and 1950s, scientists developed key principles of population genetics to try and explain the aging process. Almost a century later, these aging theories, including antagonistic pleiotropy and mutation accumulation, have been experimentally validated in animals. Although the theories have been much harder to test in humans despite research dating back to the 1970s, recent research is closing this evidence gap. Here we examine the strength of evidence for antagonistic pleiotropy in humans, one of the leading evolutionary explanations for the retention of genetic risk variation for non-communicable diseases. We discuss the analytical tools and types of data that are used to test for patterns of antagonistic pleiotropy and provide a primer of evolutionary theory on types of selection as a guide for understanding this mechanism and how it may manifest in other diseases. We find an abundance of non-experimental evidence for antagonistic pleiotropy in many diseases. In some cases, several studies have independently found corroborating evidence for this mechanism in the same or related sets of diseases including cancer and neurodegenerative diseases. Recent studies also suggest antagonistic pleiotropy may be involved in cardiovascular disease and diabetes. There are also compelling examples of disease risk variants that confer fitness benefits ranging from resistance to other diseases or survival in extreme environments. This provides increasingly strong support for the theory that antagonistic pleiotropic variants have enabled improved fitness but have been traded for higher burden of disease later in life. Future research in this field is required to better understand how this mechanism influences contemporary disease and possible consequences for their treatment.

Positive Selection of Squalene Synthase in Cucurbitaceae Plants

Triterpenoid saponins are secondary metabolites synthesized through isoprenoid pathways in plants. Cucurbitaceae represent an important plant family in which many species contain cucurbitacins as secondary metabolites synthesized through isoprenoid and triterpenoid pathways. Squalene synthase (SQS) is required for the biosynthesis of isoprenoids, but the forces driving the evolution of SQS remain undetermined. In this study, 10 SQS cDNA sequences cloned from 10 species of Cucurbitaceae and 49 sequences of SQS downloaded from GenBank and UniProt databases were analyzed in a phylogenetic framework to identify the evolutionary forces for functional divergence. Through phylogenetic construction and positive selection analysis, we found that SQS sequences are under positive selection. The sites of positive selection map to functional and transmembrane domains. 180L, 189S, 194S, 196S, 265I, 289P, 389P, 390T, 407S, 408A, 410R, and 414N were identified as sites of positive selection that are important during terpenoid synthesis and map to transmembrane domains. 196S and 407S are phosphorylated and influence SQS catalysis and triterpenoid accumulation. These results reveal that positive selection is an important evolutionary force for SQS in plants. This provides new information into the molecular evolution of SQS within the Cucurbitaceae family.

Selection Acting on Genomes

Populations evolve as mutations arise in individual organisms and, through hereditary transmission, may become "fixed" (shared by all individuals) in the population. Most mutations are lethal or have negative fitness consequences for the organism. Others have essentially no effect on organismal fitness and can become fixed through the neutral stochastic process known as random drift. However, mutations may also produce a selective advantage that boosts their chances of reaching fixation. Regions of genomes where new mutations are beneficial, rather than neutral or deleterious, tend to evolve more rapidly due to positive selection. Genes involved in immunity and defense are a well-known example; rapid evolution in these genes presumably occurs because new mutations help organisms to prevail in evolutionary "arms races" with pathogens. In recent years genome-wide scans for selection have enlarged our understanding of the genome evolution of various species. In this chapter, we will focus on methods to detect selection on the genome. In particular, we will discuss probabilistic models and how they have changed with the advent of new genome-wide data now available.

Adaptive molecular evolution of MC1R gene reveals the evidence for positive diversifying selection in indigenous goat populations

Detecting signatures of selection can provide a new insight into the mechanism of contemporary breeding and artificial selection and further reveal the causal genes associated to the phenotypic variation. However, the signatures of selection on genes entailing for profitable traits between Chinese commercial and indigenous goats have been poorly interpreted. We noticed footprints of positive selection at MC1R gene containing SNPs genotyped in five Chinese native goat breeds. An experimental distribution of F_ST was built based on approximations of F_ST for each SNP across five breeds. We identified selection using the high F_ST outlier method and found that MC1R candidate gene show evidence of positive selection. Furthermore, adaptive selection pressure on specific codons was determined using different codon based on maximum‐likelihood methods; signature of positive selection in mammalian MC1R was explored in individual codons. Evolutionary analyses were inferred under maximum likelihood models, the HyPhy package implemented in the DATAMONKEY Web Server. The results of codon selection displayed positive diversifying selection at the sites were mainly involved in development of genetic variations in coat color in various mammalian species. Positive diversifying selection inferred with recent evolutionary changes in domesticated goat MC1R provides new insights that the gene evolution may have been modulated by domestication events in goats.

Evolution of Brain Active Gene Promoters in Human Lineage Towards the Increased Plasticity of Gene Regulation

Adaptability to a variety of environmental conditions is a prominent feature of Homo sapiens. We hypothesize that this feature can be explained by evolutionary changes in gene promoters active in the brain prefrontal cortex leading to a more flexible gene regulation network. The genotype-dependent range of gene expression can be broader in humans than in other higher primates. Thus, we searched for specific signatures of evolutionary changes in promoter architectures of multiple hominid genes, including the genes active in human cortical neurons that may indicate an increase of variability of gene expression rather than just changes in the level of expression, such as downregulation or upregulation of the genes. We performed a whole-genome search for genetic-based alterations that may impact gene regulation "flexibility" in a process of hominids evolution, such as (i) CpG dinucleotide content, (ii) predicted nucleosome-DNA dissociation constant, and (iii) predicted affinities for TATA-binding protein (TBP) in gene promoters. We tested all putative promoter regions across the human genome and especially gene promoters in active chromatin state in neurons of prefrontal cortex, the brain region critical for abstract thinking and social and behavioral adaptation. Our data imply that the origin of modern man has been associated with an increase of flexibility of promoter-driven gene regulation in brain. In contrast, after splitting from the ancestral lineages of H. sapiens, the evolution of ape species is characterized by reduced flexibility of gene promoter functioning, underlying reduced variability of the gene expression.

Positive selection and functional divergence of farnesyl pyrophosphate synthase genes in plants

Background

Farnesyl pyrophosphate synthase (FPS) belongs to the short-chain prenyltransferase family, and it performs a conserved and essential role in the terpenoid biosynthesis pathway. However, its classification, evolutionary history, and the forces driving the evolution of FPS genes in plants remain poorly understood.

Results

Phylogeny and positive selection analysis was used to identify the evolutionary forces that led to the functional divergence of FPS in plants, and recombinant detection was undertaken using the Genetic Algorithm for Recombination Detection (GARD) method. The dataset included 68 FPS variation pattern sequences (2 gymnosperms, 10 monocotyledons, 54 dicotyledons, and 2 outgroups). This study revealed that the FPS gene was under positive selection in plants. No recombinant within the FPS gene was found. Therefore, it was inferred that the positive selection of FPS had not been influenced by a recombinant episode. The positively selected sites were mainly located in the catalytic center and functional areas, which indicated that the 98S and 234D were important positively selected sites for plant FPS in the terpenoid biosynthesis pathway. They were located in the FPS conserved domain of the catalytic site. We inferred that the diversification of FPS genes was associated with functional divergence and could be driven by positive selection.

Conclusions

It was clear that protein sequence evolution via positive selection was able to drive adaptive diversification in plant FPS proteins. This study provides information on the classification and positive selection of plant FPS genes, and the results could be useful for further research on the regulation of triterpenoid biosynthesis.

Electronic supplementary material

The online version of this article (doi:10.1186/s12867-017-0081-4) contains supplementary material, which is available to authorized users.

Genetic Variation of Goat Interferon Regulatory Factor 3 Gene and Its Implication in Goat Evolution

The immune systems are fundamentally vital for evolution and survival of species; as such, selection patterns in innate immune loci are of special interest in molecular evolutionary research. The interferon regulatory factor (IRF) gene family control many different aspects of the innate and adaptive immune responses in vertebrates. Among these, IRF3 is known to take active part in very many biological processes. We assembled and evaluated 1356 base pairs of the IRF3 gene coding region in domesticated goats from Africa (Nigeria, Ethiopia and South Africa) and Asia (Iran and China) and the wild goat (Capra aegagrus). Five segregating sites with θ value of 0.0009 for this gene demonstrated a low diversity across the goats’ populations. Fu and Li tests were significantly positive but Tajima’s D test was significantly negative, suggesting its deviation from neutrality. Neighbor joining tree of IRF3 gene in domesticated goats, wild goat and sheep showed that all domesticated goats have a closer relationship than with the wild goat and sheep. Maximum likelihood tree of the gene showed that different domesticated goats share a common ancestor and suggest single origin. Four unique haplotypes were observed across all the sequences, of which, one was particularly common to African goats (MOCH-K14-0425, Poitou and WAD). In assessing the evolution mode of the gene, we found that the codon model d_N/d_S ratio for all goats was greater than one. Phylogenetic Analysis by Maximum Likelihood (PAML) gave a ω₀ (d_N/d_S) value of 0.067 with LnL value of -6900.3 for the first Model (M1) while ω2 = 1.667 in model M2 with LnL value of -6900.3 with positive selection inferred in 3 codon sites. Mechanistic empirical combination (MEC) model for evaluating adaptive selection pressure on particular codons also confirmed adaptive selection pressure in three codons (207, 358 and 408) in IRF3 gene. Positive diversifying selection inferred with recent evolutionary changes in domesticated goat IRF3 led us to conclude that the gene evolution may have been influenced by domestication processes in goats.

Apolipoprotein E and Sex Bias in Cerebrovascular Aging of Men and Mice

Alzheimer disease (AD) research has mainly focused on neurodegenerative processes associated with the classic neuropathologic markers of senile plaques and neurofibrillary tangles. Additionally, cerebrovascular contributions to dementia are increasingly recognized, particularly from cerebral small vessel disease (SVD). Remarkably, in AD brains, the apolipoprotein E (ApoE) ɛ4 allele shows male excess for cerebral microbleeds (CMBs), a marker of SVD, which is opposite to the female excess of plaques and tangles. Mouse transgenic models add further complexities to sex-ApoE ɛ4 allele interactions, with female excess of both CMBs and brain amyloid. We conclude that brain aging and AD pathogenesis cannot be understood in humans without addressing major gaps in the extent of sex differences in cerebrovascular pathology.

Functional test of PCDHB11, the most human-specific neuronal surface protein

Background

Brain-expressed proteins that have undergone functional change during human evolution may contribute to human cognitive capacities, and may also leave us vulnerable to specifically human diseases, such as schizophrenia, autism or Alzheimer’s disease. In order to search systematically for those proteins that have changed the most during human evolution and that might contribute to brain function and pathology, all proteins with orthologs in chimpanzee, orangutan and rhesus macaque and annotated as being expressed on the surface of cells in the human central nervous system were ordered by the number of human-specific amino acid differences that are fixed in modern populations.

Results

PCDHB11, a beta-protocadherin homologous to murine cell adhesion proteins, stood out with 12 substitutions and maintained its lead after normalizing for protein size and applying weights for amino acid exchange probabilities. Human PCDHB11 was found to cause homophilic cell adhesion, but at lower levels than shown for other clustered protocadherins. Homophilic adhesion caused by a PCDHB11 with reversion of human-specific changes was as low as for modern human PCDHB11; while neither human nor reverted PCDHB11 adhered to controls, they did adhere to each other. A loss of function in PCDHB11 is unlikely because intra-human variability did not increase relative to the other human beta-protocadherins.

Conclusions

The brain-expressed protein with the highest number of human-specific substitutions is PCDHB11. In spite of its fast evolution and low intra-human variability, cell-based tests on the only proposed function for PCDHB11 did not indicate a functional change.

The APOE4 allele shows opposite sex bias in microbleeds and Alzheimer's disease of humans and mice

The apolipoprotein APOE4 allele confers greater risk of Alzheimer's disease (AD) for women than men, in conjunction with greater clinical deficits per unit of AD neuropathology (plaques, tangles). Cerebral microbleeds, which contribute to cognitive dysfunctions during AD, also show APOE4 excess, but sex-APOE allele interactions are not described. We report that elderly men diagnosed for mild cognitive impairment and AD showed a higher risk of cerebral cortex microbleeds with APOE4 allele dose effect in 2 clinical cohorts (ADNI and KIDS). Sex-APOE interactions were further analyzed in EFAD mice carrying human APOE alleles and familial AD genes (5XFAD (+/-) /human APOE(+/+)). At 7 months, E4FAD mice had cerebral cortex microbleeds with female excess, in contrast to humans. Cerebral amyloid angiopathy, plaques, and soluble Aβ also showed female excess. Both the cerebral microbleeds and cerebral amyloid angiopathy increased in proportion to individual Aβ load. In humans, the opposite sex bias of APOE4 allele for microbleeds versus the plaques and tangles is the first example of organ-specific, sex-linked APOE allele effects, and further shows AD as a uniquely human condition.

Genomic Patterns of Positive Selection at the Origin of Rust Fungi

Understanding the origin and evolution of pathogenicity and biotrophic life-style of rust fungi has remained a conundrum for decades. Research on the molecular mechanisms responsible for rust fungi evolution has been hampered by their biotrophic life-style until the sequencing of some rust fungi genomes. With the availability of multiple whole genomes and EST data for this group, it is now possible to employ genome-wide surveys and investigate how natural selection shaped their evolution. In this work, we employed a phylogenomics approach to search for positive selection and genes undergoing accelerated evolution at the origin of rust fungi on an assembly of single copy genes conserved across a broad range of basidiomycetes. Up to 985 genes were screened for positive selection on the phylogenetic branch leading to rusts, revealing a pervasive signal of positive selection throughout the data set with the proportion of positively selected genes ranging between 19.6–33.3%. Additionally, 30 genes were found to be under accelerated evolution at the origin of rust fungi, probably due to a mixture of positive selection and relaxation of purifying selection. Functional annotation of the positively selected genes revealed an enrichment in genes involved in the biosynthesis of secondary metabolites and several metabolism and transporter classes.

The transcriptome of the bowhead whale Balaena mysticetus reveals adaptations of the longest-lived mammal

Mammals vary dramatically in lifespan, by at least two-orders of magnitude, but the molecular basis for this difference remains largely unknown. The bowhead whale Balaena mysticetus is the longest-lived mammal known, with an estimated maximal lifespan in excess of two hundred years. It is also one of the two largest animals and the most cold-adapted baleen whale species. Here, we report the first genome-wide gene expression analyses of the bowhead whale, based on the de novo assembly of its transcriptome. Bowhead whale or cetacean-specific changes in gene expression were identified in the liver, kidney and heart, and complemented with analyses of positively selected genes. Changes associated with altered insulin signaling and other gene expression patterns could help explain the remarkable longevity of bowhead whales as well as their adaptation to a lipid-rich diet. The data also reveal parallels in candidate longevity adaptations of the bowhead whale, naked mole rat and Brandt's bat. The bowhead whale transcriptome is a valuable resource for the study of this remarkable animal, including the evolution of longevity and its important correlates such as resistance to cancer and other diseases.

Darwin and Fisher meet at biotech: on the potential of computational molecular evolution in industry

Background

Today computational molecular evolution is a vibrant research field that benefits from the availability of large and complex new generation sequencing data – ranging from full genomes and proteomes to microbiomes, metabolomes and epigenomes. The grounds for this progress were established long before the discovery of the DNA structure. Specifically, Darwin’s theory of evolution by means of natural selection not only remains relevant today, but also provides a solid basis for computational research with a variety of applications. But a long-term progress in biology was ensured by the mathematical sciences, as exemplified by Sir R. Fisher in early 20th century. Now this is true more than ever: The data size and its complexity require biologists to work in close collaboration with experts in computational sciences, modeling and statistics.

Results

Natural selection drives function conservation and adaptation to emerging pathogens or new environments; selection plays key role in immune and resistance systems. Here I focus on computational methods for evaluating selection in molecular sequences, and argue that they have a high potential for applications. Pharma and biotech industries can successfully use this potential, and should take the initiative to enhance their research and development with state of the art bioinformatics approaches.

Conclusions

This review provides a quick guide to the current computational approaches that apply the evolutionary principles of natural selection to real life problems – from drug target validation, vaccine design and protein engineering to applications in agriculture, ecology and conservation.

Characterization of the avian Trojan gene family reveals contrasting evolutionary constraints

"Trojan" is a leukocyte-specific, cell surface protein originally identified in the chicken. Its molecular function has been hypothesized to be related to anti-apoptosis and the proliferation of immune cells. The Trojan gene has been localized onto the Z sex chromosome. The adjacent two genes also show significant homology to Trojan, suggesting the existence of a novel gene/protein family. Here, we characterize this Trojan family, identify homologues in other species and predict evolutionary constraints on these genes. The two Trojan-related proteins in chicken were predicted as a receptor-type tyrosine phosphatase and a transmembrane protein, bearing a cytoplasmic immuno-receptor tyrosine-based activation motif. We identified the Trojan gene family in ten other bird species and found related genes in three reptiles and a fish species. The phylogenetic analysis of the homologues revealed a gradual diversification among the family members. Evolutionary analyzes of the avian genes predicted that the extracellular regions of the proteins have been subjected to positive selection. Such selection was possibly a response to evolving interacting partners or to pathogen challenges. We also observed an almost complete lack of intracellular positively selected sites, suggesting a conserved signaling mechanism of the molecules. Therefore, the contrasting patterns of selection likely correlate with the interaction and signaling potential of the molecules.

The molecular evolutionary dynamics of the vomeronasal receptor (class 1) genes in primates: a gene family on the verge of a functional breakdown

Olfaction plays a critical role in both survival of the individual and in the propagation of species. Studies from across the mammalian clade have found a remarkable correlation between organismal lifestyle and molecular evolutionary properties of receptor genes in both the main olfactory system (MOS) and the vomeronasal system (VNS). When a large proportion of intact (and putatively functional) copies is observed, the inference is made that a particular mode of chemoreception is critical for an organism’s fit to its environment and is thus under strong positive selection. Conversely, when the receptors in question show a disproportionately large number of pseudogene copies, this contraction is interpreted as evidence of relaxed selection potentially leading to gene family extinction. Notably, it appears that a risk factor for gene family extinction is a high rate of nonsynonymous substitution. A survey of intact vs. pseudogene copies among primate vomeronasal receptor Class one genes (V1Rs) appears to substantiate this hypothesis. Molecular evolutionary complexities in the V1R gene family combine rapid rates of gene duplication, gene conversion, lineage-specific expansions, deletions, and/or pseudogenization. An intricate mix of phylogenetic footprints and current adaptive landscapes have left their mark on primate V1Rs suggesting that the primate clade offers an ideal model system for exploring the molecular evolutionary and functional properties of the VNS of mammals. Primate V1Rs tell a story of ancestral function and divergent selection as species have moved into ever diversifying adaptive regimes. The sensitivity to functional collapse in these genes, consequent to their precariously high rates of nonsynonymous substitution, confer a remarkable capacity to reveal the lifestyles of the genomes that they presently occupy as well as those of their ancestors.

Molecular evolution in the CREB1 signal pathway and a rare haplotype in CREB1 with genetic predisposition to schizophrenia

CREB1 is a cAMP responsive transcriptional factor which plays a key role in neural development. CREB1 signal pathway (CSP) has been implicated repeatedly in studies of predisposition for schizophrenia. We speculated that CSP has undergone positive selection during evolution of modern human and some genes that have undergone natural selection in the past may predispose to schizophrenia (SCZ) in modern time. Positive selection and association analysis were employed to explore the molecular evolution of CSP and association with schizophrenia. Our results showed a pan-ethnic selection event on NRG1 and CREB1, as confirmed in all 14 ethnic populations studied, which also suggested a selection process occurred before the "Out of Africa" scenario. Analysis of 62 SNPs covering 6 CSP genes in 2019 Han Chinese (976 SCZ patients and 1043 healthy individuals) showed an association of two SNPs (rs4379857, P = 0.009, OR [95% CI]: 1.200 [1.379-1.046]; rs2238751, P = 0.023, OR [95% CI]: 1.253 [1.522-1.032]) with SCZ. However, none of these significances survived after multiple testing corrections. Nonetheless, we observed an association of a rare CREB1 haplotype CCGGC (Bonferroni corrected P = 1.74 × 10(-5)) with SCZ. Our study showed that there was substantial population heterogeneity in genetic predisposition to SCZ, and different genes in the CSP pathway may predispose to SCZ in different populations.

Patterns of positive selection in seven ant genomes

The evolution of ants is marked by remarkable adaptations that allowed the development of very complex social systems. To identify how ant-specific adaptations are associated with patterns of molecular evolution, we searched for signs of positive selection on amino-acid changes in proteins. We identified 24 functional categories of genes which were enriched for positively selected genes in the ant lineage. We also reanalyzed genome-wide data sets in bees and flies with the same methodology to check whether positive selection was specific to ants or also present in other insects. Notably, genes implicated in immunity were enriched for positively selected genes in the three lineages, ruling out the hypothesis that the evolution of hygienic behaviors in social insects caused a major relaxation of selective pressure on immune genes. Our scan also indicated that genes implicated in neurogenesis and olfaction started to undergo increased positive selection before the evolution of sociality in Hymenoptera. Finally, the comparison between these three lineages allowed us to pinpoint molecular evolution patterns that were specific to the ant lineage. In particular, there was ant-specific recurrent positive selection on genes with mitochondrial functions, suggesting that mitochondrial activity was improved during the evolution of this lineage. This might have been an important step toward the evolution of extreme lifespan that is a hallmark of ants.

Medical relevance of UK-funded non-human primate research published from January 1997 to July 2012

In 2012, the Bateson Review of research using non-human primates (NHPs) recommended the commissioning of a working group to identify and follow-up the results of UK-funded NHP research of potential benefit for human health (Recommendation 4), but the Medical Research Council (MRC) has postponed implementation of the recommendation. Information on results and potential benefits of NHP research therefore remains unavailable. To fill this gap in knowledge, this study identified all published NHP research studies funded by the MRC, Wellcome Trust and Biotechnology and Biological Sciences Research Council (BBSRC) from January 1997 to July 2012 and assessed full texts for medical relevance. In total, 284 papers were identified, of which 51 (18%) involved invasive NHP research, compared to 176 (61%) which used NHP tissue and cell lines, indicating a shift in research emphasis from invasive whole animal to cell-based research. Of these studies, 98 (35%) were medically relevant, of which 22 had potential therapeutic or public health applications. The relatively low proportion of medical studies together with the small number of applied studies raises questions over the level of investment in medical research and the effectiveness of knowledge transfer from basic to applied research. Implementation of the Bateson Review's Recommendation 4 would address these questions.

Translational neuropharmacology: the use of human isolated gastrointestinal tissues

Translational sciences increasingly emphasize the measurement of functions in native human tissues. However, such studies must confront variations in patient age, gender, genetic background and disease. Here, these are discussed with reference to neuromuscular and neurosecretory functions of the human gastrointestinal (GI) tract. Tissues are obtained after informed consent, in collaboration with surgeons (surgical techniques help minimize variables) and pathologists. Given the difficulties of directly recording from human myenteric neurones (embedded between muscle layers), enteric motor nerve functions are studied by measuring muscle contractions/relaxations evoked by electrical stimulation of intrinsic nerves; responses are regionally dependent, often involving cholinergic and nitrergic phenotypes. Enteric sensory functions can be studied by evoking the peristaltic reflex, involving enteric sensory and motor nerves, but this has rarely been achieved. As submucosal neurones are more accessible (after removing the mucosa), direct neuronal recordings are possible. Neurosecretory functions are studied by measuring changes in short-circuit current across the mucosa. For all experiments, basic questions must be addressed. Because tissues are from patients, what are the controls and the influence of disease? How long does it take before function fully recovers? What is the impact of age- and gender-related differences? What is the optimal sample size? Addressing these and other questions minimizes variability and raises the scientific credibility of human tissue research. Such studies also reduce animal use. Further, the many differences between animal and human GI functions also means that human tissue research must question the ethical validity of using strains of animals with unproved translational significance.

Functional characterisation of bovine TLR5 indicates species-specific recognition of flagellin

Mammalian toll-like receptor 5 (TLR5) senses flagellin of several bacterial species and has been described to activate the innate immune system. To assess the role of bovine TLR5 (boTLR5) in the cattle system, we cloned and successfully expressed boTLR5 in human embryonic kidney (HEK) 293 cells, as indicated by quantitative PCR and confocal microscopy. However, in contrast to huTLR5-transfected cells, exposure of boTLR5-transfected cells to flagellin neither activated nuclear factor kappa-light-chain-enhancer of activated B cells (NF-κB) nor CXCL8 production. Subsequent comparison of the flagellin response induced in human and bovine primary macrophages revealed that flagellin did not lead to phosphorylation of major signalling molecules. Furthermore, the CXCL8 and TNFα response of primary bovine macrophages stimulated with flagellin was very low compared to that observed in human primary macrophages. Our results indicate that cattle express a functional TLR5 albeit with different flagellin sensing qualities compared to human TLR5. However, boTLR5 seemed to play a different role in the bovine system compared to the human system in recognizing flagellin, and other potentially intracellular expressed receptors may play a more important role in the bovine system to detect flagellin.

Immune-related functions of the Hivep gene family in East African cichlid fishes

Immune-related genes are often characterized by adaptive protein evolution. Selection on immune genes can be particularly strong when hosts encounter novel parasites, for instance, after the colonization of a new habitat or upon the exploitation of vacant ecological niches in an adaptive radiation. We examined a set of new candidate immune genes in East African cichlid fishes. More specifically, we studied the signatures of selection in five paralogs of the human immunodeficiency virus type I enhancer-binding protein (Hivep) gene family, tested their involvement in the immune defense, and related our results to explosive speciation and adaptive radiation events in cichlids. We found signatures of long-term positive selection in four Hivep paralogs and lineage-specific positive selection in Hivep3b in two radiating cichlid lineages. Exposure of the cichlid Astatotilapia burtoni to a vaccination with Vibrio anguillarum bacteria resulted in a positive correlation between immune response parameters and expression levels of three Hivep loci. This work provides the first evidence for a role of Hivep paralogs in teleost immune defense and links the signatures of positive selection to host-pathogen interactions within an adaptive radiation.

Expression of heat shock protein 90 genes during early development and infection in Megalobrama amblycephala and evidence for adaptive evolution in teleost

Heat shock protein 90 (HSP90) are highly conserved molecular chaperones, playing a pivotal role in cellular progress. In this study, we reported the characterization of the Hsp90α and Hsp90β genes in Megalobrama amblycephala, the expression profiling during early development, in various healthy tissues and in response to bacterial infection, and the assessment of their adaptive evolution. The Hsp90α cDNA contains an open reading frame (ORF) of 2193 bp encoding 731 amino acids and the Hsp90β cDNA has an ORF of 2184 bp encoding 728 amino acids. Using quantitative real-time PCR (qRT-PCR) analysis, the mRNA of both Hsp90α and Hsp90β reached the highest level at 15th day post-hatch. Using qRT-PCR and Western blot, both Hsp90α and Hsp90β were widely expressed in various healthy tissues and significantly higher in blood than in other tissues. Expression of both Hsp90α and Hsp90β were up-regulated upon bacterial infection and reached the peak level at 4 h post infection. Site model analysis indicated that one positive selection site (T717) in Hsp90α was found, while no positive selection site was observed in Hsp90β. Branch-site model test showed that there were adaptively evolutionary evidences in the branches of Salmoniformes and Gasterosteiformes for Hsp90α gene, and in the branch of Salmoniformes for Hsp90β gene.

Functional proteomics, human genetics and cancer biology of GIPC family members

GIPC1, GIPC2 and GIPC3 consist of GIPC homology 1 (GH1) domain, PDZ domain and GH2 domain. The regions around the GH1 and GH2 domains of GIPC1 are involved in dimerization and interaction with myosin VI (MYO6), respectively. The PDZ domain of GIPC1 is involved in interactions with transmembrane proteins [IGF1R, NTRK1, ADRB1, DRD2, TGFβR3 (transforming growth factorβ receptor type III), SDC4, SEMA4C, LRP1, NRP1, GLUT1, integrin α5 and VANGL2], cytosolic signaling regulators (APPL1 and RGS19) and viral proteins (HBc and HPV-18 E6). GIPC1 is an adaptor protein with dimerizing ability that loads PDZ ligands as cargoes for MYO6-dependent endosomal trafficking. GIPC1 is required for cell-surface expression of IGF1R and TGFβR3. GIPC1 is also required for integrin recycling during cell migration, angiogenesis and cytokinesis. On early endosomes, GIPC1 assembles receptor tyrosine kinases (RTKs) and APPL1 for activation of PI3K-AKT signaling, and G protein-coupled receptors (GPCRs) and RGS19 for attenuation of inhibitory Gα signaling. GIPC1 upregulation in breast, ovarian and pancreatic cancers promotes tumor proliferation and invasion, whereas GIPC1 downregulation in cervical cancer with human papillomavirus type 18 infection leads to resistance to cytostatic transforming growth factorβ signaling. GIPC2 is downregulated in acute lymphocytic leukemia owing to epigenetic silencing, while Gipc2 is upregulated in estrogen-induced mammary tumors. Somatic mutations of GIPC2 occur in malignant melanoma, and colorectal and ovarian cancers. Germ-line mutations of the GIPC3 or MYO6 gene cause nonsyndromic hearing loss. As GIPC proteins are involved in trafficking, signaling and recycling of RTKs, GPCRs, integrins and other transmembrane proteins, dysregulation of GIPCs results in human pathologies, such as cancer and hereditary deafness.

Natural selection and functional diversification of the epidermal growth factor receptor EGFR family in vertebrates

BACKGROUND

Genes that have been subject to adaptive evolution can produce varying degrees of pathology or differing symptomatology. ErbB family receptor activation will initiate a number of downstream signaling pathways, such as mitogen-activated protein kinase (MAPK), activator of transcription (STAT), the modulation of calcium channels, and so on, all of which lead to aggressive tumor behavior. However, the evolutionary mechanisms operating in the retention of ErbB family genes and the changes in selection pressures are not clear.

RESULTS

Sixty-two full-length cDNA sequences from 27 vertebrate species were extracted from the UniProt protein database, NCBI's GenBank and the Ensembl database. The result of phylogenetic analysis showed that the four ErbB family members in vertebrates might be formed by gene duplication. In order to determine the mode of evolution in vertebrates, selection analysis and functional divergence analysis were combined to explain the relationship of the site-specific evolution and functional divergence in the vertebrate ErbB family. Our results indicate that the acceleration of asymmetric evolutionary rates and purifying selection together were the main force for the production of ErbBs, and positive selections were detected in the ErbB family.

CONCLUSION

An evolutional phylogeny of 27 vertebrates was presented in our study; the tree showed that the genes have evolved through duplications followed by purifying selection, except for seven sites, which evolved by positive selection. There was one common site with positive selection and functional divergence. In the process of functional differentiation evolving through gene duplication, relaxed selection may play an important part.

Primate aging in the mammalian scheme: the puzzle of extreme variation in brain aging

At later ages, humans have high risk of developing Alzheimer disease (AD) which may afflict up to 50% by 90 years. While prosimians and monkeys show more substantial changes, the great apes brains examined show mild neurodegenerative changes. Compared with rodents, primates develop and reproduce slowly and are long lived. The New World primates contain some of the shortest as well as some of the longest-lived monkey species, while the prosimians develop the most rapidly and are the shortest lived. Great apes have the largest brains, slowest development, and longest lives among the primates. All primates share some level of slowly progressive, age-related neurodegenerative changes. However, no species besides humans has yet shown regular drastic neuron loss or cognitive decline approaching clinical grade AD. Several primates accumulate extensive deposits of diffuse amyloid-beta protein (Aβ) but only a prosimian-the gray mouse lemur-regularly develops a tauopathy approaching the neurofibrillary tangles of AD. Compared with monkeys, nonhuman great apes display even milder brain-aging changes, a deeply puzzling observation. The genetic basis for these major species differences in brain aging remains obscure but does not involve the Aβ coding sequence which is identical in nonhuman primates and humans. While chimpanzees merit more study, we note the value of smaller, shorter-lived species such as marmosets and small lemurs for aging studies. A continuing concern for all aging studies employing primates is that relative to laboratory rodents, primate husbandry is in a relatively primitive state, and better husbandry to control infections and obesity is needed for brain aging research.

Adaptive evolution of Toll-like receptor 5 in domesticated mammals

Background

Previous studies have proposed that mammalian toll like receptors (TLRs) have evolved under diversifying selection due to their role in pathogen detection. To determine if this is the case, we examined the extent of adaptive evolution in the TLR5 gene in both individual species and defined clades of the mammalia.

Results

In support of previous studies, we find evidence of adaptive evolution of mammalian TLR5. However, we also show that TLR5 genes of domestic livestock have a concentration of single nucleotide polymorphisms suggesting a specific signature of adaptation. Using codon models of evolution we have identified a concentration of rapidly evolving codons within the TLR5 extracellular domain a site of interaction between host and the bacterial surface protein flagellin.

Conclusions

The results suggest that interactions between pathogen and host may be driving adaptive change in TLR5 by competition between species. In support of this, we have identified single nucleotide polymorphisms (SNP) in sheep and cattle TLR5 genes that are co-localised and co-incident with the predicted adaptive codons suggesting that adaptation in this region of the TLR5 gene is on-going in domestic species.

Variation in CCL3L1 copy number in rhesus macaques (Macaca mulatta)

We used real-time quantitative PCR (qPCR) methodology to examine copy number variation (CNV) of the CCL3L1 gene among pure Indian-origin, pure Chinese-origin, and hybrid Indian-Chinese rhesus macaques (Macaca mulatta). CNV among purebred macaques fell within expected ranges, with Indian macaques having lower copy numbers than those of Chinese macaques. Compared with the purebred macaques, Indian-Chinese hybrid rhesus macaques showed much greater variance in copy number and an intermediate average copy number. Copy numbers of CCL3L1 in rhesus macaque trios (sire, dam, and offspring) were consistent with Mendelian inheritance.

Alzheimer disease susceptibility loci: evidence for a protein network under natural selection

Recent genome-wide association studies have identified a number of susceptibility loci for Alzheimer disease (AD). To understand the functional consequences and potential interactions of the associated loci, we explored large-scale data sets interrogating the human genome for evidence of positive natural selection. Our findings provide significant evidence for signatures of recent positive selection acting on several haplotypes carrying AD susceptibility alleles; interestingly, the genes found in these selected haplotypes can be assembled, independently, into a molecular complex via a protein-protein interaction (PPI) network approach. These results suggest a possible coevolution of genes encoding physically-interacting proteins that underlie AD susceptibility and are coexpressed in different tissues. In particular, PICALM, BIN1, CD2AP, and EPHA1 are interconnected through multiple interacting proteins and appear to have coordinated evidence of selection in the same human population, suggesting that they may be involved in the execution of a shared molecular function. This observation may be AD-specific, as the 12 loci associated with Parkinson disease do not demonstrate excess evidence of natural selection. The context for selection is probably unrelated to AD itself; it is likely that these genes interact in another context, such as in immune cells, where we observe cis-regulatory effects at several of the selected AD loci.

Lessons from chimpanzee-based research on human disease: the implications of genetic differences

Assertions that the use of chimpanzees to investigate human diseases is valid scientifically are frequently based on a reported 98-99% genetic similarity between the species. Critical analyses of the relevance of chimpanzee studies to human biology, however, indicate that this genetic similarity does not result in sufficient physiological similarity for the chimpanzee to constitute a good model for research, and furthermore, that chimpanzee data do not translate well to progress in clinical practice for humans. Leading examples include the minimal citations of chimpanzee research that is relevant to human medicine, the highly different pathology of HIV/AIDS and hepatitis C virus infection in the two species, the lack of correlation in the efficacy of vaccines and treatments between chimpanzees and humans, and the fact that chimpanzees are not useful for research on human cancer. The major molecular differences underlying these inter-species phenotypic disparities have been revealed by comparative genomics and molecular biology - there are key differences in all aspects of gene expression and protein function, from chromosome and chromatin structure to post-translational modification. The collective effects of these differences are striking, extensive and widespread, and they show that the superficial similarity between human and chimpanzee genetic sequences is of little consequence for biomedical research. The extrapolation of biomedical data from the chimpanzee to the human is therefore highly unreliable, and the use of the chimpanzee must be considered of little value, particularly given the breadth and potential of alternative methods of enquiry that are currently available to science.

Selection on the protein-coding genome

Populations evolve as mutations arise in individual organisms and, through hereditary transmission, may become "fixed" (shared by all individuals) in the population. Most mutations are lethal or have negative fitness consequences for the organism. Others have essentially no effect on organismal fitness and can become fixed through the neutral stochastic process known as random drift. However, mutations may also produce a selective advantage that boosts their chances of reaching fixation. Regions of genes where new mutations are beneficial, rather than neutral or deleterious, tend to evolve more rapidly due to positive selection. Genes involved in immunity and defense are a well-known example; rapid evolution in these genes presumably occurs because new mutations help organisms to prevail in evolutionary "arms races" with pathogens. In recent years, genome-wide scans for selection have enlarged our understanding of the evolution of the protein-coding regions of the various species. In this chapter, we focus on the methods to detect selection in protein-coding genes. In particular, we discuss probabilistic models and how they have changed with the advent of new genome-wide data now available.

Overestimation of nonsynonymous/synonymous rate ratio by reverse-translation of aligned amino acid sequences

In the analysis of protein-coding nucleotide sequences, the ratio of the number of nonsynonymous substitutions to that of synonymous substitutions (d(N)/d(S)) is used as an indicator for the direction and magnitude of natural selection operating at the amino acid sequence level. The d(S) and d(N) values are estimated based on the comparison of homologous codons, which are often identified by converting (reverse-translating) aligned amino acid sequences into codon sequences. In this method, however, homologous codons may be mis-identified when frame-shifts occurred or amino acid sequences were mis-aligned, which may lead to overestimation of the d(N)/d(S) ratio. Here the effect of reverse-translating aligned amino acid sequences on the estimation of d(N)/d(S) ratio was examined through a large-scale analysis of protein-coding nucleotide sequences from vertebrate species. Apparently, 1-9% of codon sites that were identified as homologous with reverse-translation contained non-homologous codons, where the d(N)/d(S) ratio was unduly high. By correcting the d(N)/d(S) ratio for these codon sites, it was inferred that the ratio was 5-43% overestimated with reverse-translation. These results suggest that caution should be exerted in the study of natural selection using the d(N)/d(S) ratio by reverse-translating aligned amino acid sequences.

[Does Alzheimer's disease exist in all primates? Alzheimer pathology in non-human primates and its pathophysiological implications (II)]

INTRODUCTION

In the ageing process there are some species of non-human primates which can show some of the defining characteristics of the Alzheimer's disease (AD) of man, both in neuropathological changes and cognitive-behavioural symptoms. The study of these species is of prime importance to understand AD and develop therapies to combat this neurodegenerative disease.

DEVELOPMENT

In this second part of the study, these AD features are discussed in the most important non-experimental AD models (Mouse Lemur -Microcebus murinus, Caribbean vervet -Chlorocebus aethiops, and the Rhesus and stump-tailed macaque -Macaca mulatta and M. arctoides) and experimental models (lesional, neurotoxic, pharmacological, immunological, etc.) non-human primates. In all these models cerebral amyloid neuropathology can occur in senility, although with different levels of incidence (100% in vervets;<30% in macaques). The differences between normal and pathological (Alzheimer's) senility in these species are difficult to establish due to the lack of cognitive-behavioural studies in the many groups analysed, as well as the controversy in the results of these studies when they were carried out. However, in some macaques, a correlation between a high degree of functional brain impairment and a large number of neuropathological changes ("possible AD") has been found.

CONCLUSIONS

In some non-human primates, such as the macaque, the existence of a possible continuum between "normal" ageing process, "normal" ageing with no deep neuropathological and cognitive-behavioural changes, and "pathological ageing" (or "Alzheimer type ageing"), may be considered. In other cases, such as the Caribbean vervet, neuropathological changes are constant and quite marked, but its impact on cognition and behaviour does not seem to be very important. This does assume the possible existence in the human senile physiological regression of a stable phase without dementia even if neuropathological changes appeared.

High-throughput single-nucleotide polymorphism discovery and the search for candidate genes for long-term SIVmac nonprogression in Chinese rhesus macaques (Macaca mulatta)

BACKGROUND

Genetic differences between Indian and Chinese rhesus macaques contribute to the phenotypic variance of clinical trials, including infection with SIVmac. The completion of the rhesus genome has facilitated the discovery of several thousand markers.

METHODS

We developed a genome-wide SNP map for rhesus macaques containing 3869 validated markers with an average distance of 0.88 Mb and used the program VarLD to identify genomic areas with significant differences in linkage disequilibrium (LD) between Indian-derived and Chinese rhesus macaques.

RESULTS

Forty-one statistically significant differences in LD between Chinese and Indian-origin rhesus were detected on chromosomes 1, 4, 5 and 11. The region of greatest LD difference was located on the proximal end of chromosome one, which also contained the genes ELAVL4, MAST2 and HIVEP3.

CONCLUSION

These genomic areas provide entry to more detailed studies of gene function. This method is also applicable to the study of differences in biomarkers between regional populations of other species.