The Effects of Background and Interference Selection on Patterns of Genetic Variation in Subdivided Populations

It is well known that most new mutations that affect fitness exert deleterious effects and that natural populations are often composed of subpopulations (demes) connected by gene flow. To gain a better understanding of the joint effects of purifying selection and population structure, we focus on a scenario where an ancestral population splits into multiple demes and study neutral diversity patterns in regions linked to selected sites. In the background selection regime of strong selection, we first derive analytic equations for pairwise coalescent times and F_ST as a function of time after the ancestral population splits into two demes and then construct a flexible coalescent simulator that can generate samples under complex models such as those involving multiple demes or nonconservative migration. We have carried out extensive forward simulations to show that the new methods can accurately predict diversity patterns both in the nonequilibrium phase following the split of the ancestral population and in the equilibrium between mutation, migration, drift, and selection. In the interference selection regime of many tightly linked selected sites, forward simulations provide evidence that neutral diversity patterns obtained from both the nonequilibrium and equilibrium phases may be virtually indistinguishable for models that have identical variance in fitness, but are nonetheless different with respect to the number of selected sites and the strength of purifying selection. This equivalence in neutral diversity patterns suggests that data collected from subdivided populations may have limited power for differentiating among the selective pressures to which closely linked selected sites are subject.     

Quantifying Population Genetic Differentiation from Next-Generation Sequencing Data

Over the past few years, new high-throughput DNA sequencing technologies have dramatically increased speed and reduced sequencing costs. However, the use of these sequencing technologies is often challenged by errors and biases associated with the bioinformatical methods used for analyzing the data. In particular, the use of naïve methods to identify polymorphic sites and infer genotypes can inflate downstream analyses. Recently, explicit modeling of genotype probability distributions has been proposed as a method for taking genotype call uncertainty into account. Based on this idea, we propose a novel method for quantifying population genetic differentiation from next-generation sequencing data. In addition, we present a strategy for investigating population structure via principal components analysis. Through extensive simulations, we compare the new method herein proposed to approaches based on genotype calling and demonstrate a marked improvement in estimation accuracy for a wide range of conditions. We apply the method to a large-scale genomic data set of domesticated and wild silkworms sequenced at low coverage. We find that we can infer the fine-scale genetic structure of the sampled individuals, suggesting that employing this new method is useful for investigating the genetic relationships of populations sampled at low coverage.     

Population genomics of parallel phenotypic evolution in stickleback across stream-lake ecological transitions

Understanding the genetics of adaptation is a central focus in evolutionary biology. Here, we use a population genomics approach to examine striking parallel morphological divergences of parapatric stream-lake ecotypes of threespine stickleback fish in three watersheds on the Haida Gwaii archipelago, western Canada. Genome-wide variation at greater than 1000 single nucleotide polymorphism loci indicate separate origin of giant lake and small-bodied stream fish within each watershed (mean F(ST) between watersheds = 0.244 and within = 0.114). Genome scans within watersheds identified a total of 21 genomic regions that are highly differentiated between ecotypes and are probably subject to directional selection. Most outliers were watershed-specific, but genomic regions undergoing parallel genetic changes in multiple watersheds were also identified. Interestingly, several of the stream-lake outlier regions match those previously identified in marine-freshwater and benthic-limnetic genome scans, indicating reuse of the same genetic loci in different adaptive scenarios. We also identified multiple new outlier loci, which may contribute to unique aspects of differentiation in stream-lake environments. Overall, our data emphasize the important role of ecological boundaries in driving both local and broadly occurring parallel genetic changes during adaptation.     

Heritability and quantitative genetic divergence of serotiny,a fire-persistence plant trait

Background and Aims

Although it is well known that fire acts as a selective pressure shaping plant phenotypes, there are no quantitative estimates of the heritability of any trait related to plant persistence under recurrent fires, such as serotiny. In this study, the heritability of serotiny in Pinus halepensis is calculated, and an evaluation is made as to whether fire has left a selection signature on the level of serotiny among populations by comparing the genetic divergence of serotiny with the expected divergence of neutral molecular markers (Q_ST–F_ST comparison).

Methods

A common garden of P. halepensis was used, located in inland Spain and composed of 145 open-pollinated families from 29 provenances covering the entire natural range of P. halepensis in the Iberian Peninsula and Balearic Islands. Narrow-sense heritability (h²) and quantitative genetic differentiation among populations for serotiny (Q_ST) were estimated by means of an ‘animal model’ fitted by Bayesian inference. In order to determine whether genetic differentiation for serotiny is the result of differential natural selection, Q_ST estimates for serotiny were compared with F_ST estimates obtained from allozyme data. Finally, a test was made of whether levels of serotiny in the different provenances were related to different fire regimes, using summer rainfall as a proxy for fire regime in each provenance.

Key Results

Serotiny showed a significant narrow-sense heritability (h²) of 0·20 (credible interval 0·09–0·40). Quantitative genetic differentiation among provenances for serotiny (Q_ST = 0·44) was significantly higher than expected under a neutral process (F_ST = 0·12), suggesting adaptive differentiation. A significant negative relationship was found between the serotiny level of trees in the common garden and summer rainfall of their provenance sites.

Conclusions

Serotiny is a heritable trait in P. halepensis, and selection acts on it, giving rise to contrasting serotiny levels among populations depending on the fire regime, and supporting the role of fire in generating genetic divergence for adaptive traits.     

Molecular and quantitative trait variation within and among small fragmented populations of the endangered plant species Psilopeganum sinense

Background and Aims

Natural selection and genetic drift are important evolutionary forces in determining genetic and phenotypic differentiation in plant populations. The extent to which these two distinct evolutionary forces affect locally adaptive quantitative traits has been well studied in common plant and animal species. However, we know less about how quantitative traits respond to selection pressures and drift in endangered species that have small population sizes and fragmented distributions. To address this question, this study assessed the relative strengths of selection and genetic drift in shaping population differentiation of phenotypic traits in Psilopeganum sinense, a naturally rare and recently endangered plant species.

Methods

Population differentiation at five quantitative traits (Q_ST) obtained from a common garden experiment was compared with differentiation at putatively neutral microsatellite markers (F_ST) in seven populations of P. sinense. Q_ST estimates were derived using a Bayesian hierarchical variance component method.

Key Results

Trait-specific Q_ST values were equal to or lower than F_ST. Neutral genetic diversity was not correlated with quantitative genetic variation within the populations of P. sinense.

Conclusions

Despite the prevalent empirical evidence for Q_ST > F_ST, the results instead suggest a definitive role of stabilizing selection and drift leading to phenotypic differentiation among small populations. Three traits exhibited a significantly lower Q_ST relative to F_ST, suggesting that populations of P. sinense might have experienced stabilizing selection for the same optimal phenotypes despite large geographical distances between populations and habitat fragmentation. For the other two traits, Q_ST estimates were of the same magnitude as F_ST, indicating that divergence in these traits could have been achieved by genetic drift alone. The lack of correlation between molecular marker and quantitative genetic variation suggests that sophisticated considerations are required for the inference of conservation measures of P. sinense from neutral genetic markers.     

Allozymic Differentiation Among Geographically Distant Populations of Patella vulgata (Mollusca, Patellogastropoda)

Patella vulgata is a boreal cold temperate species and is the dominant limpet in northern Europe. Few works have focussed on the population genetics of this species. Therefore, the aim of this work was to assess the degree of genetic and morphological differentiation of P. vulgata on a macroscale by using 20 allozyme loci and 6 morphological variables. Samples were taken from the following locations: Dingle Peninsula (Southwest Ireland), Port Erin (Southwest Isle of Man), St. Bees Head (north Cumbria, England), St. Agnes Head (north Cornwall, England), Cellar Beach (south Devon, England), Whitley Bay (north Newcastle-Upon-Tyne, England), Sines (Portugal), and Pointe de Chanchardon, La Rochelle (Bay of Biscay, France). Morphological variables were analysed by the multivariate Canonical discriminant analysis. Genetic variation was assessed by diversity measures such as polymorphism and heterozygosity; genetic subdivision of P. vulgata population was determined by the estimator θ of F _ST, and the genetic similarity between populations was measured by Nei’s genetic identity. No significant morphological differentiation was observed among samples. Moderate genetic population subdivision was observed (θ = 0.137±0.074) despite great geographic distances. The minimum genetic identity observed was between Ireland and France (I = 0.942) and maximum was observed between Portugal and north-east England (I=0.998). Two main groups were shown by UPGMA cluster analysis (I = 0.965). One formed by Irish, Manx, north Cumbria, and curiously, south Devon samples, while the second includes Portuguese, French, north-Newcastle-upon-thyne, and north Cornwall samples. No association (g = 0.956; p>0.050) was found between pair-wise genetic divergence and geographic distance separating subpopulations, mainly due to an unexpected pattern of genetic heterogeneity found in Southwest England.     

Detecting the footprints of divergent selection in oaks with linked markers

Genome scans are increasingly used to study ecological speciation, providing a useful genome-wide perspective on divergent selection in the presence of gene flow. Here, we compare current approaches to detect footprints of divergent selection in closely related species. We analyzed 192 individuals from two interfertile European temperate oak species using 30 nuclear microsatellites from eight linkage groups. These markers present little intraspecific differentiation and can be used in combination to assign individual genotypes to species. We first show that different outlier detection tests give somewhat different results, possibly due to model constraints. Second, using linkage information for these markers, we further characterize the signature of divergent selection in the presence of gene flow. In particular, we show that recombination estimates for regions with outlier markers are lower than those for a control region, in line with a prediction from ecological speciation theory. Most importantly, we show that analyses at the haplotype level can distinguish between truly divergent (bi-directional) selection and positive selection in one of the two species, offering a new and improved method for characterizing the speciation process.     

The influence of stochastic and selective forces in the population divergence of female colour polymorphism in damselflies of the genus Ischnura

Disentangling the relative importance and potential interactions of selection and genetic drift in driving phenotypic divergence of species is a classical research topic in population genetics and evolutionary biology. Here, we evaluate the role of stochastic and selective forces on population divergence of a colour polymorphism in seven damselfly species of the genus Ischnura, with a particular focus on I. elegans and I. graellsii. Colour-morph frequencies in Spanish I. elegans populations varied greatly, even at a local scale, whereas more similar frequencies were found among populations in eastern Europe. In contrast, I. graellsii and the other five Ischnura species showed little variation in colour-morph frequencies between populations. F(ST)-outlier analyses revealed that the colour locus deviated strongly from neutral expectations in Spanish populations of I. elegans, contrasting the pattern found in eastern European populations, and in I. graellsii, where no such discrepancy between morph divergence and neutral divergence could be detected. This suggests that divergent selection has been operating on the colour locus in Spanish populations of I. elegans, whereas processes such as genetic drift, possibly in combination with other forms of selection (such as negative frequency-dependent selection), appear to have been present in other regions, such as eastern Europe. Overall, the results indicate that both selective and stochastic processes operate on these colour polymorphisms, and suggest that the relative importance of factors varies between geographical regions.     

Divergent selection along climatic gradients in a rare central European endemic species,Saxifraga sponhemica

Background and Aims The effects of habitat fragmentation on quantitative genetic variation in plant populations are still poorly known. Saxifraga sponhemica is a rare endemic of Central Europe with a disjunct distribution, and a stable and specialized habitat of treeless screes and cliffs. This study therefore used S. sponhemica as a model species to compare quantitative and molecular variation in order to explore (1) the relative importance of drift and selection in shaping the distribution of quantitative genetic variation along climatic gradients; (2) the relationship between plant fitness, quantitative genetic variation, molecular genetic variation and population size; and (3) the relationship between the differentiation of a trait among populations and its evolvability.Methods Genetic variation within and among 22 populations from the whole distribution area of S. sponhemica was studied using RAPD (random amplified polymorphic DNA) markers, and climatic variables were obtained for each site. Seeds were collected from each population and germinated, and seedlings were transplanted into a common garden for determination of variation in plant traits.Key Results In contrast to previous results from rare plant species, strong evidence was found for divergent selection. Most population trait means of S. sponhemica were significantly related to climate gradients, indicating adaptation. Quantitative genetic differentiation increased with geographical distance, even when neutral molecular divergence was controlled for, and Q_ST exceeded F_ST for some traits. The evolvability of traits was negatively correlated with the degree of differentiation among populations (Q_ST), i.e. traits under strong selection showed little genetic variation within populations. The evolutionary potential of a population was not related to its size, the performance of the population or its neutral genetic diversity. However, performance in the common garden was lower for plants from populations with reduced molecular genetic variation, suggesting inbreeding depression due to genetic erosion.Conclusions The findings suggest that studies of molecular and quantitative genetic variation may provide complementary insights important for the conservation of rare species. The strong differentiation of quantitative traits among populations shows that selection can be an important force for structuring variation in evolutionarily important traits even for rare endemic species restricted to very specific habitats.     

Biased morph ratios and skewed mating success contribute to loss of genetic diversity in the distylous Pulmonaria officinalis

Background and Aims

In heterostylous plant species, skewed morph ratios are not uncommon and may arise from a range of factors. Despite the recognized importance of skewed morph ratios on overall reproductive success within populations, little is known about the impact of skewed morph ratios on population genetic diversity and differentiation in heterostylous species. This study specifically aimed to clarify the effect of population size and morph bias on population genetic diversity and differentiation in the temperate forest herb Pulmonaria officinalis. This species is characterized by a distylous breeding system and shows morph-specific differences in reproductive success.

Methods

Genetic diversity was determined for 27 P. officinalis populations in northern Belgium by using eight recently developed microsatellite markers. Multiple regressions were used to assess the relationship between genetic diversity, morph bias and population size, and F_ST-values were calculated for short- and long-styled morphs separately to study genetic differentiation as a function of morph type.

Key Results

For all genetic measures used, morph bias was more important in explaining patterns of genetic diversity than population size, and in all cases patterns of population genetic diversity followed a quadratic function, which showed a symmetrical decrease in genetic diversity with increasing morph bias. However, probably due to the reproductive advantage of L-morphs relative to S-morphs, maximum genetic diversity was found in populations showing an excess of L-morphs (60·7 % L-morph). On the other hand, no significant difference in pairwise genetic distances between populations was observed between L- (0·107) and S-morphs (0·106).

Conclusions

Our results indicate that significant deviations from equal morph ratios not only affect plant reproductive success but also population genetic diversity of heterostylous plant species. Hence, when defining conservation measures for populations of heterostylous plant species, morph ratios should be considered as an important trait affecting their long-term population viability.     

Novel measures of linkage disequilibrium that correct the bias due to population structure and relatedness

Among the several linkage disequilibrium measures known to capture different features of the non-independence between alleles at different loci, the most commonly used for diallelic loci is the r(2) measure. In the present study, we tackled the problem of the bias of r(2) estimate, which results from the sample structure and/or the relatedness between genotyped individuals. We derived two novel linkage disequilibrium measures for diallelic loci that are both extensions of the usual r(2) measure. The first one, r(S)(2), uses the population structure matrix, which consists of information about the origins of each individual and the admixture proportions of each individual genome. The second one, r(V)(2), includes the kinship matrix into the calculation. These two corrections can be applied together in order to correct for both biases and are defined either on phased or unphased genotypes.We proved that these novel measures are linked to the power of association tests under the mixed linear model including structure and kinship corrections. We validated them on simulated data and applied them to real data sets collected on Vitis vinifera plants. Our results clearly showed the usefulness of the two corrected r(2) measures, which actually captured 'true' linkage disequilibrium unlike the usual r(2) measure.     

Divergent selection in trailing- versus leading-edge populations of Biscutella laevigata

Background and Aims

Knowledge on how climate-induced range shifts might affect natural selection is crucial to understand the evolution of species ranges.

Methods

Using historical demographic perspectives gathered from regional-scale phylogeography on the alpine herb Biscutella laevigata, indirect inferences on gene flow and signature of selection based on AFLP genotyping were compared between local populations persisting at the trailing edge and expanding at the leading edge.

Key Results

Spatial autocorrelation revealed that gene flow was two times more restricted at the trailing edge and genome scans indicated divergent selection in this persisting population. In contrast, no pattern of selection emerged in the expanding population at the leading edge.

Conclusions

Historical effects may determine different architecture of genetic variation and selective patterns within local populations, what is arguably important to understand evolutionary processes acting across the species ranges.     

Invasion Genetics of New World Medflies: Testing Alternative Colonization Scenarios

The Mediterranean fruit fly (Ceratitis capitata) is an invasive agricultural pest with a wide host range and a nearly global distribution. Efforts to forgo the medfly's spread into the United States are dependent on an understanding of population dynamics in newly established populations elsewhere. To explore the potential influence of demographic and historical parameters in six medfly populations distributed from Mexico to Peru, we created population genetic null models using Monte Carlo simulations. Null expectations for genetic differentiation (F _ST) were compared with actual sequence variation from four highly polymorphic nuclear loci. Four colonization scenarios that were modeled led to unique genetic signatures that could be used to interpret empirical data. Unless current gene flow across Latin America was assumed to be very high, we could reject colonizations consisting of multiple introductions, each of low genetic diversity. Further, if simulated populations were small (N _e = 5 × 10² individuals per population), small invasions from a single source consistently produced F _ST values comparable to those currently observed in Latin America. In contrast, only large invasions from diverse sources were compatible with the observed data for large populations (N _e 5 × 10³). This study demonstrates that alternative population genetic hypotheses can be tested empirically even when departures from equilibrium are extreme, and that population genetic theory can be used to explore the processes that underlie biological invasions.     

Positive selection on the Plasmodium falciparum sporozoite threonine-asparagine-rich protein: analysis of isolates mainly from low endemic areas

The sporozoite threonine-asparagine-rich protein (STARP) of Plasmodium falciparum is an attractive target for a pre-erythrocytic stage malaria vaccine because both naturally acquired and experimentally induced anti-STARP antibodies can block sporozoite invasion of hepatocytes. To explore the extent of sequence variation, we surveyed nucleotide polymorphism across the entire gene, encompassing 2 exons and an intron, of 124 P. falciparum-infected blood samples from Thailand and 10 from 4 other endemic areas. In total 24 haplotypes were identified despite low-level nucleotide diversity at this locus. The mean number of nonsynonymous substitutions per nonsynonymous site (d(N)) significantly exceeded that of synonymous substitutions per synonymous site (d(S)), suggesting that the STARP gene has evolved under positive selection, probably from host immune pressure. The preponderance of conservative amino acid exchanges and a strongly biased T-nucleotide toward the third position of codons in repeat arrays have reflected simultaneous constraints on this molecule, probably from its respective unknown function and nucleotide composition. Sequence conservation in the STARP locus among clinical isolates from different disease endemic areas would not compromise vaccine incorporation.     

Genomic diversity and affinities in population groups of North West India: An analysis of Alu insertion and a single nucleotide polymorphism

The North West region of India is extremely important to understand the peopling of India, as it acted as a corridor to the foreign invaders from Eurasia and Central Asia. A series of these invasions along with multiple migrations led to intermixture of variable populations, strongly contributing to genetic variations. The present investigation was designed to explore the genetic diversities and affinities among the five major ethnic groups from North West India; Brahmin, Jat Sikh, Bania, Rajput and Gujjar. A total of 327 individuals of the abovementioned ethnic groups were analyzed for 4 Alu insertion marker loci (ACE, PV92, APO and D1) and a Single Nucleotide Polymorphism (SNP) rs2234693 in the intronic region of the ESR1 gene. Statistical analysis was performed to interpret the genetic structure and diversity of the population groups. Genotypes for ACE, APO, ESR1 and PV92 loci were found to be in Hardy–Weinberg equilibrium in all the ethnic groups, while significant departures were observed at the D1 locus in every investigated population after Bonferroni's correction. The average heterozygosity for all the loci in these ethnic groups was fairly substantial ranging from 0.3927 ± 0.1877 to 0.4333 ± 0.1416. Inbreeding coefficient indicated an overall 10% decrease in heterozygosity in these North West Indian populations. The gene differentiation among the populations was observed to be of the order of 0.013. Genetic distance estimates revealed that Gujjars were close to Banias and Jat Sikhs were close to Rajputs. Overall the study favored the recent division of the populations of North West India into largely endogamous groups. It was observed that the populations of North West India represent a more or less homogenous genetic entity, owing to their common ancestral history as well as geographical proximity.     

Microsatellite markers uncover cryptic species of Odontotermes (Termitoidae: Termitidae) from Peninsular Malaysia

Termites from the genus Odontotermes are known to contain numerous species complexes that are difficult to tell apart morphologically or with mitochondrial DNA sequences. We developed markers for one such cryptic species complex, that is, Odontotermes srinakarinensis sp. nov. from Maxwell Hill Forest Reserve (Perak, Malaysia), and characterised them using a sample of 41 termite workers from three voucher samples from the same area. We then genotyped 150 termite individuals from 23 voucher samples/colonies of this species complex from several sites in Peninsular Malaysia. We analysed their population by constructing dendograms from the proportion of shared-alleles between individuals and genetic distances between colonies; additionally, we examined the Bayesian clustering pattern of their genotype data. All methods of analysis indicated that there were two distinct clusters within our data set. After the morphologies of specimens from each cluster were reexamined, we were able to separate the two species morphologically and found that a single diagnostic character found on the mandibles of its soldiers could be used to separate the two species quite accurately. The additional species in the clade was identified as Odontotermes denticulatus after it was matched to type specimens at the NHM London and Cambridge Museum of Zoology.     

Detecting novel SNPs and breed-specific haplotypes at calpastatin gene in Iranian fat- and thin-tailed sheep breeds and their effects on protein structure

Calpastatin has been introduced as a potential candidate gene for growth and meat quality traits. In this study, genetic variability was investigated in the exon 6 and its intron boundaries of ovine CAST gene by PCR-SSCP analysis and DNA sequencing. Also a protein sequence and structural analysis were performed to predict the possible impact of amino acid substitutions on physicochemical properties and structure of the CAST protein. A total of 487 animals belonging to four ancient Iranian sheep breeds with different fat metabolisms, Lori-Bakhtiari and Chall (fat-tailed), Zel-Atabay cross-bred (medium fat-tailed) and Zel (thin-tailed), were analyzed. Eight unique SSCP patterns, representing eight different sequences or haplotypes, CAST-1, CAST-2 and CAST-6 to CAST-11, were identified. Haplotypes CAST-1 and CAST-2 were most common with frequency of 0.365 and 0.295. The novel haplotype CAST-8 had considerable frequency in Iranian sheep breeds (0.129). All the consensus sequences showed 98–99%, 94–98%, 92–93% and 82–83% similarity to the published ovine, caprine, bovine and porcine CAST locus sequences, respectively. Sequence analysis revealed four SNPs in intron 5 (C24T, G62A, G65T and T69-) and three SNPs in exon 6 (c.197A > T, c.282G > T and c.296C > G). All three SNPs in exon 6 were missense mutations which would result in p.Gln 66 Leu, p.Glu 94 Asp and p.Pro 99 Arg substitutions, respectively, in CAST protein. All three amino acid substitutions affected the physicochemical properties of ovine CAST protein including hydrophobicity, amphiphilicity and net charge and subsequently might influence its structure and effect on the activity of Ca2 + channels; hence, they might regulate calpain activity and afterwards meat tenderness and growth rate. The Lori-Bakhtiari population showed the highest heterozygosity in the ovine CAST locus (0.802). Frequency difference of haplotypes CAST-10 and CAST-8 between Lori-Bakhtiari (fat-tailed) and Zel (thin-tailed) breeds was highly significant (P < 0.001), indicating that these two haplotypes might be breed-specific haplotypes that distinguish between fat-tailed and thin-tailed sheep breeds.     

Population genetic structure of the Plasmodium vivax circumsporozoite protein (Pvcsp) in Sri Lanka

Molecular methods elucidate evolutionary and ecological processes in parasites, where interaction between hosts and parasites enlighten the evolution of parasite lifestyles and host defenses. Population genetics of Plasmodium vivax parasites accurately describe transmission dynamics of the parasites and evaluation of malaria control measures. As a first generation vaccine candidate against malaria, the Circumsporozoite Protein (CSP) has demonstrated significant potential in P. falciparum. Extensive polymorphism hinders the development of a potent malaria vaccine. Hence, the genetic diversity of Pvcsp was investigated for the first time in 60 Sri Lankan clinical isolates by obtaining the nucleotide sequence of the central repeat (CR) domain and examining the polymorphism of the peptide repeat motifs (PRMs), the genetic diversity indices and phylogenetic relationships. PCR amplicons determined size polymorphism of 610, 700 and 710 bp in Pvcsp of Sri Lanka where all amino acid sequences obtained were of the VK210 variant, consisting variable repeats of 4 different PRMs. The two most abundant PRMs of the CR domain, GDRADGQPA and GDRAAGQPA consisted ~ 2-4 repeats, while GNRAAGQPA was unique to the island. Though, different nucleotide sequences termed repeat allotypes (RATs) were observed for each PRM, these were synonymous contributing to a less polymorphic CR domain. The genetic diversity of Pvcsp in Sri Lanka was due to the number of repetitive peptide repeat motifs, point mutations, and intragenic recombination. The 19 amino acid haplotypes defined were exclusive to Sri Lanka, whereas the 194 Pvcsp sequences of global isolates generated 57 more distinct a.a. haplotypes of the VK210 variant. Strikingly, the CR domain of both VK210 and VK247 variants was under purifying selection interpreting the scarcity of CSP non-synonymous polymorphisms. Insights to the distribution of RATs in the CR region with geographic clustering of the P. vivax VK210 variant were revealed. The cladogram reiterated this unique geographic clustering of local (VK210) and global isolates (VK210 and VK247), which was further validated by the elevated fixation index values of the VK210 variant.     

Mutation of a vitelline membrane protein, BmEP80, is responsible for the silkworm “Ming” lethal egg mutant

The egg stage is an important stage in the silkworm (Bombyx mori) life cycle. Normal silkworm eggs are usually short, elliptical, and laterally flattened, with a sometimes hollowed surface on the lateral side. However, the eggs laid by homozygous recessive “Ming” lethal egg mutants (l-e^m) lose water and become concaved around 1 h, ultimately exhibiting a triangular shape on the egg surfaces. We performed positional cloning, and narrowed down the region containing the gene responsible for the l-e^m mutant to 360 kb on chromosome 10 using 2287 F₂ individuals. Using expression analysis and RNA interference, the best l-e^m candidate gene was shown to be BmEP80. The results of the inverse polymerase chain reaction showed that an ~ 1.9 kb region from the 3′ untranslated region of BmVMP23 to the forepart of BmEP80 was replaced by a > 100 kb DNA fragment in the l-e^m mutant. Several eggs laid by the normal moths injected with BmEP80 small interfering RNAs were evidently depressed and exhibited a triangular shape on the surface. The phenotype exhibited was consistent with the eggs laid by the l-e^m mutant. Moreover, two-dimensional gel electrophoresis showed that the BmEP80 protein was expressed in the ovary from the 9th day of the pupa stage to eclosion in the wild-type silkworm, but was absent in the l-e^m mutant. These results indicate that BmEP80 is responsible for the l-e^m mutation.     

Protective effect of the KIR2DS1 gene in atopic dermatitis

Atopic dermatitis (AD) is a common skin disease of complex etiology including affected humoral and cellular immune responses. The role of NK cells in development of this disease has been recently postulated, but is still poorly documented. The current study was undertaken to determine the impact of genes for the most polymorphic NK cell receptors, known as killer cell immunoglobulin-like receptors (KIRs), on the development of AD.