Contrasting patterns of selection between MHC I and II across populations of Humboldt and Magellanic penguins

The evolutionary and adaptive potential of populations or species facing an emerging infectious disease depends on their genetic diversity in genes, such as the major histocompatibility complex (MHC). In birds, MHC class I deals predominantly with intracellular infections (e.g., viruses) and MHC class II with extracellular infections (e.g., bacteria). Therefore, patterns of MHC I and II diversity may differ between species and across populations of species depending on the relative effect of local and global environmental selective pressures, genetic drift, and gene flow. We hypothesize that high gene flow among populations of Humboldt and Magellanic penguins limits local adaptation in MHC I and MHC II, and signatures of selection differ between markers, locations, and species. We evaluated the MHC I and II diversity using 454 next‐generation sequencing of 100 Humboldt and 75 Magellanic penguins from seven different breeding colonies. Higher genetic diversity was observed in MHC I than MHC II for both species, explained by more than one MHC I loci identified. Large population sizes, high gene flow, and/or similar selection pressures maintain diversity but limit local adaptation in MHC I. A pattern of isolation by distance was observed for MHC II for Humboldt penguin suggesting local adaptation, mainly on the northernmost studied locality. Furthermore, trans‐species alleles were found due to a recent speciation for the genus or convergent evolution. High MHC I and MHC II gene diversity described is extremely advantageous for the long‐term survival of the species.     

Multilocus patterns of genetic variation across Cryptosporidium species suggest balancing selection at the gp60 locus

Cryptosporidium is an apicomplexan protozoan that lives in most vertebrates, including humans. Its gp60 gene is functionally involved in its attachment to host cells, and its high level of genetic variation has made it the reference marker for sample typing in epidemiological studies. To understand the origin of such high diversity and to determine the extent to which this classification applies to the rest of the genome, we analysed the patterns of variation at gp60 and nine other nuclear loci in isolates of three Cryptosporidium species. Most loci showed low genetic polymorphism (π_S <1%) and similar levels of between‐species divergence. Contrastingly, gp60 exhibited very different characteristics: (i) it was nearly ten times more variable than the other loci; (ii) it displayed a significant excess of polymorphisms relative to between‐species differences in a maximum‐likelihood Hudson–Kreitman–Aguadé test; (iii) gp60 subtypes turned out to be much older than the species they were found in; and (iv) showed a significant excess of polymorphic variants shared across species from random expectations. These observations suggest that this locus evolves under balancing selection and specifically under negative frequency‐dependent selection (FDS). Interestingly, genetic variation at the other loci clusters very well within the groups of isolates defined by gp60 subtypes, which may provide new tools to understand the genome‐wide patterns of genetic variation of the parasite in the wild. These results suggest that gp60 plays an active and essential role in the life cycle of the parasite and that genetic variation at this locus might be essential for the parasite's long‐term success.     

Experimental evidence for accelerated adaptation to desiccation through sexual selection on males

The impact of sexual selection on the adaptive process remains unclear. On the one hand, sexual selection might hinder adaptation by favouring costly traits and preferences that reduce nonsexual fitness. On the other hand, condition dependence of success in sexual selection may accelerate adaptation. Here, we used replicate populations of Drosophila melanogaster to artificially select on male desiccation resistance while manipulating the opportunity for precopulatory sexual selection in a factorial design. Following five generations of artificial selection, we measured the desiccation resistance of males and females to test whether the addition of sexual selection accelerated adaptation. We found a significant interaction between the effects of natural selection and sexual selection: desiccation resistance was highest in populations where sexual selection was allowed to operate. Despite only selecting on males, we also found a correlated response in females. These results provide empirical support for the idea that sexual selection can accelerate the rate of adaptation.     

Mito‐nuclear discordance with evidence of shared ancestral polymorphism and selection in cactophilic species of Drosophila

The Drosophila serido haplogroup is a monophyletic group composed of the following four cryptic and cactophilic species that are endemic to eastern Brazil: D. borborema, D. gouveai, D. seriema and D. serido. Here, we investigate the mito‐nuclear discordance in these species found among the cytochrome c oxidase subunit I (COI) mitochondrial gene, the autosomal alpha‐Esterase‐5 (α‐Est5) and the X‐linked period gene (per). Our analysis indicates that shared polymorphisms in these three molecular markers may be explained by the maintenance of ancestral polymorphisms rather than introgressive hybridization. The primary structures of COI, per and α‐Est5 genes evolve primarily under purifying selection, but we detected some sites that evolved under positive selection in α‐Est5. Considering the high variability of cacti species in eastern Brazil and the role attributed to Drosophila esterases in digestion metabolism and/or the detoxification of several compounds found in cactus tissues, we conjecture about the role of natural selection triggered by host shifts as an important factor in the intraspecific diversification of the D. serido haplogroup.     

Gene flow,divergent selection and resistance to introgression in two species of morning glories (Ipomoea)

Gene flow is thought to impede genetic divergence and speciation by homogenizing genomes. Recent theory and research suggest that sufficiently strong divergent selection can overpower gene flow, leading to loci that are highly differentiated compared to others. However, there are also alternative explanations for this pattern. Independent evidence that loci in highly differentiated regions are under divergent selection would allow these explanations to be distinguished, but such evidence is scarce. Here, we present multiple lines of evidence that many of the highly divergent SNPs in a pair of sister morning glory species, Ipomoea cordatotriloba and I. lacunosa, are the result of divergent selection in the face of gene flow. We analysed a SNP data set across the genome to assess the amount of gene flow, resistance to introgression and patterns of selection on loci resistant to introgression. We show that differentiation between the two species is much lower in sympatry than in allopatry, consistent with interspecific gene flow in sympatry. Gene flow appears to be substantially greater from I. lacunosa to I. cordatotriloba than in the reverse direction, resulting in sympatric and allopatric I. cordatotriloba being substantially more different than sympatric and allopatric I. lacunosa. Many SNPs highly differentiated in allopatry have experienced divergent selection, and, despite gene flow in sympatry, resist homogenization in sympatry. Finally, five out of eight floral and inflorescence characteristics measured exhibit asymmetric convergence in sympatry. Consistent with the pattern of gene flow, I. cordatotriloba traits become much more like those of I. lacunosa than the reverse. Our investigation reveals the complex interplay between selection and gene flow that can occur during the early stages of speciation.     

QTL sequencing strategy to map genomic regions associated with resistance to ascochyta blight in chickpea

Whole‐genome sequencing‐based bulked segregant analysis (BSA) for mapping quantitative trait loci (QTL) provides an efficient alternative approach to conventional QTL analysis as it significantly reduces the scale and cost of analysis with comparable power to QTL detection using full mapping population. We tested the application of next‐generation sequencing (NGS)‐based BSA approach for mapping QTLs for ascochyta blight resistance in chickpea using two recombinant inbred line populations CPR‐01 and CPR‐02. Eleven QTLs in CPR‐01 and six QTLs in CPR‐02 populations were mapped on chromosomes Ca1, Ca2, Ca4, Ca6 and Ca7. The QTLs identified in CPR‐01 using conventional biparental mapping approach were used to compare the efficiency of NGS‐based BSA in detecting QTLs for ascochyta blight resistance. The QTLs on chromosomes Ca1, Ca4, Ca6 and Ca7 overlapped with the QTLs previously detected in CPR‐01 using conventional QTL mapping method. The QTLs on chromosome Ca4 were detected in both populations and overlapped with the previously reported QTLs indicating conserved region for ascochyta blight resistance across different chickpea genotypes. Six candidate genes in the QTL regions identified using NGS‐based BSA on chromosomes Ca2 and Ca4 were validated for their association with ascochyta blight resistance in the CPR‐02 population. This study demonstrated the efficiency of NGS‐based BSA as a rapid and cost‐effective method to identify QTLs associated with ascochyta blight in chickpea.     

Evolution and selection of Rhg1, a copy‐number variant nematode‐resistance locus

The soybean cyst nematode (SCN) resistance locus Rhg1 is a tandem repeat of a 31.2 kb unit of the soybean genome. Each 31.2‐kb unit contains four genes. One allele of Rhg1, Rhg1‐b, is responsible for protecting most US soybean production from SCN. Whole‐genome sequencing was performed, and PCR assays were developed to investigate allelic variation in sequence and copy number of the Rhg1 locus across a population of soybean germplasm accessions. Four distinct sequences of the 31.2‐kb repeat unit were identified, and some Rhg1 alleles carry up to three different types of repeat unit. The total number of copies of the repeat varies from 1 to 10 per haploid genome. Both copy number and sequence of the repeat correlate with the resistance phenotype, and the Rhg1 locus shows strong signatures of selection. Significant linkage disequilibrium in the genome outside the boundaries of the repeat allowed the Rhg1 genotype to be inferred using high‐density single nucleotide polymorphism genotyping of 15 996 accessions. Over 860 germplasm accessions were found likely to possess Rhg1 alleles. The regions surrounding the repeat show indications of non‐neutral evolution and high genetic variability in populations from different geographic locations, but without evidence of fixation of the resistant genotype. A compelling explanation of these results is that balancing selection is in operation at Rhg1.     

Positive selection on the mitochondrial ATP synthase 6 and the NADH dehydrogenase 2 genes across 22 hare species (genus Lepus)

Amino acid changes in mitochondrial (mt) oxidative phosphorylation (OXPHOS) genes have been suggested as a key adaptation to environmental variation. Here, we analyzed 416 sequences of ATPase synthase 6 (MT‐ATP6) and NADH dehydrogenase 2 (MT‐ND2) in 22 different hare (Lepus) species from across a wide range of habitats and climates. We used site‐ and branch‐based methods to test for positive selection on specific codons and lineages. We found four codons in MT‐ATP6 and five in MT‐ND2 under positive selection, affecting several species lineages. We investigated the association of protein variants at each locus with climate zone, using multinomial generalized linear models (glm), including species, regions, historical introgression events, and the co‐occurring protein variant at the other locus as additional explanatory variables. A significant climate effect as based on the “Köppen climate classification” was observed for MT‐ND2 protein variants as translated from our nucleotide sequences. Moreover, MT‐ND2 protein variants were significantly affected by the co‐occurring MT‐ATP6 protein variant in the same mtDNA molecule. Contrary to the expectation for non‐recombining mitochondrial DNA molecules, the presence of an evolutionarily relatively ancestral protein variant at one locus was associated with a relatively derived protein at the other locus in the same mitochondrial molecule, respectively. The relative evolutionary status of a protein variant was evaluated according to its positions relative to the respective out‐group protein variant in a network analysis of nucleotide sequences. All our results suggest a complex effect of various climatic parameters acting on multiple mtOXPHOS genes in a co‐adaptive way, favoring combinations of ancestral and derived variants.     

Artificial selection for increased dispersal results in lower fitness

Dispersal often covaries with other traits, and this covariation was shown to have a genetic basis. Here, we wanted to explore to what extent genetic constraints and correlational selection can explain patterns of covariation between dispersal and key life‐history traits—lifespan and reproduction. A prediction from the fitness‐associated dispersal hypothesis was that lower genetic quality is associated with higher dispersal propensity as driven by the benefits of genetic mixing. We wanted to contrast it with a prediction from a different model that individuals putting more emphasis on current rather than future reproduction disperse more, as they are expected to be more risk‐prone and exploratory. However, if dispersal has inherent costs, this will also result in a negative genetic correlation between higher rates of dispersal and some aspects of performance. To explore this issue, we used the dioecious nematode Caenorhabditis remanei and selected for increased and decreased dispersal propensity for 10 generations, followed by five generations of relaxed selection. Dispersal propensity responded to selection, and females from high‐dispersal lines dispersed more than females from low‐dispersal lines. Females selected for increased dispersal propensity produced fewer offspring and were more likely to die from matricide, which is associated with a low physiological condition in Caenorhabditis nematodes. There was no evidence for differences in age‐specific reproductive effort between high‐ and low‐dispersal females. Rather, reproductive output of high‐dispersal females was consistently reduced. We argue that our data provide support for the fitness‐associated dispersal hypothesis.     

Genome‐wide detection of signatures of selection in Korean Hanwoo cattle

The Korean Hanwoo cattle have been intensively selected for production traits, especially high intramuscular fat content. It is believed that ancient crossings between different breeds contributed to forming the Hanwoo, but little is known about the genomic differences and similarities between other cattle breeds and the Hanwoo. In this work, cattle breeds were grouped by origin into four types and used for comparisons: the Europeans (represented by six breeds), zebu (Nelore), African taurine (N'Dama) and Hanwoo. All animals had genotypes for around 680 000 SNPs after quality control of genotypes. Average heterozygosity was lower in Nelore and N'Dama (0.22 and 0.21 respectively) than in Europeans (0.26–0.31, with Shorthorn as outlier at 0.24) and Hanwoo (0.29). Pairwise F_ST analyses demonstrated that Hanwoo are more related to European cattle than to Nelore, with N'Dama in an intermediate position. This finding was corroborated by principal components and unsupervised hierarchical clustering. Using genome‐wide smoothed F_ST, 55 genomic regions potentially under positive selection in Hanwoo were identified. Among these, 29 were regions also detected in previous studies. Twenty‐four regions were exclusive to Hanwoo, and a number of other regions were shared with one or two of the other groups. These regions overlap a number of genes that are related to immune, reproduction and fatty acid metabolism pathways. Further analyses are needed to better characterize the ancestry of the Hanwoo cattle and to define the genes responsible to the identified selection peaks.     

Analysis of ROH patterns in the Noriker horse breed reveals signatures of selection for coat color and body size

Overlapping runs of homozygosity (ROH islands) shared by the majority of a population are hypothesized to be the result of selection around a target locus. In this study we investigated the impact of selection for coat color within the Noriker horse on autozygosity and ROH patterns. We analyzed overlapping homozygous regions (ROH islands) for gene content in fragments shared by more than 50% of horses. Long‐term assortative mating of chestnut horses and the small effective population size of leopard spotted and tobiano horses resulted in higher mean genome‐wide ROH coverage (S_ROH) within the range of 237.4–284.2 Mb, whereas for bay, black and roan horses, where rotation mating is commonly applied, lower autozygosity (S_ROH from 176.4–180.0 Mb) was determined. We identified seven common ROH islands considering all Noriker horses from our dataset. Specific islands were documented for chestnut, leopard spotted, roan and bay horses. The ROH islands contained, among others, genes associated with body size (ZFAT, LASP1 and LCORL/NCAPG), coat color (MC1R in chestnut and the factor PATN1 in leopard spotted horses) and morphogenesis (HOXB cluster in all color strains except leopard spotted horses). This study demonstrates that within a closed population sharing the same founders and ancestors, selection on a single phenotypic trait, in this case coat color, can result in genetic fragmentation affecting levels of autozygosity and distribution of ROH islands and enclosed gene content.     

Response to joint selection on germination and flowering phenology depends on the direction of selection

Flowering and germination time are components of phenology, a complex phenotype that incorporates a number of traits. In natural populations, selection is likely to occur on multiple components of phenology at once. However, we have little knowledge of how joint selection on several phenological traits influences evolutionary response. We conducted one generation of artificial selection for all combinations of early and late germination and flowering on replicated lines within two independent base populations in the herb Campanula americana. We then measured response to selection and realized heritability for each trait. Response to selection and heritability were greater for flowering time than germination time, indicating greater evolutionary potential of this trait. Selection for earlier phenology, both flowering and germination, did not depend on the direction of selection on the other trait, whereas response to selection to delay germination and flowering was greater when selection on the other trait was in the opposite direction (e.g., early germination and late flowering), indicating a negative genetic correlation between the traits. Therefore, the extent to which correlations shaped response to selection depended on the direction of selection. Furthermore, the genetic correlation between timing of germination and flowering varies across the trait distributions. The negative correlation between germination and flowering time found when selecting for delayed phenology follows theoretical predictions of constraint for traits that jointly determine life history schedule. In contrast, the lack of constraint found when selecting for an accelerated phenology suggests a reduction of the covariance due to strong selection favoring earlier flowering and a shorter life cycle. This genetic architecture, in turn, will facilitate further evolution of the early phenology often favored in warm climates.     

An efficient method for measuring copy number variation applied to improvement of nematode resistance in soybean

Copy number variation (CNV) is implicated in important traits in multiple crop plants, but can be challenging to genotype using conventional methods. The Rhg1 locus of soybean, which confers resistance to soybean cyst nematode (SCN), is a CNV of multiple 31.2‐kb genomic units each containing four genes. Reliable, high‐throughput methods to quantify Rhg1 and other CNVs for selective breeding were developed. The CNV genotyping assay described here uses a homeologous gene copy within the paleopolyploid soybean genome to provide the internal control for a single‐tube TaqMan copy number assay. Using this assay, CNV in breeding populations can be tracked with high precision. We also show that extensive CNV exists within Fayette, a released, inbred SCN‐resistant soybean cultivar with a high copy number at Rhg1 derived from a single donor parent. Copy number at Rhg1 is therefore unstable within a released variety over a relatively small number of generations. Using this assay to select for individuals with altered copy number, plants were obtained with both increased copy number and increased SCN resistance relative to control plants. Thus, CNV genotyping technologies can be used as a new type of marker‐assisted selection to select for desirable traits in breeding populations, and to control for undesirable variation within cultivars.     

Nonlinear frequency‐dependent selection promotes long‐term coexistence between bacteria species

Negative frequency‐dependent selection (NFDS) is an important mechanism for species coexistence and for the maintenance of genetic polymorphism. Long‐term coexistence nevertheless requires NFDS interactions to be resilient to further evolution of the interacting species or genotypes. For closely related genotypes, NFDS interactions have been shown to be preserved through successive rounds of evolution in coexisting lineages. On the contrary, the evolution of NFDS interactions between distantly related species has received less attention. Here, we tracked the co‐evolution of Escherichia coli and Citrobacter freundii that initially differ in their ecological characteristics. We showed that these two bacterial species engaged in an NFDS interaction particularly resilient to further evolution: despite a very strong asymmetric rate of adaptation, their coexistence was maintained owing to an NFDS pattern where fitness increases steeply as the frequency decreases towards zero. Using a model, we showed how and why such NFDS pattern can emerge. These findings provide a robust explanation for the long‐term maintenance of species at very low frequencies.     

Comparative study of sequence aligners for detecting antibiotic resistance in bacterial metagenomes

We aim to compare the performance of Bowtie2 , bwa‐mem , blastn and blastx when aligning bacterial metagenomes against the Comprehensive Antibiotic Resistance Database (CARD). Simulated reads were used to evaluate the performance of each aligner under the following four performance criteria: correctly mapped, false positives, multi‐reads and partials. The optimal alignment approach was applied to samples from two wastewater treatment plants to detect antibiotic resistance genes using next generation sequencing. blastn mapped with greater accuracy among the four sequence alignment approaches considered followed by Bowtie2 . blastx generated the greatest number of false positives and multi‐reads when aligned against the CARD. The performance of each alignment tool was also investigated using error‐free reads. Although each aligner mapped a greater number of error‐free reads as compared to Illumina‐error reads, in general, the introduction of sequencing errors had little effect on alignment results when aligning against the CARD. Given each performance criteria, blastn was found to be the most favourable alignment tool and was therefore used to assess resistance genes in sewage samples. Beta‐lactam and aminoglycoside were found to be the most abundant classes of antibiotic resistance genes in each sample.

Significance and Impact of the Study

Antibiotic resistance genes (ARGs) are pollutants known to persist in wastewater treatment plants among other environments, thus methods for detecting these genes have become increasingly relevant. Next generation sequencing has brought about a host of sequence alignment tools that provide a comprehensive look into antimicrobial resistance in environmental samples. However, standardizing practices in ARG metagenomic studies is challenging since results produced from alignment tools can vary significantly. Our study provides sequence alignment results of synthetic, and authentic bacterial metagenomes mapped against an ARG database using multiple alignment tools, and the best practice for detecting ARGs in environmental samples.     

Ecology and multilevel selection explain aggression in spider colonies

Progress in sociobiology continues to be hindered by abstract debates over methodology and the relative importance of within‐group vs. between‐group selection. We need concrete biological examples to ground discussions in empirical data. Recent work argued that the levels of aggression in social spider colonies are explained by group‐level adaptation. Here, we examine this conclusion using models that incorporate ecological detail while remaining consistent with kin‐ and multilevel selection frameworks. We show that although levels of aggression are driven, in part, by between‐group selection, incorporating universal within‐group competition provides a striking fit to the data that is inconsistent with pure group‐level adaptation. Instead, our analyses suggest that aggression is favoured primarily as a selfish strategy to compete for resources, despite causing lower group foraging efficiency or higher risk of group extinction. We argue that sociobiology will benefit from a pluralistic approach and stronger links between ecologically informed models and data.