Molecular Inversion Probes for targeted resequencing in non-model organisms

Applications that require resequencing of hundreds or thousands of predefined genomic regions in numerous samples are common in studies of non-model organisms. However few approaches at the scale intermediate between multiplex PCR and sequence capture methods are available. Here we explored the utility of Molecular Inversion Probes (MIPs) for the medium-scale targeted resequencing in a non-model system. Markers targeting 112 bp of exonic sequence were designed from transcriptome of Lissotriton newts. We assessed performance of 248 MIP markers in a sample of 85 individuals. Among the 234 (94.4%) successfully amplified markers 80% had median coverage within one order of magnitude, indicating relatively uniform performance; coverage uniformity across individuals was also high. In the analysis of polymorphism and segregation within family, 77% of 248 tested MIPs were confirmed as single copy Mendelian markers. Genotyping concordance assessed using replicate samples exceeded 99%. MIP markers for targeted resequencing have a number of advantages: high specificity, high multiplexing level, low sample requirement, straightforward laboratory protocol, no need for preparation of genomic libraries and no ascertainment bias. We conclude that MIP markers provide an effective solution for resequencing targets of tens or hundreds of kb in any organism and in a large number of samples

Linkage map of Lissotriton newts provides insight into the genetic basis of reproductive isolation

Linkage maps are widely used to investigate structure, function, and evolution of genomes. In speciation research, maps facilitate the study of the genetic architecture of reproductive isolation by allowing identification of genomic regions underlying reduced fitness of hybrids. Here we present a linkage map for European newts of the Lissotriton vulgaris species complex, constructed using two families of F2 L. montandoni × L. vulgaris hybrids. The map consists of 1146 protein-coding genes on 12 linkage groups, equal to the haploid chromosome number, with a total length of 1484 cM (1.29 cM per marker). It is notably shorter than two other maps available for salamanders, but the differences in map length are consistent with cytogenetic estimates of the number of chiasmata per chromosomal arm. Thus, large salamander genomes do not necessarily translate into long linkage maps, as previously suggested. Consequently, salamanders are an excellent model to study evolutionary consequences of recombination rate variation in taxa with large genomes and a similar number of chromosomes. A complex pattern of transmission ratio distortion (TRD) was detected: TRD occurred mostly in one family, in one breeding season, and was clustered in two genomic segments. This is consistent with environment-dependent mortality of individuals carrying L. montandoni alleles in these two segments and suggests a role of TRD blocks in reproductive isolation. The reported linkage map will empower studies on the genomic architecture of divergence and interactions between the genomes of hybridizing newts

Constraint and adaptation in newt Toll-like receptor genes

Acute die-offs of amphibian populations worldwide have been linked to the emergence of viral and fungal diseases. Inter and intraspecific immunogenetic differences may influence the outcome of infection. Toll-like receptors (TLRs) are an essential component of innate immunity and also prime acquired defenses. We report the first comprehensive assessment of TLR gene variation for urodele amphibians. The Lissotriton newt TLR repertoire includes representatives of 13 families and is compositionally most similar to that of the anuran Xenopus. Both ancient and recent gene duplications have occurred in urodeles, bringing the total number of TLR genes to at least 21. Purifying selection has predominated the evolution of newt TLRs in both long (∼70 Ma) and medium (∼18 Ma) timescales. However, we find evidence for both purifying and positive selection acting on TLRs in two recently diverged (2–5 Ma) allopatric evolutionary lineages (Lissotriton montandoni and L. vulgaris graecus). Overall, both forms of selection have been stronger in L. v. graecus, while constraint on most TLR genes in L. montandoni appears relaxed. The differences in selection regimes are unlikely to be biased by demographic effects because these were controlled by means of a historical demographic model derived from an independent data set of 62 loci. We infer that TLR genes undergo distinct trajectories of adaptive evolution in closely related amphibian lineages, highlight the potential of TLRs to capture the signatures of different assemblages of pathogenic microorganisms, and suggest differences between lineages in the relative roles of innate and acquired immunity

The distributions of the six species constituting the smooth newt species complex (Lissotriton vulgaris sensu lato and L. montandoni) – an addition to the New Atlas of Amphibians and Reptiles of Europe

The ‘smooth newt’, the taxon traditionally referred to as Lissotriton vulgaris, consists of multiple morphologically distinct taxa. Given the uncertainty concerning the validity and rank of these taxa, L. vulgaris sensu lato has often been treated as a single, polytypic species. A recent study, driven by genetic data, proposed to recognize five species, L. graecus, L. kosswigi, L. lantzi, L. schmidtleri and a more restricted L. vulgaris. The Carpathian newt L. montandoni was confirmed to be a closely related sister species. We propose to refer to this collective of six Lissotriton species as the smooth newt or Lissotriton vulgaris species complex. Guided by comprehensive genomic data from throughout the range of the smooth newt species complex we 1) delineate the distribution ranges, 2) provide a distribution database, and 3) produce distribution maps according to the format of the New Atlas of Amphibians and Reptiles of Europe, for the six constituent species. This allows us to 4) highlight regions where more research is needed to determine the position of contact zones

Genetic and antigenic variation of the bovine tick-borne pathogen Theileria parva in the Great Lakes region of Central Africa

BACKGROUND : Theileria parva causes East Coast fever (ECF), one of the most economically important tick-borne diseases of cattle in sub-Saharan Africa. A live immunisation approach using the infection and treatment method (ITM) provides a strong long-term strain-restricted immunity. However, it typically induces a tick-transmissible carrier state in cattle and may lead to spread of antigenically distinct parasites. Thus, understanding the genetic composition of T. parva is needed prior to the use of the ITM vaccine in new areas. This study examined the sequence diversity and the evolutionary and biogeographical dynamics of T. parva within the African Great Lakes region to better understand the epidemiology of ECF and to assure vaccine safety. Genetic analyses were performed using sequences of two antigencoding genes, Tp1 and Tp2, generated among 119 T. parva samples collected from cattle in four agro-ecological zones of DRC and Burundi. RESULTS : The results provided evidence of nucleotide and amino acid polymorphisms in both antigens, resulting in 11 and 10 distinct nucleotide alleles, that predicted 6 and 9 protein variants in Tp1 and Tp2, respectively. Theileria parva samples showed high variation within populations and a moderate biogeographical sub-structuring due to the widespread major genotypes. The diversity was greater in samples from lowlands and midlands areas compared to those from highlands and other African countries. The evolutionary dynamics modelling revealed a signal of selective evolution which was not preferentially detected within the epitope-coding regions, suggesting that the observed polymorphism could be more related to gene flow rather than recent host immune-based selection. Most alleles isolated in the Great Lakes region were closely related to the components of the trivalent Muguga vaccine. CONCLUSIONS : Our findings suggest that the extensive sequence diversity of T. parva and its biogeographical distribution mainly depend on host migration and agro-ecological conditions driving tick population dynamics. Such patterns are likely to contribute to the epidemic and unstable endemic situations of ECF in the region. However, the fact that ubiquitous alleles are genetically similar to the components of the Muguga vaccine together with the limited geographical clustering may justify testing the existing trivalent vaccine for cross-immunity in the region.Additional file 1: Table S1. Cattle blood sample distribution across agroecological zones.Additional file 2: Table S2. Nucleotide and amino acid sequences of Tp1 and Tp2 antigen epitopes from T. parva Muguga reference sequence.Additional file 3: Table S3. Characteristics of 119 T. parva samples obtained from cattle in different agro-ecological zones (AEZs) of The Democratic Republic of Congo and Burundi.Additional file 4: Figure S1. Multiple sequence alignment of the 11 Tp1 gene alleles obtained in this study.Additional file 5: Table S4. Estimates of evolutionary divergence between gene alleles for Tp1 and Tp2, using proportion nucleotide distance.Additional file 6: Table S5. Tp1 and Tp2 genes alleles with their corresponding antigen variants.Additional file 7: Table S6. Amino acid variants of Tp1 and Tp2 CD8+ T cell target epitopes of T. parva from DRC and Burundi.Additional file 8: Figure S2. Multiple sequence alignment of the 10 Tp2 gene alleles obtained in this study.Additional file 9: Table S7. Distribution of Tp1 gene alleles of T. parva from cattle and buffalo in the sub-Saharan region of Africa.Additional file 10: Table S8. Distribution of Tp2 gene alleles of T. parva from cattle and buffalo in the sub-Saharan region of Africa.Additional file 11: Figure S3. Neighbor-joining tree showing phylogenetic relationships among 48 Tp1 gene alleles described in Africa.Additional file 12: Figure S4. Phylogenetic tree showing the relationships among concatenated Tp1 and Tp2 nucleotide sequences of 93 T. parva samples from cattle in DRC and Burundi.This study is part of the PhD work supported by the University of Namur (UNamur, Belgium) through the UNamur-CERUNA institutional PhD grant awarded to GSA for bioinformatic analyses, interpretation of data and manuscript write up in Belgium. The laboratory aspects (molecular biology analysis) of the project were supported by the BecA-ILRI Hub through the Africa Biosciences Challenge Fund (ABCF) programme. The ABCF Programme is funded by the Australian Department for Foreign Affairs and Trade (DFAT) through the BecA-CSIRO partnership; the Syngenta Foundation for Sustainable Agriculture (SFSA); the Bill & Melinda Gates Foundation (BMGF); the UK Department for International Development (DFID); and the Swedish International Development Cooperation Agency (Sida). The ABCF Fellowship awarded to GAS was funded by BMGF grant (OPP1075938). Sample collection, field equipment and preliminary sample processing were supported through the “Theileria” project co-funded to the Université Evangélique en Afrique (UEA) by the Agence Universitaire de la Francophonie (AUF) and the Communauté Economique des Pays des Grands Lacs (CEPGL). The International Foundation for Science (IFS, Stockholm, Sweden) supported the individual scholarship awarded to GSA (grant no. IFS-92890CA3) for field work and part of field equipment to the “Theileria” project.http://www.parasitesandvectors.comam2020Veterinary Tropical Disease

Restriction associated DNA-genotyping at multiple spatial scales in Arabidopsis lyrata reveals signatures of pathogen-mediated selection

Background: Genome scans based on outlier analyses have revolutionized detection of genes involved in adaptive processes, but reports of some forms of selection, such as balancing selection, are still limited. It is unclear whether high throughput genotyping approaches for identification of single nucleotide polymorphisms have sufficient power to detect modes of selection expected to result in reduced genetic differentiation among populations. In this study, we used Arabidopsis lyrata to investigate whether signatures of balancing selection can be detected based on genomic smoothing of Restriction Associated DNA sequencing (RAD-seq) data. We compared how different sampling approaches (both within and between subspecies) and different background levels of polymorphism (inbreeding or outcrossing populations) affected the ability to detect genomic regions showing key signatures of balancing selection, specifically elevated polymorphism, reduced differentiation and shifts towards intermediate allele frequencies. We then tested whether candidate genes associated with disease resistance (R-gene analogs) were detected more frequently in these regions compared to other regions of the genome. Results: We found that genomic regions showing elevated polymorphism contained a significantly higher density of R-gene analogs predicted to be under pathogen-mediated selection than regions of non-elevated polymorphism, and that many of these also showed evidence for an intermediate site-frequency spectrum based on Tajima’s D. However, we found few genomic regions that showed both elevated polymorphism and reduced FST among populations, despite strong background levels of genetic differentiation among populations. This suggests either insufficient power to detect the reduced population structure predicted for genes under balancing selection using sparsely distributed RAD markers, or that other forms of diversifying selection are more common for the R-gene analogs tested. Conclusions: Genome scans based on a small number of individuals sampled from a wide range of populations were sufficient to confirm the relative scarcity of signatures of balancing selection across the genome, but also identified new potential disease resistance candidates within genomic regions showing signatures of balancing selection that would be strong candidates for further sequencing efforts

Para-allopatry in hybridizing fire-bellied toads (Bombina bombina and B. variegata):Inference from transcriptome-wide coalescence analyses

Ancient origins, profound ecological divergence, and extensive hybridization make the fire‐bellied toads Bombina bombina and B. variegata (Anura: Bombinatoridae) an intriguing test case of ecological speciation. Previous modeling has proposed that the narrow Bombina hybrid zones represent strong barriers to neutral introgression. We test this prediction by inferring the rate of gene exchange between pure populations on either side of the intensively studied Kraków transect. We developed a method to extract high confidence sets of orthologous genes from de novo transcriptome assemblies, fitted a range of divergence models to these data and assessed their relative support with analytic likelihood calculations. There was clear evidence for postdivergence gene flow, but, as expected, no perceptible signal of recent introgression via the nearby hybrid zone. The analysis of two additional Bombina taxa (B. v. scabra and B. orientalis) validated our parameter estimates against a larger set of prior expectations. Despite substantial cumulative introgression over millions of years, adaptive divergence of the hybridizing taxa is essentially unaffected by their lack of reproductive isolation. Extended distribution ranges also buffer them against small‐scale environmental perturbations that have been shown to reverse the speciation process in other, more recent ecotypes

Selective landscapes in newt immune genes inferred from patterns of nucleotide variation

Host–pathogen interactions may result in either directional selection or in pressure for the maintenance of polymorphism at the molecular level. Hence signatures of both positive and balancing selection are expected in immune genes. Because both overall selective pressure and specific targets may differ between species, large-scale population genomic studies are useful in detecting functionally important immune genes and comparing selective landscapes between taxa. Such studies are of particular interest in amphibians, a group threatened worldwide by emerging infectious diseases. Here, we present an analysis of polymorphism and divergence of 634 immune genes in two lineages of Lissotriton newts: L. montandoni and L. vulgaris graecus. Variation in newt immune genes has been shaped predominantly by widespread purifying selection and strong evolutionary constraint, implying long-term importance of these genes for functioning of the immune system. The two evolutionary lineages differ in the overall strength of purifying selection which can partially be explained by demographic history but may also signal differences in long-term pathogen pressure. The prevalent constraint notwithstanding, 23 putative targets of positive selection and 11 putative targets of balancing selection were identified. The latter were detected by composite tests involving the demographic model and further validated in independent population samples. Putative targets of balancing selection encode proteins which may interact closely with pathogens but include also regulators of immune response. The identified candidates will be useful for testing whether genes affected by balancing selection are more prone to interspecific introgression than other genes in the genome