Genome-wide association study for yield and quality of granulated cassava processed product.

The starchy storage roots of cassava are commonly processed into a variety of products, including cassava granulated processed products (gari). The commercial value of cassava roots depends on the yield and quality of processed products, directly influencing the acceptance of new varieties by farmers, processors, and consumers. This study aims to estimate genetic advance through phenotypic selection and identify genomic regions associated and candidate genes linked with gari yield and quality. Higher single nucleotide polymorphism (SNP)-based heritability estimates compared to broad-sense heritability estimates were observed for most traits highlighting the influence of genetic factors on observed variation. Using genome-wide association analysis of 188 clones, genotyped using 53,150 genome-wide SNPs, nine SNPs located on seven chromosomes were significantly associated with peel loss, gari yield, color parameters for gari and eba, bulk density, swelling index, and textural properties of eba. Future research will focus on validating and understanding the functions of identified genes and their influence on gari yield and quality traits.

Q&A: Methods for estimating genetic gain in sub-Saharan Africa and achieving improved gains.

Regular measurement of realized genetic gain allows plant breeders to assess and review the effectiveness of their strategies, allocate resources efficiently, and make informed decisions throughout the breeding process. Realized genetic gain estimation requires separating genetic trends from nongenetic trends using the linear mixed model (LMM) on historical multi-environment trial data. The LMM, accounting for the year effect, experimental designs, and heterogeneous residual variances, estimates best linear unbiased estimators of genotypes and regresses them on their years of origin. An illustrative example of estimating realized genetic gain was provided by analyzing historical data on fresh cassava (Manihot esculenta Crantz) yield in West Africa (https://github.com/Biometrics-IITA/Estimating-Realized-Genetic-Gain). This approach can serve as a model applicable to other crops and regions. Modernization of breeding programs is necessary to maximize the rate of genetic gain. This can be achieved by adopting genomics to enable faster breeding, accurate selection, and improved traits through genomic selection and gene editing. Tracking operational costs, establishing robust, digitalized data management and analytics systems, and developing effective varietal selection processes based on customer insights are also crucial for success. Capacity building and collaboration of breeding programs and institutions also play a significant role in accelerating genetic gains.

Targeted knockout of early nodulin-like 3 (MusaENODL3) gene in banana reveals its function in resistance to Xanthomonas wilt disease.

Nodulins and nodulin-like proteins play an essential role in the symbiotic associations between legumes and Rhizobium bacteria. Their role extends beyond the leguminous species, as numerous nodulin-like proteins, including early nodulin-like proteins (ENODL), have been identified in various non-leguminous plants, implying their involvement in functions beyond nodulation, such as nutrient transport and growth modulation. Some ENODL proteins have been associated with plant defense against pathogens, as evident in banana infected with Xanthomonas campestris pv. musacearum (Xcm) causing banana Xanthomonas wilt (BXW) disease. Nonetheless, the specific role of ENODL in plant defense remains to be fully elucidated. The MusaENODL3 gene was found to be repressed in BXW-resistant banana progenitor 'Musa balbisiana' and 20-fold upregulated in BXW-susceptible cultivar 'Gonja Manjaya' upon early infection with Xcm. To further unravel the role of the ENODL gene in disease resistance, the CRISPR/Cas9 system was employed to disrupt the MusaENODL3 gene in 'Gonja Manjaya' precisely. Analysis of the enodl3 edited events confirmed the accurate manipulation of the MusaENODL3 gene. Disease resistance and gene expression analysis demonstrated that editing the MusaENODL3 gene resulted in resistance to BXW disease, with 50% of the edited plants remaining asymptomatic. The identification and manipulation of the MusaENODL3 gene highlight its potential as a critical player in plant-pathogen interactions, offering new opportunities for enhancing disease resistance in crops like banana, an important staple food crop and source of income for resource-poor farmers in the tropics. This study provides the first evidence of the direct role of the ENODL3 gene in developing disease-resistant plants.

Haplotype-resolved genome of heterozygous African cassava cultivar TMEB117 (Manihot esculenta).

Cassava (Manihot esculenta Crantz) is a vital tropical root crop providing essential dietary energy to over 800 million people in tropical and subtropical regions. As a climate-resilient crop, its significance grows as the human population expands. However, yield improvement faces challenges from biotic and abiotic stress and limited breeding. Advanced sequencing and assembly techniques enabled the generation of a highly accurate, nearly complete, haplotype-resolved genome of the African cassava cultivar TMEB117. It is the most accurate cassava genome sequence to date with a base-level accuracy of QV > 64, N50 > 35 Mbp, and 98.9% BUSCO completeness. Over 60% of the genome comprises repetitive elements. We predicted over 45,000 gene models for both haplotypes. This achievement offers valuable insights into the heterozygosity genome organization of the cassava genome, with improved accuracy, completeness, and phased genomes. Due to its high susceptibility to African Cassava Mosaic Virus (ACMV) infections compared to other cassava varieties, TMEB117 provides an ideal reference for studying virus resistance mechanisms, including epigenetic variations and smallRNA expressions.

Candidate genes for field resistance to cassava brown streak disease revealed through the analysis of multiple data sources.

Cassava (Manihot esculenta Crantz) is a food and industrial storage root crop with substantial potential to contribute to managing risk associated with climate change due to its inherent resilience and in providing a biodegradable option in manufacturing. In Africa, cassava production is challenged by two viral diseases, cassava brown streak disease (CBSD) and cassava mosaic disease. Here we detect quantitative trait loci (QTL) associated with CBSD in a biparental mapping population of a Tanzanian landrace, Nachinyaya and AR37-80, phenotyped in two locations over three years. The purpose was to use the information to ultimately facilitate either marker-assisted selection or adjust weightings in genomic selection to increase the efficiency of breeding. Results from this study were considered in relation to those from four other biparental populations, of similar genetic backgrounds, that were phenotyped and genotyped simultaneously. Further, we investigated the co-localization of QTL for CBSD resistance across populations and the genetic relationships of parents based on whole genome sequence information. Two QTL on chromosome 4 for resistance to CBSD foliar symptoms and one on each of chromosomes 11 and 18 for root necrosis were of interest. Of significance within the candidate genes underlying the QTL on chromosome 4 are Phenylalanine ammonia-lyase (PAL) and Cinnamoyl-CoA reductase (CCR) genes and three PEPR1-related kinases associated with the lignin pathway. In addition, a CCR gene was also underlying the root necrosis-resistant QTL on chromosome 11. Upregulation of key genes in the cassava lignification pathway from an earlier transcriptome study, including PAL and CCR, in a CBSD-resistant landrace compared to a susceptible landrace suggests a higher level of basal lignin deposition in the CBSD-resistant landrace. Earlier RNAscope® in situ hybridisation imaging experiments demonstrate that cassava brown streak virus (CBSV) is restricted to phloem vessels in CBSV-resistant varieties, and phloem unloading for replication in mesophyll cells is prevented. The results provide evidence for the involvement of the lignin pathway. In addition, five eukaryotic initiation factor (eIF) genes associated with plant virus resistance were found within the priority QTL regions.

The Progression in Developing Genomic Resources for Crop Improvement.

Sequencing technologies have rapidly evolved over the past two decades, and new technologies are being continually developed and commercialized. The emerging sequencing technologies target generating more data with fewer inputs and at lower costs. This has also translated to an increase in the number and type of corresponding applications in genomics besides enhanced computational capacities (both hardware and software). Alongside the evolving DNA sequencing landscape, bioinformatics research teams have also evolved to accommodate the increasingly demanding techniques used to combine and interpret data, leading to many researchers moving from the lab to the computer. The rich history of DNA sequencing has paved the way for new insights and the development of new analysis methods. Understanding and learning from past technologies can help with the progress of future applications. This review focuses on the evolution of sequencing technologies, their significant enabling role in generating plant genome assemblies and downstream applications, and the parallel development of bioinformatics tools and skills, filling the gap in data analysis techniques.

The banana genome hub: a community database for genomics in the Musaceae.

The Banana Genome Hub provides centralized access for genome assemblies, annotations, and the extensive related omics resources available for bananas and banana relatives. A series of tools and unique interfaces are implemented to harness the potential of genomics in bananas, leveraging the power of comparative analysis, while recognizing the differences between datasets. Besides effective genomic tools like BLAST and the JBrowse genome browser, additional interfaces enable advanced gene search and gene family analyses including multiple alignments and phylogenies. A synteny viewer enables the comparison of genome structures between chromosome-scale assemblies. Interfaces for differential expression analyses, metabolic pathways and GO enrichment were also added. A catalogue of variants spanning the banana diversity is made available for exploration, filtering, and export to a wide variety of software. Furthermore, we implemented new ways to graphically explore gene presence-absence in pangenomes as well as genome ancestry mosaics for cultivated bananas. Besides, to guide the community in future sequencing efforts, we provide recommendations for nomenclature of locus tags and a curated list of public genomic resources (assemblies, resequencing, high density genotyping) and upcoming resources-planned, ongoing or not yet public. The Banana Genome Hub aims at supporting the banana scientific community for basic, translational, and applied research and can be accessed at https://banana-genome-hub.southgreen.fr.

Breedbase: a digital ecosystem for modern plant breeding.

Modern breeding methods integrate next-generation sequencing and phenomics to identify plants with the best characteristics and greatest genetic merit for use as parents in subsequent breeding cycles to ultimately create improved cultivars able to sustain high adoption rates by farmers. This data-driven approach hinges on strong foundations in data management, quality control, and analytics. Of crucial importance is a central database able to (1) track breeding materials, (2) store experimental evaluations, (3) record phenotypic measurements using consistent ontologies, (4) store genotypic information, and (5) implement algorithms for analysis, prediction, and selection decisions. Because of the complexity of the breeding process, breeding databases also tend to be complex, difficult, and expensive to implement and maintain. Here, we present a breeding database system, Breedbase (https://breedbase.org/, last accessed 4/18/2022). Originally initiated as Cassavabase (https://cassavabase.org/, last accessed 4/18/2022) with the NextGen Cassava project (https://www.nextgencassava.org/, last accessed 4/18/2022), and later developed into a crop-agnostic system, it is presently used by dozens of different crops and projects. The system is web based and is available as open source software. It is available on GitHub (https://github.com/solgenomics/, last accessed 4/18/2022) and packaged in a Docker image for deployment (https://hub.docker.com/u/breedbase, last accessed 4/18/2022). The Breedbase system enables breeding programs to better manage and leverage their data for decision making within a fully integrated digital ecosystem.

Geographical Distribution and Genetic Diversity of the Banana Fusarium Wilt Fungus in Laos and Vietnam.

Fusarium wilt, caused by the fungus Fusarium oxysporum f. sp. cubense (Foc), poses a major threat to global banana production. The tropical race 4 (TR4) variant of Foc is a highly virulent form with a large host range, and severely affects Cavendish bananas. Foc TR4 was recently observed within the Greater Mekong Subregion, after Chinese private companies expanded Cavendish production to the region. In this study, extensive surveys conducted across Laos and Vietnam show that Foc TR4 is still mainly constricted to the northern regions of these countries and is limited to Cavendish cultivation settings. In Laos, Foc TR4 is associated with large-scale Cavendish plantations owned by or involved with Chinese companies through which infected planting material could have been imported. In Vietnam, mostly small-holder Cavendish farmers and backyard gardens were affected by Foc TR4. In Vietnam, no direct link is found with Chinese growers, and it is expected the pathogen mainly spreads through local and regional movement of infected planting materials. Foc TR4 was not recorded on banana cultivars other than Cavendish. The extensively cultivated 'Pisang Awak' cultivar was solely infected by VCGs belonging to Foc race 1 and 2, with a high occurrence of VCG 0123 across Laos, and of VCG 0124/5 in Vietnam. Substantial diversity of Foc VCGs was recorded (VCGs 0123, 0124/5, 01218 and 01221) from northern to southern regions in both countries, suggesting that Fusarium wilt is well established in the region. Interviews with farmers indicated that the local knowledge of Fusarium wilt epidemiology and options for disease management was limited. Clear communication efforts on disease epidemiology and management with emphasis on biosecurity practices need to be improved in order to prevent further spread of Foc TR4 to mixed variety smallholder settings.

Continuous Mapping Identifies Loci Associated With Weevil Resistance [Cosmopolites sordidus (Germar)] in a Triploid Banana Population.

The first step toward marker-assisted selection is linking the phenotypes to molecular markers through quantitative trait loci (QTL) analysis. While the process is straightforward in self-pollinating diploid (2x) species, QTL analysis in polyploids requires unconventional methods. In this study, we have identified markers associated with weevil Cosmopolites sordidus (Germar) resistance in bananas using 138 triploid (2n = 3x) hybrids derived from a cross between a tetraploid "Monyet" (2n = 4x) and a 2x "Kokopo" (2n = 2x) banana genotypes. The population was genotyped by Diversity Arrays Technology Sequencing (DArTSeq), resulting in 18,009 polymorphic single nucleotide polymorphisms (SNPs) between the two parents. Marker-trait association was carried out by continuous mapping where the adjusted trait means for the corm peripheral damage (PD) and total cross-section damage (TXD), both on the logit scale, were regressed on the marker allele frequencies. Forty-four SNPs that were associated with corm PD were identified on the chromosomes 5, 6, and 8, with 41 of them located on chromosome 6 and segregated in "Kokopo." Eleven SNPs associated with corm total TXD were identified on chromosome 6 and segregated in "Monyet." The additive effect of replacing one reference allele with the alternative allele was determined at each marker position. The PD QTL was confirmed using conventional QTL linkage analysis in the simplex markers segregating in "Kokopo" (AAAA × RA). We also identified 43 putative genes in the vicinity of the markers significantly associated with the two traits. The identified loci associated with resistance to weevil damage will be used in the efforts of developing molecular tools for marker-assisted breeding in bananas.

Genetic analysis and QTL mapping for multiple biotic stress resistance in cassava.

In Sub-Saharan Africa cassava (Manihot esculenta Crantz) is one of the most important food crops where more than 40% of the population relies on it as their staple carbohydrate source. Biotic constraints such as viral diseases, mainly Cassava Mosaic Disease (CMD) and Cassava Brown Streak Disease (CBSD), and arthropod pests, particularly Cassava Green Mite (CGM), are major constraints to the realization of cassava's full production potential in Africa. To address these problems, we aimed to map the quantitative trait loci (QTL) associated with resistance to CBSD foliar and root necrosis symptoms, foliar CMD and CGM symptoms in a full-sib mapping population derived from the genotypes AR40-6 and Albert. A high-density linkage map was constructed with 2,125 SNP markers using a genotyping-by-sequencing approach. For phenotyping, clonal evaluation trials were conducted with 120 F1 individuals for two consecutive field seasons using an alpha-lattice design at Chambezi and Naliendele, Tanzania. Previously identified QTL for resistance to CBSD foliar symptoms were corroborated, and a new putative QTL for CBSD root necrosis identified (qCBSDRNc14AR) from AR40-6. Two QTL were identified within the region of the previously recognized CMD2 locus from this population in which both parents are thought to possess the CMD2 locus. Interestingly, a minor but consistent QTL, qCGM18AR, for CGM resistance at 3 months after planting stage was also detected and co-localized with a previously identified SSR marker, NS346, linked with CGM resistance. Markers underlying these QTL may be used to increase efficiencies in cassava breeding programs.

A global overview of cassava genetic diversity.

Although numerous studies of diversity have been conducted in cassava, there is no comprehensive assessment of global genetic diversity. Here we draw on previous studies and breeders' knowledge to select diversity sets from the International Institute of Tropical Agriculture (IITA) and the International Center for Tropical Agriculture (CIAT) genebanks and breeders' germplasm, as well as elite germplasm and landraces from eastern, southern and central (ESC) Africa to make a global assessment of diversity in cassava, using a SNP based GoldenGate (Illumina Inc.) assay. A synthesis of results from genetic distance and ADMIXTURE analysis essentially revealed four populations (i) South American germplasm characterised by relatively higher genetic diversity with hypothetical ancestral founder genotypes from Brazil, (ii) a smaller group of African introduction germplasm which is more distantly related to all other germplasm, (iii) West Africa germplasm dominated by IITA breeding lines, containing sources of cassava mosaic disease resistance, and IITA genebank accessions from West Africa, both characterised by slightly lower diversity, and (iv) a less cohesive group of African germplasm, termed 'Other', with moderate levels of diversity and a majority of germplasm from ESC Africa. This study highlights opportunities for heterosis breeding, purging of duplicates in genebanks and the need for conservation of ESC Africa landraces.

Molecular Basis of Disease Resistance in Banana Progenitor Musa balbisiana against Xanthomonas campestris pv. musacearum.

Banana Xanthomonas wilt disease, caused by Xanthomonas campestris pv. musacearum (Xcm), is a major threat to banana production in east Africa. All cultivated varieties of banana are susceptible to Xcm and only the progenitor species Musa balbisiana was found to be resistant. The molecular basis of susceptibility and resistance of banana genotypes to Xcm is currently unknown. Transcriptome analysis of disease resistant genotype Musa balbisiana and highly susceptible banana cultivar Pisang Awak challenged with Xcm was performed to understand the disease response. The number of differentially expressed genes (DEGs) was higher in Musa balbisiana in comparison to Pisang Awak. Genes associated with response to biotic stress were up-regulated in Musa balbisiana. The DEGs were further mapped to the biotic stress pathways. Our results suggested activation of both PAMP-triggered basal defense and disease resistance (R) protein-mediated defense in Musa balbisiana as early response to Xcm infection. This study reports the first comparative transcriptome profile of the susceptible and resistant genotype of banana during early infection with Xcm and provide insights on the defense mechanism in Musa balbisiana, which can be used for genetic improvement of commonly cultivated banana varieties.

Agromorphologic, genetic and methylation profiling of Dioscorea and Musa species multiplied under three micropropagation systems.

Plant in vitro vegetative propagation using classical semi-solid culture medium is limited due to the low degree of automation, suboptimal nutrient availability and induced physiological stress which often reduce its efficiency. Temporary Immersion System (TIS) emerged as an innovative approach to optimize and eliminate the drawbacks associated with the conventional system of micropropagation. In this study, both Dioscorea and Musa spp. were subjected to conventional semi-solid culture media, complete immersion in shaking liquid culture media and TIS using RITA bioreactor. In vitro grown plantlets were screened for possible vegetative changes using agro-morphological descriptors while genetic and methylation differences were assessed using amplified fragment length polymorphism (AFLP) and methylation-sensitive amplification polymorphism (MSAP). In vitro results showed that the number of shoots produced in Musa spp. varied significantly (P≤0.001) with the type of culture system. The highest mean shoot produced was observed with TIS (28.40) and the least using semi-solid culture medium (1.13). For Dioscorea spp., there was no significant interaction between the hormone combination and the culture system. However, the lowest mean shoot value (1.55) was observed in the semi-solid culture medium. Genetic analysis via AFLP using 15 primer pair combinations revealed that the 3 culture systems maintained genetic variation for Musa and Dioscorea spp. under in vitro and field conditions. Results showed 99% and 91% of the total bands were polymorphic under in vitro and field conditions respectively for Musa and 100% polymorphism for Dioscorea under in vitro and field conditions. Methylation investigation via MSAP using 12 primer pair combinations showed 25% and 46% polymorphic methylated-sensitive loci, 100% and 78% of non-methylated loci of the total bands generated under in vitro and field conditions respectively. Unmethylated (HPA+/MSP+) levels were highest in TIS (0.0842) as compared to CI (0.0227) and SS (0.0161) while full methylation or absence of target (HPA-/MSP-) was lowest in TIS (0.5890) and highest in SS (0.7138). For Dioscorea, 52% and 53% methylated sensitive loci and 100% non-methylated loci were polymorphic under in vitro and field conditions respectively. Although in vitro plant tissue culture techniques led to methylation at some loci of both species, there were no observable changes in the phenotype of both crops under field conditions. This also confirmed that not all methylation events lead to phenotypic changes.

Genome-Wide Identification and Analysis of Arabidopsis Sodium Proton Antiporter (NHX) and Human Sodium Proton Exchanger (NHE) Homologs in Sorghum bicolor.

Na⁺ transporters play an important role during salt stress and development. The present study is aimed at genome-wide identification, in silico analysis of sodium-proton antiporter (NHX) and sodium-proton exchanger (NHE)-type transporters in Sorghum bicolor and their expression patterns under varied abiotic stress conditions. In Sorghum, seven NHX and nine NHE homologs were identified. Amiloride (a known inhibitor of Na⁺/H⁺ exchanger activity) binding motif was noticed in both types of the transporters. Chromosome 2 was found to be a hotspot region with five sodium transporters. Phylogenetic analysis inferred six ortholog and three paralog groups. To gain an insight into functional divergence of SbNHX/NHE transporters, real-time gene expression was performed under salt, drought, heat, and cold stresses in embryo, root, stem, and leaf tissues. Expression patterns revealed that both SbNHXs and SbNHEs are responsive either to single or multiple abiotic stresses. The predicted protein⁻protein interaction networks revealed that only SbNHX7 is involved in the calcineurin B-like proteins (CBL)- CBL interacting protein kinases (CIPK) pathway. The study provides insights into the functional divergence of SbNHX/NHE transporter genes with tissue specific expressions in Sorghum under different abiotic stress conditions.

QTL Mapping for Pest and Disease Resistance in Cassava and Coincidence of Some QTL with Introgression Regions Derived from Manihot glaziovii.

Genetic mapping of quantitative trait loci (QTL) for resistance to cassava brown streak disease (CBSD), cassava mosaic disease (CMD), and cassava green mite (CGM) was performed using an F1 cross developed between the Tanzanian landrace, Kiroba, and a breeding line, AR37-80. The population was evaluated for two consecutive years in two sites in Tanzania. A genetic linkage map was derived from 106 F1 progeny and 1,974 SNP markers and spanned 18 chromosomes covering a distance of 1,698 cM. Fifteen significant QTL were identified; two are associated with CBSD root necrosis only, and were detected on chromosomes V and XII, while seven were associated with CBSD foliar symptoms only and were detected on chromosomes IV, VI, XVII, and XVIII. QTL on chromosomes 11 and 15 were associated with both CBSD foliar and root necrosis symptoms. Two QTL were found to be associated with CMD and were detected on chromosomes XII and XIV, while two were associated with CGM and were identified on chromosomes V and X. There are large Manihot glaziovii introgression regions in Kiroba on chromosomes I, XVII, and XVIII. The introgression segments on chromosomes XVII and XVIII overlap with QTL associated with CBSD foliar symptoms. The introgression region on chromosome I is of a different haplotype to the characteristic "Amani haplotype" found in the landrace Namikonga and others, and unlike some other genotypes, Kiroba does not have a large introgression block on chromosome IV. Kiroba is closely related to a sampled Tanzanian "tree cassava." This supports the observation that some of the QTL associated with CBSD resistance in Kiroba are different to those observed in another variety, Namikonga.

Trait Associations in Diversity Panels of the Two Common Bean (Phaseolus vulgaris L.) Gene Pools Grown under Well-watered and Water-Stress Conditions.

Common beans are a warm-season, food legume cultivated in areas prone to water limitation throughout their growing season. This study assessed the magnitude and pattern of trait associations for a total of 202 common bean genotypes divided into panels of 81 Andean and 121 Mesoamerican gene pool accessions grown under contrasting treatments of well-watered, non-stress, and water-limited, terminal drought-stress conditions. Linear correlation, complex path coefficient, and genetic divergence analyses were used to dissect the relationship dynamics between traits and the relative contribution of adaptive traits to differentiation among gene pools and genotypes based on drought stress. Drought severity level for the trial was high and created the ideal condition to reveal genotypic differences, as seen by the differential response of the genotypes for the various traits measured. The value for phenotypic coefficients of variation for all traits was higher than the corresponding genotypic values. Seed yield had positive and strong genotypic and phenotypic correlation with pods per plant across gene pools and stress levels. The overall amount of genetic correlation was greater than the corresponding phenotypic correlation matrix for all the traits within the gene pool and across stress levels. Moreover, the results depicted the phenotypic correlation as equal or better than its genotypic counterpart in estimating drought tolerance in common bean plants. Clustering analysis with Mahanalobis's coefficient of generalized distance grouped genotypes with a differential level of drought adaptation into different classes within each panel. This indicates drought tolerance involves different mechanisms of plant response and is present separately in each gene pool panel. Pods per plant, seed weight, pod partitioning index, and harvest index are useful selection objectives to improve drought adaptation in common bean, but must be differentially weighted in each gene pool. The analysis of genetic variation and association between adaptive traits on the two panels provided useful insights on which traits could be used to improve common bean adaptation to low water availability during the growth season.

Functional mechanisms of drought tolerance in subtropical maize (Zea mays L.) identified using genome-wide association mapping.

BACKGROUND: Earlier studies were focused on the genetics of temperate and tropical maize under drought. We identified genetic loci and their association with functional mechanisms in 240 accessions of subtropical maize using a high-density marker set under water stress. RESULTS: Out of 61 significant SNPs (11 were false-discovery-rate-corrected associations), identified across agronomic traits, models, and locations by subjecting the accessions to water stress at flowering stage, 48% were associated with drought-tolerant genes. Maize gene models revealed that SNPs mapped for agronomic traits were in fact associated with number of functional traits as follows: stomatal closure, 28; flowering, 15; root development, 5; detoxification, 4; and reduced water potential, 2. Interactions of these SNPS through the functional traits could lead to drought tolerance. The SNPs associated with ABA-dependent signalling pathways played a major role in the plant's response to stress by regulating a series of functions including flowering, root development, auxin metabolism, guard cell functions, and scavenging reactive oxygen species (ROS). ABA signalling genes regulate flowering through epigenetic changes in stress-responsive genes. ROS generated by ABA signalling are reduced by the interplay between ethylene, ABA, and detoxification signalling transductions. Integration of ABA-signalling genes with auxin-inducible genes regulates root development which in turn, maintains the water balance by regulating electrochemical gradient in plant. CONCLUSIONS: Several genes are directly or indirectly involved in the functioning of agronomic traits related to water stress. Genes involved in these crucial biological functions interacted significantly in order to maintain the primary as well as exclusive functions related to coping with water stress. SNPs associated with drought-tolerant genes involved in strategic biological functions will be useful to understand the mechanisms of drought tolerance in subtropical maize.

An integrated SNP mining and utilization (ISMU) pipeline for next generation sequencing data.

Open source single nucleotide polymorphism (SNP) discovery pipelines for next generation sequencing data commonly requires working knowledge of command line interface, massive computational resources and expertise which is a daunting task for biologists. Further, the SNP information generated may not be readily used for downstream processes such as genotyping. Hence, a comprehensive pipeline has been developed by integrating several open source next generation sequencing (NGS) tools along with a graphical user interface called Integrated SNP Mining and Utilization (ISMU) for SNP discovery and their utilization by developing genotyping assays. The pipeline features functionalities such as pre-processing of raw data, integration of open source alignment tools (Bowtie2, BWA, Maq, NovoAlign and SOAP2), SNP prediction (SAMtools/SOAPsnp/CNS2snp and CbCC) methods and interfaces for developing genotyping assays. The pipeline outputs a list of high quality SNPs between all pairwise combinations of genotypes analyzed, in addition to the reference genome/sequence. Visualization tools (Tablet and Flapjack) integrated into the pipeline enable inspection of the alignment and errors, if any. The pipeline also provides a confidence score or polymorphism information content value with flanking sequences for identified SNPs in standard format required for developing marker genotyping (KASP and Golden Gate) assays. The pipeline enables users to process a range of NGS datasets such as whole genome re-sequencing, restriction site associated DNA sequencing and transcriptome sequencing data at a fast speed. The pipeline is very useful for plant genetics and breeding community with no computational expertise in order to discover SNPs and utilize in genomics, genetics and breeding studies. The pipeline has been parallelized to process huge datasets of next generation sequencing. It has been developed in Java language and is available at http://hpc.icrisat.cgiar.org/ISMU as a standalone free software.