Differentiation of Venoms of Predatory Marine Gastropods: Divergence of Orthologous Toxin Genes of Closely Related Conus Species with Different Dietary Specializations

Venoms of Conus are remarkably diverse among species and the genes that encode conotoxins show high rates of evolution. Yet no prior studies have specifically explored how conotoxin gene evolution contributes to the differentiation of venoms of closely related Conus species. Previous investigations of four-loop conotoxin expression patterns of six closely related Conus species identified 12 sets of putative orthologous loci from these species, including eight pairs of loci that are coexpressed by two of these six species, C. abbreviatus and C. miliaris. Here I analyze the molecular evolution of orthologous conotoxin loci of these species and specifically examine the divergence of the eight orthologous counterparts of C. abbreviatus and C. miliaris. Tree and maximum likelihood-based analyses of these sequences reveal that positive selection promotes the divergence of orthologous genes among species and that the evolution of orthologues of C. abbreviatus and C. miliaris is asymmetric among species. The asymmetric evolution of conotoxin loci among species may result from lineage-specific dietary shifts or interspecific differences in the impact of selection from predator-prey interactions on conotoxin loci.     

Extensive and continuous duplication facilitates rapid evolution and diversification of gene families

The origin of novel gene functions through gene duplication, mutation, and natural selection represents one of the mechanisms by which organisms diversify and one of the possible paths leading to adaptation. Nonetheless, the extent, role, and consequences of duplications in the origins of ecological adaptations, especially in the context of species interactions, remain unclear. To explore the evolution of a gene family that is likely linked to species associations, we investigated the evolutionary history of the A-superfamily of conotoxin genes of predatory marine cone snails (Conus species). Members of this gene family are expressed in the venoms of Conus species and are presumably involved in predator-prey associations because of their utility in prey capture. We recovered sequences of this gene family from genomic DNA of four closely related species of Conus and reconstructed the evolutionary history of these genes. Our study is the first to directly recover conotoxin genes from Conus genomes to investigate the evolution of conotoxin gene families. Our results revealed a phenomenon of rapid and continuous gene turnover that is coupled with heightened rates of evolution. This continuous duplication pattern has not been observed previously, and the rate of gene turnover is at least two times higher than estimates from other multigene families. Conotoxin genes are among the most rapidly evolving protein-coding genes in metazoans, a phenomenon that may be facilitated by extensive gene duplications and have driven changes in conotoxin functions through neofunctionalization. Together these mechanisms led to dramatically divergent arrangements of A-superfamily conotoxin genes among closely related species of Conus. Our findings suggest that extensive and continuous gene duplication facilitates rapid evolution and drastic divergence in venom compositions among species, processes that may be associated with evolutionary responses to predator-prey interactions.     

Rapid divergence of two classes of Haemophilus ducreyi

Haemophilus ducreyi, the etiologic agent of chancroid, expresses variants of several key virulence factors. While previous reports suggested that H. ducreyi strains formed two clonal populations, the differences between, and diversity within, these populations were unclear. To assess their variability, we examined sequence diversity at 11 H. ducreyi loci, including virulence and housekeeping genes, augmenting published data sets with PCR-amplified genes to acquire data for at least 10 strains at each locus. While sequences from all 11 loci place strains into two distinct groups, there was very little variation within each group. The difference between alleles of the two groups was variable and large at 3 loci encoding surface-exposed proteins (0.4 < K(S) < 1.3, where K(S) is divergence at synonymous sites) but consistently small at genes encoding cytoplasmic or periplasmic proteins (K(S) < 0.09). The data suggest that the two classes have recently diverged, that recombination has introduced variant alleles into at least 3 distinct loci, and that these alleles have been confined to one of the two classes. In addition, recombination is evident among alleles within, but not between, classes. Rather than clones of the same species, these properties indicate that the two classes may form distinct species.     

Isolation and characterization of microsatellite loci in Schizolobium parahyba (Leguminosae)

Primer sequences for 16 microsatellite loci were developed from Schizolobium parahyba, a tropical tree species. Twelve loci were found to be polymorphic after screening diversity in individuals from Belize. A total of 39 alleles were found at nine loci. The markers are invaluable tools for studying the population genetics and mating system of the species.     

Molecular Evolution and Diversity of Conus Peptide Toxins,as Revealed by Gene Structure and Intron Sequence Analyses

Cone snails, which are predatory marine gastropods, produce a cocktail of venoms used for predation, defense and competition. The major venom component, conotoxin, has received significant attention because it is useful in neuroscience research, drug development and molecular diversity studies. In this study, we report the genomic characterization of nine conotoxin gene superfamilies from 18 Conus species and investigate the relationships among conotoxin gene structure, molecular evolution and diversity. The I1, I2, M, O2, O3, P, S, and T superfamily precursors all contain three exons and two introns, while A superfamily members contain two exons and one intron. The introns are conserved within a certain gene superfamily, and also conserved across different Conus species, but divergent among different superfamilies. The intronic sequences contain many simple repeat sequences and regulatory elements that may influence conotoxin gene expression. Furthermore, due to the unique gene structure of conotoxins, the base substitution rates and the number of positively selected sites vary greatly among exons. Many more point mutations and trinucleotide indels were observed in the mature peptide exon than in the other exons. In addition, the first example of alternative splicing in conotoxin genes was found. These results suggest that the diversity of conotoxin genes has been shaped by point mutations and indels, as well as rare gene recombination or alternative splicing events, and that the unique gene structures could have made a contribution to the evolution of conotoxin genes.     

Personal receptor repertoires: olfaction as a model

ABSTRACT: BACKGROUND: Information on nucleotide diversity along completely sequenced human genomes has increased tremendously over the last few years. This makes it possible to reassess the diversity status of distinct receptor proteins in different human individuals. To this end, we focused on the complete inventory of human olfactory receptor coding regions as a model for personal receptor repertoires. RESULTS: By performing data-mining from public and private sources we scored genetic variations in 413 intact OR loci, for which one or more individuals had an intact open reading frame. Using 1000 Genomes Project haplotypes, we identified a total of 4069 full-length polypeptide variants encoded by these OR loci, average of ~10 per locus, constituting a lower limit for the effective human OR repertoire. Each individual is found to harbor as many as 600 OR allelic variants, ~50% higher than the locus count. Because OR neuronal expression is allelically excluded, this has direct effect on smell perception diversity of the species. We further identified 244 OR segregating pseudogenes (SPGs), loci showing both intact and pseudogene forms in the population, twenty-six of which are annotatively "resurrected" from a pseudogene status in the reference genome. Using a custom SNP microarray we validated 150 SPGs in a cohort of 468 individuals, with every individual genome averaging 36 disrupted sequence variations, 15 in homozygote form. Finally, we generated a multi-source compendium of 63 OR loci harboring deletion Copy Number Variations (CNVs). Our combined data suggest that 271 of the 413 intact OR loci (66%) are affected by nonfunctional SNPs/indels and/or CNVs. CONCLUSIONS: These results portray a case of unusually high genetic diversity, and suggest that individual humans have a highly personalized inventory of functional olfactory receptors, a conclusion that might apply to other receptor multigene families.     

Length polymorphism and allele structure of trinucleotide microsatellites in natural accessions of Arabidopsis thaliana

 The objective of this work was to assess the degree of trinucleotide microsatellite length polymorphism in the selfing species Arabidopsis thaliana. PCR amplifications of 12 microsatellite loci among 49 natural populations revealed between one to eight length variants (alleles) for each locus. The average number of alleles per locus was four and the average genetic diversity index was 0.43. Divergence between length variants was investigated at the nucleotide level. Several observations emerge from the sequence data: (1) for most loci, length polymorphism results only from variations in the number of trinucleotide repeats; (2) for a few others, some variability was noted in the flanking sequences; (3) for compound and interrupted loci containing two arrays of trinucleotide repeats, length variations preferentially affect the longest one. Five of the Arabidopsis thaliana accessions were clearly composed of two sublines. In 2 other accessions, some heterozygous individual plants, probably resulting from recent outcrosses, were found. A phylogenetic tree constructed on the basis of trinucleotide microsatellite allelic diversity shows that genetic relationships among the accessions are not correlated with their geographic origin. Received: 4 November 1997 / Accepted: 3 March 1998     

Transition stages of molecular drive in multiple-copy DNA families in Drosophila

Multigene and non-genic DNA families are in a state of turnover and hence are continually being replaced throughout a population by new variant repeats. To quantify such molecular processes, in the absence of selection, it is necessary to find and compare stages of transistion during the homogenization of at least two non-genic families evolving in parallel in a closely related group of species. Detailed sequence analysis of patterns of variation, at each nucleotide position considered independently, amongst repeats of two tandem DNA families from seven related Drosophila species, reveals all stages of transition during the spread of randomly produced variant repeats. Variant repeats are found at different stages of homogenization and fixation in a population, irrespective of the loci, chromosomes or individuals from which they were cloned. Differences between the families in the relatively small number of variants at each transition stage and the greater number of fully homogenized and fixed variants between species of greater divergence indicate that the process of spread (molecular drive) is rapid relative to the mutation rate and occurs at seemingly different constant rates for each family. Occasional gene conversions, in addition to unequal exchanges, have contributed to family turnover. The significance of these results to the evolution of functional multigene families and divergence and conservation of sequences is discussed.     

Microsatellite loci for the North American tortoises (genus Gopherus) and their applicability to other turtle species

We describe the development of nine microsatellite loci from the gopher tortoise (Gopherus polyphemus). Screening 270 individuals, we found these loci to be highly variable (from two to 15 alleles per locus) and thus likely to be applicable in parentage and population‐level analyses for G. polyphemus. All nine loci amplified readily in the other three species of Gopherus as well. This observation, together with successful cross amplification of several loci in two additional families (Cheloniidae, Kinosternidae), underscores their potential utility among turtles in general.     

中华猕猴桃和美味猕猴桃自然居群遗传结构及其种间杂交渐渗

 利用9对SSR引物对中华猕猴桃（Actinidia chinensis）和美味猕猴桃（A. deliciosa）两近缘种的5个同域分布复合体和各自1个非同域分布居群进行了遗传多样性、居群遗传结构的分析以及种间杂交渐渗的探讨。结果表明：1）两物种共有等位基因比例高达81.13%,物种特有等位基因较少（中华猕猴桃：13.27%,美味猕猴桃：5.61%）,但共享等位基因表型频率在两近缘种间存在差异,而且与各同域复合体中两物种样本的交错程度或间距存在关联;2）两种猕猴桃均具有极高遗传多样性,美味猕猴桃的遗传多样性（H_o＝0 .749, PIC＝0.818）都略高于中华猕猴桃（H_o＝0.686,PIC＝0.799）;3）两猕 猴桃物种均具有较低的Nei’s居群遗传分化度,但AMOVA分析结果揭示种内异域居群间（F_ST=0.091 5）和同域复合体种间（F_ST=0.111 5）均存在一定程度的遗传分化;中华猕猴桃居群遗传分化（G_ST=0.086; F_ST=0.212 1）高于美味猕猴桃（G_ST= 0.080;F_ST=0.142 0）;4）同域分布复合体两物种间的遗传分化（G_ST=0.020）低于物种内异域居群间的遗传分化（中华猕猴桃：G_ST=0.086; 美味猕猴桃：G_ST=0.080）,同域复合体物种间的基因流（N_m＝7.89 -29.75）远远高于 同种异域居群间（中华猕猴桃：N_m =2.663; 美味猕猴桃：N_m=2.880）; 5）居群UPGMA聚类揭示在同一地域的居群优先聚类,个体聚类结果显示多数个体聚在各自居群组内,但各地理居群并不按地理距离的远近聚类,这与Mantel相关性检测所揭示的居群间遗传距离与地理距离没有显著性相关的结果一致。进一步分析表明两种猕猴桃的遗传多样性和居群遗传结构不仅受其广域分布、远交、晚期分化等生活史特性的影响,同时还与猕猴桃的染色体基数高 (x=29)、倍性复杂和种间杂交等因素密切相关,其中两种猕猴桃的共享祖先多态性和同域分布种间杂交基因渗透对两猕猴桃的居群遗传结构产生了重要影响。     

MHC class I of saltwater crocodiles (Crocodylus porosus): polymorphism and balancing selection

Saltwater crocodiles are in high demand for the production of luxury fashion items. However, their susceptibility to disease incurs substantial losses and it is hoped to be able to genetically select these animals for disease resistance. So far, this has only been enabled by phenotypic selection. Investigating the major histocompatibility complex (MHC) could provide insight into the ability of an individual to respond to pathogens acting as a selective pressure on the host. Here, we assessed genetic diversity and a role of selection in shaping the diversity of MHC class I exon 3 among 42 saltwater crocodiles from nine river basins in the Northern Territory, Australia. We generated 640 sequences using cloning and sequencing methods and identified 43 MHC variants among them. Phylogenetic analyses clustered these variants into two major clades, which may suggest two gene lineages. We found the number of variants within an individual varying between one and seven, indicating that there are at least four gene loci in this species. Selection detection analyses revealed an elevated ratio of nonsynonymous to synonymous substitutions (mean?=?1.152 per codon), suggesting balancing selection. Population differentiation analyses revealed that the MHC did not show structuring among the river basins, and there were some shared variants among them. This may be a result of possible gene flow and/or similar selection pressures among populations. These findings provide background knowledge to identify potential MHC markers, which could be used for selecting genetically variable individuals for future disease associations. All MHC class I exon 3 sequences reported in this paper were submitted to the GenBank database with following accession numbers: HQ008785–HQ008789, HQ008791–HQ008798, HQ008808–HQ008815, HQ008824, HQ008826–HQ008830, HQ008835, HQ008839, HQ008842–HQ008850, and JX023536–JX023540.     

A locus-wide approach to assessing variation in the avian MHC: the B-locus of the wild turkey

Studies of major histocompatibility complex (MHC) diversity in non-model vertebrates typically focus on structure and sequence variation in the antigen-presenting loci: the highly variable and polymorphic class I and class IIB genes. Although these studies provide estimates of the number of genes and alleles/locus, they often overlook variation in functionally related and co-inherited genes important in the immune response. This study utilizes the sequence of the MHC B-locus derived from a commercial turkey to investigate MHC variation in wild birds. Sequences were obtained for nine interspersed MHC amplicons (non-class I/II) from each of 40 birds representing 3 subspecies of wild turkey (Meleagris gallopavo). Analysis of aligned sequences identified 238 single-nucleotide variants approximately one-third of which had minor allele frequencies >0.2 in the sampled birds. PHASE analysis identified 70 prospective MHC haplotypes in the wild turkeys, whereas a combined analysis with commercial birds identified almost 100 haplotypes in the species. Denaturing gradient gel electrophoresis (DGGE) of the class IIB loci was used to test the efficacy of single-nucleotide polymorphism (SNP) haplotyping to capture locus-wide variation. Diversity in SNP haplotypes and haplotype sharing among individuals was directly reflected in the DGGE patterns. Utilization of a reference haplotype to sequence interspersed regions of the MHC has significant advantages over other methods of surveying diversity while identifying high-frequency SNPs for genotyping. SNP haplotyping provides a means to identify both divergent haplotypes and homozygous individuals for assessment of immunological variation in wild and domestic populations.     

Coevolution of immunoglobulin heavy- and light-chain variable-region gene families

The gene families encoding the immunoglobulin variable regions of heavy (VH) and light (VL) chains in vertebrates are composed of many genes. However, the gene number and the extent of diversity among VH and VL gene copies vary with species. To examine the causes of this variation and the evolutionary forces for these multigene families, we conducted a phylogenetic analysis of VH and VL genes from the species of amniotes. The results of our analysis showed that for each species, VH and VL genes have the same pattern of clustering in the trees, and, according to this clustering pattern, the species can be divided into two groups. In the first group of species (humans and mice), VH and VL genes were extensively intermingled with genes from other organisms; in the second group of species (chickens, rabbits, cattle, sheep, swine, and horses), the genes tended to form clusters within the same group of organisms. These results suggest that the VH and VL multigene families have evolved in the same fashion: they have undergone coordinated contraction and expansion of gene repertoires such that each group of organisms is characterized by a certain level of diversity of VH and VL genes. The extent of diversity among copies of VH and VL genes in each species is related to the mechanism of generation of antibody variety. In humans and mice, DNA rearrangement of immunoglobulin variable, diversity, and joining-segment genes is a main source of antibody diversity, whereas in chickens, rabbits, cattle, sheep, swine, and horses, somatic hypermutation and somatic gene conversion play important roles. The evolutionary pattern of VH and VL multigene families is consistent with the birth-and-death model of evolution, yet different levels of diversifying selection seem to operate in the VH and VL genes of these two groups of species.      

Length Polymorphisms in Human Proline-Rich Protein Genes Generated by Intragenic Unequal Crossing over

Southern blot hybridization analysis of genomic DNAs from 44 unrelated individuals revealed extensive insertion/deletion polymorphisms within the BstNI-type loci (PRB1, PRB2, PRB3 and PRB4) of the human proline-rich protein (PRP) multigene family. Ten length variants were cloned, including alleles at each of the four PRB loci, and in every case the region of length difference was localized to the tandemly repetitious third exon. DNA sequences covering the region of length variation were determined for seven of the alleles. The data indicate (1) that the PRB loci can be divided into two subtypes, PRB1 plus PRB2, and PRB3 plus PRB4, and (2) that the length differences result from different numbers of tandem repeats in the third exons. Variant chromosomes were also identified with different numbers of PRP loci resulting from homologous but unequal exchange between the PRB1 and PRB2 loci. The overall data are compatible with the observed length variants having been generated via homologous but unequal intragenic exchange. The results also indicate that these crossover events are sensitive to the amount of homology shared between the interacting DNA strands. Allelic length variants have arisen independently at least 20 times at the PRB loci, but only one has been detected at a PRH locus. Comparison of the detailed structures of the repetitious regions in PRB and PRH loci shows that the repeats in PRB genes are very similar to each other in sequence and in length. The PRH genes contain fewer repeats, which differ considerably in their individual lengths. These differences suggest that the larger number of length variants in PRB genes is related to their greater ease of homologous but unequal pairing compared to PRH genes.     

Population and genetic status of a critically endangered species in Korea, Euchresta japonica (Leguminosae), and their implications for conservation

The conservation status of Euchresta japonica Hook. f. ex Regel in Korea was investigated, with an emphasis on its genetic diversity. From field surveys, we obtained the only locality record for a wild population in Korea, which contained eight individuals. Genotyping was performed using nine microsatellite markers for all 20 remaining individuals, including those in ex situ collections. Among nine microsatellite loci that amplified within this group, five showed polymorphism with low hererozygosities, and a total of 12 multilocus genotypes were detected. Wild-specific alleles were detected in two individuals, and ex situ-specific alleles were detected in six individuals. Five individuals proved to have individual-specific alleles. The Korean population was also distinguished from the previously reported Japanese population by different alleles and higher diversity. To conserve this species more effectively in Korea, we recommend the following: (1) fencing the remaining wild population; (2) no relocation of wild individuals, as nine ex situ plants are already available; (3) complete ex situ conservation of all genetic diversity via clonal propagation of wild individuals; and (4) continuous protection and monitoring of the wild population.     

Microsatellite analysis of the inbreeding grass weed Barren Brome (Anisantha sterilis) reveals genetic diversity at the within- and between-farm scales

Nine microsatellites were used to screen 131 samples of Barren Brome (Anisantha sterilis: synonym Bromus sterilis) collected from within the fields of three English farms [from Oxfordshire (Oxon), Leicestershire (Leics) and Wiltshire (Wilts)] and eight seeds taken from samples of each of 10 farms across England, UK. Most individuals (approximately 97%) were homozygous. Polymorphism occurred at all nine loci in all three farms sampled at the field scale, and at most loci for nine of the other 10 farm samples. Between three and 11 alleles were found per locus. Gene diversity (D = 1 - summation operator p(i)2) ranged from 0.088 to 0.760. Polymorphism occurred among individuals within and among fields, and farms. Some alleles were found in only one farm. On the basis of the alleles at all nine loci in the 211 sampled plants, a total of 92 (44%) different genotypes was identified. Clustering analysis using the unweighted pair group method with arithmetic averages (UPGMA) for the combined Oxon, Wilts and Leics samples did not cluster them into their respective farms. Similarly, a phenogram of samples from all 10 farms showed considerable mixing of individuals with respect to farm origins. Identification of genotypes on field plans showed evidence of both spatial localization and mixing. Previous reports have suggested that A. sterilis is strictly inbreeding with little intrapopulation variation at the genetic level. Our data reveal that A. sterilis exists as numerous separate and genetically different lines, which are maintained by inbreeding but which very occasionally outcross. Possible explanations for this pattern of high genetic diversity are discussed.