The utility of single nucleotide polymorphisms in inferences of population history

Genetic diversity and demographic history of the largest remaining migratory population of brindled wildebeest (Connochaetes taurinus taurinus) in southern Africa

The blue wildebeest (Connochaetes taurinus) is a keystone species in the savannahs of southern Africa, where it maintains shortgrass plains and serves as an important prey source for large carnivores. Despite being the second-largest migratory wildebeest population, the brindled wildebeest (C. t. taurinus) of the Greater Liuwa Ecosystem (GLE) of western Zambia have remained largely unstudied, until recently. While studies have increased understanding of recent demography, migration, and population limiting factors, the level of genetic diversity, patterns of gene flow, and long-term demographic history of brindled wildebeest in the GLE remains unknown. Most genetic studies of wildebeest have focused on small, heavily-managed populations, rather than large, migratory populations of high conservation significance. We used restriction-site associated DNA sequencing (RAD-seq) to assess genetic diversity, population structure, and demographic history of brindled wildebeest in the GLE. Using SNPs from 1,730 loci genotyped across 75 individuals, we found moderate levels of genetic diversity in GLE brindled wildebeest (He =  0.210), very low levels of inbreeding (FIS =  0.033), and an effective population size of about one tenth the estimated population size. No genetic population structure was evident within the GLE. Analyses of the site frequency spectrum found signatures of expansion during the Middle Pleistocene followed by population decline in the Late Pleistocene and early Holocene, a pattern previously observed in other African ungulates. These results will supplement field studies in developing effective conservation plans for wildebeest as they face continued and increasing threats of habitat loss, poaching, and other human impacts across their remaining range.

Typing of Echinococcus multilocularis by Region-Specific Extraction and Next-Generation Sequencing of the mitogenome

Background

Infection by the fox tapeworm Echinococcus multilocularis may lead to a severe zoonosis in humans, alveolar echinococcosis, which may be fatal if left untreated. Typing is important to understand the epidemiology of this parasite, yet there is limited knowledge on the microdiversity of E. multilocularis on the local scale, since the typing resolution of established methods is restricted.

Methods

The mitogenome of E. multilocularis was used as the target regions to modify, apply and validate the Region-Specific Extraction (RSE) method in combination with Next-Generation Sequencing (NGS). Single Nucleotide Polymorphisms (SNPs) were detected in the mitochondrial DNA (mtDNA) and analysed bioinformatically. To validate the success and the accuracy of the RSE protocol, the mitogenomes of some E. multilocularis isolates were also analysed by the Whole-Genome Sequencing (WGS).

Results

With the chosen combination of methods, the entire mitogenome (~13 kb) of E. multilocularis could be captured and amplified. The read depth (median ≥ 156X) was sufficient to detect existing SNPs. The comparison of mitogenome sequences extracted by RSE with mitogenome sequences obtained by WGS showed that the accuracy of the RSE method was consistently comparable to direct Whole-Genome Sequencing.

Conclusion

The results demonstrate that the RSE method in combination with NGS is suitable to analyse the microdiversity of E. multilocularis at the whole mitogenome level. For the capture and sequencing of large (several kb) genomic regions of E. multilocularis and other applications, this method can be very helpful.

Short INDELs and SNPs as markers of evolutionary processes in hybrid zones

Polymorphic short insertions and deletions (INDELs ≤ 50 bp) are abundant, although less common than single nucleotide polymorphisms (SNPs). Evidence from model organisms shows INDELs to be more strongly influenced by purifying selection than SNPs. Partly for this reason, INDELs are rarely used as markers for demographic processes or to detect divergent selection. Here, we compared INDELs and SNPs in the intertidal snail Littorina saxatilis, focussing on hybrid zones between ecotypes, in order to test the utility of INDELs in the detection of divergent selection. We computed INDEL and SNP site frequency spectra using capture sequencing data. We assessed the impact of divergent selection by analyzing allele frequency clines across habitat boundaries. We also examined the influence of GC-biased gene conversion because it may be confounded with signatures of selection. We show evidence that short INDELs are affected more by purifying selection than SNPs, but part of the observed site frequency spectra difference can be attributed to GC-biased gene conversion. We did not find a difference in the impact of divergent selection between short INDELs and SNPs. Short INDELs and SNPs were similarly distributed across the genome and so are likely to respond to indirect selection in the same way. A few regions likely affected by divergent selection were revealed by INDELs and not by SNPs. Short INDELs can be useful (additional) genetic markers helping to identify genomic regions important for adaptation and population divergence.

Non-Invasive Prenatal Diagnosis of Chromosomal and Monogenic Disease by a Novel Bioinspired Micro-Nanochip for Isolating Fetal Nucleated Red Blood Cells

Purpose

Fetal nucleated red blood cells (fNRBCs) in the peripheral blood of pregnant women contain comprehensive fetal genetic information, making them an ideal target for non-invasive prenatal diagnosis (NIPD). However, challenges in identifying, enriching, and detecting fNRBCs limit their diagnostic potential.

Methods

To overcome these obstacles, we developed a novel biomimetic chip, replicating the micro-nano structure of red rose petals on polydimethylsiloxane (PDMS). The surface was modified with gelatin nanoparticles (GNPs) and affinity antibodies to enhance cell adhesion and facilitate specific cell identification. We subsequently investigated the chip's characteristics, along with its in vitro capture and release system, and conducted further experiments using peripheral blood samples from pregnant women.

Results

In the cell line capture and release assay, the chip achieved a cell capture efficiency of 90.4%. Following metalloproteinase-9 (MMP-9) enzymatic degradation, the release efficiency was 84.08%, with cell viability at 85.97%. Notably, fNRBCs can be captured from the peripheral blood of pregnant women as early as 7 weeks of gestation. We used these fNRBCs to diagnose a case of single-gene disease and instances of chromosomal aneuploidies, yielding results consistent with those obtained from amniotic fluid punctures.

Conclusion

This novel chip not only enables efficient enrichment of fNRBCs for NIPD but also extends the diagnostic window for genetic and developmental disorders to as early as 7 weeks of gestation, potentially allowing for earlier interventions. By improving the accuracy and reliability of NIPD, this technology could reduce reliance on invasive diagnostic techniques, offering a new pathway for diagnosing fetal genetic conditions in clinical practice.

Identification of Insertion and Deletion (InDel) Markers for Chickpea (Cicer arietinum L.) Based on Double-Digest Restriction Site-Associated DNA Sequencing

Enhancing the marker repository and the development of breeder-friendly markers in chickpeas is important in relation to chickpea genomics-assisted breeding applications. Insertion-deletion (InDel) markers are widely distributed across genomes and easily observed with specifically designed primers, leading to less time, cost, and labor requirements. In light of this, the present study focused on the identification and development of InDel markers through the use of double-digest restriction site-associated DNA sequencing (ddRADSeq) data from 20 chickpea accessions. Bioinformatic analysis identified 20,700 InDel sites, including 15,031 (72.61%) deletions and 5669 (27.39%) insertions, among the chickpea accessions. The InDel markers ranged from 1 to 25 bp in length, while single-nucleotide-length InDel markers were found to represent the majority of the InDel sites and account for 79% of the total InDel markers. However, we focused on InDel markers wherein the length was greater than a single nucleotide to avoid any read or alignment errors. Among all of the InDel markers, 96.1% were less than 10 bp, 3.6% were between 10 and 20 bp, and 0.3% were more than 20 bp in length. We examined the InDel markers that were 10 bp and longer for the development of InDel markers based on a consideration of the genomic distribution and low-cost genotyping with agarose gels. A total of 29 InDel regions were selected, and primers were successfully designed to evaluate their efficiency. Annotation analysis of the InDel markers revealed them to be found with the highest frequency in the intergenic regions (82.76%), followed by the introns (6.90%), coding sequences (6.90%), and exons (3.45%). Genetic diversity analysis demonstrated that the polymorphic information content of the markers varied from 0.09 to 0.37, with an average of 0.20. Taken together, these results showed the efficiency of InDel marker development for chickpea genetic and genomic studies using the ddRADSeq method. The identified markers might prove valuable for chickpea breeders.

Draft Genome of Akame (Lates Japonicus) Reveals Possible Genetic Mechanisms for Long-Term Persistence and Adaptive Evolution with Low Genetic Diversity

It is known that some endangered species have persisted for thousands of years despite their very small effective population sizes and low levels of genetic polymorphisms. To understand the genetic mechanisms of long-term persistence in threatened species, we determined the whole genome sequences of akame (Lates japonicus), which has survived for a long time with extremely low genetic variations. Genome-wide heterozygosity in akame was estimated to be 3.3 to 3.4 × 10-4/bp, one of the smallest values in teleost fishes. Analysis of demographic history revealed that the effective population size in akame was around 1,000 from 30,000 years ago to the recent past. The relatively high ratio of nonsynonymous to synonymous heterozygosity in akame indicated an increased genetic load. However, a detailed analysis of genetic diversity in the akame genome revealed that multiple genomic regions, including genes involved in immunity, synaptic development, and olfactory sensory systems, have retained relatively high nucleotide polymorphisms. This implies that the akame genome has preserved the functional genetic variations by balancing selection, to avoid a reduction in viability and loss of adaptive potential. Analysis of synonymous and nonsynonymous nucleotide substitution rates has detected signs of positive selection in many akame genes, suggesting adaptive evolution to temperate waters after the speciation of akame and its close relative, barramundi (Lates calcarifer). Our results indicate that the functional genetic diversity likely contributed to the long-term persistence of this species by avoiding the harmful effects of the population size reduction.

Whole genome assembly and annotation of the King Angelfish (Holacanthus passer) gives insight into the evolution of marine fishes of the Tropical Eastern Pacific

Holacanthus angelfishes are some of the most iconic marine fishes of the Tropical Eastern Pacific (TEP). However, very limited genomic resources currently exist for the genus. In this study we: (i) assembled and annotated the nuclear genome of the King Angelfish (Holacanthus passer), and (ii) examined the demographic history of H. passer in the TEP. We generated 43.8 Gb of ONT and 97.3 Gb Illumina reads representing 75× and 167× coverage, respectively. The final genome assembly size was 583 Mb with a contig N50 of 5.7 Mb, which captured 97.5% of the complete Actinoterygii Benchmarking Universal Single-Copy Orthologs (BUSCOs). Repetitive elements accounted for 5.09% of the genome, and 33,889 protein-coding genes were predicted, of which 22,984 were functionally annotated. Our demographic analysis suggests that population expansions of H. passer occurred prior to the last glacial maximum (LGM) and were more likely shaped by events associated with the closure of the Isthmus of Panama. This result is surprising, given that most rapid population expansions in both freshwater and marine organisms have been reported to occur globally after the LGM. Overall, this annotated genome assembly provides a novel molecular resource to study the evolution of Holacanthus angelfishes, while facilitating research into local adaptation, speciation, and introgression in marine fishes.

Modern day breeding approaches for improvement of castor

Castor (Ricinus communis L.) is an industrially important oil producing crop belongs to Euphorbiaceae family. Castor oil has unique chemical properties make it industrially important crop. It is a member of monotypic genus even though it has ample amount of variability. Using this variability, conventionally many varieties and hybrids have been developed. But, like other crops, the modern and unconventional methods of crop improvement has not fully explored in castor. This article discusses the use of polyploidy induction, distant/wide hybridization and mutation breeding as tools for generating variety. Modern approaches accelerate the speed of crop breeding as an alternative tool. To achieve this goal, molecular markers are employed in breeding to capture the genetic variability through molecular analysis and population structuring. Allele mining is used to trace the evolution of alleles, identify new haplotypes and produce allele specific markers for use in marker aided selection using Genome wide association studies (GWAS) and quantitative trait loci (QTL) mapping. Plant genetic transformation is a rapid and effective mode of castor improvement is also discussed here. The efforts towards developing stable regeneration protocol provide a wide range of utility like embryo rescue in distant crosses, development of somaclonal variation, haploid development using anther culture and callus development for stable genetic transformation has reviewed in this article. Omics has provided intuitions to the molecular mechanisms of (a)biotic stress management in castor along with dissected out the possible genes for improving the yield. Relating genes to traits offers additional scientific inevitability leading to enhancement and sympathetic mechanisms of yield improvement and several stress tolerance.

Genetic Diversity and Virulence Variation of Metarhizium rileyi from Infected Spodoptera frugiperda in Corn Fields

Metarhizium rileyi is an entomopathogenic fungus that naturally infects the larvae of Spodoptera frugiperda, and has biocontrol potential. To explore more natural entomopathogenic fungi resources, a total of 31 strains were isolated from 13 prefectures in Yunnan Province. All the strains were identified using morphology and molecular biology. The genetic diversity of the 31 isolates of M. rileyi was analyzed using inter-simple sequence repeat (ISSR) techniques. Seven primers with good polymorphism were selected, and fifty-four distinct amplification sites were obtained by polymerase chain reaction amplification. Among them, 50 were polymorphic sites, and the percentage of polymorphic sites was 94.44%. The thirty-one strains were divided into eight subpopulations according to the regions. The Nei's gene diversity was 0.2945, and the Shannon information index was 0.4574, indicating that M. rileyi had rich genetic diversity. The average total genetic diversity of the subpopulations in the different regions was 0.2962, the gene diversity within the populations was 0.1931, the genetic differentiation coefficient was 0.3482 (>0.25), and the gene flow was 0.9360 (<1). The individual cluster analysis showed that there was no obvious correlation between the genetic diversity of the strains and their geographical origin, which also indicated that the virulence of the strains was not related to their phylogeny. Thus, the genetic distance of the different populations of M. rileyi in Yunnan Province was not related to the geographical distance. The virulence of those 32 strains against the 3rd-instar larvae of S. frugiperda were varied with the differences in geographical locations. On the 10th day of inoculation, seventeen strains had an insect mortality rate of 70.0%, and seven strains had an insect mortality rate of 100%. The half-lethal times of the M. rileyi SZCY201010, XSBN200920, and MDXZ200803 strains against the S. frugiperda larvae were less than 4 d. Thus, they have the potential to be developed into fungal insecticidal agents.

A computational analysis to evaluate deleterious SNPs of GSK3β, a multifunctional and regulatory protein, for metabolism, wound healing, and migratory processes

This study focused on GSK-3β, a critical serine/threonine kinase with diverse cellular functions. However, there is limited understanding of the impact of non-synonymous single nucleotide polymorphisms (nsSNPs) on its structure and function. Through an exhaustive in-silico investigation 12 harmful nsSNPs were predicted from a pool of 172 acquired from the NCBI dbSNP database using 12 established tools that detects deleterious SNPs. Consistently, these nsSNPs were discovered in locations with high levels of conservation. Notably, the three harmful nsSNPs F67C, A83T, and T138I were situated in the active/binding site of GSK-3β, which may affect the protein's capacity to bind to substrates and other proteins. Molecular dynamics simulations revealed that the F67C and T138I mutants had stable structures, indicating rigidness, whereas the A83T mutant was unstable. Analysis of secondary structures revealed different modifications in all mutant forms, which may affect the stability, functioning, and interactions of the protein. These mutations appear to alter the structural dynamics of GSK-3β, which may have functional ramifications, such as the formation of novel secondary structures and variations in coil-to-helix transitions. In conclusion, this study illuminates the possible structural and functional ramifications of these GSK-3 nsSNPs, revealing how protein compactness, stiffness, and interactions may affect biological activities.

Genetic basis and origin of resistance to acetolactate synthase inhibitors in Amaranthus palmeri from Spain and Italy

BACKGROUND

Amaranthus palmeri is an aggressive annual weed native to the United States, which has become invasive in some European countries. Populations resistant to acetolactate synthase (ALS) inhibitors have been recorded in Spain and Italy, but the evolutionary origin of the resistance traits remains unknown. Bioassays were conducted to identify cross-resistance to ALS inhibitors and a haplotype-based genetic approach was used to elucidate the origin and distribution of resistance in both countries.

RESULTS

Amaranthus palmeri populations were resistant to thifensulfuron-methyl and imazamox, and the 574-Leu mutant ALS allele was found to be the main cause of resistance among them. In two Spanish populations, 376-Glu and 197-Thr mutant ALS alleles were also found. The haplotype analyses revealed the presence of two and four distinct 574-Leu mutant haplotypes in the Italian and Spanish populations, respectively. None was common to both countries, but some mutant haplotypes were shared between geographically close populations or between populations more than 100 km apart. Wide genetic diversity was found in two very close Spanish populations.

CONCLUSION

ALS-resistant A. palmeri populations were introduced to Italy and Spain from outside Europe. Populations from both countries have different evolutionary histories and originate from independent introduction events. ALS resistance then spread over short and long distances by seed dispersal. The higher number and genetic diversity among mutant haplotypes from the Spanish populations indicated recurrent invasions. The implementation of control tactics to limit seed dispersal and the establishment of A. palmeri is recommended in both countries. © 2023 The Authors. Pest Management Science published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

Comparative Analysis of Microsatellite and SNP Markers for Genetic Management of Red Deer

The analysis of population genetic structure and individual multilocus heterozygosity are crucial for wildlife management and conservation. Microsatellite markers have traditionally been used to assess these genetic parameters. However, single-nucleotide polymorphisms (SNPs) are becoming increasingly popular. Our goal here was to determine to what extent SNPs can provide better insights than microsatellites into the overall genetic status and population genetic processes in the species. To this end, we genotyped 210 red deer (Cervus elaphus) in the Spanish wild population with both 11 microsatellites and 31,712 SNPs. We compared parameters related to population genetic structure and individual multilocus heterozygosity obtained with both types of markers. Our results showed correlations between parameters measured using both microsatellites and SNPs, particularly those related to the level of genetic diversity and genetic differentiation. However, we found notably lower precision of microsatellites in measuring the distribution of genetic diversity among individuals. We conclude that microsatellites can be used to monitor the overall genetic status and detect broad patterns in red deer populations. Nevertheless, the greater precision of SNPs in inferring genetic structure and multilocus heterozygosity leads us to encourage scientists and wildlife managers to prioritize their use whenever possible.

Evaluation of Genetic Diversity in Dog Breeds Using Pedigree and Molecular Analysis: A Review

Domestic dogs are important for many economic and social reasons, and they have become a well-known model species for human disease. According to research, dog breeds exhibit significant levels of inbreeding and genetic diversity loss, decreasing the population’s ability to adapt in certain conditions, and indicating the need of conservation strategies. Before the development of molecular markers, pedigree information was used for genetic diversity management. In recent years, genomic tools are frequently applied for accurate estimation of genetic diversity and improved genetic conservation due to incomplete pedigrees and pedigree errors. The most frequently used molecular markers include PCR-based microsatellite markers (STRs) and DNA sequencing-based single-nucleotide polymorphism markers (SNP). The aim of this review was to highlight genetic diversity studies on dog breeds conducted using pedigree and molecular markers, as well as the importance of genetic diversity conservation in increasing the adaptability and survival of dog breed populations.

Associations of MDM2 rs2279744 and TP53 rs1042522 polymorphisms with cervical cancer risk: A meta-analysis and systematic review

Background

Although the association between MDM2 rs2279744 and TP53 rs1042522 polymorphisms and cervical cancer has been reported, the results of its correlation were contradictory. Thus, we conducted a meta-analysis to precisely verify the relationships between MDM2 rs2279744 and TP53 rs1042522 polymorphisms and cervical cancer.

Methods

We thoroughly searched the PubMed, Web of Science, Embase, and Scopus databases for all potential articles from inception to June 2022 and used R Version 4.1.2 and STATA software 12.0 for the meta-analysis. The odds ratios (ORs), 95% confidence intervals (CIs) and 95% prediction intervals (PIs) were calculated to evaluate the associations. Subgroup analyses stratified by ethnicity, source of control, quality score and adjustment were further conducted to assess the relationship between MDM2 rs2279744 and TP53 rs1042522 polymorphisms and cervical cancer.

Results

A total of 30 case-control studies involving 5025 cases and 6680 controls were included. All the included studies were population-based or hospital-based studies. The overall analysis showed that MDM2 rs2279744 polymorphism was closely related to an increased risk of cervical cancer in the recessive model (GG vs GT + TT: OR = 1.602, 95% CI: 1.077-2.383, P = 0.020) and homozygote model (GG vs TT: OR = 1.469, 95% CI: 1.031-2.095, P = 0.033, 95% PI: 0.516-4.184). A significant correlation between TP53 rs1042522 polymorphism and cervical cancer was observed in two models (CC + CG vs GG: OR = 1.759, 95% CI: 1.192-2.596, P = 0.004, 95% PI: 0.474-6.533; GG vs CC: OR = 2.442, 95% CI: 1.433-4.162, P = 0.001, 95% PI: 0.456-13.071).

Conclusions

This meta-analysis revealed that MDM2 SNP309T&gt;G and TP53 rs1042522 C&gt;G polymorphisms were associated with the increased risk of cervical cancer.

The Genome of the Mimosoid Legume Prosopis cineraria, a Desert Tree

The mimosoid legumes are a clade of ~40 genera in the Caesalpinioideae subfamily of the Fabaceae that grow in tropical and subtropical regions. Unlike the better studied Papilionoideae, there are few genomic resources within this legume group. The tree Prosopis cineraria is native to the Near East and Indian subcontinent, where it thrives in very hot desert environments. To develop a tool to better understand desert plant adaptation mechanisms, we sequenced the P. cineraria genome to near-chromosomal assembly, with a total sequence length of ~691 Mb. We predicted 77,579 gene models (76,554 CDS, 361 rRNAs and 664 tRNAs) from the assembled genome, among them 55,325 (~72%) protein-coding genes that were functionally annotated. This genome was found to consist of over 58% repeat sequences, primarily long terminal repeats (LTR-)-retrotransposons. We find an expansion of terpenoid metabolism genes in P. cineraria and its relative Prosopis alba, but not in other legumes. We also observed an amplification of NBS-LRR disease-resistance genes correlated with LTR-associated retrotransposition, and identified 410 retrogenes with an active burst of chimeric retrogene creation that approximately occurred at the same time of divergence of P. cineraria from a common lineage with P. alba~23 Mya. These retrogenes include many biotic defense responses and abiotic stress stimulus responses, as well as the early Nodulin 93 gene. Nodulin 93 gene amplification is consistent with an adaptive response of the species to the low nitrogen in arid desert soil. Consistent with these results, our differentially expressed genes show a tissue specific expression of isoprenoid pathways in shoots, but not in roots, as well as important genes involved in abiotic salt stress in both tissues. Overall, the genome sequence of P. cineraria enriches our understanding of the genomic mechanisms of its disease resistance and abiotic stress tolerance. Thus, it is a very important step in crop and legume improvement.

Prioritization and Meta-analysis of regulatory SNPs identified IL6, TGFB1, TLR9 and MMP7 as significantly associated with cervical cancer

Cervical cancer is a leading women cancer globally with respect to both incidence and mortality. Its increased risk has been linked with HPV infection and genetic variations like single nucleotide polymorphisms (SNPs). Although, studies have been published which evaluates the effect of SNPs in a few candidate genes, however the role of number of regulatory SNPs (rSNPs) in cervical cancer is not available. As literature evidence has shown that non-coding rSNPs are related with increasing cervical cancer risk, we undertook this study to prioritize the important rSNPs and elucidate their role. A search was conducted in PubMed up to December 2020, which led to the identification of 263 articles and 969 SNPs in the non-coding region. These 969 SNPs were analysed through rSNPBase and RegulomeDB, leading to identification of 105 rSNPs. Afterwards, a regulatory module was constructed using protein-protein interaction data and a hub of highly interacting 23 target genes (corresponding to 34 rSNPs) was identified using MCODE. To further understand the mechanism of action of the 34 rSNPs, their transcription factor information with respect to cervical cancer was retrieved. To evaluate the pooled effect of these prioritized polymorphisms in cervical cancer patients, a meta-analysis was performed on 10,537 cases and 11,252 controls from 30 studies corresponding to 8 rSNPs. It led to identification of polymorphisms in IL6 (rs2069837), TGFB1 (rs1800469), TLR9 (rs187084) and MMP7 (rs11568818) which are significantly (p < 0.05) associated with increased cervical cancer risk at the population level. Overall, the study demonstrates that rSNPs targeting immune and inflammatory genes (IL1B, IL6, IL10, IL18, TGFB1, CCR5, CD40, TLR9, and MMP7) are associated with cervical cancer.

Genome-wide SNPs detect fine-scale genetic structure in threatened populations of squirrel glider Petaurus norfolcensis

Australian arboreal mammals are experiencing significant population declines, particularly due to land clearing and resulting habitat fragmentation. The squirrel glider, Petaurus norfolcensis, is a threatened species in New South Wales, with a stronghold population in the Lake Macquarie Local Government Area (LGA) where fragmentation due to urbanization is an ongoing problem for the species conservation. Here we report on the use of squirrel glider mitochondrial (385 bp cytochrome b gene, 70 individuals) and nuclear DNA (6,834 SNPs, 87 individuals) markers to assess their population genetic structure and connectivity across 14 locations sampled in the Lake Macquarie LGA. The mitochondrial DNA sequences detected evidence of a historical genetic bottleneck, while the genome-wide SNPs detected significant population structure in the Lake Macquarie squirrel glider populations at scales as fine as one kilometer. There was no evidence of inbreeding within patches, however there were clear effects of habitat fragmentation and biogeographical barriers on gene flow. A least cost path analysis identified thin linear corridors that have high priority for conservation. These areas should be protected to avoid further isolation of squirrel glider populations and the loss of genetic diversity through genetic drift.

Population Genetic Structure and Biodiversity Conservation of a Relict and Medicinal Subshrub Capparis spinosa in Arid Central Asia

As a Tertiary Tethyan relict, Capparis spinosa is a typical wind-preventing and sand-fixing deciduous subshrub in arid central Asia. Due to its medicinal and energy value, this species is at risk of potential threat from human overexploitation, habitat destruction and resource depletion. In this study, our purpose was to evaluate the conservation strategies of C. spinosa according to its genetic structure characteristics and genetic diversity pattern among 37 natural distributional populations. Based on genomic SNP data generated from dd-RAD sequencing, genetic diversity analysis, principal component analysis, maximum likelihood phylogenetic trees and ADMIXTURE clustering, the significant population structure and differentiation were explored. The results showed the following: (1) Six distinct lineages were identified corresponding to geographic locations, and various levels of genetic diversity existed among the lineages for the natural habitat heterogeneity or human interferences; (2) The lineage divergences were influenced by isolation by distances, vicariance and restricted gene flow under complex topographic and climatic conditions. Finally, for the preservation of the genetic integrity of C. spinosa, we suggest that conservation units should be established corresponding to different geographic groups, and that attention should be paid to isolated and peripheral populations that are experiencing biodiversity loss. Simultaneously, monitoring and reducing anthropogenic disturbances in addition to rationally and sustainably utilizing wild resources would be beneficial to guarantee population resilience and evolutionary potential of this xerophyte in response to future environmental changes.

Population genetic structure of the carrot weevil (Listronotus oregonensis) in North America

Population genetic studies of insect pests enhance our ability to anticipate problems in agroecosystems, such as pest outbreaks, insecticide resistance, or expansions of the host range. This study focuses on geographic distance and host plant selection as potential determinants of genetic differentiation of the carrot weevil Listronotus oregonensis, a major pest of several apiaceous crops in North America. To undertake genetic studies on this species, we assembled the first complete genome sequence for L. oregonensis. Then, we used both haplotype discrimination with mitochondrial DNA (mtDNA) and a genotyping-by-sequencing (GBS) approach to characterize the genetic population structure. A total of 220 individuals were sampled from 17 localities in the provinces of Québec, Ontario, Nova Scotia (Canada), and the state of Ohio (USA). Our results showed significant genetic differences between distant populations across North America, indicating that geographic distance represents an important factor of differentiation for the carrot weevil. Furthermore, the GBS analysis revealed more different clusters than COI analysis between Québec and Nova Scotia populations, suggesting a recent differentiation in the latter province. In contrast, we found no clear evidence of population structure associated with the four cultivated apiaceous plants tested (carrot, parsley, celery, and celeriac) using populations from Québec. This first characterization of the genetic structure of the carrot weevil contributes to a better understanding of the gene flow of the species and helps to adapt local pest management measures to better control this agricultural pest.

Integrating genomic and epidemiologic data to accelerate progress toward schistosomiasis elimination

The global community has adopted ambitious goals to eliminate schistosomiasis as a public health problem, and new tools are needed to achieve them. Mass drug administration programs, for example, have reduced the burden of schistosomiasis, but the identification of hotspots of persistent and reemergent transmission threaten progress toward elimination and underscore the need to couple treatment with interventions that reduce transmission. Recent advances in DNA sequencing technologies make whole-genome sequencing a valuable and increasingly feasible option for population-based studies of complex parasites such as schistosomes. Here, we focus on leveraging genomic data to tailor interventions to distinct social and ecological circumstances. We consider two priority questions that can be addressed by integrating epidemiological, ecological, and genomic information: (1) how often do non-human host species contribute to human schistosome infection? and (2) what is the importance of locally acquired versus imported infections in driving transmission at different stages of elimination? These questions address processes that can undermine control programs, especially those that rely heavily on treatment with praziquantel. Until recently, these questions were difficult to answer with sufficient precision to inform public health decision-making. We review the literature related to these questions and discuss how whole-genome approaches can identify the geographic and taxonomic sources of infection, and how such information can inform context-specific efforts that advance schistosomiasis control efforts and minimize the risk of reemergence.

Suitability of GWAS as a Tool to Discover SNPs Associated with Tick Resistance in Cattle: A Review

Understanding the biological mechanisms underlying tick resistance in cattle holds the potential to facilitate genetic improvement through selective breeding. Genome wide association studies (GWAS) are popular in research on unraveling genetic determinants underlying complex traits such as tick resistance. To date, various studies have been published on single nucleotide polymorphisms (SNPs) associated with tick resistance in cattle. The discovery of SNPs related to tick resistance has led to the mapping of associated candidate genes. Despite the success of these studies, information on genetic determinants associated with tick resistance in cattle is still limited. This warrants the need for more studies to be conducted. In Africa, the cost of genotyping is still relatively expensive; thus, conducting GWAS is a challenge, as the minimum number of animals recommended cannot be genotyped. These population size and genotype cost challenges may be overcome through the establishment of collaborations. Thus, the current review discusses GWAS as a tool to uncover SNPs associated with tick resistance, by focusing on the study design, association analysis, factors influencing the success of GWAS, and the progress on cattle tick resistance studies.

Relict inland mangrove ecosystem reveals Last Interglacial sea levels

With geological sea-level fluctuations driven by climate change, the distribution of mangrove forests has expanded and contracted through time. We studied an inland, isolated mangrove forest located 170 km away from the nearest coastline in the interior of the rainforests of the Yucatan Peninsula (Mexico). Combining multiple lines of evidence, we demonstrate that this extant forest is a relict from a past, warmer world when relative sea levels were 6 to 9 m higher than at the present. Our finding highlights the extensive landscape impacts of past climate change on the world’s coastline and opens opportunities to better understand future scenarios of relative sea level rise.

Climatic oscillations during the Pleistocene played a major role in shaping the spatial distribution and demographic dynamics of Earth's biota, including our own species. The Last Interglacial (LIG) or Eemian Period (ca. 130 to 115 thousand years B.P.) was particularly influential because this period of peak warmth led to the retreat of all ice sheets with concomitant changes in global sea level. The impact of these strong environmental changes on the spatial distribution of marine and terrestrial ecosystems was severe as revealed by fossil data and paleogeographic modeling. Here, we report the occurrence of an extant, inland mangrove ecosystem and demonstrate that it is a relict of the LIG. This ecosystem is currently confined to the banks of the freshwater San Pedro Mártir River in the interior of the Mexico–Guatemala El Petén rainforests, 170 km away from the nearest ocean coast but showing the plant composition and physiognomy typical of a coastal lagoon ecosystem. Integrating genomic, geologic, and floristic data with sea level modeling, we present evidence that this inland ecosystem reached its current location during the LIG and has persisted there in isolation ever since the oceans receded during the Wisconsin glaciation. Our study provides a snapshot of the Pleistocene peak warmth and reveals biotic evidence that sea levels substantially influenced landscapes and species ranges in the tropics during this period.

Effects of single nucleotide polymorphism ascertainment on population structure inferences

Single nucleotide polymorphism (SNP) data are widely used in research on natural populations. Although they are useful, SNP genotyping data are known to contain bias, normally referred to as ascertainment bias, because they are conditioned by already confirmed variants. This bias is introduced during the genotyping process, including the selection of populations for novel SNP discovery and the number of individuals involved in the discovery panel and selection of SNP markers. It is widely recognized that ascertainment bias can cause inaccurate inferences in population genetics and several methods to address these bias issues have been proposed. However, especially in natural populations, it is not always possible to apply an ideal ascertainment scheme because natural populations tend to have complex structures and histories. In addition, it was not fully assessed if ascertainment bias has the same effect on different types of population structure. Here, we examine the effects of bias produced during the selection of population for SNP discovery and consequent SNP marker selection processes under three demographic models: the island, stepping-stone, and population split models. Results show that site frequency spectra and summary statistics contain biases that depend on the joint effect of population structure and ascertainment schemes. Additionally, population structure inferences are also affected by ascertainment bias. Based on these results, it is recommended to evaluate the validity of the ascertainment strategy prior to the actual typing process because the direction and extent of ascertainment bias vary depending on several factors.

500 years of breeding in the Carthusian Strain of Pura Raza Español horse: An evolutional analysis using genealogical and genomic data

The Carthusian horse is a Pura Raza Español (PRE) strain (CS), bred as a closed population since its creation more than 500 years ago. The aim of this study was to analyse for the first time its population structure and situation of variability combining both genealogical (GEL) and genomic (GEN) data. The GEL data comprised 348,429 pedigree records (56,105 CS horses), while the GEN analysis included the high-density genotypes (670,804 SNPs) of 287 horses. Pedigree completeness demonstrated its accuracy, showing a good correlation of GEL (F) and GEN (F_ROH ) inbreeding coefficient in the case of PRE subpopulations partially related and non-related to Carthusian strain (0.68) but a lower value in the 100% Carthusian horses (0.42), due to the high weight of founders not detected by GEL analysis. GEN (PCA, AMOVA, and Admixture) and GEL analysis showed a good differentiation of subpopulations, but also a high level of introgression of the CS in the breed during past decades. A recent change in this trend was noteworthy, with a considerable reduction in CS variability and a genetic bottleneck (effective population sizes of 31.57 and 30.20 in GEL and GEN analysis, respectively, in last generation). The PRE has maintained its variability, and a considerable difference in estimated N_e by GEL (60.77) and GEN (188.0) data was observed. Using two sources of complementary information, it was found the existence of an ancient PRE strain with a unique genetic landmark, practically free from the influence of other equine populations.

Genome-wide SNP analysis to assess the genetic population structure and diversity of Acrocomia species

Acrocomia (Arecaceae) is a genus widely distributed in tropical and subtropical America that has been achieving economic interest due to the great potential of oil production of some of its species. In particular A. aculeata, due to its vocation to supply oil with the same productive capacity as the oil palm (Elaeis guineenses) even in areas with water deficit. Although eight species are recognized in the genus, the taxonomic classification based on morphology and geographic distribution is still controversial. Knowledge about the genetic diversity and population structure of the species is limited, which has limited the understanding of the genetic relationships and the orientation of management, conservation, and genetic improvement activities of species of the genus. In the present study, we analyzed the genomic diversity and population structure of Acrocomia genus, including 172 samples from seven species, with a focus on A. aculeata with 117 samples covering a wide geographical area of occurrence of the species, using Single Nucleotide Polymorphism (SNP) markers originated from Genotyping By Sequencing (GBS).The genetic structure of the Acrocomia species were partially congruent with the current taxonomic classification based on morphological characters, recovering the separation of the species A. aculeata, A. totai, A. crispa and A. intumescens as distinct taxonomic groups. However, the species A. media was attributed to the cluster of A. aculeata while A. hassleri and A. glauscescens were grouped together with A. totai. The species that showed the highest and lowest genetic diversity were A. totai and A. media, respectively. When analyzed separately, the species A. aculeata showed a strong genetic structure, forming two genetic groups, the first represented mainly by genotypes from Brazil and the second by accessions from Central and North American countries. Greater genetic diversity was found in Brazil when compared to the other countries. Our results on the genetic diversity of the genus are unprecedented, as is also establishes new insights on the genomic relationships between Acrocomia species. It is also the first study to provide a more global view of the genomic diversity of A. aculeata. We also highlight the applicability of genomic data as a reference for future studies on genetic diversity, taxonomy, evolution and phylogeny of the Acrocomia genus, as well as to support strategies for the conservation, exploration and breeding of Acrocomia species and in particular A. aculeata.

Parentage Analysis in Giant Grouper (Epinephelus lanceolatus) Using Microsatellite and SNP Markers from Genotyping-by-Sequencing Data

Pedigree information is necessary for the maintenance of diversity for wild and captive populations. Accurate pedigree is determined by molecular marker-based parentage analysis, which may be influenced by the polymorphism and number of markers, integrity of samples, relatedness of parents, or different analysis programs. Here, we described the first development of 208 single nucleotide polymorphisms (SNPs) and 11 microsatellites for giant grouper (Epinephelus lanceolatus) taking advantage of Genotyping-by-sequencing (GBS), and compared the power of SNPs and microsatellites for parentage and relatedness analysis, based on a mixed family composed of 4 candidate females, 4 candidate males and 289 offspring. CERVUS, PAPA and COLONY were used for mutually verification. We found that SNPs had a better potential for relatedness estimation, exclusion of non-parentage and individual identification than microsatellites, and > 98% accuracy of parentage assignment could be achieved by 100 polymorphic SNPs (MAF cut-off < 0.4) or 10 polymorphic microsatellites (mean H_o = 0.821, mean PIC = 0.651). This study provides a reference for the development of molecular markers for parentage analysis taking advantage of next-generation sequencing, and contributes to the molecular breeding, fishery management and population conservation.

Novel Structural Variation and Evolutionary Characteristics of Chloroplast tRNA in Gossypium Plants

Cotton is one of the most important fiber and oil crops in the world. Chloroplast genomes harbor their own genetic materials and are considered to be highly conserved. Transfer RNAs (tRNAs) act as "bridges" in protein synthesis by carrying amino acids. Currently, the variation and evolutionary characteristics of tRNAs in the cotton chloroplast genome are poorly understood. Here, we analyzed the structural variation and evolution of chloroplast tRNA (cp tRNA) based on eight diploid and two allotetraploid cotton species. We also investigated the nucleotide evolution of chloroplast genomes in cotton species. We found that cp tRNAs in cotton encoded 36 or 37 tRNAs, and 28 or 29 anti-codon types with lengths ranging from 60 to 93 nucleotides. Cotton chloroplast tRNA sequences possessed specific conservation and, in particular, the Ψ-loop contained the conserved U-U-C-X3-U. The cp tRNAs of Gossypium L. contained introns, and cp tRNA^Ile contained the anti-codon (C-A-U), which was generally the anti-codon of tRNA^Met. The transition and transversion analyses showed that cp tRNAs in cotton species were iso-acceptor specific and had undergone unequal rates of evolution. The intergenic region was more variable than coding regions, and non-synonymous mutations have been fixed in cotton cp genomes. On the other hand, phylogeny analyses indicated that cp tRNAs of cotton were derived from several inferred ancestors with greater gene duplications. This study provides new insights into the structural variation and evolution of chloroplast tRNAs in cotton plants. Our findings could contribute to understanding the detailed characteristics and evolutionary variation of the tRNA family.

Unboxing mutations: Connecting mutation types with evolutionary consequences

A key step in understanding the genetic basis of different evolutionary outcomes (e.g., adaptation) is to determine the roles played by different mutation types (e.g., SNPs, translocations and inversions). To do this we must simultaneously consider different mutation types in an evolutionary framework. Here, we propose a research framework that directly utilizes the most important characteristics of mutations, their population genetic effects, to determine their relative evolutionary significance in a given scenario. We review known population genetic effects of different mutation types and show how these may be connected to different evolutionary outcomes. We provide examples of how to implement this framework and pinpoint areas where more data, theory and synthesis are needed. Linking experimental and theoretical approaches to examine different mutation types simultaneously is a critical step towards understanding their evolutionary significance.

Soybean anthracnose caused by Colletotrichum species: Current status and future prospects

Soybean (Glycine max) is one of the most important cultivated plants worldwide as a source of protein-rich foods and animal feeds. Anthracnose, caused by different lineages of the hemibiotrophic fungus Colletotrichum, is one of the main limiting factors to soybean production. Losses due to anthracnose have been neglected, but their impact may threaten up to 50% of the grain production.

TAXONOMY

While C. truncatum is considered the main species associated with soybean anthracnose, recently other species have been reported as pathogenic on this host. Until now, it has not been clear whether the association of new Colletotrichum species with the disease is related to emerging species or whether it is due to the undergoing changes in the taxonomy of the genus.

DISEASE SYMPTOMS

Typical anthracnose symptoms are pre- and postemergence damping-off; dark, depressed, and irregular spots on cotyledons, stems, petioles, and pods; and necrotic laminar veins on leaves that can result in premature defoliation. Symptoms may evolve to pod rot, immature opening of pods, and premature germination of grains.

CHALLENGES

As accurate species identification of the causal agent is decisive for disease control and prevention, in this work we review the taxonomic designation of Colletotrichum isolated from soybean to understand which lineages are pathogenic on this host. We also present a comprehensive literature review of soybean anthracnose, focusing on distribution, symptomatology, epidemiology, disease management, identification, and diagnosis. We consider the knowledge emerging from population studies and comparative genomics of Colletotrichum spp. associated with soybean providing future perspectives in the identification of molecular factors involved in the pathogenicity process.

USEFUL WEBSITE

Updates on Colletotrichum can be found at http://www.colletotrichum.org/. All available Colletotrichum genomes on GenBank can be viewed at http://www.colletotrichum.org/genomics/.

Development of Breeder-Friendly KASP Markers for Low Concentration of Kunitz Trypsin Inhibitor in Soybean Seeds

Trypsin inhibitors (TI), a common anti-nutritional factor in soybean, prevent animals' protein digestibility reducing animal growth performance. No commercial soybean cultivars with low or null concentration of TI are available. The availability of a high throughput genotyping assay will be beneficial to incorporate the low TI trait into elite breeding lines. The aim of this study is to develop and validate a breeder friendly Kompetitive Allele Specific PCR (KASP) assay linked to low Kunitz trypsin inhibitor (KTI) in soybean seeds. A total of 200 F3:5 lines derived from PI 547656 (low KTI) X Glenn (normal KTI) were genotyped using the BARCSoySNP6K_v2 Beadchip. F3:4 and F3:5 lines were grown in Blacksburg and Orange, Virginia in three years, respectively, and were measured for KTI content using a quantitative HPLC method. We identified three SNP markers tightly linked to the major QTL associated to low KTI in the mapping population. Based on these SNPs, we developed and validated the KASP assays in a set of 93 diverse germplasm accessions. The marker Gm08_44814503 has 86% selection efficiency for the accessions with low KTI and could be used in marker assisted breeding to facilitate the incorporation of low KTI content in soybean seeds.

Diversity Patterns of Bermuda Grass along Latitudinal Gradients at Different Temperatures in Southeastern China

Cynodon dactylon (L.) Pers. (common Bermuda grass) has a limited capacity to grow at low temperatures, which limits its geographical range. Exploring its evolutionary relationship across different environmental gradients is necessary to understand the effects of temperature change on the genetics of common Bermuda grass. In this study, high-throughput transcriptome sequencing was performed on 137 samples of C. dactylon from 16 latitudinal gradients to explore the differential molecular markers and analyze genetic diversity and structure along latitudinal gradients at different temperatures. We primarily sampled more high-quality single nucleotide polymorphisms (SNPs) from populations at lower and middle latitudes. Greater intraspecific genetic variation at each level of temperature treatment could be due to factors such as wind pollination and asexual breeding. Populations of C. dactylon at high latitudes differed from populations at middle and low latitudes, which was supported by a principal component analysis (PCA) and genetic structure analysis, performed at different temperatures. We observed more genetic variation for low-latitude populations at 5 °C, according to an analysis of three phylogenetic trees at different temperature levels, suggesting that low temperatures affected samples with low cold resistance. Based on the results of phylogenetic analysis, we found that samples from high latitudes evolved earlier than most samples at low latitudes. The results provide a comprehensive understanding of the evolutionary phenomenon of landscape genetics, laying the groundwork for future structural and comparative genomic studies of C. dactylon.

Modern Techniques in Colorado Potato Beetle (Leptinotarsa decemlineata Say) Control and Resistance Management: History Review and Future Perspectives

Simple Summary

The Colorado potato beetle (CPB) is one of the most important potato pest worldwide. It is native to U.S. but during the 20th century it has dispersed through Europe, Asia and western China. It continues to expand in an east and southeast direction. Damages are caused by larvae and adults. Their feeding on potato plant leaves can cause complete defoliation and lead to a large yield loss. After the long period of using only chemical control measures, the emergence of resistance increased and some new and different methods come to the fore. The main focus of this review is on new approaches to the old CPB control problem. We describe the use of Bacillus thuringiensis and RNA interference (RNAi) as possible solutions for the future in CPB management. RNAi has proven successful in controlling many pests and shows great potential for CPB control. Better understanding of the mechanisms that affect efficiency will enable the development of this technology and boost potential of RNAi to become part of integrated plant protection in the future. We described also the possibility of using single nucleotide polymorphisms (SNPs) as a way to go deeper into our understanding of resistance and how it influences genotypes.

Abstract

Colorado potato beetle, CPB (Leptinotarsa decemlineata Say), is one of the most important pests of the potato globally. Larvae and adults can cause complete defoliation of potato plant leaves and can lead to a large yield loss. The insect has been successfully suppressed by insecticides; however, over time, has developed resistance to insecticides from various chemical groups, and its once successful control has diminished. The number of available active chemical control substances is decreasing with the process of testing, and registering new products on the market are time-consuming and expensive, with the possibility of resistance ever present. All of these concerns have led to the search for new methods to control CPB and efficient tools to assist with the detection of resistant variants and monitoring of resistant populations. Current strategies that may aid in slowing resistance include gene silencing by RNA interference (RNAi). RNAi, besides providing an efficient tool for gene functional studies, represents a safe, efficient, and eco-friendly strategy for CPB control. Genetically modified (GM) crops that produce the toxins of Bacillus thuringiensis (Bt) have many advantages over agro-technical, mechanical, biological, and chemical measures. However, pest resistance that may occur and public acceptance of GM modified food crops are the main problems associated with Bt crops. Recent developments in the speed, cost, and accuracy of next generation sequencing are revolutionizing the discovery of single nucleotide polymorphisms (SNPs) and field of population genomics. There is a need for effective resistance monitoring programs that are capable of the early detection of resistance and successful implementation of integrated resistance management (IRM). The main focus of this review is on new technologies for CPB control (RNAi) and tools (SNPs) for detection of resistant CPB populations.

Genetic Differentiation of a New World Screwworm Fly Population from Uruguay Detected by SNPs, Mitochondrial DNA and Microsatellites in Two Consecutive Years

The New World screwworm (NWS) fly, Cochliomyia hominivorax (Diptera: Calliphoridae), is an economically important ectoparasite currently distributed in South America and in the Caribbean basin. The successful eradication of this species in USA, Mexico and continental Central America was achieved by a control program based on the sterile insect technique (SIT). In order to implement a genetic control strategy over the NWS fly's current area of occurrence, first, it is necessary to understand the species dynamics and population structure. In order to address this objective, the spatial genetic structure of the NWS fly was previously reported in South America based on different genetic markers; however, to date, no study has investigated temporal changes in the genetic composition of its populations. In the current study, the temporal genetic structure of a NWS fly population from Uruguay was investigated through two consecutive samplings from the same locality over an interval of approximately 18 generations. The genetic structure was accessed with neutral and under selection SNPs obtained with genotyping-by-sequencing. The results gathered with these data were compared to estimates achieved with mitochondrial DNA sequences and eight microsatellite markers. Temporal changes in the genetic composition were revealed by all three molecular markers, which may be attributed to seasonal changes in the NWS fly's southern distribution. SNPs were employed for the first time for estimating the genetic structure in a NWS fly population; these results provide new clues and perspectives on its population genetic structure. This approach could have significant implications for the planning and implementation of management programs.

Analysis of Comprehensive Pharmacogenomic Profiling of VIP Variants Among the Genetically Isolated Chechen Subpopulation from Jordan

Background

Profiling rare variants in isolated populations can significantly clarify and understand the development of a clinically relevant process. Therefore, leading to a better identifying novel targeted treatment.

Objective

This study aimed to determine the allele frequencies of 56 single nucleotide polymorphisms (SNPs) within several important pharmacogenes.

Methods

This study consisted of 166 unrelated subjects from a genetically isolated group (Chechen) who were living in Jordan. In this study, the distribution of the variants among Chechen was compared to other ethnic groups available at two databases (Genome 1000 and (ExAC)). The frequency of genotypes and alleles was calculated and tested using the chi-square test and the Hardy–Weinberg equilibrium equation (HWE).

Results

Our results revealed that the distribution of allele frequencies within different pharmacogenes among Chechen showed different similarities with other populations. The CEU and TSI showed the highest resemblance with the Chechen population (75% similarity), in contrast to LWK which had the lowest similarity (30%).

Conclusion

This study sheds light on clinically relevant SNPs to enhance medical research and apply pharmacogenomics in clinical settings.

A GT-seq panel for walleye (Sander vitreus) provides important insights for efficient development and implementation of amplicon panels in non-model organisms

Targeted amplicon sequencing methods, such as genotyping-in-thousands by sequencing (GT-seq), facilitate rapid, accurate, and cost-effective analysis of hundreds of genetic loci in thousands of individuals. Development of GT-seq panels is nontrivial, but studies describing trade-offs associated with different steps of GT-seq panel development are rare. Here, we construct a dual-purpose GT-seq panel for walleye (Sander vitreus), discuss trade-offs associated with different development and genotyping approaches, and provide suggestions for researchers constructing their own GT-seq panels. Our GT-seq panel was developed using an ascertainment set consisting of restriction site-associated DNA data from 954 individuals sampled from 23 populations in Minnesota and Wisconsin, USA. We conducted simulations to test the utility of all loci for parentage analysis and genetic stock identification and designed 600 primer pairs to maximize joint accuracy for these analyses. We then performed three rounds of primer optimization to remove loci that overamplified and our final panel consisted of 436 loci. We also explored different approaches for DNA extraction, multiplexed polymerase chain reaction (PCR) amplification, and cleanup steps during the GT-seq process and discovered the following: (i) inexpensive Chelex extractions performed well for genotyping; (ii) the exonuclease I and shrimp alkaline phosphatase (ExoSAP) procedure included in some current protocols did not improve results substantially and was probably unnecessary; and (iii) it was possible to PCR amplify panels separately and combine them prior to adapter ligation. Well-optimized GT-seq panels are valuable resources for conservation genetics and our findings and suggestions should aid in their construction in myriad taxa.

Sexual reproduction and genetic polymorphism within the cosmopolitan marine diatom Pseudo-nitzschia pungens

Different clades belonging to the cosmopolitan marine diatom Pseudo-nitzschia pungens appear to be present in different oceanic environments, however, a ‘hybrid zone’, where populations of different clades interbreed, has also been reported. Many studies have investigated the sexual reproduction of P. pungens, focused on morphology and life cycle, rather than the role of sexual reproduction in mixing the genomes of their parents. We carried out crossing experiments to determine the sexual compatibility/incompatibility between different clades of P. pungens, and examined the genetic polymorphism in the ITS2 region. Sexual reproduction did not occur only between clades II and III under any of experimental temperature conditions. Four offspring strains were established between clade I and III successfully. Strains established from offspring were found interbreed with other offspring strains as well as viable with their parental strains. We confirmed the hybrid sequence patterns between clades I and III and found novel sequence types including polymorphic single nucleotide polymorphisms (SNPs) in the offspring strains. Our results implicate that gene exchange and mixing between different clades are still possible, and that sexual reproduction is a significant ecological strategy to maintain the genetic diversity within this diatom species.

Coping with Pleistocene climatic fluctuations: Demographic responses in remote endemic reef fishes

Elucidating demographic history during the settlement of ecological communities is crucial for properly inferring the mechanisms that shape patterns of species diversity and their persistence through time. Here, we used genomic data and coalescent-based approaches to elucidate for the first time the demographic dynamics associated with the settlement by endemic reef fish fauna of one of the most remote peripheral islands of the Pacific Ocean, Rapa Nui (Easter Island). We compared the demographic history of nine endemic species in order to explore their demographic responses to Pleistocene climatic fluctuations. We found that species endemic to Rapa Nui share a common demographic history, as signatures of population expansions were retrieved for almost all of the species studied here, and synchronous demographic expansions initiated during the last glacial period were recovered for more than half of the studied species. These results suggest that eustatic fluctuations associated with Milankovitch cycles have played a central role in species demographic histories and in the final stage of the community assembly of many Rapa Nui reef fishes. Specifically, sea level lowstands resulted in the maximum reef habitat extension for Rapa Nui endemic species; we discuss the potential role of seamounts in allowing endemic species to cope with Pleistocene climatic fluctuations, and we highlight the importance of local historical processes over regional ones. Overall, our results shed light on the mechanisms by which endemism arises and is maintained in peripheral reef fish fauna.

An empirical comparison of population genetic analyses using microsatellite and SNP data for a species of conservation concern

BACKGROUND

Use of genomic tools to characterize wildlife populations has increased in recent years. In the past, genetic characterization has been accomplished with more traditional genetic tools (e.g., microsatellites). The explosion of genomic methods and the subsequent creation of large SNP datasets has led to the promise of increased precision in population genetic parameter estimates and identification of demographically and evolutionarily independent groups, as well as questions about the future usefulness of the more traditional genetic tools. At present, few empirical comparisons of population genetic parameters and clustering analyses performed with microsatellites and SNPs have been conducted.

RESULTS

Here we used microsatellite and SNP data generated from Gunnison sage-grouse (Centrocercus minimus) samples to evaluate concordance of the results obtained from each dataset for common metrics of genetic diversity (HO, HE, FIS, AR) and differentiation (FST, GST, DJost). Additionally, we evaluated clustering of individuals using putatively neutral (SNPs and microsatellites), putatively adaptive, and a combined dataset of putatively neutral and adaptive loci. We took particular interest in the conservation implications of any differences. Generally, we found high concordance between microsatellites and SNPs for HE, FIS, AR, and all differentiation estimates. Although there was strong correlation between metrics from SNPs and microsatellites, the magnitude of the diversity and differentiation metrics were quite different in some cases. Clustering analyses also showed similar patterns, though SNP data was able to cluster individuals into more distinct groups. Importantly, clustering analyses with SNP data suggest strong demographic independence among the six distinct populations of Gunnison sage-grouse with some indication of evolutionary independence in two or three populations; a finding that was not revealed by microsatellite data.

CONCLUSION

We demonstrate that SNPs have three main advantages over microsatellites: more precise estimates of population-level diversity, higher power to identify groups in clustering methods, and the ability to consider local adaptation. This study adds to a growing body of work comparing the use of SNPs and microsatellites to evaluate genetic diversity and differentiation for a species of conservation concern with relatively high population structure and using the most common method of obtaining SNP genotypes for non-model organisms.

Sea Turtle Population Genomic Discovery: Global and Locus-Specific Signatures of Polymorphism, Selection, and Adaptive Potential

In the era of genomics, single-nucleotide polymorphisms (SNPs) have become a preferred molecular marker to study signatures of selection and population structure and to enable improved population monitoring and conservation of vulnerable populations. We apply a SNP calling pipeline to assess population differentiation, visualize linkage disequilibrium, and identify loci with sex-specific genotypes of 45 loggerhead sea turtles (Caretta caretta) sampled from the southeastern coast of the United States, including 42 individuals experimentally confirmed for gonadal sex. By performing reference-based SNP calling in independent runs of Stacks, 3,901–6,998 SNPs and up to 30 potentially sex-specific genotypes were identified. Up to 68 pairs of loci were found to be in complete linkage disequilibrium, potentially indicating regions of natural selection and adaptive evolution. This study provides a valuable SNP diagnostic workflow and a large body of new biomarkers for guiding targeted studies of sea turtle genome evolution and for managing legally protected nonmodel iconic species that have high economic and ecological importance but limited genomic resources.

Pest Management Challenges and Control Practices in Codling Moth: A Review

The codling moth, Cydia pomonella L., is a serious insect pest in pome fruit production worldwide with a preference for apple. The pest is known for having developed resistance to several chemical groups of insecticides, making its control difficult. The control and management of the codling moth is often hindered by a lack of understanding about its biology and ecology, including aspects of its population genetics. This review summarizes the information about the origin and biology of the codling moth, describes the mechanisms of resistance in this pest, and provides an overview of current research of resistant pest populations and genetic research both in Europe and globally. The main focus of this review is on non-pesticide control measures and anti-resistance strategies which help to reduce the number of chemical pesticides used and their residues on food and the local environment. Regular monitoring for insecticide resistance is essential for proactive management to mitigate potential insecticide resistance. Here we describe techniques for the detection of resistant variants and possibilities for monitoring resistance populations. Also, we present our present work on developing new methods to maintain effective control using appropriate integrated resistance management (IRM) strategies for this economically important perennial pest.

Detecting sequence variants in clinically important protozoan parasites

Second and third generation sequencing methods are crucial for population genetic studies, and variant detection is a popular approach for exploiting this sequence data. While mini- and microsatellites are historically useful markers for studying important Protozoa such as Toxoplasma and Plasmodium spp., detecting non-repetitive variants such as those found in genes can be fundamental to investigating a pathogen's biology. These variants, namely single nucleotide polymorphisms and insertions and deletions, can help elucidate the genetic basis of an organism's pathogenicity, identify selective pressures, and resolve phylogenetic relationships. They also have the added benefit of possessing a comparatively low mutation rate, which contributes to their stability. However, there is a plethora of variant analysis tools with nuanced pipelines and conflicting recommendations for best practise, which can be confounding. This lack of standardisation means that variant analysis requires careful parameter optimisation, an understanding of its limitations, and the availability of high quality data. This review explores the value of variant detection when applied to non-model organisms such as clinically important protozoan pathogens. The limitations of current methods are discussed, including special considerations that require the end-users' attention to ensure that the results generated are reproducible, and the biological conclusions drawn are valid.

Development of a large SNPs resource and a low-density SNP array for brown trout (Salmo trutta) population genetics

Background

The brown trout (Salmo trutta) is an economically and ecologically important species for which population genetic monitoring is frequently performed. The most commonly used genetic markers for this species are microsatellites and mitochondrial markers that lack replicability among laboratories, and a large genome coverage. An alternative that may be particularly efficient and universal is the development of small to large panels of Single Nucleotide Polymorphism markers (SNPs). Here, we used Restriction site Associated DNA sequences (RADs) markers to identify a set of 12,204 informative SNPs positioned on the brown trout linkage map and suitable for population genetics studies. Then, we used this novel resource to develop a cost-effective array of 192 SNPs (96 × 2) evenly spread on this map. This array was tested for genotyping success in five independent rivers occupied by two main brown trout evolutionary lineages (Atlantic -AT- and Mediterranean -ME-) on a total of 1862 individuals. Moreover, inference of admixture rate with domestic strains and population differentiation were assessed for a small river system (the Taurion River, 190 individuals) and results were compared to a panel of 13 microsatellites.

Results

A high genotyping success was observed for all rivers (< 1% of non-genotyped loci per individual), although some initially used SNP failed to be amplified, probably because of mutations in primers, and were replaced. These SNPs permitted to identify patterns of isolation-by-distance for some rivers. Finally, we found that microsatellite and SNP markers yielded very similar patterns for population differentiation and admixture assessments, with SNPs having better ability to detect introgression and differentiation.

Conclusions

The novel resources provided here opens new perspectives for universality and genome-wide studies in brown trout populations.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5958-9) contains supplementary material, which is available to authorized users.

Genotypic Characterization of Torymus sinensis (Hymenoptera: Torymidae) After Its Introduction in Tuscany (Italy) for the Biological Control of Dryocosmus kuriphilus (Hymenoptera: Cynipidae)

Torymus sinensis Kamijo (Hymenoptera: Torymidae) is an alien parasitoid that is used in many areas of the world for biological control the Asian chestnut gall wasp, Dryocosmus kuriphilus Yasumatsu (Hymenoptera: Cynipidae). In Italy, this parasitoid was imported from Japan in 2003 and subsequently multiplied and released throughout the country. In this study, a phylogenetic investigation was carried out on insects from three different sites in northern Tuscany (Italy). Moreover, the possible hybridization between T. sinensis and some native Torymus species was evaluated. The conserved region 18S rRNA gene and the hypervariable ITS2 (Internal Transcribed Spacer 2) region of the ribosomal cistrone were selected as molecular markers. Sequencing the amplified products, after cloning, ruled out any hybridization between T. sinensis and the native Torymus species, and also confirmed the presence of two haplotypes for the Tuscan population of T. sinensis both for the region of the 18S rRNA gene as well as for the ITS2 region. These results confirm that the environmental impact of the alien parasitoid T. sinensis in the study site is acceptable, although an extensive and repeated monitoring would be desirable.

Alternated selection mechanisms maintain adaptive diversity in different demographic scenarios of a large carnivore

Background

Different population trajectories are expected to impact the signature of neutral and adaptive processes at multiple levels, challenging the assessment of the relative roles of different microevolutionary forces. Here, we integrate adaptive and neutral variability patterns to disentangle how adaptive diversity is driven under different demographic scenarios within the Iberian wolf (Canis lupus) range. We studied the persistent, the expanding and a small, isolated group within the Iberian wolf population, using 3 MHC class II genes (DRB1, DQA1, and DQB1), which diversity was compared with 39 microsatellite loci.

Results

Both the persistent and the expanding groups show evidence of balancing selection, revealed by a significant departure from neutrality at MHC loci, significant higher observed and expected heterozygosity and lower differentiation at MHC than at neutral loci, and signs of positive selection. However, despite exhibiting a significantly higher genetic diversity than the isolated group, the persistent group did not show significant excess of MHC heterozygotes. The expanding group, while showing a similar level of genetic diversity than the persistent group, displays by contrast a significant excess of MHC heterozygotes, which is compatible with the heterozygote advantage mechanism. Results are not clear regarding the role of drift and selection in the isolated group due to the small size of this population. Although diversity indices of MHC loci correspond to neutral expectations in the isolated group, accelerated MHC divergence, revealed by a higher differentiation at MHC than neutral loci, may indicate diversifying selection.

Conclusion

Different selective pressures were observed in the three different demographic scenarios, which are possibly driven by different selection mechanisms to maintain adaptive diversity.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1420-5) contains supplementary material, which is available to authorized users.

Female choice for related males in wild red-backed toadlets (Pseudophryne coriacea)

Mate choice for genetic benefits is assumed to be widespread in nature, yet very few studies have comprehensively examined relationships between female mate choice and male genetic quality in wild populations. Here, we use exhaustive sampling and single nucleotide polymorphisms to provide a partial test of the “good genes as heterozygosity” hypothesis and the “genetic compatibility” hypothesis in an entire population of terrestrial breeding red-backed toadlets, Pseudophryne coriacea. We found that successful males did not display higher heterozygosity, despite a positive relationship between male heterozygosity and offspring heterozygosity. Rather, in the larger of 2 breeding events, we found that successful males were more genetically similar to their mate than expected under random mating, indicating that females can use pre- or post-copulatory mate choice mechanisms to bias paternity toward more related males. These findings provide no support for the good genes as heterozygosity hypothesis but lend support to the genetic compatibility hypothesis. A complete test of this hypothesis will now require evaluating how parental genetic similarity impacts offspring fitness. Terrestrial toadlets show a high degree of site fidelity, high levels of genetic structuring between populations, and frequently hybridize with sister species. As such, female mate choice for related males may be an adaptive strategy to reduce outbreeding depression. Our findings provide the first population-wide evidence for non-random preferential inbreeding in a wild amphibian. We argue that such reproductive patterns may be common in amphibians because extreme genetic differentiation within meta-populations creates an inherently high risk of outbreeding depression.