Tolerance of nonsynonymous variation is closely correlated between human and mouse orthologues

ABSTRACTGenic constraint describes how tolerant a gene is of nonsynonymous variation before it is removed from the population by negative selection. Here, we provide the first estimates of intraspecific constraint for mouse genes genome-wide, and show constraint is positively correlated between human and mouse orthologues (r = 0.806). We assess the relationships between mouse gene constraint and knockout phenotypes, showing gene constraint is positively associated with pleiotropy (ie an increased number of phenotype associations (R2 = 0.65)), in addition to an enrichment in lethal, developmental, and craniofacial knockout phenotypes amongst the most constrained genes. Finally, we show mouse constraint can be used to predict human genes associated with Mendelian disease, and is positively correlated with an increase in the number of known pathogenic variants in the human orthologue (R2 = 0.23). Our metrics of mouse and human constraint are available to inform future research using mouse models.

Resequencing DCDC5 in the Flanking Region of an LD-SNP Derived from a Kidney-Yang Deficiency Syndrome Family

Objective. To explore the genetic traits of Kidney-yang deficiency syndrome (KDS).Design. Twelve KDS subjects and three spouses from a typical KDS family were recruited. Their genomic DNA samples were genotyped by Affymetrix 100K single-nucleotide polymorphism (SNP) arrays. The linkage disequilibrium (LD) SNPs were generated using GeneChip DNA analysis software (GDAS, Affymetrix). Genes located within 100 bp of the flanks of LD SNPs were mined via GeneView. 29 exons of the doublecortin domain containing 5 (DCDC5), a representative gene within the flank of an LD SNP, were resequenced.Results. Five LD SNPs display midrange linkage with KDS. Two genes with established functions, DCDC5 and Leucyl-tRNA synthetase, were mined in the flanks of LD SNPs. Resequencing of DCDC5 revealed a nonsynonymous variation, in which 3764T/A was replaced by C/G. Accordingly, the Ser1172was substituted by Pro1172. The S1172P substitution effect was evaluated as “possibly damaging” by PolyPhen.Conclusion. We have identified a genomic variation of DCDC5 based on the LD SNPs derived from a KDS family. DCDC5 and other genes surrounding these SNPs display some relationships with key symptoms of KDS.

Patterns of Intra- and Interhost Nonsynonymous Variation Reveal Strong Purifying Selection in Dengue Virus

ABSTRACT Considerable uncertainty surrounds the evolutionary rates of and selection pressures acting on arthropod-borne RNA viruses (arboviruses). In particular, it is unclear why arboviruses such as dengue virus show substantial genetic variation within individual humans and mosquitoes yet low long-term rates of amino acid substitution. To address this question, I compared patterns of nonsynonymous variation in populations of dengue virus sampled at different levels of evolutionary divergence. Although nonsynonymous variation was abundant in viral populations within individual humans, there was a marked reduction in the frequency of nonsynonymous mutations in interhost comparisons. Moreover, intrahost genetic variation corresponded to a random pattern of mutation, and most of the sites that exhibited nonsynonymous variation within hosts were invariant at deeper phylogenetic levels. This loss of long-term nonsynonymous variation is the signature of extensive purifying selection such that more than 90% of all nonsynonymous mutations are deleterious. Consequently, although arboviruses are able to successfully adapt to diverse cell types, they are characterized by a high rate of deleterious mutation.

Different Forces Drive the Evolution of the Acp26Aa and Acp26Ab Accessory Gland Genes in the Drosophila melanogaster Species Complex

The Acp26Aa and Acp26Ab genes that code for male accessory gland proteins are tandemly arranged in the species of the Drosophila melanogaster complex. An ∼1.6-kb region encompassing both genes has been sequenced in 10, 24, and 18 lines from Spain, Ivory Coast, and Malawi, respectively; the previously studied 10 lines from North Carolina have also been included in the analyses. A total of 110 nucleotide and 4 length polymorphisms were detected. Silent variation for the whole Acp26A region was slightly higher in African than in non-African populations, while for both genes nonsynonymous variation was similar in all populations studied. Based on Fst estimates no major genetic differentiation was detected between East and West Africa, while in general non-African populations were strongly differentiated from both African populations. Comparison of polymorphism and divergence at synonymous and nonsynonymous sites revealed that directional selection acting on amino acid replacement changes has driven the evolution of the Acp26Aa protein in the last 2.5 myr.

Shortening of the Symptom-Free Period in Rhesus Macaques Is Associated with Decreasing Nonsynonymous Variation in theenv Variable Regions of Simian Immunodeficiency Virus SIVsm during Passage

ABSTRACT During six blood passages of simian immunodeficiency virus SIVsm in rhesus macaques, the asymptomatic period shortened from 18 months to 1 month. To study SIVsm envelope gene (env) evolution during passage in rhesus macaques, the C1 to CD4 binding regions of multiple clones were sequenced at seroconversion and again at death. Theenv variation found during adaptation was almost completely confined to the variable regions. Intrasample sequence variation among clones at seroconversion was lower than the variation among clones at death. Intrasample variation among clones from a single time point as well as intersample variation decreased during the passage. In the variable regions, the mean number of intrasample nonsynonymous nucleotide substitutions decreased from the first passage (5.26 × 10−2 ± 0.6 × 10−2 per site) to the fifth passage (2.24 × 10−2 ± 0.4 × 10−2 per site), whereas in the constant regions, the mean number of intrasample nonsynonymous nucleotide substitutions differed less between the first and fifth passages (1.14 × 10−2 ± 0.27 × 10−2 and 0.80 × 10−2 ± 0.24 × 10−2 per site). Shortening of the asymptomatic period coincided with a rise in theKs/Ka ratio (ratio between the number of synonymous [Ks] and the number of nonsynonymous [Ka] substitutions) from 1.080 in passage one to 1.428 in passage five and mimicked the difference seen in the intrahost evolution between asymptomatic and fast-progressing individuals infected with human immunodeficiency virus type 1. The distribution of nonsynonymous substitutions was biphasic, with most of the adaptation ofenv variable regions occurring in the first three passages. This phase, in which the symptom-free period fell to 4 months, was followed by a plateau phase of apparently reduced adaptation. Analysis of codon usage revealed decreased codon redundancy in the variable regions. Overall, the results suggested a biphasic pattern of adaptation and evolution, with extremely rapid selection in the first three passages followed by an equilibrium or stabilization of the variation between env clones at different time points in passages four to six.