Next‐generation transgenic cotton: pyramiding RNAi and Bt counters insect resistance

Transgenic crops producing insecticidal proteins from the bacterium Bacillus thuringiensis (Bt) are extensively cultivated worldwide. To counter rapidly increasing pest resistance to crops that produce single Bt toxins, transgenic plant ‘pyramids’ producing two or more Bt toxins that kill the same pest have been widely adopted. However, cross‐resistance and antagonism between Bt toxins limit the sustainability of this approach. Here we describe development and testing of the first pyramids of cotton combining protection from a Bt toxin and RNA interference (RNAi). We developed two types of transgenic cotton plants producing double‐stranded RNA (dsRNA) from the global lepidopteran pest Helicoverpa armigera designed to interfere with its metabolism of juvenile hormone (JH). We focused on suppression of JH acid methyltransferase (JHAMT), which is crucial for JH synthesis, and JH‐binding protein (JHBP), which transports JH to organs. In 2015 and 2016, we tested larvae from a Bt‐resistant strain and a related susceptible strain of H. armigera on seven types of cotton: two controls, Bt cotton, two types of RNAi cotton (targeting JHAMT or JHBP) and two pyramids (Bt cotton plus each type of RNAi). Both types of RNAi cotton were effective against Bt‐resistant insects. Bt cotton and RNAi acted independently against the susceptible strain. In computer simulations of conditions in northern China, where millions of farmers grow Bt cotton as well as abundant non‐transgenic host plants of H. armigera, pyramided cotton combining a Bt toxin and RNAi substantially delayed resistance relative to using Bt cotton alone.     

Predation of cotton bollworm by green lacewings in the presence of cotton aphid as alternative prey on transgenic Bt and conventional cotton

Experiments were conducted in small arenas and on whole plants to explore the effect of cotton aphids, Aphis gossypii Glover (Hemiptera: Aphididae), as alternative prey on the predation of Helicoverpa armigera Hübner (Lepidoptera: Noctuidae) larvae by green lacewing larvae, Mallada signatus Schneider (Neuroptera: Chrysopidae). Transgenic Bt (Bollgard II^®) and conventional cotton plants were included to explore potential differences in the predator's performance on these cotton types. In small arenas, the presence of 20 aphids reduced predation on H. armigera larvae by 22% (from 5.5 to 3.3 of 10) by a single lacewing larva over a 24‐h period. The presence of H. armigera reduced predation on aphids by ca. 29% (from 16.8 to 11.0 of 20) over 24 h. On whole plants, the presence of alternative prey had no effect on the number of H. armigera larvae or aphids remaining after 3 days. The presence of H. armigera larvae alone, without the predator, caused a 24% reduction in the numbers of aphids on conventional, but not on Bt cotton plants. The combination of Bt cotton and lacewing larvae caused a 96.6% removal of early‐stage H. armigera larvae, a statistically significant increase over the addition of the proportions (91.6%) removed by each factor measured separately, providing evidence of synergism. These studies suggest that the presence of aphids as alternative prey would not necessarily disrupt the predation by green lacewing on larvae of H. armigera, especially on Bt cotton.     

High degree of genetic differentiation in marine three‐spined sticklebacks (Gasterosteus aculeatus)

Populations of widespread marine organisms are typically characterized by a low degree of genetic differentiation in neutral genetic markers, but much less is known about differentiation in genes whose functional roles are associated with specific selection regimes. To uncover possible adaptive population divergence and heterogeneous genomic differentiation in marine three‐spined sticklebacks (Gasterosteus aculeatus), we used a candidate gene‐based genome‐scan approach to analyse variability in 138 microsatellite loci located within/close to (<6 kb) functionally important genes in samples collected from ten geographic locations. The degree of genetic differentiation in markers classified as neutral or under balancing selection—as determined with several outlier detection methods—was low (F_ST = 0.033 or 0.011, respectively), whereas average F_ST for directionally selected markers was significantly higher (F_ST = 0.097). Clustering analyses provided support for genomic and geographic heterogeneity in selection: six genetic clusters were identified based on allele frequency differences in the directionally selected loci, whereas four were identified with the neutral loci. Allelic variation in several loci exhibited significant associations with environmental variables, supporting the conjecture that temperature and salinity, but not optic conditions, are important drivers of adaptive divergence among populations. In general, these results suggest that in spite of the high degree of physical connectivity and gene flow as inferred from neutral marker genes, marine stickleback populations are strongly genetically structured in loci associated with functionally relevant genes.     

Fitness of Bt‐resistant cabbage loopers on Bt cotton plants

Development of resistance to the insecticidal toxins from Bacillus thuringiensis (Bt) in insects is the major threat to the continued success of transgenic Bt crops in agriculture. The fitness of Bt‐resistant insects on Bt and non‐Bt plants is a key parameter that determines the development of Bt resistance in insect populations. In this study, a comprehensive analysis of the fitness of Bt‐resistant Trichoplusia ni strains on Bt cotton leaves was conducted. The Bt‐resistant T. ni strains carried two genetically independent mechanisms of resistance to Bt toxins Cry1Ac and Cry2Ab. The effects of the two resistance mechanisms, individually and in combination, on the fitness of the T. ni strains on conventional non‐Bt cotton and on transgenic Bt cotton leaves expressing a single‐toxin Cry1Ac (Bollgard I) or two Bt toxins Cry1Ac and Cry2Ab (Bollgard II) were examined. The presence of Bt toxins in plants reduced the fitness of resistant insects, indicated by decreased net reproductive rate (R₀) and intrinsic rate of increase (r). The reduction in fitness in resistant T. ni on Bollgard II leaves was greater than that on Bollgard I leaves. A 12.4‐day asynchrony of adult emergence between the susceptible T. ni grown on non‐Bt cotton leaves and the dual‐toxin‐resistant T. ni on Bollgard II leaves was observed. Therefore, multitoxin Bt plants not only reduce the probability for T. ni to develop resistance but also strongly reduce the fitness of resistant insects feeding on the plants.     

Bacillus thuringiensis plants expressing Cry1Ac,Cry2Ab and Cry1F are not toxic to the assassin bug,Zelus renardii

Cotton‐ and maize‐producing insecticidal crystal (Cry) proteins from the bacterium, Bacillus thuringiensis (Bt), have been commercialized since 1996. Bt plants are subjected to environmental risk assessments for non‐target organisms, including natural enemies that suppress pest populations. Here, we used Cry1F‐resistant Spodoptera frugiperda (J.E. Smith) and Cry1Ac and Cry2Ab‐resistant Trichoplusia ni (Hübner) as prey for the assassin bug, Zelus renardii (Kolenati), a common predator in maize and cotton fields. In tritrophic studies, we assessed several fitness parameters of Z. renardii when it fed on resistant S. frugiperda that had fed on Bt maize expressing Cry1F or on resistant T. ni that had fed on Bt cotton expressing Cry1Ac and Cry2Ab. Survival, nymphal duration, adult weight, adult longevity and female fecundity of Z. renardii were not different when they were fed resistant‐prey larvae (S. frugiperda or T. ni) reared on either a Bt crop or respective non‐Bt crops. ELISA tests demonstrated that the Cry proteins were present in the plant at the highest levels, at lower levels in the prey and at the lowest levels in the predator. While Z. renardii was exposed to Cry1F and Cry1Ac and Cry2Ab when it fed on hosts that consumed Bt‐transgenic plants, the proteins did not affect important fitness parameters in this common and important predator.     

Consumptive and non‐consumptive effects of wolf spiders on cotton bollworms

Larvae of the cotton bollworm, Helicoverpa armigera (Hübner) (Lepidoptera: Noctuidae) that survive on genetically modified Bt cotton (Gossypium hirsutum L., Malvaceae) contribute to the risk of widespread resistance to Bt toxins. Current resistance management techniques include pupae busting, which involves deep tilling of the soil to kill overwintering pupae. Unfortunately, pupae busting runs counter to soil and water conserving techniques, such as minimum tillage. This problem could be relieved with biological control methods, whereby predators attack either larvae going to ground to pupate or moths emerging from the ground. We found that the wolf spider Tasmanicosa leuckartii (Thorell) (Araneae: Lycosidae), a common inhabitant of Australian cotton agroecosystems, is an effective predator of H. armigera, attacking and killing most larvae (66%) and emerging moths (77%) in simple laboratory arenas. Tasmanicosa leuckartii also reduced the number of emerging moths by 66% on average in more structurally complex glasshouse arenas. Males, females, and late‐instar juveniles of T. leuckartii were similarly effective. Tasmanicosa leuckartii also imposed non‐consumptive effects on H. armigera, as when a spider was present larvae in the laboratory areas spent less time on the cotton boll and more time on the soil and more mass was lost from the cotton boll. Increased loss of boll mass likely reflects changes in H. armigera foraging behavior induced by the presence of spiders (indirect non‐consumptive effects). Helicoverpa armigera spent more time as pupae when the spider was present in simple laboratory arenas, but not in more complex glasshouse enclosures. Overall, results indicate that T. leuckartii spiders can be effective predators of H. armigera late instars and moths but also suggest that, under some conditions, the presence of spiders could increase the damage to individual cotton bolls.     

Genome scans for the contemporary response to selection in quantitative traits

Genome scans have been an important approach for discovering historical signatures of selection in both model and nonmodel species. An exciting new experimental design for genome scans is to measure the change in allele frequency before and after contemporary selection within a generation, from a single population. The most widely‐used methods, however, have two major limitations: they are based on testing one locus at a time, and they only have power to uncover loci that have evolved under relatively strong selection. On the other hand, complex quantitative traits are common in nature and are caused by several loci of small effect. Selection on a quantitative trait at the phenotypic level is predicted to be accompanied by subtle allele frequency changes in many loci that covary (a polygenic soft sweep), rather than a large, single‐effect allele (a selective sweep). In this issue of Molecular Ecology, Bourret et al. (2014) measure the contemporary response to natural selection across the genome in multiple cohorts of Atlantic salmon during their first year at sea. They introduce a multilocus framework based on groups of markers that covary in their genotypic distribution. While the traditional, single‐locus approach did not find evidence for repeated patterns of selection, the multivariate approach found that a group of covarying SNPs was selected for in different cohorts at one site. Their multilocus framework has potential to be a more fruitful approach for uncovering the genomic basis of adaptation in quantitative traits, although caution should be applied as the framework has yet to be validated with simulated data.     

Genome‐wide analyses suggest parallel selection for universal traits may eclipse local environmental selection in a highly mobile carnivore

Ecological and environmental heterogeneity can produce genetic differentiation in highly mobile species. Accordingly, local adaptation may be expected across comparatively short distances in the presence of marked environmental gradients. Within the European continent, wolves (Canis lupus) exhibit distinct north–south population differentiation. We investigated more than 67‐K single nucleotide polymorphism (SNP) loci for signatures of local adaptation in 59 unrelated wolves from four previously identified population clusters (northcentral Europe n = 32, Carpathian Mountains n = 7, Dinaric‐Balkan n = 9, Ukrainian Steppe n = 11). Our analyses combined identification of outlier loci with findings from genome‐wide association study of individual genomic profiles and 12 environmental variables. We identified 353 candidate SNP loci. We examined the SNP position and neighboring megabase (1 Mb, one million bases) regions in the dog (C. lupus familiaris) genome for genes potentially under selection, including homologue genes in other vertebrates. These regions included functional genes for, for example, temperature regulation that may indicate local adaptation and genes controlling for functions universally important for wolves, including olfaction, hearing, vision, and cognitive functions. We also observed strong outliers not associated with any of the investigated variables, which could suggest selective pressures associated with other unmeasured environmental variables and/or demographic factors. These patterns are further supported by the examination of spatial distributions of the SNPs associated with universally important traits, which typically show marked differences in allele frequencies among population clusters. Accordingly, parallel selection for features important to all wolves may eclipse local environmental selection and implies long‐term separation among population clusters.     

Survival of Helicoverpa armigera larvae on and Bt toxin expression in various parts of transgenic Bt cotton (Bollgard II) plants

There is no conclusive evidence that Helicoverpa spp. (Lepidoptera: Noctuidae) in Australia have evolved significant levels of resistance to Bollgard II^® cotton (which expresses two Bt toxin genes, cry1Ac and cry2Ab). However, there is evidence of surviving larvae on Bollgard II cotton in the field. The distribution and survival of early‐instar Helicoverpa armigera (Hübner) (Lepidoptera: Noctuidae) larvae were examined on whole Bollgard II and non‐Bt cotton plants in greenhouse bioassays. The expression of Cry toxins in various parts of Bollgard II plants was compared to the survival of larvae in those locations. Only 1% of larvae survived after 6 days on greenhouse‐grown Bollgard II plants compared to 31% on non‐Bt cotton plants. Overall, and across all time intervals, more larvae survived on reproductive parts (squares, flowers, and bolls) than on vegetative parts (leaves, stems, and petioles) on Bollgard II plants. The concentration of Cry1Ac toxin did not differ between plant structures, whereas Cry2Ab toxin differed significantly, but there was no relationship between the level of expression and the location of larvae. This study provides no evidence that lower expression of Cry toxins in the reproductive parts of plants explains the survival of H. armigera larvae on Bollgard II cotton.     

Communication interference in sympatrically occurring moth species

In moth species, females emit a species‐specific sex pheromone that is perceived over long distance by conspecific males. The species‐specificity in the chemical communication channel is achieved by a combination of unique components in specific ratios and sometimes also by interspecific behavioural antagonists to deter sympatrically occurring heterospecific males. In this study, we determined possible antagonistic effects in Helicoverpa gelotopoeon Dyar (Lepidoptera: Noctuidae) males to the major sex pheromone component of sympatrically occurring heliothine moths, Z11‐16:Ald, as well as to the sex pheromone of the sympatrically occurring Heliothis virescens (Fabricius) (Lepidoptera: Noctuidae) (Z11‐16:Ald and Z9‐14:Ald). We also explored whether other co‐occurring species are attracted to these pheromone blends. Our field experiments showed that the addition of Z11‐16:Ald alone or in combination with Z9‐14:Ald inhibited trap catches of H. gelotopoeon males and that this inhibition depended on the concentration of these compounds. In addition, other moth species were attracted to the blends. Together, our results confirm the antagonistic effect of heterospecific sex pheromone compounds of H. virescens to H. gelotopoeon.     

Haplotype‐based genotyping‐by‐sequencing in oat genome research

In a de novo genotyping‐by‐sequencing (GBS) analysis of short, 64‐base tag‐level haplotypes in 4657 accessions of cultivated oat, we discovered 164741 tag‐level (TL) genetic variants containing 241224 SNPs. From this, the marker density of an oat consensus map was increased by the addition of more than 70000 loci. The mapped TL genotypes of a 635‐line diversity panel were used to infer chromosome‐level (CL) haplotype maps. These maps revealed differences in the number and size of haplotype blocks, as well as differences in haplotype diversity between chromosomes and subsets of the diversity panel. We then explored potential benefits of SNP vs. TL vs. CL GBS variants for mapping, high‐resolution genome analysis and genomic selection in oats. A combined genome‐wide association study (GWAS) of heading date from multiple locations using both TL haplotypes and individual SNP markers identified 184 significant associations. A comparative GWAS using TL haplotypes, CL haplotype blocks and their combinations demonstrated the superiority of using TL haplotype markers. Using a principal component‐based genome‐wide scan, genomic regions containing signatures of selection were identified. These regions may contain genes that are responsible for the local adaptation of oats to Northern American conditions. Genomic selection for heading date using TL haplotypes or SNP markers gave comparable and promising prediction accuracies of up to r = 0.74. Genomic selection carried out in an independent calibration and test population for heading date gave promising prediction accuracies that ranged between r = 0.42 and 0.67. In conclusion, TL haplotype GBS‐derived markers facilitate genome analysis and genomic selection in oat.     

Evidence for adaptation from standing genetic variation on an antimicrobial peptide gene in the mussel Mytilus edulis

Genome scans of population differentiation identify candidate loci for adaptation but provide little information on how selection has influenced the genetic structure of these loci. Following a genome scan, we investigated the nature of the selection responsible for the outlying differentiation observed between populations of the marine mussel Mytilus edulis at a leucine/arginine polymorphism (L31R) in the antimicrobial peptide MGD2. We analysed DNA sequence polymorphisms, allele frequencies and population differentiation of polymorphisms closely linked to L31R, and pairwise and third‐order linkage disequilibria. An outlying level of population differentiation was observed at L31R only, while no departure from panmixia was observed at linked loci surrounding L31R, as in most of the genome. Selection therefore seems to affect L31R directly. Three hypotheses can explain the lack of differentiation in the chromosomal region close to L31R: (i) hitchhiking has occurred but migration and recombination subsequently erased the signal, (ii) selection was weak enough and recombination strong enough to limit the hitchhiking effect to a very small chromosomal region or (iii) selection acted on a pre‐existing polymorphism (i.e. standing variation) at linkage equilibrium with its background. Linkage equilibrium was observed between L31R and linked polymorphisms in every population analysed, as expected under the three hypotheses. However, linkage disequilibrium was observed in some populations between pairs of loci located upstream and downstream to L31R, generating a complex pattern of third‐order linkage disequilibria which is best explained by the hypothesis of selection on a pre‐existing polymorphism. We hypothesise that selection could be either balanced, maintaining alleles at different frequencies depending on the pathogen community encountered locally by mussels, or intermittent, resulting in sporadic fluctuations in allele frequency.     

Development of a novel‐type transgenic cotton plant for control of cotton bollworm

The transgenic Bt cotton plant has been widely planted throughout the world for the control of cotton budworm Helicoverpa armigera (Hubner). However, a shift towards insect tolerance of Bt cotton is now apparent. In this study, the gene encoding neuropeptide F (NPF) was cloned from cotton budworm H. armigera, an important agricultural pest. The npf gene produces two splicing mRNA variants—npf1 and npf2 (with a 120‐bp segment inserted into the npf1 sequence). These are predicted to form the mature NPF1 and NPF2 peptides, and they were found to regulate feeding behaviour. Knock down of larval npf with dsNPF in vitro resulted in decreases of food consumption and body weight, and dsNPF also caused a decrease of glycogen and an increase of trehalose. Moreover, we produced transgenic tobacco plants transiently expressing dsNPF and transgenic cotton plants with stably expressed dsNPF. Results showed that H. armigera larvae fed on these transgenic plants or leaves had lower food consumption, body size and body weight compared to controls. These results indicate that NPF is important in the control of feeding of H. armigera and valuable for production of potential transgenic cotton.     

Water availability drives signatures of local adaptation in whitebark pine (Pinus albicaulis Engelm.) across fine spatial scales of the Lake Tahoe Basin,USA

Patterns of local adaptation at fine spatial scales are central to understanding how evolution proceeds, and are essential to the effective management of economically and ecologically important forest tree species. Here, we employ single and multilocus analyses of genetic data (n = 116 231 SNPs) to describe signatures of fine‐scale adaptation within eight whitebark pine (Pinus albicaulis Engelm.) populations across the local extent of the environmentally heterogeneous Lake Tahoe Basin, USA. We show that despite highly shared genetic variation (F_ST = 0.0069), there is strong evidence for adaptation to the rain shadow experienced across the eastern Sierra Nevada. Specifically, we build upon evidence from a common garden study and find that allele frequencies of loci associated with four phenotypes (mean = 236 SNPs), 18 environmental variables (mean = 99 SNPs), and those detected through genetic differentiation (n = 110 SNPs) exhibit significantly higher signals of selection (covariance of allele frequencies) than could be expected to arise, given the data. We also provide evidence that this covariance tracks environmental measures related to soil water availability through subtle allele frequency shifts across populations. Our results replicate empirical support for theoretical expectations of local adaptation for populations exhibiting strong gene flow and high selective pressures and suggest that ongoing adaptation of many P. albicaulis populations within the Lake Tahoe Basin will not be constrained by the lack of genetic variation. Even so, some populations exhibit low levels of heritability for the traits presumed to be related to fitness. These instances could be used to prioritize management to maintain adaptive potential. Overall, we suggest that established practices regarding whitebark pine conservation be maintained, with the additional context of fine‐scale adaptation.     

SNPs in SNCA,MCCC1, DLG2, GBF1 and MBNL2 are associated with Parkinson's disease in southern Chinese population

Numerous single nucleotide polymorphisms (SNPs), which have been identified as susceptibility factors for Parkinson's disease (PD) as per genome‐wide association studies, have not been fully characterized for PD patients in China. This study aimed to replicate the relationship between 12 novel SNPs of 12 genes and PD risk in southern Chinese population. Twelve SNPs of 12 genes were detected in 231 PD patients and 249 controls, using the SNaPshot technique. Meta‐analysis was used to assess heterogeneity of effect sizes between this study and published data. The impact of SNPs on gene expression was investigated by analysing the SNP‐gene association in the expression quantitative trait loci (eQTL) data sets. rs8180209 of SNCA (allele model: P = .047, OR = 0.77; additive model: P = .047, OR = 0.77), rs2270968 of MCCC1 (dominant model: P = .024, OR = 1.52), rs7479949 of DLG2 (recessive model; P = .019, OR = 1.52), rs10748818 of GBF1 (additive model: P < .001, OR = 0.37), and rs4771268 of MBNL2 (recessive model: P = .003, OR = 0.48) were replicated to be significantly associated with the increased risk of PD. Noteworthy, a meta‐analysis of previous studies suggested rs8180209, rs2270968, rs7479949 and rs4771268 were in line with those of our cohort. Our study replicated five novel functional SNPs in SNCA, MCCC1, DLG2, GBF1 and MBNL2 could be associated with increased risk of PD in southern Chinese population.     

Molecular diversity and landscape genomics of the crop wild relative Triticum urartu across the Fertile Crescent

Modern plant breeding can benefit from the allelic variation that exists in natural populations of crop wild relatives that evolved under natural selection in varying pedoclimatic conditions. In this study, next‐generation sequencing was used to generate 1.3 million genome‐wide single nucleotide polymorphisms (SNPs) on ex situ collections of Triticum urartu L., the wild donor of the A^u subgenome of modern wheat. A set of 75 511 high‐quality SNPs were retained to describe 298 T. urartu accessions collected throughout the Fertile Crescent. Triticum urartu showed a complex pattern of genetic diversity, with two main genetic groups distributed sequentially from west to east. The incorporation of geographical information on sampling points showed that genetic diversity was correlated to the geographical distance (R² = 0.19) separating samples from Jordan and Lebanon, from Syria and southern Turkey, and from eastern Turkey, Iran and Iraq. The wild emmer genome was used to derive the physical positions of SNPs on the seven chromosomes of the A^u subgenome, allowing us to describe a relatively slow decay of linkage disequilibrium in the collection. Outlier loci were described on the basis of the geographic distribution of the T. urartu accessions, identifying a hotspot of directional selection on chromosome 4A. Bioclimatic variation was derived from grid data and related to allelic variation using a genome‐wide association approach, identifying several marker–environment associations (MEAs). Fifty‐seven MEAs were associated with altitude and temperature measures while 358 were associated with rainfall measures. The most significant MEAs and outlier loci were used to identify genomic loci with adaptive potential (some already reported in wheat), including dormancy and frost resistance loci. We advocate the application of genomics and landscape genomics on ex situ collections of crop wild relatives to efficiently identify promising alleles and genetic materials for incorporation into modern crop breeding.     

Sublethal effects of Bt toxin and chlorpyrifos on various Spodoptera exigua populations

Bt cotton (Cry1Ac) has been commercially grown in China since 1997, saving China's cotton production from attack by Bt‐target pests and also tremendously reducing pesticide usage. In recent years, however, Bt cotton, with 4.2 million ha of cultivation, has suffered from a secondary target pest, Spodoptera exigua (Hübner) (Lepidoptera: Noctuidae). In China, growers have even had to re‐adopt conventional pesticides to control the pest, and this practice has already caused serious pesticide residue. In order to clarify the sublethal effects of chemical pesticide, the responses of a Bt‐susceptible and a Bt‐tolerant (Bt10) S. exigua strain to three treatment combinations were examined, including Bt toxin, sublethal chlorpyrifos, and Bt + sublethal chlorpyrifos. The susceptible and the Bt10 strain responded differently to dual pressure. Bt toxin + sublethal chlorpyrifos treatment lowered larval mortality and stimulated population increase of the susceptible S. exigua, whereas it delayed growth and development of the Bt10 strain. Under dual pressure, although larvae of the Bt10 strain developed faster than larvae of the susceptible strain, the Bt10 population experienced higher larval mortality, prolonged pupal duration, decreased pupal weight, decreased emergence rate, and shortened adult longevity. Compared with the susceptible strain, the Bt10 strain was deleteriously affected by sublethal chlorpyrifos. The Bt‐tolerant/resistant S. exigua population was more vulnerable to chemical pesticides like chlorpyrifos regardless of whether it was exposed to Bt toxin or not. Our study provides a reference for increasing the efficacy of control of S. exigua in Bt‐cotton planting areas.     

Population genomics reveals a fine‐scale recombination landscape for genetic improvement of cotton

Recombination breaks up ancestral linkage disequilibrium, creates combinations of alleles, affects the efficiency of natural selection, and plays a major role in crop domestication and improvement. However, there is little knowledge regarding the variation in the population‐scaled recombination rate in cotton. We constructed recombination maps and characterized the difference in the genomic landscape of the population‐scaled recombination rate between Gossypium hirsutum and G. arboreum and sub‐genomes based on the 381 sequenced G. hirsutum and 215 G. arboreum accessions. Comparative genomics identified large structural variations and syntenic genes in the recombination regions, suggesting that recombination was related to structural variation and occurred preferentially in the distal chromosomal regions. Correlation analysis indicated that recombination was only slightly affected by geographical distribution and breeding period. A genome‐wide association study (GWAS) was performed with 15 agronomic traits using 267 cotton accessions and identified 163 quantitative trait loci (QTL) and an important candidate gene (Ghir_COL2) for early maturity traits. Comparative analysis of recombination and a GWAS revealed that the QTL of fibre quality traits tended to be more common in high‐recombination regions than were those of yield and early maturity traits. These results provide insights into the population‐scaled recombination landscape, suggesting that recombination contributed to the domestication and improvement of cotton, which provides a useful reference for studying recombination in other species.     

Footprints of selection in wild populations of Bicyclus anynana along a latitudinal cline

One of the major questions in ecology and evolutionary biology is how variation in the genome enables species to adapt to divergent environments. Here, we study footprints of thermal selection in candidate genes in six wild populations of the afrotropical butterfly Bicyclus anynana sampled along a c. 3000 km latitudinal cline. We sequenced coding regions of 31 selected genes with known functions in metabolism, pigment production, development and heat shock responses. These include genes for which we expect a priori a role in thermal adaptation and, thus, varying selection pressures along a latitudinal cline, and genes we do not expect to vary clinally and can be used as controls. We identified amino acid substitution polymorphisms in 13 genes and tested these for clinal variation by correlation analysis of allele frequencies with latitude. In addition, we used two F_ST‐based outlier methods to identify loci with higher population differentiation than expected under neutral evolution, while accounting for potentially confounding effects of population structure and demographic history. Two metabolic enzymes of the glycolytic pathway, UGP and Treh, showed clinal variation. The same loci showed elevated population differentiation and were identified as significant outliers. We found no evidence of clines in the pigmentation genes, heat shock proteins and developmental genes. However, we identified outlier loci in more localized parts of the range in the pigmentation genes yellow and black. We discuss that the observed clinal variation and elevated population divergence in UGP and Treh may reflect adaptation to a geographic thermal gradient.     

Microevolution in time and space: SNP analysis of historical DNA reveals dynamic signatures of selection in Atlantic cod

Little is known about how quickly natural populations adapt to changes in their environment and how temporal and spatial variation in selection pressures interact to shape patterns of genetic diversity. We here address these issues with a series of genome scans in four overfished populations of Atlantic cod (Gadus morhua) studied over an 80‐year period. Screening of >1000 gene‐associated single‐nucleotide polymorphisms (SNPs) identified 77 loci that showed highly elevated levels of differentiation, likely as an effect of directional selection, in either time, space or both. Exploratory analysis suggested that temporal allele frequency shifts at certain loci may correlate with local temperature variation and with life history changes suggested to be fisheries induced. Interestingly, however, largely nonoverlapping sets of loci were temporal outliers in the different populations and outliers from the 1928 to 1960 period showed almost complete stability during later decades. The contrasting microevolutionary trajectories among populations resulted in sequential shifts in spatial outliers, with no locus maintaining elevated spatial differentiation throughout the study period. Simulations of migration coupled with observations of temporally stable spatial structure at neutral loci suggest that population replacement or gene flow alone could not explain all the observed allele frequency variation. Thus, the genetic changes are likely to at least partly be driven by highly dynamic temporally and spatially varying selection. These findings have important implications for our understanding of local adaptation and evolutionary potential in high gene flow organisms and underscore the need to carefully consider all dimensions of biocomplexity for evolutionarily sustainable management.