Nonrandom Intragenic Variations in Patterns of Codon Bias Implicate a Sequential Interplay Between Transitional Genetic Drift and Functional Amino Acid Selection

Although most codon third bases appear to be functionless, the synonymous codons so defined exhibit a strikingly nonrandom distribution (codon bias) within human and other genes. To examine this phenomenon further, we generated a database of DNA sequences encoding human transmembrane cell-surface receptor proteins. Using this database we show here that the guanine and cytosine content of codon third bases (GC3) varies intragenically with the nature of the specified receptor domains (transmembrane > extracellular > intracellular domains; p < 0.001), the phenotype of the encoded amino acids (hydrophobic > hydrophilic > neutral amino acids; p < 0.001), and the receptor affiliation of the transmembrane domain superfamily (G-protein- coupled receptors > receptor tyrosine kinases; p < 0.001). Within gene regions specifying transmembrane domains, GC3 declines as domain functionality becomes redundant with increasing hydrophobicity (p < 0.001). Codons containing the second-base cytosine (XCZ, which encodes neutral amino acids) are selectively depleted of third-base adenine content (A3: XCA codons) when encoding transmembrane domain residues, consistent with positive selection for transitional mutation of XCG to XTG (which encodes hydrophobic amino acids) rather than to the synonymous XCA. Supporting this XCG XTG mechanism of codon bias, the G3:A3 ratio of codons specifying the transmembrane amino acid glycine (GGZ) is intermediate between that of its functional homolog alanine (GCZ) and that of hydrophobic valine (GTZ), even though the C3:T3 ratios are similar. Conversely, nearest-neighbor analysis of third bases 5 to codons specifying valine and leucine (CTZ) confirms a significant difference in C3:T3 but not G3:A3 ratios (i.e., C3/G1 T3/G1 > C3/A1; p < 0.001), consistent with the functionally advantageous retention of hydrophobic residues. These data raise the possibility that patterns of intragenic codon bias reflect a balance between negative and positive selection, suggesting in turn that analysis of codon third-base usage may help to predict the functional significance of encoded products. Supplementary information: Current address: (K. Lin) College of Life Sciences, Beijing Normal University, Beijing 100875, China     

How to infer reliable diploid genotypes from NGS or traditional sequence data: from basic probability to experimental optimization

The use of diploid sequence markers is still challenging despite the good quality of the information they provide. There is a common problem to all sequencing approaches [traditional cloning and sequencing of PCR amplicons as well as next-generation sequencing (NGS)]: when no variation is found within the sequences from a given individual, homozygozity can never be asserted with certainty. As a consequence, sequence data from diploid markers are mostly analysed at the population (not the individual level) particularly in animal studies. This study aims at contributing to solve this. Using the Bayes theorem and the binomial law, useful results are derived, among which: (i) the number of sequence reads per individual (or sequencing depth) which is required to ensure, at a given probability threshold, that some heterozygotes are not considered erroneously as homozygotes, as a function of the observed heterozygozity (H(o) ) of the locus in the population; (ii) a way of estimating H(o) from low coverage NGS data; (iii) a way of testing the null hypothesis that a genetic marker corresponds to a single and diploid locus, in the absence of data from controlled crosses; (iv) strategies for characterizing sequence genotypes in populations minimizing the average number of sequence reads per individual; (v) a rationale to decide which are the variations that one needs to consider along the sequence, as a function of the sequencing depth affordable, the level of polymorphism desired and the risk of sequencing error. For traditional sequencing technology, optimal strategies appear surprisingly different from the usual empirical ones. The average number of sequence reads required to obtain 99% of fully determined genotypes never exceeds six, this value corresponding to the worst situation when H(o) equals 0.6. This threshold value of H(o) is strikingly stable when the tolerated proportion of nonfully resolved genotypes varies in a reasonable range. These results do not rely on the Hardy-Weinberg equilibrium assumption or on diallelism of nucleotidic sites.     

新一代生物防治及其在天敌瓢虫中的研究进展

生物防治是一种利用自然天敌控制有害生物的手段.当前生物防治的发展面临着使用成本高和效果不稳定的两大难题.这要求我们更精准地选育和释放天敌.随着测序技术的进步,学界提出了新一代生物防治(Next Generation Biological Control),旨在利用遗传学和基因组学来认识天敌昆虫种下水平的多样性和变异规律...     

Staphylococcus aureus SH1000 and 8325‐4: comparative genome sequences of key laboratory strains in staphylococcal research

Aims: To provide comparative genome sequence data for two related model strains of Staphylococcus aureus (SH1000 and 8325‐4) that are used extensively in laboratory research. Methods and Results: Comparative genome sequencing was used to identify genetic differences between Staph. aureus SH1000 and the fully genome‐sequenced ancestral strain, Staph. aureus NCTC 8325. PCR amplification and DNA sequencing were employed to determine which of the genetic polymorphisms identified were also present in Staph. aureus 8325‐4, a direct derivative of 8325 and the parent strain of SH1000. Aside from known genetic differences between these strains, Staph. aureus SH1000 harboured 15 single‐nucleotide polymorphisms compared with 8325 (of which 12 were also found in 8325‐4), and a 63‐bp deletion upstream of the spa gene not present in either 8325 or 8325‐4. Conclusions: Staphylococcus aureus SH1000 and 8325‐4 contain a number of genetic polymorphisms relative to the progenitor strain of the lineage (8325) and to each other. Significance and Impact of the Study: The comparative genome sequences of SH1000 and 8325‐4 presented here define the genotypes of two key strains in staphylococcal laboratory research and reveal genetic polymorphisms that may impact their phenotypic properties.     

From where did the Western honeybee (Apis mellifera) originate?

The native range of the honeybee Apis mellifera encompasses Europe, Africa, and the Middle East, whereas the nine other species of Apis are found exclusively in Asia. It is therefore commonly assumed that A. mellifera arose in Asia and expanded into Europe and Africa. However, other hypotheses for the origin of A. mellifera have also been proposed based on phylogenetic trees constructed from genetic markers. In particular, an analysis based on >1000 single-nucleotide polymorphism markers placed the root of the tree of A. mellifera subspecies among samples from Africa, suggestive of an out-of-Africa expansion. Here, we re-evaluate the evidence for this and other hypotheses by testing the robustness of the tree topology to different tree-building methods and by removing specimens with a potentially hybrid background. These analyses do not unequivocally place the root of the tree of A. mellifera subspecies within Africa, and are potentially consistent with a variety of hypotheses for honeybee evolution, including an expansion out of Asia. Our analyses also support high divergence between western and eastern European populations of A. mellifera, suggesting they are likely derived from two distinct colonization routes, although the sources of these expansions are still unclear.     

Do genetic drift and accumulation of deleterious mutations preclude adaptation? Empirical investigation using RADseq in a northern lacustrine fish

Understanding genomic signatures of divergent selection underlying long‐term adaptation in populations located in heterogeneous environments is a key goal in evolutionary biology. In this study, we investigated neutral, adaptive and deleterious genetic variation using 7,192 SNPs in 31 Lake Trout (Salvelinus namaycush) populations (n = 673) from Québec, Canada. Average genetic diversity was low, weakly shared among lakes, and positively correlated with lake size, indicating a major role for genetic drift subsequent to lake isolation. Putatively deleterious mutations were on average at lower frequencies than the other SNPs, and their abundance relative to the entire polymorphism in each population was positively correlated with inbreeding, suggesting that the effectiveness of purifying selection was negatively correlated with inbreeding, as predicted from theory. Despite evidence for pronounced genetic drift and inbreeding, several outlier loci were associated with temperature and found in or close to genes with biologically relevant functions notably related to heat stress and immune responses. Outcomes of gene–temperature associations were influenced by the inclusion of the most inbred populations, in which allele frequencies deviated the most from model predictions. This result illustrates challenge in identifying gene–environment associations in cases of high genetic drift and restricted gene flow and suggests limited adaptation in populations experiencing higher inbreeding. We discuss the relevance of these findings for the conservation and management, notably regarding stocking and genetic rescue, of Lake Trout populations and other species inhabiting highly fragmented habitats.