Academics' Beliefs Regarding the Importance of Curriculum Internationalization in the Ethiopian Research Universities.

This study set out to investigate the beliefs of academics regarding the importance of curriculum internationalization in Ethiopian research universities. A convergent-parallel mixed-method design was used to achieve this purpose. Quantitative data was gathered from 415 randomly chosen academics taken from a sample frame of 6808 through a questionnaire. Using a semi-structured interview, 16 purposefully chosen individuals provided qualitative data. The questionnaire's construct, face, content, and pilot test validity were examined. Moreover, triangulation, use of multiple methods, external audits, and low inference descriptors were employed for the qualitative instrument. Descriptive statistics (frequency, percent, mean, and standard deviation) and inferential statics (Chi-square, Cramer's-V test, and one-way ANOVA) were used to examine the quantitative data. Bell's "Spectrum of Acceptance of Internationalizing Curriculum," integrated with Ellingboe's "Great Divide," was used to identify the positions of the academics in their beliefs. The qualitative data was analyzed using descriptive statements and presented thematically based on research questions. Hence, major findings revealed that academics' representations based on their demographic variables were not as expected; their variation was practically not significant among universities; the majority of academics had crossed Ellingboe's "Great Divide"; significant differences among disciplines and/or universities were not observed in this study. Finally, urging academics to actively integrate international perspectives into their curricula, elevating curriculum internationalization of higher education to a national priority, and undertaking extensive national research emerged as key recommendations for the Ethiopian Research Universities, Ethiopian Ministry of Education, and future research endeavors, respectively.

Experiences and lessons from structural interventions against COVID-19 in Addis Ababa, Ethiopia.

Introduction: Fighting pandemics like COVID-19 requires implementing successful structural and behavioral interventions that attempt to change the social and political environments to increase adherence to preventive behavior among community members. However, studying structural interventions implemented during pandemics and their challenges remains to be uncharted territory in developing implemented countries. Objectives: Given this, we documented the experiences of implementing such interventions in Ethiopia with the aim of drawing lessons for future efforts to fight similar outbreaks in resource limited and low-income settings. Methods: We conducted a qualitative study between September and October 2021. Data were collected through face to face and telephone interviews from purposefully selected stakeholders from government and private sectors engaged in social interventions to prevent COVID-19. The systematization and the analysis of the data were conducted with MAXQDA 2020 software. Results: Ethiopia implemented structural and social interventions to respond to the COVID-19 pandemic. This included: developing national policy and guidelines, mainstreaming COVID-19 interventions to local organizations, implementing capacity development programs, and developing strategies to engage the community, through traditional institutions, in intervention activities. In addition, a mass communication approach was used to deliver risk messages. This yielded a promising result in slowing down the spread of COVID-19 in the capital of Ethiopia-Addis Ababa. On the other hand, competing interests, misconceptions, capacity constraints among professionals and organizations, limited capacity to enforce legislation and lack of motivation for change from the community side affected the implementation and the outcomes of interventions. Conclusion: Going forward, these challenges need to be taken into consideration when designing and implementing structural interventions to contain disease outbreaks effectively. The study highlighted that attempts to withstand pandemic in low- and middle-income settings shall successfully utilize local resources, act swiftly when pandemics outbreak and adjust themselves to the dynamic challenges and limitations of structural interventions.

Neisseria meningitidis carriage rate, antibiotic susceptibility profile, and associated factors among prisoners at Jimma zonal correction facility in Jimma Town, Southwestern Ethiopia: a cross-sectional study.

BACKGROUND: Neisseria meningitidis causes severe life-threatening meningococcal disease with a case fatality rate of 10-15% even with proper treatment. In Ethiopia, particularly in our study area, inadequate information is found on meningococcal disease. So, this study aimed to assess N. meningitidis carriage rate, antibiotic susceptibility profile, and associated factors among prisoners in Jimma Town, Southwestern Ethiopia. METHODS: A cross-sectional study was conducted in Jimma town, Southwest Ethiopia, from May to October 2019. A stratified sampling technique was used and proportional allocation was done. A total of 550 oropharyngeal swabs were collected, processed, isolated, and identified N. meningitidis using standard microbiological techniques. Antibiotics susceptibility test was done for isolates using the disk diffusion method. Data on demographic and associated factors for carriage were collected using a structured questionnaire. Data were summarized using frequency, percentage, graph, and table. A logistic regression model was used to see the association between the dependent and independent variables. Variables with a p-value < 0.25 during bivariate analysis were included in multivariate analysis to identify factors significantly associated with the meningococcal carriage and, a p-value < 0.05 was considered statistically significant. RESULT: Out of the 550 study participants, 76(13.8%) with (CI: 7.20-18.20) were found carriers of N meningitidis. The predominant isolates were non-serogroupable 26(34.2%) and serogroup W/Y 22(28.9%), respectively. N. meningitidis isolates showed highest sensitivity to chloramphenicol 74(97.4%). Meningococcal carriage rate was significantly associated with being age group of 16-20 years; having respiratory symptoms within 3 months and active cigarette smoking within 3 months. CONCLUSIONS: The majority of participants harbor most of the serogroups responsible for invasive cases of meningococcal disease. Respiratory symptoms, active cigarette smoking, and age group of 16-20 years increased the risk of N. meningitidis pharyngeal carriage rate. This study suggests providing better health education to control respiratory symptoms, smoking, and providing antibiotic prophylaxis for prisoners.

Identification of α-Glucosidase Inhibitors from Leaf Extract of Pepper (Capsicum spp.) through Metabolomic Analysis.

Metabolomics and in vitro α-glucosidase inhibitory (AGI) activities of pepper leaves were used to identify bioactive compounds and select genotypes for the management of type 2 diabetes mellitus (T2DM). Targeted metabolite analysis using UPLC-DAD-QToF-MS was employed and identified compounds that belong to flavone and hydroxycinnamic acid derivatives from extracts of pepper leaves. A total of 21 metabolites were detected from 155 samples and identified based on MS fragmentations, retention time, UV absorbance, and previous reports. Apigenin-O-(malonyl) hexoside, luteolin-O-(malonyl) hexoside, and chrysoeriol-O-(malonyl) hexoside were identified for the first time from pepper leaves. Pepper genotypes showed a huge variation in their inhibitory activity against α-glucosidase enzyme(AGE) ranging from 17% to 79%. Genotype GP38 with inhibitory activity of 79% was found to be more potent than the positive control acarbose (70.8%.). Orthogonal partial least square (OPLS) analyses were conducted for the prediction of the AGI activities of pepper leaves based on their metabolite composition. Compounds that contributed the most to the bioactivity prediction model (VIP >1.5), showed a strong inhibitory potency. Caffeoyl-putrescine was found to show a stronger inhibitory potency (IC50 = 145 µM) compared to acarbose (IC50 = 197 µM). The chemometric procedure combined with high-throughput AGI screening was effective in selecting polyphenols of pepper leaf for T2DM management.

Alpha Glucosidase Inhibitory Activities of Plants with Focus on Common Vegetables.

Type-2 diabetes mellitus is one of the most prevalent metabolic diseases in the world, and is characterized by hyperglycemia (i.e., high levels of glucose in the blood). Alpha-glucosidases are enzymes in the digestive tract that hydrolyze carbohydrates into glucose. One strategy that has been developed to treat type-2 diabetes is inhibition of the activity of alpha-glucosidases using synthetic drugs. However, these inhibitors are usually associated with gastrointestinal side effects. Therefore, the development of inhibitors from natural products offers an alternative option for the control of hyperglycemia. In recent years, various studies have been conducted to identify alpha-glucosidases inhibitors from natural sources such as plants, and many candidates have transpired to be secondary metabolites including alkaloids, flavonoids, phenols, and terpenoids. In this review, we focus on the alpha-glucosidases inhibitors found in common vegetable crops and the major classes of phytochemicals responsible for the inhibitory activity, and also as potential/natural drug candidates for the treatment of type-2 diabetes mellitus. In addition, possible breeding strategies for production of improved vegetable crops with higher content of the inhibitors are also described.

Population genetic structure and adaptation of malaria parasites on the edge of endemic distribution.

To determine whether the major human malaria parasite Plasmodium falciparum exhibits fragmented population structure or local adaptation at the northern limit of its African distribution where the dry Sahel zone meets the Sahara, samples were collected from diverse locations within Mauritania over a range of ~1000 km. Microsatellite genotypes were obtained for 203 clinical infection samples from eight locations, and Illumina paired-end sequences were obtained to yield high coverage genomewide single nucleotide polymorphism (SNP) data for 65 clinical infection samples from four locations. Most infections contained single parasite genotypes, reflecting low rates of transmission and superinfection locally, in contrast to the situation seen in population samples from countries further south. A minority of infections shared related or identical genotypes locally, indicating some repeated transmission of parasite clones without recombination. This caused some multilocus linkage disequilibrium and local divergence, but aside from the effect of repeated genotypes there was minimal differentiation between locations. Several chromosomal regions had elevated integrated haplotype scores (|iHS|) indicating recent selection, including those containing drug resistance genes. A genomewide FST scan comparison with previous sequence data from an area in West Africa with higher infection endemicity indicates that regional gene flow prevents genetic isolation, but revealed allele frequency differentiation at three drug resistance loci and an erythrocyte invasion ligand gene. Contrast of extended haplotype signatures revealed none to be unique to Mauritania. Discrete foci of infection on the edge of the Sahara are genetically highly connected to the wider continental parasite population, and local elimination would be difficult to achieve without very substantial reduction in malaria throughout the region.

Culture adaptation of malaria parasites selects for convergent loss-of-function mutants.

Cultured human pathogens may differ significantly from source populations. To investigate the genetic basis of laboratory adaptation in malaria parasites, clinical Plasmodium falciparum isolates were sampled from patients and cultured in vitro for up to three months. Genome sequence analysis was performed on multiple culture time point samples from six monoclonal isolates, and single nucleotide polymorphism (SNP) variants emerging over time were detected. Out of a total of five positively selected SNPs, four represented nonsense mutations resulting in stop codons, three of these in a single ApiAP2 transcription factor gene, and one in SRPK1. To survey further for nonsense mutants associated with culture, genome sequences of eleven long-term laboratory-adapted parasite strains were examined, revealing four independently acquired nonsense mutations in two other ApiAP2 genes, and five in Epac. No mutants of these genes exist in a large database of parasite sequences from uncultured clinical samples. This implicates putative master regulator genes in which multiple independent stop codon mutations have convergently led to culture adaptation, affecting most laboratory lines of P. falciparum. Understanding the adaptive processes should guide development of experimental models, which could include targeted gene disruption to adapt fastidious malaria parasite species to culture.

Microsatellite genotyping and genome-wide single nucleotide polymorphism-based indices of Plasmodium falciparum diversity within clinical infections.

BACKGROUND: In regions where malaria is endemic, individuals are often infected with multiple distinct parasite genotypes, a situation that may impact on evolution of parasite virulence and drug resistance. Most approaches to studying genotypic diversity have involved analysis of a modest number of polymorphic loci, although whole genome sequencing enables a broader characterisation of samples. METHODS: PCR-based microsatellite typing of a panel of ten loci was performed on Plasmodium falciparum in 95 clinical isolates from a highly endemic area in the Republic of Guinea, to characterize within-isolate genetic diversity. Separately, single nucleotide polymorphism (SNP) data from genome-wide short-read sequences of the same samples were used to derive within-isolate fixation indices (F ws), an inverse measure of diversity within each isolate compared to overall local genetic diversity. The latter indices were compared with the microsatellite results, and also with indices derived by randomly sampling modest numbers of SNPs. RESULTS: As expected, the number of microsatellite loci with more than one allele in each isolate was highly significantly inversely correlated with the genome-wide F ws fixation index (r = -0.88, P < 0.001). However, the microsatellite analysis revealed that most isolates contained mixed genotypes, even those that had no detectable genome sequence heterogeneity. Random sampling of different numbers of SNPs showed that an F ws index derived from ten or more SNPs with minor allele frequencies of >10 % had high correlation (r > 0.90) with the index derived using all SNPs. CONCLUSIONS: Different types of data give highly correlated indices of within-infection diversity, although PCR-based analysis detects low-level minority genotypes not apparent in bulk sequence analysis. When whole-genome data are not obtainable, quantitative assay of ten or more SNPs can yield a reasonably accurate estimate of the within-infection fixation index (F ws).

Genomic variation in two gametocyte non-producing Plasmodium falciparum clonal lines.

BACKGROUND: Transmission of the malaria parasite Plasmodium falciparum from humans to the mosquito vector requires differentiation of a sub-population of asexual forms replicating within red blood cells into non-dividing male and female gametocytes. The nature of the molecular mechanism underlying this key differentiation event required for malaria transmission is not fully understood. METHODS: Whole genome sequencing was used to examine the genomic diversity of the gametocyte non-producing 3D7-derived lines F12 and A4. These lines were used in the recent detection of the PF3D7_1222600 locus (encoding PfAP2-G), which acts as a genetic master switch that triggers gametocyte development. RESULTS: The evolutionary changes from the 3D7 parental strain through its derivatives F12 (culture-passage derived cloned line) and A4 (transgenic cloned line) were identified. The genetic differences including the formation of chimeric var genes are presented. CONCLUSION: A genomics resource is provided for the further study of gametocytogenesis or other phenotypes using these parasite lines.

Population genomic structure and adaptation in the zoonotic malaria parasite Plasmodium knowlesi.

Malaria cases caused by the zoonotic parasite Plasmodium knowlesi are being increasingly reported throughout Southeast Asia and in travelers returning from the region. To test for evidence of signatures of selection or unusual population structure in this parasite, we surveyed genome sequence diversity in 48 clinical isolates recently sampled from Malaysian Borneo and in five lines maintained in laboratory rhesus macaques after isolation in the 1960s from Peninsular Malaysia and the Philippines. Overall genomewide nucleotide diversity (π = 6.03 × 10(-3)) was much higher than has been seen in worldwide samples of either of the major endemic malaria parasite species Plasmodium falciparum and Plasmodium vivax. A remarkable substructure is revealed within P. knowlesi, consisting of two major sympatric clusters of the clinical isolates and a third cluster comprising the laboratory isolates. There was deep differentiation between the two clusters of clinical isolates [mean genomewide fixation index (FST) = 0.21, with 9,293 SNPs having fixed differences of FST = 1.0]. This differentiation showed marked heterogeneity across the genome, with mean FST values of different chromosomes ranging from 0.08 to 0.34 and with further significant variation across regions within several chromosomes. Analysis of the largest cluster (cluster 1, 38 isolates) indicated long-term population growth, with negatively skewed allele frequency distributions (genomewide average Tajima's D = -1.35). Against this background there was evidence of balancing selection on particular genes, including the circumsporozoite protein (csp) gene, which had the top Tajima's D value (1.57), and scans of haplotype homozygosity implicate several genomic regions as being under recent positive selection.

Comparison of genomic signatures of selection on Plasmodium falciparum between different regions of a country with high malaria endemicity.

BACKGROUND: Genome wide sequence analyses of malaria parasites from widely separated areas of the world have identified contrasting population structures and signatures of selection. To compare relatively closely situated but ecologically contrasting regions within an endemic African country, population samples of Plasmodium falciparum clinical isolates were collected in Ghana from Kintampo in the central forest-savannah area, and Navrongo in a drier savannah area ~350 km to the north with more seasonally-restricted transmission. Parasite DNA was sequenced and paired-end reads mapped to the P. falciparum reference genome. RESULTS: High coverage genome wide sequence data for 85 different clinical isolates enabled analysis of 121,712 single nucleotide polymorphisms (SNPs). The local populations had similar proportions of mixed genotype infections, similar SNP allele frequency distributions, and eleven chromosomal regions had elevated integrated haplotype scores (|iHS|) in both. A between-population Rsb metric comparing extended haplotype homozygosity indicated a stronger signal within Kintampo for one of these regions (on chromosome 14) and in Navrongo for two of these regions (on chromosomes 10 and 13). At least one gene in each of these identified regions is a potential target of locally varying selection. The candidates include genes involved in parasite development in mosquitoes, members of variant-expressed multigene families, and a leading vaccine-candidate target of immunity. CONCLUSIONS: Against a background of very similar population structure and selection signatures in the P. falciparum populations of Ghana, three narrow genomic regions showed evidence indicating local differences in historical timing or intensity of selection. Sampling of closely situated populations across heterogeneous environments has potential to refine the mapping of important loci under temporally or spatially varying selection.

Admixture in Humans of Two Divergent Plasmodium knowlesi Populations Associated with Different Macaque Host Species.

Human malaria parasite species were originally acquired from other primate hosts and subsequently became endemic, then spread throughout large parts of the world. A major zoonosis is now occurring with Plasmodium knowlesi from macaques in Southeast Asia, with a recent acceleration in numbers of reported cases particularly in Malaysia. To investigate the parasite population genetics, we developed sensitive and species-specific microsatellite genotyping protocols and applied these to analysis of samples from 10 sites covering a range of >1,600 km within which most cases have occurred. Genotypic analyses of 599 P. knowlesi infections (552 in humans and 47 in wild macaques) at 10 highly polymorphic loci provide radical new insights on the emergence. Parasites from sympatric long-tailed macaques (Macaca fascicularis) and pig-tailed macaques (M. nemestrina) were very highly differentiated (FST = 0.22, and K-means clustering confirmed two host-associated subpopulations). Approximately two thirds of human P. knowlesi infections were of the long-tailed macaque type (Cluster 1), and one third were of the pig-tailed-macaque type (Cluster 2), with relative proportions varying across the different sites. Among the samples from humans, there was significant indication of genetic isolation by geographical distance overall and within Cluster 1 alone. Across the different sites, the level of multi-locus linkage disequilibrium correlated with the degree of local admixture of the two different clusters. The widespread occurrence of both types of P. knowlesi in humans enhances the potential for parasite adaptation in this zoonotic system.

Whole-genome scans provide evidence of adaptive evolution in Malawian Plasmodium falciparum isolates.

BACKGROUND: Selection by host immunity and antimalarial drugs has driven extensive adaptive evolution in Plasmodium falciparum and continues to produce ever-changing landscapes of genetic variation. METHODS: We performed whole-genome sequencing of 69 P. falciparum isolates from Malawi and used population genetics approaches to investigate genetic diversity and population structure and identify loci under selection. RESULTS: High genetic diversity (π = 2.4 × 10(-4)), moderately high multiplicity of infection (2.7), and low linkage disequilibrium (500-bp) were observed in Chikhwawa District, Malawi, an area of high malaria transmission. Allele frequency-based tests provided evidence of recent population growth in Malawi and detected potential targets of host immunity and candidate vaccine antigens. Comparison of the sequence variation between isolates from Malawi and those from 5 geographically dispersed countries (Kenya, Burkina Faso, Mali, Cambodia, and Thailand) detected population genetic differences between Africa and Asia, within Southeast Asia, and within Africa. Haplotype-based tests of selection to sequence data from all 6 populations identified signals of directional selection at known drug-resistance loci, including pfcrt, pfdhps, pfmdr1, and pfgch1. CONCLUSIONS: The sequence variations observed at drug-resistance loci reflect differences in each country's historical use of antimalarial drugs and may be useful in formulating local malaria treatment guidelines.

A barcode of organellar genome polymorphisms identifies the geographic origin of Plasmodium falciparum strains.

Malaria is a major public health problem that is actively being addressed in a global eradication campaign. Increased population mobility through international air travel has elevated the risk of re-introducing parasites to elimination areas and dispersing drug-resistant parasites to new regions. A simple genetic marker that quickly and accurately identifies the geographic origin of infections would be a valuable public health tool for locating the source of imported outbreaks. Here we analyse the mitochondrion and apicoplast genomes of 711 Plasmodium falciparum isolates from 14 countries, and find evidence that they are non-recombining and co-inherited. The high degree of linkage produces a panel of relatively few single-nucleotide polymorphisms (SNPs) that is geographically informative. We design a 23-SNP barcode that is highly predictive (~92%) and easily adapted to aid case management in the field and survey parasite migration worldwide.

estMOI: estimating multiplicity of infection using parasite deep sequencing data.

Individuals living in endemic areas generally harbour multiple parasite strains. Multiplicity of infection (MOI) can be an indicator of immune status and transmission intensity. It has a potentially confounding effect on a number of population genetic analyses, which often assume isolates are clonal. Polymerase chain reaction-based approaches to estimate MOI can lack sensitivity. For example, in the human malaria parasite Plasmodium falciparum, genotyping of the merozoite surface protein (MSP1/2) genes is a standard method for assessing MOI, despite the apparent problem of underestimation. The availability of deep coverage data from massively parallizable sequencing technologies means that MOI can be detected genome wide by considering the abundance of heterozygous genotypes. Here, we present a method to estimate MOI, which considers unique combinations of polymorphisms from sequence reads. The method is implemented within the estMOI software. When applied to clinical P.falciparum isolates from three continents, we find that multiple infections are common, especially in regions with high transmission.

PlasmoView: a web-based resource to visualise global Plasmodium falciparum genomic variation.

Malaria is a global public health challenge, with drug resistance a major barrier to disease control and elimination. To meet the urgent need for better treatments and vaccines, a deeper knowledge of Plasmodium biology and malaria epidemiology is required. An improved understanding of the genomic variation of malaria parasites, especially the most virulent Plasmodium falciparum (Pf) species, has the potential to yield new insights in these areas. High-throughput sequencing and genotyping is generating large amounts of genomic data across multiple parasite populations. The resulting ability to identify informative variants, particularly single-nucleotide polymorphisms (SNPs), will lead to the discovery of intra- and inter-population differences and thus enable the development of genetic barcodes for diagnostic assays and clinical studies. Knowledge of genetic variability underlying drug resistance and other differential phenotypes will also facilitate the identification of novel mutations and contribute to surveillance and stratified medicine applications. The PlasmoView interactive web-browsing tool enables the research community to visualise genomic variation and annotation (eg, biological function) in a geographic setting. The first release contains over 600,000 high-quality SNPs in 631 Pf isolates from laboratory strains and four malaria-endemic regions (West Africa, East Africa, Southeast Asia and Oceania).

Elucidating emergence and transmission of multidrug-resistant tuberculosis in treatment experienced patients by whole genome sequencing.

BACKGROUND: Understanding the emergence and spread of multidrug-resistant tuberculosis (MDR-TB) is crucial for its control. MDR-TB in previously treated patients is generally attributed to the selection of drug resistant mutants during inadequate therapy rather than transmission of a resistant strain. Traditional genotyping methods are not sufficient to distinguish strains in populations with a high burden of tuberculosis and it has previously been difficult to assess the degree of transmission in these settings. We have used whole genome analysis to investigate M. tuberculosis strains isolated from treatment experienced patients with MDR-TB in Uganda over a period of four years. METHODS AND FINDINGS: We used high throughput genome sequencing technology to investigate small polymorphisms and large deletions in 51 Mycobacterium tuberculosis samples from 41 treatment-experienced TB patients attending a TB referral and treatment clinic in Kampala. This was a convenience sample representing 69% of MDR-TB cases identified over the four year period. Low polymorphism was observed in longitudinal samples from individual patients (2-15 SNPs). Clusters of samples with less than 50 SNPs variation were examined. Three clusters comprising a total of 8 patients were found with almost identical genetic profiles, including mutations predictive for resistance to rifampicin and isoniazid, suggesting transmission of MDR-TB. Two patients with previous drug susceptible disease were found to have acquired MDR strains, one of which shared its genotype with an isolate from another patient in the cohort. CONCLUSIONS: Whole genome sequence analysis identified MDR-TB strains that were shared by more than one patient. The transmission of multidrug-resistant disease in this cohort of retreatment patients emphasises the importance of early detection and need for infection control. Consideration should be given to rapid testing for drug resistance in patients undergoing treatment to monitor the emergence of resistance and permit early intervention to avoid onward transmission.

Genome-wide screen identifies new candidate genes associated with artemisinin susceptibility in Plasmodium falciparum in Kenya.

Early identification of causal genetic variants underlying antimalarial drug resistance could provide robust epidemiological tools for timely public health interventions. Using a novel natural genetics strategy for mapping novel candidate genes we analyzed >75,000 high quality single nucleotide polymorphisms selected from high-resolution whole-genome sequencing data in 27 isolates of Plasmodium falciparum. We identified genetic variants associated with susceptibility to dihydroartemisinin that implicate one region on chromosome 13, a candidate gene on chromosome 1 (PFA0220w, a UBP1 ortholog) and others (PFB0560w, PFB0630c, PFF0445w) with putative roles in protein homeostasis and stress response. There was a strong signal for positive selection on PFA0220w, but not the other candidate loci. Our results demonstrate the power of full-genome sequencing-based association studies for uncovering candidate genes that determine parasite sensitivity to artemisinins. Our study provides a unique reference for the interpretation of results from resistant infections.

A Poisson hierarchical modelling approach to detecting copy number variation in sequence coverage data.

BACKGROUND: The advent of next generation sequencing technology has accelerated efforts to map and catalogue copy number variation (CNV) in genomes of important micro-organisms for public health. A typical analysis of the sequence data involves mapping reads onto a reference genome, calculating the respective coverage, and detecting regions with too-low or too-high coverage (deletions and amplifications, respectively). Current CNV detection methods rely on statistical assumptions (e.g., a Poisson model) that may not hold in general, or require fine-tuning the underlying algorithms to detect known hits. We propose a new CNV detection methodology based on two Poisson hierarchical models, the Poisson-Gamma and Poisson-Lognormal, with the advantage of being sufficiently flexible to describe different data patterns, whilst robust against deviations from the often assumed Poisson model. RESULTS: Using sequence coverage data of 7 Plasmodium falciparum malaria genomes (3D7 reference strain, HB3, DD2, 7G8, GB4, OX005, and OX006), we showed that empirical coverage distributions are intrinsically asymmetric and overdispersed in relation to the Poisson model. We also demonstrated a low baseline false positive rate for the proposed methodology using 3D7 resequencing data and simulation. When applied to the non-reference isolate data, our approach detected known CNV hits, including an amplification of the PfMDR1 locus in DD2 and a large deletion in the CLAG3.2 gene in GB4, and putative novel CNV regions. When compared to the recently available FREEC and cn.MOPS approaches, our findings were more concordant with putative hits from the highest quality array data for the 7G8 and GB4 isolates. CONCLUSIONS: In summary, the proposed methodology brings an increase in flexibility, robustness, accuracy and statistical rigour to CNV detection using sequence coverage data.

VarB: a variation browsing and analysis tool for variants derived from next-generation sequencing data.

SUMMARY: There is an immediate need for tools to both analyse and visualize in real-time single-nucleotide polymorphisms, insertions and deletions, and other structural variants from new sequence file formats. We have developed VarB software that can be used to visualize variant call format files in real time, as well as identify regions under balancing selection and informative markers to differentiate user-defined groups (e.g. populations). We demonstrate its utility using sequence data from 50 Plasmodium falciparum isolates comprising two different continents and confirm known signals from genomic regions that contain important antigenic and anti-malarial drug-resistance genes. AVAILABILITY AND IMPLEMENTATION: The C++-based software VarB and user manual are available from www.pathogenseq.org/varb. CONTACT: taane.clark@lshtm.ac.uk