Engineering the genomes of wild insect populations: Challenges, and opportunities provided by synthetic Medea selfish genetic elements

Advances in insect transgenesis and our knowledge of insect physiology and genomics are making it possible to create transgenic populations of beneficial or pest insects that express novel traits. There are contexts in which we may want the transgenes responsible for these traits to spread so that all individuals within a wild population carry them, a process known as population replacement. Transgenes of interest are unlikely to confer an overall fitness benefit on those who carry them. Therefore, an essential component of any population replacement strategy is the presence of a drive mechanism that will ensure the spread of linked transgenes. We discuss contexts in which population replacement might be desirable and the requirements a drive system must satisfy to be both effective and safe. We then describe the creation of synthetic Medea elements, the first selfish genetic elements synthesized de novo, with the capability of driving population replacement, in this case in Drosophila. The strategy used to create DrosophilaMedea is applicable to a number of other insect species and the Medea system satisfies key requirements for scientific and social acceptance. Finally, we highlight several challenges to implementing population replacement in the wild.     

The distribution and spread of naturally occurring Medea selfish genetic elements in the United States

Selfish genetic elements (SGEs) are DNA sequences that are transmitted to viable offspring in greater than Mendelian frequencies. Medea SGEs occur naturally in some populations of red flour beetle (Tribolium castaneum) and are expected to increase in frequency within populations and spread among populations. The large‐scale U.S. distributions of Medea‐4 (M⁴) had been mapped based on samples from 1993 to 1995. We sampled beetles in 2011–2014 and show that the distribution of M⁴ in the United States is dynamic and has shifted southward. By using a genetic marker of Medea‐1 (M¹), we found five unique geographic clusters with high and low M¹ frequencies in a pattern not predicted by microsatellite‐based analysis of population structure. Our results indicate the absence of rigid barriers to Medea spread in the United States, so assessment of what factors have limited its current distribution requires further investigation. There is great interest in using synthetic SGEs, including synthetic Medea, to alter or suppress pest populations, but there is concern about unpredicted spread of these SGEs and potential for populations to become resistant to them. The finding of patchy distributions of Medea elements suggests that released synthetic SGEs cannot always be expected to spread uniformly, especially in target species with limited dispersal.     

Endosymbiotic Wolbachia bacteria as biological control tools of disease vectors and pests

Wolbachia bacteria are common cytoplasmic symbionts of insects, mites and filarial nematodes. They can alter the reproduction of their hosts. The symbiont could be eliminated, transferred or used through genetic alteration to take advantage or remove their possible influences on pests and/or natural enemies. Their extensive effects on reproduction and host fitness have made Wolbachia the subject of growing attention as a potential biocontrol agent. Here, we summarize the relations of Wolbachia in the control of disease vectors and pests. Furthermore, the drawbacks of these bacteria are also discussed.     

Site-specific selfish genes as tools for the control and genetic engineering of natural populations

Site-specific selfish genes exploit host functions to copy themselves into a defined target DNA sequence, and include homing endonuclease genes, group II introns and some LINE-like transposable elements. If such genes can be engineered to target new host sequences, then they can be used to manipulate natural populations, even if the number of individuals released is a small fraction of the entire population. For example, a genetic load sufficient to eradicate a population can be imposed in fewer than 20 generations, if the target is an essential host gene, the knockout is recessive and the selfish gene has an appropriate promoter. There will be selection for resistance, but several strategies are available for reducing the likelihood of it evolving. These genes may also be used to genetically engineer natural populations, by means of population-wide gene knockouts, gene replacements and genetic transformations. By targeting sex-linked loci just prior to meiosis one may skew the population sex ratio, and by changing the promoter one may limit the spread of the gene to neighbouring populations. The proposed constructs are evolutionarily stable in the face of the mutations most likely to arise during their spread, and strategies are also available for reversing the manipulations.     

Estimating quantitative genetic parameters in wild populations: a comparison of pedigree and genomic approaches

The estimation of quantitative genetic parameters in wild populations is generally limited by the accuracy and completeness of the available pedigree information. Using relatedness at genomewide markers can potentially remove this limitation and lead to less biased and more precise estimates. We estimated heritability, maternal genetic effects and genetic correlations for body size traits in an unmanaged long‐term study population of Soay sheep on St Kilda using three increasingly complete and accurate estimates of relatedness: (i) Pedigree 1, using observation‐derived maternal links and microsatellite‐derived paternal links; (ii) Pedigree 2, using SNP‐derived assignment of both maternity and paternity; and (iii) whole‐genome relatedness at 37 037 autosomal SNPs. In initial analyses, heritability estimates were strikingly similar for all three methods, while standard errors were systematically lower in analyses based on Pedigree 2 and genomic relatedness. Genetic correlations were generally strong, differed little between the three estimates of relatedness and the standard errors declined only very slightly with improved relatedness information. When partitioning maternal effects into separate genetic and environmental components, maternal genetic effects found in juvenile traits increased substantially across the three relatedness estimates. Heritability declined compared to parallel models where only a maternal environment effect was fitted, suggesting that maternal genetic effects are confounded with direct genetic effects and that more accurate estimates of relatedness were better able to separate maternal genetic effects from direct genetic effects. We found that the heritability captured by SNP markers asymptoted at about half the SNPs available, suggesting that denser marker panels are not necessarily required for precise and unbiased heritability estimates. Finally, we present guidelines for the use of genomic relatedness in future quantitative genetics studies in natural populations.     

Alternative male life-history tactics as potential vehicles for speeding introgression of farm salmon traits into wild populations

Releases of cultured organisms, such as farm Atlantic salmon (Salmo salar L.), threaten native biodiversity and the integrity of natural communities. Salmon escaping from sea farms, however, have relatively poor reproductive success, suggesting that the rate of spread of domesticated traits may be reduced. We now compare the relative reproductive success of males that mature precociously in freshwater (parr) and find that those of farm origin have higher breeding and fertilization success than wild and hybrid individuals. Specifically, hybrid parr had 57% and wild parr 25% the success of farm parr. Early maturing males could thus be important vehicles promoting introgression of domesticated and/or non‐native traits into wild populations and ultimately have long‐term impact on the genetic integrity of native populations.     

Inverse Medea as a novel gene drive system for local population replacement: a theoretical analysis

One strategy to control mosquito-borne diseases, such as malaria and dengue fever, on a regional scale is to use gene drive systems to spread disease-refractory genes into wild mosquito populations. The development of a synthetic Medea element that has been shown to drive population replacement in laboratory Drosophila populations has provided encouragement for this strategy but has also been greeted with caution over the concern that transgenes may spread into countries without their consent. Here, we propose a novel gene drive system, inverse Medea, which is strong enough to bring about local population replacement but is unable to establish itself beyond an isolated release site. The system consists of 2 genetic components--a zygotic toxin and maternal antidote--which render heterozygous offspring of wild-type mothers unviable. Through population genetic analysis, we show that inverse Medea will only spread when it represents a majority of the alleles in a population. The element is best located on an autosome and will spread to fixation provided any associated fitness costs are dominant and to very high frequency otherwise. We suggest molecular tools that could be used to build the inverse Medea system and discuss its utility for a confined release of transgenic mosquitoes.     

Effectiveness of zooprophylaxis in malaria control: a theoretical inquiry, with a model for mosquito populations with two bloodmeal hosts

A model for a vector mosquito population with two bloodmeal hosts (man and a domestic animal) was developed to study the influences of domestic animals on the frequency of mosquito bites on man and the endemicity of human malaria. The vector population model, including blood-feeding success in the adult stage (depending on host density and biting efficiency) and density-dependent regulation in the larval stage, was combined with the Ross-Macdonald malaria transmission model. Model analyses suggested that introduction of domestic animals easily fed upon by mosquitoes increases mosquito density and, in some situations, frequency of mosquito bites on man and the infection rate of malaria through increased success of blood-feeding. Extinction of malaria was predicted only when an extremely large number of easily accessible (as compared to man) domestic animals are introduced. Limitations in the concept of zooprophylaxis and problems of livestock management in malaria control are discussed.     

Mass‐dependent reproductive strategies in wild bighorn ewes: a quantitative genetic approach

In the Ram Mountain bighorn sheep (Ovis canadensis) population, ewes differing by more than 30% in body mass weaned lambs with an average mass difference of only 3%. Variability in adult body mass was partly due to additive genetic effects, but inheritance of weaning mass was weak. Maternal effects could obscure genetic effects in the phenotypic expression of weaning mass, particularly if they reflected strategies of maternal expenditure that varied according to ewe mass. We performed a quantitative genetic analysis to assess genetic and environmental influences on ewe mass and on maternal expenditure. We used the mean daughters/mother regression method and Derivative Free Restricted Maximum Likelihood models to estimate heritability (h²) of ewe mass and indices of maternal expenditure. We found additive genetic effects on phenotypic variation in maternal mass, in lamb mass at weaning (absolute maternal expenditure) and in weaning mass relative to maternal mass at weaning (relative maternal expenditure). Heritability suggests that maternal expenditure has the potential to evolve. The genetic correlation of ewe mass and absolute maternal expenditure was weak, while ewe mass and relative maternal expenditure were strongly negatively correlated. These results suggest additive genetic effects on mass‐dependent reproductive strategies in bighorn ewes. Mass‐dependent reproductive strategies could affect lamb survival and phenotypic variation in adult mass. As population density increased and reproduction became costlier, small females reduced maternal expenditure more than large females. Constraints on reproductive strategy imposed by variations in resource availability are therefore likely to differ according to ewe mass. A general trend for a decrease in maternal expenditure relative to maternal size in mammals suggests that size‐dependent negative maternal effects may be common.     

Replicated analysis of the genetic architecture of quantitative traits in two wild great tit populations

Currently, there is much debate on the genetic architecture of quantitative traits in wild populations. Is trait variation influenced by many genes of small effect or by a few genes of major effect? Where is additive genetic variation located in the genome? Do the same loci cause similar phenotypic variation in different populations? Great tits (Parus major) have been studied extensively in long‐term studies across Europe and consequently are considered an ecological ‘model organism’. Recently, genomic resources have been developed for the great tit, including a custom SNP chip and genetic linkage map. In this study, we used a suite of approaches to investigate the genetic architecture of eight quantitative traits in two long‐term study populations of great tits—one in the Netherlands and the other in the United Kingdom. Overall, we found little evidence for the presence of genes of large effects in either population. Instead, traits appeared to be influenced by many genes of small effect, with conservative estimates of the number of contributing loci ranging from 31 to 310. Despite concordance between population‐specific heritabilities, we found no evidence for the presence of loci having similar effects in both populations. While population‐specific genetic architectures are possible, an undetected shared architecture cannot be rejected because of limited power to map loci of small and moderate effects. This study is one of few examples of genetic architecture analysis in replicated wild populations and highlights some of the challenges and limitations researchers will face when attempting similar molecular quantitative genetic studies in free‐living populations.     

Bayesian inference of genetic parameters for ultrasound scanning traits of Kivircik lambs

Ultrasound scanning traits have been adapted in selection programs in many countries to improve carcass traits for lean meat production. As the genetic parameters of the traits interested are important for breeding programs, the estimation of these parameters was aimed at the present investigation. The estimated parameters were direct and maternal heritability as well as genetic correlations between the studied traits. The traits were backfat thickness (BFT), skin+backfat thickness (SBFT), eye muscle depth (MD) and live weights at the day of scanning (LW). The breed investigated was Kivircik, which has a high quality of meat. Six different multi-trait animal models were fitted to determine the most suitable model for the data using Bayesian approach. Based on deviance information criterion, a model that includes direct additive genetic effects, maternal additive genetic effects, direct maternal genetic covariance and maternal permanent environmental effects revealed to be the most appropriate for the data, and therefore, inferences were built on the results of that model. The direct heritability estimates for BFT, SBFT, MD and LW were 0.26, 0.26, 0.23 and 0.09, whereas the maternal heritability estimates were 0.27, 0.27, 0.24 and 0.20, respectively. Negative genetic correlations were obtained between direct and maternal effects for BFT, SBFT and MD. Both direct and maternal genetic correlations between traits were favorable, whereas BFT–MD and SBFT–MD had negligible direct genetic correlation. The highest direct and maternal genetic correlations were between BFT and SBFT (0.39) and between MD and LW (0.48), respectively. Our results, in general, indicated that maternal effects should be accounted for in estimation of genetic parameters of ultrasound scanning traits in Kivircik lambs, and SBFT can be used as a selection criterion to improve BFT.     

母体遗传效应对青海细毛羊生产性能遗传参数估计的影响

为了研究母体遗传效应对青海细毛羊生长性状、产毛性状的影响,文章采用平均信息最大约束似然法应用不同混合动物模型估计青海细毛羊生产性状的遗传参数,并采用似然比检验对不同模型进行比较分析。各模型中均包括固定效应、个体直接加性遗传效应、残差效应;随机效应为:个体永久环境效应、母体遗传效应、母体永久环境效应。不同模型对随机效应作了不同考虑:模型1不考虑个体永久环境效应、母体遗传效应、母体永久环境效应;模型2考虑母体永久环境效应;模型3考虑母体遗传效应;模型4考虑母体遗传效应和母体永久环境效应;模型5考虑个体永久环境效应和母体遗传效应;模型6考虑个体永久环境效应、母体遗传效应、母体永久环境效应。各模型估计的初生重遗传力为:0.1896~0.3781;断奶重遗传力为:0.2537~0.2890;周岁重遗传力范围:0.2244~0.3225;成年羊体重遗传力范围:0.2205~0.3983;产毛量遗传力为:0.1218~0.1490;羊毛细度遗传力为:0.0983~0.4802;羊毛长度遗传力为:0.1170~0.1311。与模型1相比,模型3对于初生重、断奶重差异显著(P<0.01),对于周岁重、成年羊体重各模型与模型1的似然比检验差异不显著(P>0.05);与模型6相比,模型4、5对于羊毛细度差异显著(P<0.01),模型4对羊毛长度差异显著(P<0.05),对于产毛量各模型与模型6似然比检验差异不显著(P>0.05)。生长性状中初生重、断奶重受母体遗传效应影响显著,周岁重、成年羊体重受母体遗传效应影响不显著;产毛性状中羊毛细度、长度受母体遗传效应影响显著,产毛量受母体遗传效应影响较弱。     

Environmental DNA analysis shows high potential as a tool for estimating intraspecific genetic diversity in a wild fish population

Environmental DNA (eDNA) analysis has recently been used as a new tool for estimating intraspecific diversity. However, whether known haplotypes contained in a sample can be detected correctly using eDNA‐based methods has been examined only by an aquarium experiment. Here, we tested whether the haplotypes of Ayu fish (Plecoglossus altivelis altivelis) detected in a capture survey could also be detected from an eDNA sample derived from the field that contained various haplotypes with low concentrations and foreign substances. A water sample and Ayu specimens collected from a river on the same day were analysed by eDNA analysis and Sanger sequencing, respectively. The 10 L water sample was divided into 20 filters for each of which 15 PCR replications were performed. After high‐throughput sequencing, denoising was performed using two of the most widely used denoising packages, unoise3 and dada2 . Of the 42 haplotypes obtained from the Sanger sequencing of 96 specimens, 38 (unoise3 ) and 41 (dada2 ) haplotypes were detected by eDNA analysis. When dada2 was used, except for one haplotype, haplotypes owned by at least two specimens were detected from all the filter replications. Accordingly, although it is important to note that eDNA‐based method has some limitations and some risk of false positive and false negative, this study showed that the eDNA analysis for evaluating intraspecific genetic diversity provides comparable results for large‐scale capture‐based conventional methods. Our results suggest that eDNA‐based methods could become a more efficient survey method for investigating intraspecific genetic diversity in the field.     

附子野生资源群体遗传多样性的RAPD分析

应用RAPD标记分析了分布在附子主栽区四川、重庆、陕西及湖北16个野生乌头种群的遗传变异。24个引物共检测到643个RAPD位点,多态位点602个,总的多态位点百分率达93．5％。Shannon多样性和Nei遗传分化结果一致显示重庆酉阳种群和重庆城口种群遗传多样性最高,四川盐源种群和陕西勉县种群的遗传多样性最低。Shannon指数测出的种群内的遗传变异（57．6％）略占优势,群体间的遗传分化达到42．4％;Nei基因分化系数（GST）达40．0％;分子方差分析（AMOVA）发现群体间遗传变异仅为25．37％;种群每代迁移数Nm为0．756。Nei相似性系数的UPGMA分析结果显示该地区的野生乌头分布上有一定的地域性,特别是同为附子道地产地江油供种的北川、安县和青川种群间遗传关系密切,说明种质资源在道地药材形成中具有重要作用。研究结果表明,附子主要栽培地区的乌头野生种群之间存在较大的遗传分化,遗传多样性较高,遗传种质资源较丰富,存在一定的特异性资源,为进一步开发利用乌头（川乌、附子）提供了丰富的种质资源。     

母体遗传效应对青海细毛羊生产性能遗传参数估计的影响

为了研究母体遗传效应对青海细毛羊生长性状、产毛性状的影响, 文章采用平均信息最大约束似然法应用不同混合动物模型估计青海细毛羊生产性状的遗传参数, 并采用似然比检验对不同模型进行比较分析。各模型中均包括固定效应、个体直接加性遗传效应、残差效应; 随机效应为：个体永久环境效应、母体遗传效应、母体永久环境效应。不同模型对随机效应作了不同考虑：模型1不考虑个体永久环境效应、母体遗传效应、母体永久环境效应; 模型2考虑母体永久环境效应; 模型3考虑母体遗传效应; 模型4考虑母体遗传效应和母体永久环境效应; 模型5考虑个体永久环境效应和母体遗传效应; 模型6考虑个体永久环境效应、母体遗传效应、母体永久环境效应。各模型估计的初生重遗传力为：0.1896~0.3781; 断奶重遗传力为：0.2537~0.2890; 周岁重遗传力范围：0.2244~0.3225; 成年羊体重遗传力范围：0.2205~0.3983; 产毛量遗传力为：0.1218~0.1490; 羊毛细度遗传力为：0.0983~0.4802; 羊毛长度遗传力为：0.1170~0.1311。与模型1相比, 模型3对于初生重、断奶重差异显著(P<0.01), 对于周岁重、成年羊体重各模型与模型1的似然比检验差异不显著(P>0.05); 与模型6相比, 模型4、5对于羊毛细度差异显著(P<0.01), 模型4对羊毛长度差异显著(P<0.05), 对于产毛量各模型与模型6似然比检验差异不显著(P>0.05)。生长性状中初生重、断奶重受母体遗传效应影响显著, 周岁重、成年羊体重受母体遗传效应影响不显著; 产毛性状中羊毛细度、长度受母体遗传效应影响显著, 产毛量受母体遗传效应影响较弱。     

Fine-scale population genetic structure of a wildlife disease vector: the southern house mosquito on the island of Hawaii

The southern house mosquito, Culex quinquefasciatus, is a widespread tropical and subtropical disease vector. In the Hawaiian Islands, where it was introduced accidentally almost two centuries ago, it is considered the primary vector of avian malaria and pox. Avian malaria in particular has contributed to the extinction and endangerment of Hawaii's native avifauna, and has altered the altitudinal distribution of native bird populations. We examined the population genetic structure of Cx. quinquefasciatus on the island of Hawaii at a smaller spatial scale than has previously been attempted, with particular emphasis on the effects of elevation on population genetic structure. We found significant genetic differentiation among populations and patterns of isolation by distance within the island. Elevation per se did not have a limiting effect on gene flow; however, there was significantly lower genetic diversity among populations at mid elevations compared to those at low elevations. A recent sample taken from just above the predicted upper altitudinal distribution of Cx. quinquefasciatus on the island of Hawaii was confirmed as being a temporary summer population and appeared to consist of individuals from more than one source population. Our results indicate effects of elevation gradients on genetic structure that are consistent with known effects of elevation on population dynamics of this disease vector.     

Genomic P elements and P-M characteristics of eastern Australian populations of Drosophila melanogaster

As part of our effort to monitor changes in the clinal pattern of P element-associated traits in eastern Australian Drosophila melanogaster, we investigated the genomic P elements of 293 isofemale lines collected in the period 1991–1994 from 45 localities. P elements were present in many copies in all genomes examined, with full-size P and KP element size classes accounting for the large majority. SR elements were not present in at least 92% of the lines tested. South of about 26° south Latitude (°SLat), the ratio of KP to full-size P elements (KP/P ratio) increased, correlating weakly with the P-M phenotypes of the populations, from moderately P populations (26–29°SLat) to M populations (37–38°SLat) North of 26°SLat, in weak P populations, the KP/P ratio was higher than between 26 and 29°Slat. The KP/P ratio appears to be higher in the northern populations than it was when previous studies were done. Overall, a high KP/P ratio among lines correlated roughly with a lack of P activity, but it also correlated with reduced repressor function. In a sample of 30 lines, a maternal effect of repressor function did not show a pattern with latitude, nor with KP/P ratio, nor with presence or absence of P activity. This revised version was published online in July 2006 with corrections to the Cover Date.     

Differentiation under isolation and genetic structure of Sardinian hares as revealed by craniometric analysis,mitochondrial DNA and microsatellites

Hares (Lepus capensis Linnaeus 1758) were probably introduced into Sardinia in historical times. Previous studies indicated North Africa as the most likely source area but did not exclude the occurrence of hybridization events with continental brown hares (L. europaeus Pallas 1778) perhaps introduced for hunting purposes. We implemented both morphometric and genetic approaches to verify the genetic isolation of the Sardinian population. Specifically, we conducted a multivariate analysis of craniometric data and analysed 461 bp of the mitochondrial control region and 12 autosomal microsatellites in Sardinian hares, using North African cape hares and European brown hares as reference populations. Sardinian hares displayed a peculiar skull shape. In agreement, both nuclear and mitochondrial markers remarked the distinctiveness of this population. Observed and expected heterozygosity were 0.52 and 0.61, while haplotype and nucleotide diversity were 0.822 and 0.0129. Self‐assignment based on Bayesian cluster analysis was high (average membership 0.98), and no evident signs of introgression from continental brown hares were found. Our results support the hypothesis that the Sardinian hares have been introduced from North Africa, remained genetically isolated since the founding event and evolved independently from the source population. This long‐lasting isolation and the consequent genetic drift resulted in a differentiation, perhaps accompanied by an adaptation to local environmental conditions.     

Maintaining genetic integrity of coexisting wild and domestic populations: Genetic differentiation between wild and domestic Rangifer with long traditions of intentional interbreeding

This study investigates the genetic effect of an indigenous tradition of deliberate and controlled interbreeding between wild and domestic Rangifer . The results are interpreted in the context of conservation concerns and debates on the origin of domestic animals. The study is located in Northeastern Zaba?kal'e, Russia at approximately 57 degrees North latitude. Blood and skin samples, collected from wild and domestic Rangifer , are analyzed for their mtDNA and microsatellite signatures. Local husbandry traditions are documented ethnographically. The genetic data are analyzed with special reference to indigenous understandings of the distinctions between local domestic types and wild Rangifer . The genetic results demonstrate a strong differentiation between wild and domestic populations. Notably low levels of mtDNA haplotype sharing between wild and domestic reindeer, suggest mainly male‐mediated gene flow between the two gene pools. The nuclear microsatellite results also point to distinct differences between regional domestic clusters. Our results indicate that the Evenki herders have an effective breeding technique which, while mixing pedigrees in the short term, guards against wholesale introgression between wild and domestic populations over the long term. They support a model of domestication where wild males and domestic females are selectively interbred, without hybridizing the two populations. Our conclusions inform a debate on the origins of domestication by documenting a situation where both wild and domestic types are in constant interaction. The study further informs a debate in conservation biology by demonstrating that certain types of controlled introgression between wild and domestic types need not reduce genetic diversity.