Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “Codon usage bias”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 271 records · Page 15Linked to original sources

Variable selection pressures across lineages in Trichodesmium and related cyanobacteria based on the heterocyst differentiation protein gene hetR.

Due to the irreversible inhibition of nitrogenase by O2, N2 fixation is incompatible with the oxygenic photosynthesis of cyanobacteria. These organisms have therefore evolved various strategies for growing diazotrophically. One group of N2-fixing cyanobacteria has specialized cells, heterocysts, which contain the nitrogenase, lack the oxygenic photosystem II, and are virtually anoxic inside as the result of respiratory activity and a thick glycolipid cell wall. The hetR gene encodes a serine protease which is thought to be involved in the regulation of heterocyst development and in DNA binding. Although hetR is also present in many non-heterocystous N2-fixing cyanobacteria, its function in these organisms is unknown. In this study, hetR sequences of the N2-fixing, non-heterocystous cyanobacterium Trichodesmium spp. and related genera were examined for signatures of selection. In parsimony- or distance-based hetR phylogenies, the filamentous non-heterocystous cyanobacteria Symploca sp. and Leptolyngbya sp. were closest to Trichodesmium sp. However, accommodating molecular attributes of hetR such as nucleotide frequencies and rate heterogeneity in phylogenetic analyses suggested that many other genera could not be excluded as sister taxa of Trichodesmium. Maximum likelihood analysis of the dN/dS ratio (omega) showed that-irrespective of the use of Symploca, Leptolyngbya, or more distant taxa as an outgroup-the lineage between an outgroup and Trichodesmium (omega1=0.02-0.05) and a lineage leading to Trichodesmium erythraeum (omega1=0.02) were under much stronger purifying selection than the other lineages in Trichodesmium (omega0=0.13-0.32). Although the results from the maximum likelihood analyses are most trustworthy because of codon usage bias in Trichodesmium, the results from a simpler tree-based McDonald-Kreitman test were in general agreement. Due to their quite different assumptions, the combination of these two methods of analysis circumvents multiple testing which, in general, is problematic when using branch models. Although the causal selective forces underlying the substitution patterns in hetR have not yet been identified, these findings parallel the variety of physiological, molecular, and behavioral differences in cyanobacteria related to N2 fixation. The heterogeneity of selection pressures in Trichodesmium is more surprising, because multiple adaptation mechanisms have not been described in this genus.

Bacterial Proteins↗

Expression and purification of a trivalent pertussis toxin-diphtheria toxin-tetanus toxin fusion protein in Escherichia coli.

Pertussis toxoid, diphtheria toxoid, and tetanus toxoid are key components of diphtheria-tetanus-acellular pertussis vaccines. The efficacy of the vaccines is well documented, however, the vaccines are expensive partly because the antigens are derived from three different bacteria. In this study, a fusion protein (PDT) composed of the immunoprotective S1 fragment of pertussis toxin, the full-length non-toxic diphtheria toxin, and fragment C of tetanus toxin was constructed via genetic means. The correct fusion was verified by restriction endonuclease analysis and Western immunoblotting. Escherichia coli carrying the recombinant plasmid (pCoPDT) produced a 161kDa protein that was recognized by antibodies specific to the three toxins. The expression of the PDT protein was inducible by isopropyl-beta-d-thio-galactoside but the total amount of protein produced was relatively low. Attempts to improve the protein yield by expression in an E. coli strain (Rosetta-gami 2) that could alleviate rare-codon usage bias and by supplementation of the growth media with amino acids deemed to be a limiting factor in translation were not successful. The PDT protein remained in the insoluble fraction when the recombinant E. coli was grown at 37 degrees C but the protein became soluble when the bacteria were grown at 22 degrees C. The PDT protein was isolated via affinity chromatography on a NiCAM column. The protein was associated with five other proteins via disulfide bonds and non-covalent interactions. Following treatment with beta-mercaptoethanol, the PDT fusion was purified to homogeneity by preparative polyacrylamide gel electrophoresis with a yield of 45 microg/L of culture. Antisera generated against the purified PDT protein recognized the native toxins indicating that some, if not all, of the native epitopes were conserved.

Animals↗

Coalescent processes and relaxation of selective constraints leading to contrasting genetic diversity at paralogs AtHVA22d and AtHVA22e in Arabidopsis thaliana.

Duplicate loci offer a very powerful system for understanding the complicated genome structure and adaptive evolution of a gene family. In this study, the genetic variation at paralogs AtHVA22d and AtHVA22e, members of an ABA- and stress-inducible gene family, is examined in the selfing Arabidopsis thaliana. Population genetic analysis indicates contrasting levels of nucleotide diversity at overall exon sequence and nonsynonymous sites between AtHVA22d (pi = 0.00337, pi(rep) = 0.00158) and AtHVA22e (pi = 0.00054, pi(rep) = 0.00023). The fact of Ka/Ks ratios significantly less than 1 in all sequences indicates that both genes are functional and subjected to purifying selection. In addition, rooted at barley HVA22, accelerated evolution is detected at replacement changes in the AtHVA22d locus, indicating relaxation of purifying selection after gene duplication. However, relative rate tests reveal no deviation from the neutrality at synonymous sites between the two paralogs. Based on clock-like evolution, the rate of synonymous substitution is estimated at 1.83 x 10(-9) substitutions per site per year; and the divergence of the two paralogs is traced to 90 MYA, coinciding with a period of the diversification of angiosperms. Given no codon usage bias in both genes, natural selection alone cannot account for the 6.4-fold differences in the nucleotide variation at synonymous sites between the two paralogs. Random processes resulting in different coalescence times, 3.65 MYA at AtHVA22d vs. 1.20 MYA at AtHVA22e, may have predominantly contributed to the evident differences of the genetic diversity. Partially nonoverlapping modes of expression between the two functional paralogs suggest a subfunctionalization hypothesis for explaining the fates of duplicate loci.

Arabidopsis↗

Sulphur islands in the Escherichia coli genome: markers of the cell's architecture?

Two highly contrasted images depict genomes: at first sight, genes appear to be distributed randomly along the chromosome. In contrast, their organisation into operons (or pathogenicity islands) suggests that, at least locally, related functions are in physical proximity. Analysis of the codon usage bias in orthologous genes in the genome of bacteria which diverged a long time ago suggested that some physical (architectural) selection pressure organised the distribution of genes along the chromosome. The metabolism of highly reactive species such as sulphur-containing molecules must be compartmentalised to escape the deleterious actions of diffusible reagents such as gases or radicals. We analysed the distribution of sulphur metabolism genes in the genome of Escherichia coli and found a number of them to be clustered into statistically significant islands. Another interesting feature of these genes is that the proteins they encode are significantly deprived of cysteine and methionine residues, as compared to the bulk proteins. We speculate that this clustering is associated to the organisation of sulphur metabolism proteins into islands where the sensitive sulphur-containing molecules are protected from reacting with elements in the environment such as dioxygen, nitric oxide or radicals.

Chromosomes, Bacterial↗

Nucleotide sequence of the Azospirillum brasilense Sp7 glutamine synthetase structural gene.

The complete nucleotide sequence of the glnA gene, encoding the glutamine synthetase subunit of Azospirillum brasilense Sp7, was established. This is the first Azospirillum gene sequenced. The gene encodes a 468 residue polypeptide of MW 51,917. The similarity coefficient (SAB) between the polypeptidic sequence of Azospirillum and Anabaena 7120, which is the only other glnA sequence available, is 58%. No significant homology with E. coli canonical and ntr promoters, or with the promoter region of the Anabaena glnA gene was found. When fused to an E. coli promoter, the gene could be translated in E. coli, despite a very biased codon usage and an atypical Shine-Dalgarno sequence.

Amino Acid Sequence↗

Glyceraldehyde-3-phosphate dehydrogenase from Tetrahymena pyriformis: enzyme purification and characterization of a gapC gene with primitive eukaryotic features.

Glyceraldehyde-3-phosphate dehydrogenase (GAPDH, EC.1.2.1.12) was purified to electrophoretic homogeneity from an amicronucleated strain of the ciliate Tetrahymena pyriformis using a three-step procedure. The native enzyme is an homotetramer of 145 kDa exhibiting absolute specificity for NAD. In its catalytic properties it is similar to other glycolytic GAPDHs. Chromatofocusing analysis showed the presence of only one basic GAPDH isoform with an isoelectric point of 8.8. Western blots using a monospecific polyclonal antibody raised against the T. pyriformis GAPDH showed a single 36-kDa band corresponding to the enzyme subunit in the cytosolic protein fraction of this strain and the closely related species, both from the class Oligohymenophorea, Paramecium tetraurelia. No bands were immunodetected in the ciliate Colpoda inflata (class Colpodea) and in the diverse eukaryotes and eubacteria tested. A 0.5-kb DNA fragment which corresponds to an internal region of a gapC gene was generated by polymerase chain reaction using cDNA of T. pyriformis as template. This gene codes for a basic GAPDH protein with eukaryotic-diplomonad signatures and exhibits a codon usage biased in the manner typical for T. pyriformis genes. Southern blots performed both under homologous and heterologous conditions using this amplified cDNA fragment as a probe, indicated that it should be the only gapC gene present in the macronuclear genome of this ciliate, its expression being confirmed by Northern blot analysis. These results are discussed in connection with the peculiar genomic organization of ciliates and in the context of protist evolution.

Amino Acid Sequence↗

A test of translational selection at 'silent' sites in the human genome: base composition comparisons in alternatively spliced genes.

Natural selection appears to discriminate among synonymous codons to enhance translational efficiency in a wide range of prokaryotes and eukaryotes. Codon bias is strongly related to gene expression levels in these species. In addition, between-gene variation in silent DNA divergence is inversely correlated with codon bias. However, in mammals, between-gene comparisons are complicated by distinctive nucleotide-content bias (isochores) throughout the genome. In this study, we attempted to identify translational selection by analyzing the DNA sequences of alternatively spliced genes in humans and in Drosophila melanogaster. Among codons in an alternatively spliced gene, those in constitutively expressed exons are translated more often than those in alternatively spliced exons. Thus, translational selection should act more strongly to bias codon usage and reduce silent divergence in constitutive than in alternative exons. By controlling for regional forces affecting base-composition evolution, this within-gene comparison makes it possible to detect codon selection at synonymous sites in mammals. We found that GC-ending codons are more abundant in constitutive than alternatively spliced exons in both Drosophila and humans. Contrary to our expectation, however, silent DNA divergence between mammalian species is higher in constitutive than in alternative exons.

Alternative Splicing↗

Covariation of GC content and the silent site substitution rate in rodents: implications for methodology and for the evolution of isochores.

Many attempts to test selectionist and neutralist models employ estimates of synonymous (Ks) and non-synonymous (Ka) substitution rates of orthologous genes. For example, a stronger Ka-Ks correlation than expected under neutrality has been argued to indicate a role for selection and the absence of a Ks-GC4 correlation has been argued to be inconsistent with neutral models for isochore evolution. However, both of these results, we have shown previously, are sensitive to the method by which Ka and Ks are estimated. Using a maximum likelihood (ML) estimator (GY94) we found a positive correlation between Ks and GC4 and only a weak correlation between Ka and Ks, lower than expected under neutral expectations. This ML method is computationally slow. Recently, a new ad hoc approximation of this ML method has been provided (YN00). This is effectively an extension of Li's protocol but that also allows for codon usage bias. This method is computationally near-instantaneous and therefore potentially of great utility for analysis of large datasets. Here we ask whether this method might have such applicability. To this end we ask whether it too recovers the two unusual results. We report that when the ML and earlier ad hoc methods disagree, YN00 recovers the results described by the ML methods, i.e. a positive correlation between GC4 and Ks and only a weak correlation between Ks and Ka. If the ML method can be trusted, then YN00 can also be considered an adequately reliable method for analysis of large datasets. Assuming this to be so we also analyze further the patterns. We show, for example, that the positive correlation between GC4 and Ks is probably in part a mutational bias, there being more methyl induced CpG-->TpG mutations in GC rich regions. As regards the evolution of isochores, it seems inappropriate to use the claimed lack of a correlation between GC and Ks as definitive evidence either against or for any model. If the positive correlation is real then, we argue, this is hard to reconcile with the biased gene conversion model for isochore formation as this predicts a negative correlation.

Animals↗

An algorithm for detecting directional and non-directional positive selection, neutrality and negative selection in protein coding DNA sequences.

Positive selection or adaptive evolution is thought to be responsible, at least some of the time, for the rapid accumulation of advantageous changes in protein-coding genes. The origin of new enzymatic functions, erection of barriers to heterospecific fertilization, and evasion of host response by pathogens, among other things, are thought to be instances of adaptive evolution. Detecting positive selection in protein-coding genes is fraught with difficulties. Saturation for sequence change, codon usage bias, ephemeral selection events and differential selective pressures on amino acids all contribute to the problem. A number of solutions have been proposed with varying degrees of success, however they suffer from limitations of not being accurate enough or being prohibitively computationally intensive. We have developed a character-based method of identifying lineages that undergo positive selection. In our method we assess the possibility that for each internal branch of a phylogenetic tree an event occurred that subsequently gave rise to a greater number of replacement substitutions than might be expected. We classify these replacement substitutions into two categories - whether they subsequently became invariable or changed again in at least one descendent lineage. The former situation indicates that the new character state is under strong selection to preserve its new identity (directional selection), while the latter situation indicates that there is a persistent pressure to change identity (non-directional selection). The method is fast and accurate, easy to implement, sensitive to short-lived selection events and robust with respect to sampling density and proportion of sites under the influence of positive selection.

Algorithms↗

Cloning and heterologous expression of Entamoeba histolytica adenylate kinase and uridylate/cytidylate kinase.

We have isolated two cDNA clones encoding Entamoeba histolytica nucleotide kinases, EhAK and EhUK, expressed them in E. coli and performed functional studies of the recombinant enzymes. Nucleotide sequence analysis showed that EhAK and EhUK genes exhibited the features characteristic of E. histolytica genes, such as transcripts with relatively short 5' and 3' untranslated flanking regions containing the conserved E. histolytica transcription promoter elements located 5' to the initiation codon and a polyadenylation signal in the 3' UTR, a distinctive codon usage bias for A or T in the third position and an AT bias greater than 75% in the flanking regions of the transcripts. At the protein level, both enzymes belong to the short variant nucleoside monophosphate (NMP) kinases, which lack a 29amino acid LID region present in the long variant isoenzymes. EhAK was 30-38% identical to the members of the adenylate kinase (AK) family while EhUK was more similar (48-49% identity) to UMP/CMP kinases. Both enzymes used ATP as preferred phosphate-group donor but each one exhibited strict specificity for the acceptor NMP, EhAK for AMP and EhUK for the pyrimidine nucleoside monophosphates UMP and CMP. Biochemical characterization of the enzymes and phylogenetic reconstruction showed that EhUK is an authentic and well conserved member of the UMP/CMP kinase group while EhAK is the most divergent member known of the AK1 isoenzymes.

Adenylate Kinase↗

A compact gene cluster in Drosophila: the unrelated Cs gene is compressed between duplicated amd and Ddc.

Cs, a gene with unknown function, and amd and Ddc, which encode decarboxylases, are among the most closely spaced genes in D. melanogaster. Untranslated 3' ends of the convergently transcribed genes Cs and Ddc are known to overlap by 88bp. A number of questions arise about the organization of this tightly-packed gene region and about the evolution and function of the Cs gene. We have now investigated this three-gene cluster in Scaptodrosophila lebanonensis (which diverged from D. melanogaster 60-65 MYA), as well as in D. melanogaster and D. simulans. Gene order and direction of transcription is the same in all three species. The Cs gene codes, in Scaptodrosophila, for a polypeptide of 544 amino acids; in D. melanogaster, it consists of 504 amino acids, which is twice as long as previously suggested, which makes the gene density even more spectacular. The Cs sequences exhibit higher number of non-synonymous substitutions between species, higher ratios of non-synonymous to synonymous substitutions, and lower codon usage bias than other genes, suggesting that Cs is less functionally constrained than the other genes. This is consistent with the failure of inducing phenotypic mutations in D. melanogaster. The function of Cs remains to be identified, but a high degree of similarity indicates that it is homologous to genes coding for a corticosteroid-binding protein in yeast and a polyamine oxidase in maize.

Amino Acid Sequence↗

Neutral effect of recombination on base composition in Drosophila.

Recombination is thought to have various evolutionary effects on genome evolution. In this study, we investigated the relationship between the base composition and recombination rate in the Drosophila melanogaster genome. Because of a current debate about the accuracy of the estimates of recombination rate in Drosophila, we used eight different measures of recombination rate from recent work. We confirmed that the G + C content of large introns and flanking regions is positively correlated with recombination rate, suggesting that recombination has a neutral effect on base composition in Drosophila. We also confirmed that this neutral effect of recombination is the main determinant of the correlation between synonymous codon usage bias and recombination rate in Drosophila.

Animals↗

Single chain antibody (SCA) encoding genes: one-step construction and expression in eukaryotic cells.

We report the expression, in eukaryotic cells, of a gene encoding a single chain antibody (SCA) and a rapid method for the construction of such genes. A SCA directed against the aromatic dye fluorescein was synthesized from a gene constructed by means of the simultaneous use of four PCR primers and templates of both light and heavy chain immunoglobulin cDNAs in the form of either plasmid clones or reverse transcribed hybridoma RNA. Two of the primers were partially complementary to one another and encoded the polypeptide linker which joins the immunoglobulin light and heavy chain variable domains of the SCA polypeptide. A functional, hapten-binding product was synthesized from the gene thus constructed in both E. coli and the fission yeast, Schizosaccharomyces pombe. Our results demonstrate that gene constructs encoding single chain antigen binding proteins can be synthesized very rapidly with only limited sequence information about the pertinent light and heavy chain immunoglobulin genes, and, that neither murine codon usage bias, Thermus aquaticus DNA polymerase infidelity, nor the eukaryotic cellular environment preclude the synthesis of functional single chain antigen binding proteins in non-lymphatic, non-murine eukaryotic cells.

ATP-Binding Cassette Transporters↗

Elevated evolutionary rates in the laboratory strain of Saccharomyces cerevisiae.

By using the maximum likelihood method, we made a genome-wide comparison of the evolutionary rates in the lineages leading to the laboratory strain (S288c) and a wild strain (YJM789) of Saccharomyces cerevisiae and found that genes in the laboratory strain tend to evolve faster than in the wild strain. The pattern of elevated evolution suggests that relaxation of selection intensity is the dominant underlying reason, which is consistent with recurrent bottlenecks in the S. cerevisiae laboratory strain population. Supporting this conclusion are the following observations: (i) the increases in nonsynonymous evolutionary rate occur for genes in all functional categories; (ii) most of the synonymous evolutionary rate increases in S288c occur in genes with strong codon usage bias; (iii) genes under stronger negative selection have a larger increase in nonsynonymous evolutionary rate; and (iv) more genes with adaptive evolution were detected in the laboratory strain, but they do not account for the majority of the increased evolution. The present discoveries suggest that experimental and possible industrial manipulations of the laboratory strain of yeast could have had a strong effect on the genetic makeup of this model organism. Furthermore, they imply an evolution of laboratory model organisms away from their wild counterparts, questioning the relevancy of the models especially when extensive laboratory cultivation has occurred. In addition, these results shed light on the evolution of livestock and crop species that have been under human domestication for years.

Biological Evolution↗

Roles of selection and recombination in the evolution of type I restriction-modification systems in enterobacteria.

Restriction-modification systems can protect bacteria against viral infection. Sequences of the hsdM gene, encoding one of the three subunits of type I restriction-modification systems, have been determined for four strains of enterobacteria. Comparison with the known sequences of EcoK and EcoR124 indicates that all are homologous, though they fall into three families (exemplified by EcoK, EcoA, and EcoR124), the first two of which are apparently allelic. The extent of amino acid sequence identity between EcoK and EcoA is so low that the genes encoding them might be better termed pseudoalleles; this almost certainly reflects genetic exchange among highly divergent species. Within the EcoK family the ratio of intra- to interspecific divergence is very high. The extent of divergence between the genes from Escherichia coli K-12 and Salmonella typhimurium LT2 is similar to that for other genes with the same level of codon usage bias. In contrast, intraspecific divergence (between E. coli strains B and K-12) is extremely high and may reflect the action of frequency-dependent selection mediated by bacteriophages. There is also evidence of lateral transfer of a short sequence between E. coli and S. typhimurium.

Amino Acid Sequence↗

Rapid identification of yeast proteins on two-dimensional gels.

This work describes a rapid and sensitive technique for the identification of Saccharomyces cerevisiae proteins on two-dimensional gels based on the determination of their amino acid ratios. Specific double labeling with 3H and 14C or 35S-labeled amino acids, chosen among those that are specifically incorporated into proteins without interconversion, allowed an accurate measurement of different amino acid ratios for 200 proteins. A computer program was developed to screen a yeast data base containing 1700 protein sequences and to identify proteins matching the measured Mr, pI, and amino acid ratios. The method, tested with 45 reference proteins, allowed 79 new identifications corresponding to abundant proteins belonging to a few functional families. Some protein spots correspond to homologs of mammalian proteins or to uncharacterized open reading frames. Remarkably, among identified proteins of similar abundance, the organellar proteins have a markedly lower codon usage bias than the cytosolic ones. The double labeling technique is particularly suited to the analysis, on a single two-dimensional gel, of the influence of physiological or genetic changes on yeast protein content.

Amino Acids↗

Disulfide bond assignment in human interleukin-7 by matrix-assisted laser desorption/ionization mass spectroscopy and site-directed cysteine to serine mutational analysis.

Interleukin-7 (IL-7) is a proteinaceous biological response modifier that has a bioactive tertiary structure dependent on disulfide bond formation. Disulfide bond assignments in human (h)IL-7 are based upon the results of matrix-assisted laser desorption/ionization (MALDI) mass spectroscopy and Cys to Ser mutational analyses. A gene encoding the hIL-7 was synthesized incorporating Escherichia coli codon usage bias and was used to express biologically active protein as determined by stimulation of precursor B-cell proliferation. MALDI mass spectroscopic analysis of trypsin-digested hIL-7 was performed and compared with the anticipated results of a simulated tryptic digestion. Many of the anticipated hIL-7 tryptic fragments were detected including one with a molecular mass equivalent to the sum of two polypeptides linked through a disulfide bond formed from Cys residues (Cys3 and Cys142). Subsequently, Cys to Ser substitution mutational analyses were performed. A hIL-7 variant with all six Cys substituted with Ser was found to be biologically inactive (EC50 > 1 x 10(-7) M). In contrast, a family of single disulfide bond-forming variants of hIL-7 were constructed by reintroducing Cys pairs (Cys3-Cys142, Cys35-Cys130, and Cys48-Cys93), and each could stimulate cell proliferation with an EC50 of 4 x 10(-9), 2 x 10(-8), and 2 x 10(-9) M, respectively. In single disulfide bond-forming mutants of hIL-7, the ability to stimulate cell proliferation was abolished in the presence of 2 mM dithiothreitol. The results presented strongly suggest that only a single disulfide bond is required for hIL-7 to form a tertiary structure capable of stimulating precursor B-cell proliferation.

Amino Acid Sequence↗

The Caenorhabditis elegans unc-93 gene encodes a putative transmembrane protein that regulates muscle contraction.

unc-93 is one of a set of five interacting genes involved in the regulation or coordination of muscle contraction in Caenorhabditis elegans. Rare altered-function alleles of unc-93 result in sluggish movement and a characteristic "rubber band" uncoordinated phenotype. By contrast, null alleles cause no visibly abnormal phenotype, presumably as a consequence of the functional redundancy of unc-93. To understand better the role of unc-93 in regulating muscle contraction, we have cloned and molecularly characterized this gene. We isolated transposon-insertion alleles and used them to identify the region of DNA encoding the unc-93 protein. Two unc-93 proteins differing at their NH2 termini are potentially encoded by transcripts that differ at their 5' ends. The putative unc-93 proteins are 700 and 705 amino acids in length and have two distinct regions: the NH2 terminal portion of 240 or 245 amino acids is extremely hydrophilic, whereas the rest of the protein has multiple potential membrane-spanning domains. The unc-93 transcripts are low in abundance and the unc-93 gene displays weak codon usage bias, suggesting that the unc-93 protein is relatively rare. The unc-93 protein has no sequence similarity to proteins listed in current data-bases. Thus, unc-93 is likely to encode a novel membrane-associated muscle protein. We discuss possible roles for the unc-93 protein either as a component of an ion transport system involved in excitation-contraction coupling in muscle or in coordinating muscle contraction between muscle cells by affecting the functioning of gap junctions.

Alleles↗