Genome Properties and Prospects of Genomic Prediction of Hybrid Performance in a Breeding Program of Maize

Maize (Zea mays L.) serves as model plant for heterosis research and is the crop where hybrid breeding was pioneered. We analyzed genomic and phenotypic data of 1254 hybrids of a typical maize hybrid breeding program based on the important Dent × Flint heterotic pattern. Our main objectives were to investigate genome properties of the parental lines (e.g., allele frequencies, linkage disequilibrium, and phases) and examine the prospects of genomic prediction of hybrid performance. We found high consistency of linkage phases and large differences in allele frequencies between the Dent and Flint heterotic groups in pericentromeric regions. These results can be explained by the Hill–Robertson effect and support the hypothesis of differential fixation of alleles due to pseudo-overdominance in these regions. In pericentromeric regions we also found indications for consistent marker–QTL linkage between heterotic groups. With prediction methods GBLUP and BayesB, the cross-validation prediction accuracy ranged from 0.75 to 0.92 for grain yield and from 0.59 to 0.95 for grain moisture. The prediction accuracy of untested hybrids was highest, if both parents were parents of other hybrids in the training set, and lowest, if none of them were involved in any training set hybrid. Optimizing the composition of the training set in terms of number of lines and hybrids per line could further increase prediction accuracy. We conclude that genomic prediction facilitates a paradigm shift in hybrid breeding by focusing on the performance of experimental hybrids rather than the performance of parental lines in testcrosses.     

Predicting heterosis for egg production traits in crossbred offspring of individual White Leghorn sires using genome-wide SNP data

Background

The development of a reliable method to predict heterosis would greatly improve the efficiency of commercial crossbreeding schemes. Extending heterosis prediction from the line level to the individual sire level would take advantage of variation between sires from the same pure line, and further increase the use of heterosis in crossbreeding schemes. We aimed at deriving the theoretical expectation for heterosis due to dominance in the crossbred offspring of individual sires, and investigating how much extra variance in heterosis can be explained by predicting heterosis at the individual sire level rather than at the line level. We used 53 421 SNP (single nucleotide polymorphism) genotypes of 3427 White Leghorn sires, allele frequencies of six White Leghorn dam-lines and cage-based records on egg number and egg weight of ~210 000 crossbred hens.

Results

We derived the expected heterosis for the offspring of individual sires as the between- and within-line genome-wide heterozygosity excess in the offspring of a sire relative to the mean heterozygosity of the pure lines. Next, we predicted heterosis by regressing offspring performance on the heterozygosity excess. Predicted heterosis ranged from 7.6 to 16.7 for egg number, and from 1.1 to 2.3 grams for egg weight. Between-line differences accounted for 99.0% of the total variance in predicted heterosis, while within-line differences among sires accounted for 0.7%.

Conclusions

We show that it is possible to predict heterosis at the sire level, thus to distinguish between sires within the same pure line with offspring that show different levels of heterosis. However, based on our data, variation in genome-wide predicted heterosis between sires from the same pure line was small; most differences were observed between lines. We hypothesise that this method may work better if predictions are based on SNPs with identified dominance effects.     

Monthly to interannual variability of microbial eukaryote assemblages at four depths in the eastern North Pacific

The monthly, seasonal and interannual variability of microbial eukaryote assemblages were examined at 5 m, the deep chlorophyll maximum, 150 m and 500 m at the San Pedro Ocean Time-series station (eastern North Pacific). The depths spanned transitions in temperature, light, nutrients and oxygen, and included a persistently hypoxic environment at 500 m. Terminal restriction fragment length polymorphism was used for the analysis of 237 samples that were collected between September 2000 and December 2010. Spatiotemporal variability patterns of microeukaryote assemblages indicated the presence of distinct shallow and deep communities at the SPOT station, presumably reflecting taxa that were specifically adapted for the conditions in those environments. Community similarity values between assemblages collected 1 month apart at each depth ranged between ∼20% and ∼84% (averages were ∼50–59%). The assemblage at 5 m was temporally more dynamic than deeper assemblages and also displayed substantial interannual variability during the first ∼3 years of the study. Evidence of seasonality was detected for the microbial eukaryote assemblage at 5 m between January 2008 and December 2010 and at 150 m between September 2000 and December 2003. Seasonality was not detected for assemblages at the deep chlorophyll a maximum, which varied in depth seasonally, or at 500 m. Microbial eukaryote assemblages exhibited cyclical patterns in at least 1 year at each depth, implying an annual resetting of communities. Substantial interannual variability was detected for assemblages at all depths and represented the largest source of temporal variability in this temperate coastal ecosystem.     

Something Old and Something New: Wedding Recombinant Inbred Lines with Traditional Line Cross Analysis Increases Power to Describe Gene Interactions

In this paper we present a novel approach to quantifying genetic architecture that combines recombinant inbred lines (RIL) with line cross analysis (LCA). LCA is a method of quantifying directional genetic effects (i.e. summed effects of all loci) that differentiate two parental lines. Directional genetic effects are thought to be critical components of genetic architecture for the long term response to selection and as a cause of inbreeding depression. LCA typically begins with two inbred parental lines that are crossed to produce several generations such as F1, F2, and backcrosses to each parent. When a RIL population (founded from the same P1 and P2 as was used to found the line cross population) is added to the LCA, the sampling variance of several nonadditive genetic effect estimates is greatly reduced. Specifically, estimates of directional dominance, additive x additive, and dominance x dominance epistatic effects are reduced by 92%, 94%, and 56% respectively. The RIL population can be simultaneously used for QTL identification, thus uncovering the effects of specific loci or genomic regions as elements of genetic architecture. LCA and QTL mapping with RIL provide two qualitatively different measures of genetic architecture with the potential to overcome weaknesses of each approach alone. This approach provides cross-validation of the estimates of additive and additive x additive effects, much smaller confidence intervals on dominance, additive x additive and dominance x dominance estimates, qualitatively different measures of genetic architecture, and the potential when used together to balance the weaknesses of LCA or RIL QTL analyses when used alone.     

Survival of Desulfotomaculum spores from estuarine sediments after serial autoclaving and high-temperature exposure

Bacterial spores are widespread in marine sediments, including those of thermophilic, sulphate-reducing bacteria, which have a high minimum growth temperature making it unlikely that they grow in situ. These Desulfotomaculum spp. are thought to be from hot environments and are distributed by ocean currents. Their cells and spores upper temperature limit for survival is unknown, as is whether they can survive repeated high-temperature exposure that might occur in hydrothermal systems. This was investigated by incubating estuarine sediments significantly above (40–80 °C) maximum in situ temperatures (∼23 °C), and with and without prior triple autoclaving. Sulphate reduction occurred at 40–60 °C and at 60 °C was unaffected by autoclaving. Desulfotomaculum sp. C1A60 was isolated and was most closely related to the thermophilic D. kuznetsovii^T (∼96% 16S rRNA gene sequence identity). Cultures of Desulfotomaculum sp. C1A60, D. kuznetsovii^Tand D. geothermicum B2T survived triple autoclaving while other related Desulfotomaculum spp. did not, although they did survive pasteurisation. Desulfotomaculum sp. C1A60 and D. kuznetsovii cultures also survived more extreme autoclaving (C1A60, 130 °C for 15 min; D. kuznetsovii, 135 °C for 15 min, maximum of 154 °C reached) and high-temperature conditions in an oil bath (C1A60, 130° for 30 min, D. kuznetsovii 140 °C for 15 min). Desulfotomaculum sp. C1A60 with either spores or predominantly vegetative cells demonstrated that surviving triple autoclaving was due to spores. Spores also had very high culturability compared with vegetative cells (∼30 × higher). Combined extreme temperature survival and high culturability of some thermophilic Desulfotomaculum spp. make them very effective colonisers of hot environments, which is consistent with their presence in subsurface geothermal waters and petroleum reservoirs.     

Preservation of Differentiation and Clonogenic Potential of Human Hematopoietic Stem and Progenitor Cells during Lyophilization and Ambient Storage

Progenitor cell therapies show great promise, but their potential for clinical applications requires improved storage and transportation. Desiccated cells stored at ambient temperature would provide economic and practical advantages over approaches employing cell freezing and subzero temperature storage. The objectives of this study were to assess a method for loading the stabilizing sugar, trehalose, into hematopoietic stem and progenitor cells (HPC) and to evaluate the effects of subsequent freeze-drying and storage at ambient temperature on differentiation and clonogenic potential. HPC were isolated from human umbilical cord blood and loaded with trehalose using an endogenous cell surface receptor, termed P2Z. Solution containing trehalose-loaded HPC was placed into vials, which were transferred to a tray freeze-dryer and removed during each step of the freeze-drying process to assess differentiation and clonogenic potential. Control groups for these experiments were freshly isolated HPC. Control cells formed 1450±230 CFU-GM, 430±140 BFU-E, and 50±40 CFU-GEMM per 50 µL. Compared to the values for the control cells, there was no statistical difference observed for cells removed at the end of the freezing step or at the end of primary drying. There was a gradual decrease in the number of CFU-GM and BFU-E for cells removed at different temperatures during secondary drying; however, there were no significant differences in the number of CFU-GEMM. To determine storage stability of lyophilized HPC, cells were stored for 4 weeks at 25°C in the dark. Cells reconstituted immediately after lyophilization produced 580±90 CFU-GM (∼40%, relative to unprocessed controls p<0.0001), 170±70 BFU-E (∼40%, p<0.0001), and 41±22 CFU-GEMM (∼82%, p = 0.4171), and cells reconstituted after 28 days at room temperature produced 513±170 CFU-GM (∼35%, relative to unprocessed controls, p<0.0001), 112±68 BFU-E (∼26%, p<0.0001), and 36±17 CFU-GEMM (∼82%, p = 0.2164) These studies are the first to document high level retention of CFU-GEMM following lyophilization and storage for 4 weeks at 25°C. This type of flexible storage stability would potentially permit the ability to ship and store HPC without the need for refrigeration.     

Low levels of realized seed and pollen gene flow and strong spatial genetic structure in a small,isolated and fragmented population of the tropical tree Copaifera langsdorffii Desf

Over the past century, the Brazilian Atlantic forest has been reduced to small, isolated fragments of forest. Reproductive isolation theories predict a loss of genetic diversity and increases in inbreeding and spatial genetic structure (SGS) in such populations. We analysed eight microsatellite loci to investigate the pollen and seed dispersal patterns, genetic diversity, inbreeding and SGS of the tropical tree Copaifera langsdorffii in a small (4.8 ha), isolated population. All 112 adult trees and 128 seedlings found in the stand were sampled, mapped and genotyped. Seedlings had significantly lower levels of genetic diversity (A=16.5±0.45, mean±95% s.e.; H_e=0.838±0.006) than did adult trees (A=23.2±0.81; H_e=0.893±0.030). Parentage analysis did not indicate any seed immigration (m_seeds=0) and the pollen immigration rate was very low (m_pollen=0.047). The average distance of realized pollen dispersal within the stand was 94 m, with 81% of the pollen travelling <150 m. A significant negative correlation was found between the frequency and distance of pollen dispersal (r=−0.79, P<0.01), indicating that short-distance pollinations were more frequent. A significant SGS for both adults (∼50 m) and seedlings (∼20 m) was also found, indicating that most of the seeds were dispersed over short distances. The results suggested that the spatial isolation of populations by habitat fragmentation can restrict seed and pollen gene flow, increase SGS and affect the genetic diversity of future generations.     

Short divalent ethacrynic amides as pro-inhibitors of glutathione S-transferase isozyme Mu and potent sensitisers of cisplatin-resistant ovarian cancers

The linking of ethacrynic acid with ethylenediamine and 1,4-butanediamine gave EDEA and BDEA, respectively, as membrane-permeable divalent pro-inhibitors of glutathione S-transferase (GST). Their divalent glutathione conjugates showed subnanomolar inhibition and divalence-binding to GSTmu (GSTM) (PDB: 5HWL) at ∼0.35 min⁻¹. In cisplatin-resistant SK-OV-3, COC1, SGC7901 and A549 cells, GSTM activities probed by 15 nM BDEA or EDEA revealed 5-fold and 1.0-fold increases in cisplatin-resistant SK-OV-3 and COC1 cells, respectively, in comparison with the susceptible parental cells. Being tolerable by HEK293 and LO2 cells, BDEA at 0.2 μM sensitised resistant SK-OV-3 and COC1 cells by ∼3- and ∼5-folds, respectively, released cytochrome c and increased apoptosis; EDEA at 1.0 μM sensitised resistant SK-OV-3 and A549 cells by ∼5- and ∼7-fold, respectively. EDEA at 1.7 μg/g sensitised resistant SK-OV-3 cells to cisplatin at 3.3 μg/g in nude mouse xenograft model. BDEA and EDEA are promising leads for probing cellular GSTM and sensitising cisplatin-resistant ovarian cancers.     

Association Mapping across Numerous Traits Reveals Patterns of Functional Variation in Maize

Phenotypic variation in natural populations results from a combination of genetic effects, environmental effects, and gene-by-environment interactions. Despite the vast amount of genomic data becoming available, many pressing questions remain about the nature of genetic mutations that underlie functional variation. We present the results of combining genome-wide association analysis of 41 different phenotypes in ∼5,000 inbred maize lines to analyze patterns of high-resolution genetic association among of 28.9 million single-nucleotide polymorphisms (SNPs) and ∼800,000 copy-number variants (CNVs). We show that genic and intergenic regions have opposite patterns of enrichment, minor allele frequencies, and effect sizes, implying tradeoffs among the probability that a given polymorphism will have an effect, the detectable size of that effect, and its frequency in the population. We also find that genes tagged by GWAS are enriched for regulatory functions and are ∼50% more likely to have a paralog than expected by chance, indicating that gene regulation and gene duplication are strong drivers of phenotypic variation. These results will likely apply to many other organisms, especially ones with large and complex genomes like maize.     

Demographic factors shaped diversity in the two gene pools of wild common bean Phaseolus vulgaris L.

Wild common bean (Phaseolus vulgaris L.) is distributed throughout the Americas from Mexico to northern Argentina. Within this range, the species is divided into two gene pools (Andean and Middle American) along a latitudinal gradient. The diversity of 24 wild common bean genotypes from throughout the geographic range of the species was described by using sequence data from 13 loci. An isolation–migration model was evaluated using a coalescent analysis to estimate multiple demographic parameters. Using a Bayesian approach, Andean and Middle American subpopulations with high percentage of parentages were observed. Over all loci, the Middle American gene pool was more diverse than the Andean gene pool (π_sil=0.0089 vs 0.0068). The two subpopulations were strongly genetically differentiated over all loci (F_st=0.29). It is estimated that the two current wild gene pools diverged from a common ancestor ∼111 000 years ago. Subsequently, each gene pool underwent a bottleneck immediately after divergence and lasted ∼40 000 years. The Middle American bottleneck population size was ∼46% of the ancestral population size, whereas the Andean was 26%. Continuous asymmetric gene flow was detected between the two gene pools with a larger number of migrants entering Middle American gene pool from the Andean gene pool. These results suggest that because of the complex population structure associated with the ancestral divergence, subsequent bottlenecks in each gene pool, gene pool-specific domestication and intense selection within each gene pool by breeders; association mapping would best be practised within each common bean gene pool.     

Systematic differences in the response of genetic variation to pedigree and genome-based selection methods

Genomic selection (GS) is a DNA-based method of selecting for quantitative traits in animal and plant breeding, and offers a potentially superior alternative to traditional breeding methods that rely on pedigree and phenotype information. Using a 60 K SNP chip with markers spaced throughout the entire chicken genome, we compared the impact of GS and traditional BLUP (best linear unbiased prediction) selection methods applied side-by-side in three different lines of egg-laying chickens. Differences were demonstrated between methods, both at the level and genomic distribution of allele frequency changes. In all three lines, the average allele frequency changes were larger with GS, 0.056 0.064 and 0.066, compared with BLUP, 0.044, 0.045 and 0.036 for lines B1, B2 and W1, respectively. With BLUP, 35 selected regions (empirical P<0.05) were identified across the three lines. With GS, 70 selected regions were identified. Empirical thresholds for local allele frequency changes were determined from gene dropping, and differed considerably between GS (0.167–0.198) and BLUP (0.105–0.126). Between lines, the genomic regions with large changes in allele frequencies showed limited overlap. Our results show that GS applies selection pressure much more locally than BLUP, resulting in larger allele frequency changes. With these results, novel insights into the nature of selection on quantitative traits have been gained and important questions regarding the long-term impact of GS are raised. The rapid changes to a part of the genetic architecture, while another part may not be selected, at least in the short term, require careful consideration, especially when selection occurs before phenotypes are observed.     

Anticancer activity and toxicity of new quaternary ammonium geldanamycin derivative salts and their mixtures with potentiators

Geldanamycin (GDM) has been modified by different type neutral/acidic/basic substituents (1–7) and by quinuclidine motif (8), transformed into ammonium salts (9–13) at C(17). These compounds have been characterised by spectroscopic and x-ray methods. Derivative 8 shows better potency than GDM in MCF-7, MDA-MB-231, A549 and HeLa (IC₅₀s = 0.09–1.06 µM). Transformation of 8 into salts 9–13 reduces toxicity (by 11-fold) at attractive potency, e.g. MCF-7 cell line (IC₅₀∼2 µM). Our studies show that higher water solubility contributes to lower toxicity of salts than GDM in healthy CCD39Lu and HDF cells. The use of 13 mixtures with potentiators PEI and DOX enhanced anticancer effects from IC₅₀∼2 µM to IC₅₀∼0.5 µM in SKBR-3, SKOV-3, and PC-3 cancer cells, relative to 13. Docking studies showed that complexes between quinuclidine-bearing 8–13 and Hsp90 are stabilised by extra hydrophobic interactions between the C(17)-arms and K58 or Y61 of Hsp90.     

Prediction of Heterosis from DNA Fingerprints in Chickens

To assess the value of DNA fingerprints for the prediction of heterosis in chickens, retrospective analyses of data from three crossbreeding experiments and DNA fingerprints (DFP) of parental strains were conducted using two minisatellite and one middle-repetitive DNA probes. DFP bands were assessed on pooled DNA samples of 10-15 individuals per parental genetic group. The number of DFP bands evaluated in the experiments ranged from 81 to 139. The probes varied in their predictive value, but predictability of heterosis generally increased with multiple probes. Highly significant correlations (0.68-0.87) between band sharing ratios (SH) and heterosis were found in 25 crosses of White Leghorns in the first egg production cycle for age at sexual maturity, egg production, and mature body weight: traits with heterosis of 10% or more of the means. Regressions of SH explained 78.4% of the variation in heterosis in age at sexual maturity, 60.2% in egg production and 46.4% in mature body weight. For ``broiler' traits with heterosis of <1%, none of the correlations, based on 13 crosses, were significant. It was concluded that multilocus probe DFP of pooled DNA samples show promise as predictors of heterosis.     

Systematic Screening of Drosophila Deficiency Mutations for Embryonic Phenotypes and Orphan Receptor Ligands

This paper defines a collection of Drosophila deletion mutations (deficiencies) that can be systematically screened for embryonic phenotypes, orphan receptor ligands, and genes affecting protein localization. It reports the results of deficiency screens we have conducted that have revealed new axon guidance phenotypes in the central nervous system and neuromuscular system and permitted a quantitative assessment of the number of potential genes involved in regulating guidance of specific motor axon branches. Deficiency “kits” that cover the genome with a minimum number of lines have been established to facilitate gene mapping. These kits cannot be systematically analyzed for phenotypes, however, since embryos homozygous for many deficiencies in these kits fail to develop due to the loss of key gene products encoded within the deficiency. To create new kits that can be screened for phenotype, we have examined the development of the nervous system in embryos homozygous for more than 700 distinct deficiency mutations. A kit of ∼400 deficiency lines for which homozygotes have a recognizable nervous system and intact body walls encompasses >80% of the genome. Here we show examples of screens of this kit for orphan receptor ligands and neuronal antigen expression. It can also be used to find genes involved in expression, patterning, and subcellular localization of any protein that can be visualized by antibody staining. A subset kit of 233 deficiency lines, for which homozygotes develop relatively normally to late stage 16, covers ∼50% of the genome. We have screened it for axon guidance phenotypes, and we present examples of new phenotypes we have identified. The subset kit can be used to screen for phenotypes affecting all embryonic organs. In the future, these deficiency kits will allow Drosophila researchers to rapidly and efficiently execute genome-wide anatomical screens that require examination of individual embryos at high magnification.     

Pigment epithelium-derived factor expression prolongs survival and enhances the cytotoxicity of low-dose chemotherapy in castration-refractory prostate cancer

There is currently no cure for advanced castration-refractory prostate cancer (CRPC) despite the recent approval of several new therapeutic agents. We report here the anti-tumor effect of the angio-inhibitory pigment epithelium-derived factor (PEDF) in the metastatic LNCaP-derivative CRPC CL1 model and explore PEDF anti-neoplasic efficacy in combination with low-dose chemotherapy. Androgen-sensitive LNCaP and CRPC PC3 cell lines were examined as comparison. Using a retroviral expression system, we showed that PEDF limited the proliferation of all prostatic cell lines tested; an effect attributed to interleukin 8 (IL8)-CXCR1/IL8RA inhibition. PEDF also reduced the number and size of 3D tumor spheroids in vitro, but only induced cell differentiation in CRPC spheroids. Similarly, PEDF inhibited the migration of CRPC cells suggesting both anti-proliferative and anti-migratory functions. In vivo, PEDF decreased by 85% and 65% the growth of subcutaneous (s.c.) PC3 and CL1 tumors, respectively. In the CL1 orthotopic model, tumor intake with lethal metastases was found in all animals; nevertheless, PEDF prolonged the median survival of tumor-bearing mice (95% confidence interval: 53±0.001 to 57±1 days). Accordingly, PEDF delayed the emergence of skeletal-related event in intra-tibial xenografts. Next, we evaluated low-dose docetaxel (DTX; 5, 1, 0.5 mg/kg) or cyclophosphamide (CTX; 10–20 mg/kg) on established s.c. PC3 tumors that conditionally express PEDF anti-tumoral epitope/NT3. Although NT3–DTX-5 mg/kg combination was inefficient, NT3–DTX-1 mg/kg and -0.5 mg/kg inhibited by 95% and 87.8%, respectively, tumor growth compared with control and induced tumor stasis. Both NT3–CTX combinations were advantageous. Inversely, PEDF–DTX-5 mg/kg and PEDF–CTX-10 mg/kg delayed the most CL1 tumor growth (15, 11 and 5 days for PEDF–DTX-5 mg/kg, PEDF–CTX-10 mg/kg and single treatments, respectively) with elevated apoptosis and serum thrombospondin-1 as possible mechanism and marker, respectively. As well, both PEDF–CTX-10 mg/kg and PEDF–DTX-5 mg/kg prolonged significantly the survival of tumor-bearing mice compared with single treatments. Metastases were reduced in PEDF–DTX-5 mg/kg compared with other treatments, suggesting that PEDF–DTX delayed metastases formation. Our results advocate that PEDF/low-dose chemotherapy may represent a new therapeutic alternative for CRPC.     

DNA is a co-factor for its own replication in Xenopus egg extracts

Soluble Xenopus egg extracts efficiently replicate added plasmids using a physiological mechanism, and thus represent a powerful system to understand vertebrate DNA replication. Surprisingly, DNA replication in this system is highly sensitive to plasmid concentration, being undetectable below ∼10 pM and highly efficient above ∼75 pM. DNA replication at the high plasmid concentration does not require plasmid–plasmid contacts, since replication is not inhibited when plasmids are immobilized in agarose prior to addition of egg extract. The absence of replication at low plasmid concentration is due to a defect in the assembly of pre-replication complexes (pre-RCs). pre-RC assembly requires contact-independent communication between plasmids. Our results show that in Xenopus egg extracts, aggregation of multiple replication forks is not required for efficient replication of plasmid DNA, and they suggest that DNA functions as a co-factor for its own duplication.     

Patterns of pollen dispersal in a small population of the Canarian endemic palm (Phoenix canariensis)

The genetic diversity of small populations is greatly influenced by local dispersal patterns and genetic connectivity among populations, with pollen dispersal being the major component of gene flow in many plants species. Patterns of pollen dispersal, mating system parameters and spatial genetic structure were investigated in a small isolated population of the emblematic palm Phoenix canariensis in Gran Canaria island (Canary Islands). All adult palms present in the study population (n=182), as well as 616 seeds collected from 22 female palms, were mapped and genotyped at 8 microsatellite loci. Mating system analysis revealed an average of 5.8 effective pollen donors (N_ep) per female. There was strong variation in correlated paternity rates across maternal progenies (ranging from null to 0.9) that could not be explained by the location and density of local males around focal females. Paternity analysis revealed a mean effective pollen dispersal distance of ∼71 m, with ∼70% of effective pollen originating from a distance of <75 m, and 90% from <200 m. A spatially explicit mating model indicated a leptokurtic pollen dispersal kernel, significant pollen immigration (12%) from external palm groves and a directional pollen dispersal pattern that seems consistent with local altitudinal air movement. No evidence of inbreeding or genetic diversity erosion was found, but spatial genetic structure was detected in the small palm population. Overall, the results suggest substantial pollen dispersal over the studied population, genetic connectivity among different palm groves and some resilience to neutral genetic erosion and subsequently to fragmentation.     

Recent Rapid Rise of a Permethrin Knock Down Resistance Allele in Aedes aegypti in México

Background

Aedes aegypti, the ‘yellow fever mosquito’, is the primary vector to humans of dengue and yellow fever flaviviruses (DENV, YFV), and is a known vector of the chikungunya alphavirus (CV). Because vaccines are not yet available for DENV or CV or are inadequately distributed in developing countries (YFV), management of Ae. aegypti remains the primary option to prevent and control outbreaks of the diseases caused by these arboviruses. Permethrin is one of the most widely used active ingredients in insecticides for suppression of adult Ae. aegypti. In 2007, we documented a replacement mutation in codon 1,016 of the voltage-gated sodium channel gene (para) of Ae. aegypti that encodes an isoleucine rather than a valine and confers resistance to permethrin. Ile1,016 segregates as a recessive allele conferring knockdown resistance to homozygous mosquitoes at 5–10 µg of permethrin in bottle bioassays.

Methods and Findings

A total of 81 field collections containing 3,951 Ae. aegypti were made throughout México from 1996 to 2009. These mosquitoes were analyzed for the frequency of the Ile1,016 mutation using a melting-curve PCR assay. Dramatic increases in frequencies of Ile1,016 were recorded from the late 1990''s to 2006–2009 in several states including Nuevo León in the north, Veracruz on the central Atlantic coast, and Yucatán, Quintana Roo and Chiapas in the south. From 1996 to 2000, the overall frequency of Ile1,016 was 0.04% (95% confidence interval (CI95) = 0.12%; n = 1,359 mosquitoes examined). The earliest detection of Ile1,016 was in Nuevo Laredo on the U.S. border in 1997. By 2003–2004 the overall frequency of Ile1,016 had increased ∼100-fold to 2.7% (±0.80% CI95; n = 808). When checked again in 2006, the frequency had increased slightly to 3.9% (±1.15% CI95; n = 473). This was followed in 2007–2009 by a sudden jump in Ile1,016 frequency to 33.2% (±1.99% CI95; n = 1,074 mosquitoes). There was spatial heterogeneity in Ile1,016 frequencies among 2007–2008 collections, which ranged from 45.7% (±2.00% CI95) in the state of Veracruz to 51.2% (±4.36% CI95) in the Yucatán peninsula and 14.5% (±2.23% CI95) in and around Tapachula in the state of Chiapas. Spatial heterogeneity was also evident at smaller geographic scales. For example within the city of Chetumal, Quintana Roo, Ile1,016 frequencies varied from 38.3%–88.3%. A linear regression analysis based on seven collections from 2007 revealed that the frequency of Ile1,016 homozygotes accurately predicted knockdown rate for mosquitoes exposed to permethrin in a bioassay (R² = 0.98).

Conclusions

We have recorded a dramatic increase in the frequency of the Ile1,016 mutation in the voltage-gated sodium channel gene of Ae. aegypti in México from 1996 to 2009. This may be related to heavy use of permethrin-based insecticides in mosquito control programs. Spatial heterogeneity in Ile1,016 frequencies in 2007 and 2008 collections may reflect differences in selection pressure or in the initial frequency of Ile1,016. The rapid recent increase in Ile1,016 is predicted by a simple model of positive directional selection on a recessive allele. Unfortunately this model also predicts rapid fixation of Ile1,016 unless there is negative fitness associated with Ile1,016 in the absence of permethrin. If so, then spatial refugia of susceptible Ae. aegypti or rotational schedules of different classes of adulticides could be established to slow or prevent fixation of Ile1,016.