Genome-Wide Association Study of Growth Performance and Immune Response to Newcastle Disease Virus of Indigenous Chicken in Rwanda.

A chicken genome has several regions with quantitative trait loci (QTLs). However, replication and confirmation of QTL effects are required particularly in African chicken populations. This study identified single nucleotide polymorphisms (SNPs) and putative genes responsible for body weight (BW) and antibody response (AbR) to Newcastle disease (ND) in Rwanda indigenous chicken (IC) using genome-wide association studies (GWAS). Multiple testing was corrected using chromosomal false detection rates of 5 and 10% for significant and suggestive thresholds, respectively. BioMart data mining and variant effect predictor tools were used to annotate SNPs and candidate genes, respectively. A total of four significant SNPs (rs74098018, rs13792572, rs314702374, and rs14123335) significantly (p ≤ 7.6E-5) associated with BW were identified on chromosomes (CHRs) 8, 11, and 19. In the vicinity of these SNPs, four genes such as pre-B-cell leukaemia homeobox 1 (PBX1), GPATCH1, MPHOSPH6, and MRM1 were identified. Four other significant SNPs (rs314787954, rs13623466, rs13910430, and rs737507850) all located on chromosome 1 were strongly (p ≤ 7.6E-5) associated with chicken antibody response to ND. The closest genes to these four SNPs were cell division cycle 16 (CDC16), zinc finger, BED-type containing 1 (ZBED1), myxovirus (influenza virus) resistance 1 (MX1), and growth factor receptor bound protein 2 (GRB2) related adaptor protein 2 (GRAP2). Besides, other SNPs and genes suggestively (p ≤ 1.5E-5) associated with BW and antibody response to ND were reported. This work offers a useful entry point for the discovery of causative genes accountable for essential QTLs regulating BW and antibody response to ND traits. Results provide auspicious genes and SNP-based markers that can be used in the improvement of growth performance and ND resistance in IC populations based on gene-based and/or marker-assisted breeding selection.

Analysis of ß-amylase gene (Amyß) variation reveals allele association with low enzyme activity and increased firmness in cooked sweetpotato (Ipomoea batatas) from East Africa.

ß-amylase is a thermostable enzyme that hydrolyses starch during cooking of sweetpotato (Ipomoea batatas) storage roots, thereby influencing eating quality. Its activity is known to vary amongst genotypes but the genetic diversity of the beta-amylase gene (Amyß) is not well studied. Amyß has a highly conserved region between exon V and VI, forming part of the enzyme's active site. To determine the gene diversity, a 2.3 kb fragment, including the conserved region of the Amyß gene was sequenced from 25 sweetpotato genotypes. The effect of sequence variation on gene expression, enzyme activity, and firmness in cooked roots was determined. Six genotypes carrying several SNPs within exon V, linked with an AT or ATGATA insertion in intron V were unique and clustered together. The genotypes also shared an A336E substitution in the amino acid sequence, eight residues upstream of a substrate-binding Thr344. The genotypes carrying this allele exhibited low gene expression and low enzyme activity. Enzyme activity was negatively correlated with firmness (R = -0.42) in cooked roots. This is the first report of such an allele, associated with low enzyme activity. These results suggest that genetic variation within the AmyB locus can be utilized to develop markers for firmness in sweetpotato breeding.

Introgression of Maize Lethal Necrosis Resistance Quantitative Trait Loci Into Susceptible Maize Populations and Validation of the Resistance Under Field Conditions in Naivasha, Kenya.

Maize lethal necrosis (MLN), resulting from co-infection by maize chlorotic mottle virus (MCMV) and sugarcane mosaic virus (SCMV) can cause up to 100% yield losses in maize in Africa under serious disease conditions. Maize improvement through conventional backcross (BC) takes many generations but can significantly be shortened when molecular tools are utilized in the breeding process. We used a donor parent (KS23-6) to transfer quantitative trait loci (QTL) for resistance to MLN into nine adapted but MLN susceptible lines. Nurseries were established in Kiboko, Kenya during 2015-2017 seasons and BC3F2 progeny were developed using marker assisted backcrossing (MABC) approach. Six single nucleotide polymorphism (SNP) markers linked to QTL for resistance to MLN were used to genotype 2,400 BC3F2 lines using Kompetitive Allele Specific PCR (KASP) platform. We detected that two of the six QTL had major effects for resistance to MLN under artificial inoculation field conditions in 56 candidate BC3F2 lines. To confirm whether these two QTL are reproducible under different field conditions, the 56 BC3F2 lines including their parents were evaluated in replicated trials for two seasons under artificial MLN inoculations in Naivasha, Kenya in 2018. Strong association of genotype with phenotype was detected. Consequently, 19 superior BC3F2 lines with favorable alleles and showing improved levels of resistance to MLN under artificial field inoculation were identified. These elite lines represent superior genetic resources for improvement of maize hybrids for resistance to MLN. However, 20 BC3F2 lines were fixed for both KASP markers but were susceptible to MLN under field conditions, which could suggest weak linkage between the KASP markers and target genes. The validated two major QTL can be utilized to speed up the breeding process but additional loci need to be identified between the KASP markers and the resistance genes to strengthen the linkage.

Marker-assisted selection complements phenotypic screening at seedling stage to identify cassava mosaic disease-resistant genotypes in African cassava populations.

Cassava mosaic disease (CMD) is a serious threat to cassava production in sub-Saharan Africa. The use of genomic-assisted selection at the seedling trial stage would help to reduce the time for release, breeding cost, and resources used, hence increase selection efficiency in cassava breeding programs. Five cassava populations were screened for resistance to CMD during the seedling evaluation trial at 1, 3, and 5 months after planting using a scale of 1-5. The genotypes in the five populations were also screened using six molecular markers linked to the CMD2 gene. The correlation between the phenotypic and marker data was estimated. Based on Cassava Mosaic Disease Severity Score (CMDSS), between 53 and 82% of the progenies were resistant across the populations with an average of 70.5%. About 70% of the progenies were identified to be resistant to the disease across the populations with a range of 62-80% using the marker data. With both marker data and CMDSS combined, 40-60% of the progenies in each population, with an average of 52%, were identified to be resistant to CMD. There was a fairly significant correlation between the marker data and CMDSS in each cassava population with correlation coefficients ranging from 0.2024 to 0.3460 suggesting that novel genes not associated to the markers used might be involved in the resistance to CMD. The resistant genotypes identified in this study with potential for other desirable traits were selected for evaluation at the advanced trial stage thereby shortening the period required for the breeding program.

DArTSeq SNP-based markers revealed high genetic diversity and structured population in Ethiopian cowpea [Vigna unguiculata (L.) Walp] germplasms.

Cowpea [Vigna unguiculata (L.) Walp] is one of the important climate-resilient legume crops for food and nutrition security in sub-Saharan Africa. Ethiopia is believed to harbor high cowpea genetic diversity, but this has not yet been efficiently characterized and exploited in breeding. The objective of this study was to evaluate the extent and pattern of genetic diversity in 357 cowpea accestions comprising landraces (87%), breeding lines (11%) and released varieties (2%), using single nucleotide polymorphism markers. The overall gene diversity and heterozygosity were 0.28 and 0.12, respectively. The genetic diversity indices indicated substantial diversity in Ethiopian cowpea landraces. Analysis of molecular variance showed that most of the variation was within in the population (46%) and 44% between individuals, with only 10% of the variation being among populations. Model-based ancestry analysis, the phylogenetic tree, discriminant analysis of principal components and principal coordinate analysis classified the 357 genotypes into three well-differentiated genetic populations. Genotypes from the same region grouped into different clusters, while others from different regions fell into the same cluster. This indicates that differences in regions of origin may not be the main driver determining the genetic diversity in cowpea in Ethiopia. Therefore, differences in sources of origin, as currently distributed in Ethiopia, should not necessarily be used as indices of genetic diversity. Choice of parental lines should rather be based on a systematic assessment of genetic diversity in a specific population. The study also suggested 94 accesstions as core collection which retained 100% of the genetic diversity from the entire collection. This core set represents 26% of the entire collection pinpointing a wide distribution of the diversity within the ethiopian landraces. The outcome of this study provided new insights into the genetic diversity and population structure in Ethiopian cowpea genetic resources for designing effective collection and conservation strategies for efficient utilization in breeding.

Correction: Genetic diversity and population structure of indigenous chicken in Rwanda using microsatellite markers.

[This corrects the article DOI: 10.1371/journal.pone.0225084.].

Genome-wide genetic diversity and population structure of Garcinia kola (Heckel) in Benin using DArT-Seq technology.

Garcinia kola (Heckel) is a versatile tree indigenous to West and Central Africa. All parts of the tree have value in traditional medicine. Natural populations of the species have declined over the years due to overexploitation. Assessment of genetic diversity and population structure of G. kola is important for its management and conservation. The present study investigates the genetic diversity and population structure of G. kola populations in Benin using ultra-high-throughput diversity array technology (DArT) single nucleotide polymorphism (SNP) markers. From the 102 accessions sampled, two were excluded from the final dataset owing to poor genotyping coverage. A total of 43,736 SNPs were reported, of which 12,585 were used for analyses after screening with quality control parameters including Minor allele frequency (≥ 0.05), call rate (≥ 80%), reproducibility (≥ 95%), and polymorphic information content (≥ 1%). Analysis revealed low genetic diversity with expected heterozygosity per population ranging from 0.196 to 0.228. Pairwise F-statistics (FST) revealed low levels of genetic differentiation between populations while an Analysis of molecular variance (AMOVA) indicated that the majority of variation (97.86%) was within populations. Population structure analysis through clustering and discriminant analysis on principal component revealed two admixed clusters, implying little genetic structure. However, the model-based maximum likelihood in Admixture indicated only one genetic cluster. The present study indicated low genetic diversity of G. kola, and interventions are needed to be tailored towards its conservation.

Genomic analysis of Sweet potato feathery mottle virus from East Africa.

Sweet potato feathery mottle virus is a potyvirus that infect sweet potato. The genome of the virus was analysed to understand genetic diversity, evolution and gene flow. Motifs, nucleotide identity and a phylogenetic tree were used to determine phylogroup of the isolates. Gene flow and genetic diversity were tested using DnaSP v.5. Codons evolution were tested using three methods embedded in Datamonkey. The results indicate occurrence of an isolate of phylogroup B within East Africa. Low genetic differentiation was observed between isolates from Kenya and Uganda indicating evidence of gene flow between the two countries. Four genes were found to have positively selected codons bordering or occurring within functional motifs. A motif within P1 gene evolved differently between phylogroup A and B. The evidence of gene flow indicates frequent exchange of the virus between the two countries and P1 gene motif provide a possible marker that can be used for mapping the distribution of the phylogroups.

Genetic diversity and population structure of indigenous chicken in Rwanda using microsatellite markers.

Rwanda has about 4.5 million of indigenous chicken (IC) that are very low in productivity. To initiate any genetic improvement programme, IC needs to be accurately characterized. The key purpose of this study was to ascertain the genetic diversity of IC in Rwanda using microsatellite markers. Blood samples of IC sampled from 5 agro-ecological zones were collected from which DNA was extracted, amplified by PCR and genotyped using 28 microsatellite markers. A total of 325 (313 indigenous and 12 exotic) chickens were genotyped and revealed a total number of 305 alleles varying between 2 and 22 with a mean of 10.89 per locus. One hundred eighty-six (186) distinct alleles and 60 private alleles were also observed. The frequency of private alleles was highest in samples from the Eastern region, whereas those from the North West had the lowest. The influx of genes was lower in the Eastern agro-ecological zone than the North West. The mean observed heterozygosity was 0.6155, whereas the average expected heterozygosity was 0.688. The overall inbreeding coefficient among the population was 0.040. Divergence from the Hardy-Weinberg equilibrium was significant (p<0.05) in 90% of loci in all the populations. The analysis of molecular variance revealed that about 92% of the total variation originated from variation within populations. Additionally, the study demonstrated that IC in Rwanda could be clustered into four gene groups. In conclusion, there was considerable genetic diversity in IC in Rwanda, which represents a crucial genetic resource that can be conserved or optimized through genetic improvement.

Correction: Diversity analysis and genome-wide association studies of grain shape and eating quality traits in rice (Oryza sativa L.) using DArT markers.

[This corrects the article DOI: 10.1371/journal.pone.0198012.].

Diversity analysis and genome-wide association studies of grain shape and eating quality traits in rice (Oryza sativa L.) using DArT markers.

Microarray-based markers such as Diversity Arrays Technology (DArT) have become the genetic markers of choice for construction of high-density maps, quantitative trait loci (QTL) mapping and genetic diversity analysis based on their efficiency and low cost. More recently, the DArT technology was further developed in combination with high-throughput next-generation sequencing (NGS) technologies to generate the DArTseq platform representing a new sequencing tool of complexity-reduced representations. In this study, we used DArTseq markers to investigate genetic diversity and genome-wide association studies (GWAS) of grain quality traits in rice (Oryza sativa L.). The study was performed using 59 rice genotypes with 525 SNPs derived from DArTseq platform. Population structure analysis revealed only two distinct genetic clusters where genotypes were grouped based on environmental adaptation and pedigree information. Analysis of molecular variance indicated a low degree of differentiation among populations suggesting the need for broadening the genetic base of the current germplasm collection. GWAS revealed 22 significant associations between DArTseq-derived SNP markers and rice grain quality traits in the test genotypes. In general, 2 of the 22 significant associations were in chromosomal regions where the QTLs associated with the given traits had previously been reported, the other 20 significant SNP marker loci were indicative of the likelihood discovery of novel alleles associated with rice grain quality traits. DArTseq-derived SNP markers that include SNP12_100006178, SNP13_3052560 and SNP14_3057360 individually co-localised with two functional gene groups that were associated with QTLs for grain width and grain length to width ratio on chromosome 3, indicating trait dependency or pleiotropic-effect loci. This study demonstrated that DArTseq markers were useful genomic resources for genome-wide association studies of rice grain quality traits to accelerate varietal development and release.

Molecular genetic diversity and population structure of Ethiopian white lupin landraces: Implications for breeding and conservation.

White lupin is one of the four economically important species of the Lupinus genus and is an important grain legume in the Ethiopian farming system. However, there has been limited research effort to characterize the Ethiopian white lupin landraces. Fifteen polymorphic simple sequence repeat (SSR) markers were used to assess the genetic diversity and population structure of 212 Ethiopian white lupin (Lupinus albus) landraces and two genotypes from different species (Lupinus angustifolius and Lupinus mutabilis) were used as out-group. The SSR markers revealed 108 different alleles, 98 of them from 212 landraces and 10 from out-group genotypes, with an average of 6.5 alleles per locus. The average gene diversity was 0.31. Twenty eight landraces harbored one or more private alleles from the total of 28 private alleles identified in the 212 white lupin accessions. Seventy-seven rare alleles with a frequency of less than 5% were identified and accounted for 78.6% of the total alleles detected. Analysis of molecular variance (AMOVA) showed that 92% of allelic diversity was attributed to individual accessions within populations while only 8% was distributed among populations. At 70% similarity level, the UPGMA dendrogram resulted in the formation of 13 clusters comprised of 2 to 136 landraces, with the out-group genotypes and five landraces remaining distinct and ungrouped. Population differentiation and genetic distance were relatively high between Gondar and Ethiopian white lupin populations collected by Australians. A model-based population structure analysis divided the white lupin landraces into two populations. All Ethiopian white lupin landrace populations, except most of the landraces collected by Australians (77%) and about 44% from Awi, were grouped together with significant admixtures. The study also suggested that 34 accessions, as core collections, were sufficient to retain 100% of SSR diversity. These accessions (core G-34) represent 16% of the whole 212 Ethiopian white lupin accessions and populations from West Gojam, Awi and Australian collections contributed more accessions to the core collection.

QTL Mapping in Three Rice Populations Uncovers Major Genomic Regions Associated with African Rice Gall Midge Resistance.

African rice gall midge (AfRGM) is one of the most destructive pests of irrigated and lowland African ecologies. This study aimed to identify the quantitative trait loci (QTL) associated with AfRGM pest incidence and resistance in three independent bi-parental rice populations (ITA306xBW348-1, ITA306xTOG7106 and ITA306xTOS14519), and to conduct meta QTL (mQTL) analysis to explore whether any genomic regions are conserved across different genetic backgrounds. Composite interval mapping (CIM) conducted on the three populations independently uncovered a total of 28 QTLs associated with pest incidence (12) and pest severity (16). The number of QTLs per population associated with AfRGM resistance varied from three in the ITA306xBW348-1 population to eight in the ITA306xTOG7106 population. Each QTL individually explained 1.3 to 34.1% of the phenotypic variance. The major genomic region for AfRGM resistance had a LOD score and R2 of 60.0 and 34.1% respectively, and mapped at 111 cM on chromosome 4 (qAfrGM4) in the ITA306xTOS14519 population. The meta-analysis reduced the number of QTLs from 28 to 17 mQTLs, each explaining 1.3 to 24.5% of phenotypic variance, and narrowed the confidence intervals by 2.2 cM. There was only one minor effect mQTL on chromosome 1 that was common in the TOS14519 and TOG7106 genetic backgrounds; all other mQTLs were background specific. We are currently fine-mapping and validating the major effect genomic region on chromosome 4 (qAfRGM4). This is the first report in mapping the genomic regions associated with the AfRGM resistance, and will be highly useful for rice breeders.