Functional allele of a MATE gene selected during domestication modulates seed color in chickpea.

Seed color is one of the key target traits of domestication and artificial selection in chickpeas due to its implications on consumer preference and market value. The complex seed color trait has been well dissected in several crop species; however, the genetic mechanism underlying seed color variation in chickpea remains poorly understood. Here, we employed an integrated genomics strategy involving QTL mapping, high-density mapping, map-based cloning, association analysis, and molecular haplotyping in an inter-specific RIL mapping population, association panel, wild accessions, and introgression lines (ILs) of Cicer gene pool. This delineated a MATE gene, CaMATE23, encoding a Transparent Testa (TT) and its natural allele (8-bp insertion) and haplotype underlying a major QTL governing seed color on chickpea chromosome 4. Signatures of selective sweep and a strong purifying selection reflected that CaMATE23, especially its 8-bp insertion natural allelic variant, underwent selection during chickpea domestication. Functional investigations revealed that the 8-bp insertion containing the third cis-regulatory RY-motif element in the CaMATE23 promoter is critical for enhanced binding of CaFUSCA3 transcription factor, a key regulator of seed development and flavonoid biosynthesis, thereby affecting CaMATE23 expression and proanthocyanidin (PA) accumulation in the seed coat to impart varied seed color in chickpea. Consequently, overexpression of CaMATE23 in Arabidopsis tt12 mutant partially restored the seed color phenotype to brown pigmentation, ascertaining its functional role in PA accumulation in the seed coat. These findings shed new light on the seed color regulation and evolutionary history, and highlight the transcriptional regulation of CaMATE23 by CaFUSCA3 in modulating seed color in chickpea. The functionally relevant InDel variation, natural allele, and haplotype from CaMATE23 are vital for translational genomic research, including marker-assisted breeding, for developing chickpea cultivars with desirable seed color that appeal to consumers and meet global market demand.

Rice Pangenome Genotyping Array: an efficient genotyping solution for pangenome-based accelerated genetic improvement in rice.

The advent of the pangenome era has unraveled previously unknown genetic variation existing within diverse crop plants, including rice. This untapped genetic variation is believed to account for a major portion of phenotypic variation existing in crop plants. However, the use of conventional single reference-guided genotyping often fails to capture a large portion of this genetic variation leading to a reference bias. This makes it difficult to identify and utilize novel population/cultivar-specific genes for crop improvement. Thus, we developed a Rice Pangenome Genotyping Array (RPGA) harboring probes assaying 80K single-nucleotide polymorphisms (SNPs) and presence-absence variants spanning the entire 3K rice pangenome. This array provides a simple, user-friendly and cost-effective (60-80 USD per sample) solution for rapid pangenome-based genotyping in rice. The genome-wide association study (GWAS) conducted using RPGA-SNP genotyping data of a rice diversity panel detected a total of 42 loci, including previously known as well as novel genomic loci regulating grain size/weight traits in rice. Eight of these identified trait-associated loci (dispensable loci) could not be detected with conventional single reference genome-based GWAS. A WD repeat-containing PROTEIN 12 gene underlying one of such dispensable locus on chromosome 7 (qLWR7) along with other non-dispensable loci were subsequently detected using high-resolution quantitative trait loci mapping confirming authenticity of RPGA-led GWAS. This demonstrates the potential of RPGA-based genotyping to overcome reference bias. The application of RPGA-based genotyping for population structure analysis, hybridity testing, ultra-high-density genetic map construction and chromosome-level genome assembly, and marker-assisted selection was also demonstrated. A web application (http://www.rpgaweb.com) was further developed to provide an easy to use platform for the imputation of RPGA-based genotyping data using 3K rice reference panel and subsequent GWAS.

A superior gene allele involved in abscisic acid signaling enhances drought tolerance and yield in chickpea.

Identifying potential molecular tags for drought tolerance is essential for achieving higher crop productivity under drought stress. We employed an integrated genomics-assisted breeding and functional genomics strategy involving association mapping, fine mapping, map-based cloning, molecular haplotyping and transcript profiling in the introgression lines (ILs)- and near isogenic lines (NILs)-based association panel and mapping population of chickpea (Cicer arietinum). This combinatorial approach delineated a bHLH (basic helix-loop-helix) transcription factor, CabHLH10 (Cicer arietinum bHLH10) underlying a major QTL, along with its derived natural alleles/haplotypes governing yield traits under drought stress in chickpea. CabHLH10 binds to a cis-regulatory G-box promoter element to modulate the expression of RD22 (responsive to desiccation 22), a drought/abscisic acid (ABA)-responsive gene (via a trans-expression QTL), and two strong yield-enhancement photosynthetic efficiency (PE) genes. This, in turn, upregulates other downstream drought-responsive and ABA signaling genes, as well as yield-enhancing PE genes, thus increasing plant adaptation to drought with reduced yield penalty. We showed that a superior allele of CabHLH10 introgressed into the NILs improved root and shoot biomass and PE, thereby enhancing yield and productivity during drought without compromising agronomic performance. Furthermore, overexpression of CabHLH10 in chickpea and Arabidopsis (Arabidopsis thaliana) conferred enhanced drought tolerance by improving root and shoot agro-morphological traits. These findings facilitate translational genomics for crop improvement and the development of genetically tailored, climate-resilient, high-yielding chickpea cultivars.

A multicentre retrospective cohort study on COVID-19-related physical interventions and adult hospital admissions for ENT infections.

PURPOSE: To report changes in adult hospital admission rates for acute ENT infections following the introduction of COVID-19-related physical interventions such as hand washing, use of face masks and social distancing of 2-m in the United Kingdom. METHODS: Retrospective cohort study comparing adult admissions with acute tonsillitis, peritonsillar abscess, epiglottitis, glandular fever, peri-orbital cellulitis, acute otitis media, acute mastoiditis, retropharyngeal abscess and parapharyngeal abscess in the 1-year period after the introduction of COVID-related physical interventions (2020-2021) with a 1-year period before this (2019-2020) in three UK secondary care ENT departments. RESULTS: In total, there were significantly fewer admissions for ENT infections (n = 1073, 57.56%, p < 0.001; RR 2.36, 95% CI [2.17, 2.56]) in the 2020-2021 period than in the 2019-2020 period. There were significant reductions in admissions for tonsillitis (64.4%; p < 0.001), peritonsillar abscess (60.68%; p < 0.001), epiglottitis (66.67%; p < 0.001), glandular fever (38.79%; p = 0.001), acute otitis media (26.85%; p = 0.01) and retropharyngeal and/or parapharyngeal abscesses (45.45%; p = 0.04). CONCLUSION: Our study demonstrates a sizeable reduction in adult admissions for ENT infections since the introduction of COVID-19-related physical interventions. There is evidence to support the use of physical interventions in the prevention of viral transmission of respiratory disease. Preventing ENT infections requiring admission through simple physical interventions could be of great benefit to the quality of life of patients and economical benefit to healthcare systems.

ABC Transporter-Mediated Transport of Glutathione Conjugates Enhances Seed Yield and Quality in Chickpea.

The identification of functionally relevant molecular tags is vital for genomics-assisted crop improvement and enhancement of seed yield, quality, and productivity in chickpea (Cicer arietinum). The simultaneous improvement of yield/productivity as well as quality traits often requires pyramiding of multiple genes, which remains a major hurdle given various associated epistatic and pleotropic effects. Unfortunately, no single gene that can improve yield/productivity along with quality and other desirable agromorphological traits is known, hampering the genetic enhancement of chickpea. Using a combinatorial genomics-assisted breeding and functional genomics strategy, this study identified natural alleles and haplotypes of an ABCC3-type transporter gene that regulates seed weight, an important domestication trait, by transcriptional regulation and modulation of the transport of glutathione conjugates in seeds of desi and kabuli chickpea. The superior allele/haplotype of this gene introgressed in desi and kabuli near-isogenic lines enhances the seed weight, yield, productivity, and multiple desirable plant architecture and seed-quality traits without compromising agronomic performance. These salient findings can expedite crop improvement endeavors and the development of nutritionally enriched high-yielding cultivars in chickpea.

CLAVATA signaling pathway genes modulating flowering time and flower number in chickpea.

KEY MESSAGE: A combinatorial genomic strategy delineated functionally relevant natural allele of a CLAVATA gene and its marker (haplotype)-assisted introgression led to development of the early-flowering chickpea cultivars with high flower number and enhanced yield/productivity. Unraveling the genetic components involved in CLAVATA (CLV) signaling is crucial for modulating important shoot apical meristem (SAM) characteristics and ultimately regulating diverse SAM-regulated agromorphological traits in crop plants. A genome-wide scan identified 142 CLV1-, 28 CLV2- and 6 CLV3-like genes, and their comprehensive genomic constitution and phylogenetic relationships were deciphered in chickpea. The QTL/fine mapping and map-based cloning integrated with high-resolution association analysis identified SNP loci from CaCLV3_01 gene within a major CaqDTF1.1/CaqFN1.1 QTL associated with DTF (days to 50% flowering) and FN (flower number) traits in chickpea, which was further ascertained by quantitative expression profiling. Molecular haplotyping of CaCLV3_01 gene, expressed specifically in SAM, constituted two major haplotypes that differentiated the early-DTF and high-FN chickpea accessions from late-DTF and low-FN. Enhanced accumulation of transcripts of superior CaCLV3_01 gene haplotype and known flowering promoting genes was observed in the corresponding haplotype-introgressed early-DTF and high-FN near-isogenic lines (NILs) with narrow SAM width. The superior haplotype-introgressed NILs exhibited early-flowering, high-FN and enhanced seed yield/productivity without compromising agronomic performance. These delineated molecular signatures can regulate DTF and FN traits through SAM proliferation and differentiation and thereby will be useful for translational genomic study to develop early-flowering cultivars with enhanced yield/productivity.

Genetic dissection of plant growth habit in chickpea.

A combinatorial genomics-assisted breeding strategy encompassing association analysis, genetic mapping and expression profiling is found most promising for quantitative dissection of complex traits in crop plants. The present study employed GWAS (genome-wide association study) using 24,405 SNPs (single nucleotide polymorphisms) obtained with genotyping-by-sequencing (GBS) of 92 sequenced desi and kabuli accessions of chickpea. This identified eight significant genomic loci associated with erect (E)/semi-erect (SE) vs. spreading (S)/semi-spreading (SS)/prostrate (P) plant growth habit (PGH) trait differentiation regardless of diverse desi and kabuli genetic backgrounds of chickpea. These associated SNPs in combination explained 23.8% phenotypic variation for PGH in chickpea. Five PGH-associated genes were validated successfully in E/SE and SS/S/P PGH-bearing parental accessions and homozygous individuals of three intra- and interspecific RIL (recombinant inbred line) mapping populations as well as 12 contrasting desi and kabuli chickpea germplasm accessions by selective genotyping through Sequenom MassARRAY. The shoot apical, inflorescence and floral meristems-specific expression, including upregulation (seven-fold) of five PGH-associated genes especially in germplasm accessions and homozygous RIL mapping individuals contrasting with E/SE PGH traits was apparent. Collectively, this integrated genomic strategy delineated diverse non-synonymous SNPs from five candidate genes with strong allelic effects on PGH trait variation in chickpea. Of these, two vernalization-responsive non-synonymous SNP alleles carrying SNF2 protein-coding gene and B3 transcription factor associated with PGH traits were found to be the most promising in chickpea. The SNP allelic variants associated with E/SE/SS/S PGH trait differentiation were exclusively present in all cultivated desi and kabuli chickpea accessions while wild species/accessions belonging to primary, secondary and tertiary gene pools mostly contained prostrate PGH-associated SNP alleles. This indicates strong adaptive natural/artificial selection pressure (Tajima's D 3.15 to 4.57) on PGH-associated target genomic loci during chickpea domestication. These vital leads thus have potential to decipher complex transcriptional regulatory gene function of PGH trait differentiation and for understanding the selective sweep-based PGH trait evolution and domestication pattern in cultivated and wild chickpea accessions adapted to diverse agroclimatic conditions. Collectively, the essential inputs generated will be of profound use in marker-assisted genetic enhancement to develop cultivars with desirable plant architecture of erect growth habit types in chickpea.

Quality of life, disability scores, and distress index in fungal rhinosinusitis.

Our goal was to determine quality of life (QOL), disability, and distress in the day-to-day lives of patients suffering from fungal rhinosinusitis (FRS) based on a prospective questionnaire. The study included 125 patients suffering from FRS, excluding those with acute fulminant FRS and any other comorbid illness. The data were compared with data for 50 age- and sex-matched controls who did not have any other chronic illness and obtained outpatient services from the Department of Otolaryngology and Head and Neck Surgery of our institute. Analysis showed that patients with FRS had worse QOL, with an average score of 87.34 compared with 94.15 for the control group. QOL score further decreased to 85.31 for patients with extensive disease that included intracranial or intraorbital extension and to 71.1 in those with recurrent disease. Patients showed significant disability and had decreased work efficiency in disability parameters. Distress among these patients was also greater and further increased in those with extensive disease or recurrence. We conclude that patients with FRS have worse QOL, more disability, and more distress compared with age- and sex-matched controls. This issue needs to be addressed while treating cases of FRS.

Parameter estimation in stochastic chemical kinetic models using derivative free optimization and bootstrapping.

Recent years have seen increasing popularity of stochastic chemical kinetic models due to their ability to explain and model several critical biological phenomena. Several developments in high resolution fluorescence microscopy have enabled researchers to obtain protein and mRNA data on the single cell level. The availability of these data along with the knowledge that the system is governed by a stochastic chemical kinetic model leads to the problem of parameter estimation. This paper develops a new method of parameter estimation for stochastic chemical kinetic models. There are three components of the new method. First, we propose a new expression for likelihood of the experimental data. Second, we use sample path optimization along with UOBYQA-Fit, a variant of of Powell's unconstrained optimization by quadratic approximation, for optimization. Third, we use a variant of Efron's percentile bootstrapping method to estimate the confidence regions for the parameter estimates. We apply the parameter estimation method in an RNA dynamics model of E. coli. We test the parameter estimates obtained and the confidence regions in this model. The testing of the parameter estimation method demonstrates the efficiency, reliability, and accuracy of the new method.

Comparison of finite difference based methods to obtain sensitivities of stochastic chemical kinetic models.

Sensitivity analysis is a powerful tool in determining parameters to which the system output is most responsive, in assessing robustness of the system to extreme circumstances or unusual environmental conditions, in identifying rate limiting pathways as a candidate for drug delivery, and in parameter estimation for calculating the Hessian of the objective function. Anderson [SIAM J. Numer. Anal. 50, 2237 (2012)] shows the advantages of the newly developed coupled finite difference (CFD) estimator over the common reaction path (CRP) [M. Rathinam, P. W. Sheppard, and M. Khammash, J. Chem. Phys. 132, 034103 (2010)] estimator. In this paper, we demonstrate the superiority of the CFD estimator over the common random number (CRN) estimator in a number of scenarios not considered previously in the literature, including the sensitivity of a negative log likelihood function for parameter estimation, the sensitivity of being in a rare state, and a sensitivity with fast fluctuating species. In all examples considered, the superiority of CFD over CRN is demonstrated. We also provide an example in which the CRN method is superior to the CRP method, something not previously observed in the literature. These examples, along with Anderson's results, lead to the conclusion that CFD is currently the best estimator in the class of finite difference estimators of stochastic chemical kinetic models.

Thuja occidentalis: An Unexplored Phytomedicine with Therapeutic Applications.

BACKGROUND: The recent outbreak of SARS-CoV-2 has received global attention. Due to a lack of recommended treatment regimens, the world faced various limitations resulting in improper management of the disease. Phytomedicines have played a prominent role in the prevention of various epidemics and pandemics in the past. OBJECTIVE: Here, we attempt to focus on safe and feasible use of Thuja occidentalis to manage and alleviate the panic of viral respiratory infections, including COVID-19, by strengthening an individual's immunity. The relevant information was collected from the web-based databases PubMed, Google Scholar, and MEDLINE, as well as other internet sources to review the applicability of T. occidentalis as a phytomedicine in managing respiratory infections and strengthening immunity. CONCLUSION: As important phytomedicine, and antipsychotic, T. occidentalis possesses a plethora of immunological properties that can be used effectively in the management of viral respiratory infections and has the potential to prevent further progression of the disease. Importantly, this could be only a part of the approach for treatment during the current outbreak that should be considered along with other measures.

The stochastic quasi-steady-state assumption: reducing the model but not the noise.

Highly reactive species at small copy numbers play an important role in many biological reaction networks. We have described previously how these species can be removed from reaction networks using stochastic quasi-steady-state singular perturbation analysis (sQSPA). In this paper we apply sQSPA to three published biological models: the pap operon regulation, a biochemical oscillator, and an intracellular viral infection. These examples demonstrate three different potential benefits of sQSPA. First, rare state probabilities can be accurately estimated from simulation. Second, the method typically results in fewer and better scaled parameters that can be more readily estimated from experiments. Finally, the simulation time can be significantly reduced without sacrificing the accuracy of the solution.

The Use of Lasers in Stapes Surgery.

OBJECTIVES: Otosclerosis is a disease process that usually starts around the oval window, causing fixation of the stapes, resulting in conductive hearing loss. Treatment of the conductive hearing loss caused by otosclerosis consists of either rehabilitation with hearing aids or performing surgery. Given the risks of hearing impairment and vertigo associated with the surgery, there has been a desire to advance the practice to minimize the complications. The so-called "non-contact" or "no touch" techniques with the use of various lasers are in current practice. This review article will cover the surgical aspects, the theory behind laser and the various types used in stapes surgery. It will also review the evidence of laser versus conventional stapes surgery and the comparison of different laser types. METHODS: A literature search up to December 2019 was performed using Pubmed and a nonsystematic review of appropriate articles was undertaken. Keywords used were stapes, surgery, laser, stapedectomy, and stapedotomy. RESULTS: Overall, there is no evidence to say laser fenestration is better than conventional fenestration techniques; however, with the micro drill, there is an increased risk of footplate fracture and sensorineural hearing loss. There is an increased risk of tinnitus with the laser compared to conventional techniques. Studies have favored the CO2 laser over potassium titanyl phosphate (KTP) and erbium-doped yttrium aluminium garnet (Erbium-YAG) lasers for postoperative closure of the air-bone gap; and KTP laser has less thermal, mechanical, and sound effects compared with the thulium and carbon dioxide (CO2) lasers. There is an increased risk if inner ear complications with the thulium laser. CONCLUSIONS: It can be deduced that theoretically and practically, the thulium laser is less safe compared to the KTP and CO2 lasers. The choice of laser used depends on the surgeon's preference, as well as availability, cost, side effects profile, as well as ease of use.

Genome-wide generation and genotyping of informative SNPs to scan molecular signatures for seed yield in chickpea.

We discovered 2150 desi and 2199 kabuli accessions-derived SNPs by cultivar-wise individual assembling of sequence-reads generated through genotyping-by-sequencing of 92 chickpea accessions. Subsequent large-scale validation and genotyping of these SNPs discovered 619 desi accessions-derived (DAD) SNPs, 531 kabuli accessions-derived (KAD) SNPs, 884 multiple accessions-derived (MAD) SNPs and 1083 two accessions (desi ICC 4958 and kabuli CDC Frontier)-derived (TAD) SNPs that were mapped on eight chromosomes. These informative SNPs were annotated in coding/non-coding regulatory sequence components of genes. The MAD-SNPs were efficient to detect high intra-specific polymorphic potential and wide natural allelic diversity level including high-resolution admixed-population genetic structure and precise phylogenetic relationship among 291 desi and kabuli accessions. This signifies their effectiveness in introgression breeding and varietal improvement studies targeting useful agronomic traits of chickpea. Six trait-associated genes with SNPs including quantitative trait nucleotides (QTNs) in combination explained 27.5% phenotypic variation for seed yield per plant (SYP). A pentatricopeptide repeat (PPR) gene with a synonymous-coding SNP/QTN significantly associated with SYP trait was found most-promising in chickpea. The essential information delineated can be of immense utility in genomics-assisted breeding applications to develop high-yielding chickpea cultivars.

Measuring the response of human head and neck squamous cell carcinoma to irradiation in a microfluidic model allowing customized therapy.

Radiotherapy is the standard treatment for head and neck squamous cell carcinoma (HNSCC), however, radioresistance remains a major clinical problem despite significant improvements in treatment protocols. Therapeutic outcome could potentially be improved if a patient's tumour response to irradiation could be predicted ex vivo before clinical application. The present study employed a bespoke microfluidic device to maintain HNSCC tissue whilst subjecting it to external beam irradiation and measured the responses using a panel of cell death and proliferation markers. HNSCC biopsies from five newly-presenting patients [2 lymph node (LN); 3 primary tumour (PT)] were divided into parallel microfluidic devices and replicates of each tumour were subjected to single-dose irradiation (0, 5, 10, 15 and 20 Gy). Lactate dehydrogenase (LDH) release was measured and tissue sections were stained for cytokeratin (CK), cleaved-CK18 (cCK18), phosphorylated-H2AX (γH2AX) and Ki­67 by immunohistochemistry. In addition, fragmented DNA was detected using terminal deoxynucleotidyl transferase dUTP nick end labelling (TUNEL). Compared with non­irradiated controls, higher irradiation doses resulted in elevated CK18-labelling index in two lymph nodes [15 Gy; 34.8% on LN1 and 31.7% on LN2 (p=0.006)] and a single laryngeal primary tumour (20 Gy; 31.5%; p=0.014). Significantly higher levels of DNA fragmentation were also detected in both lymph node samples and one primary tumour but at varying doses of irradiation, i.e., LN1 (20 Gy; 27.6%; p=0.047), LN2 (15 Gy; 15.3%; p=0.038) and PT3 (10 Gy; 35.2%; p=0.01). The γH2AX expression was raised but not significantly in the majority of samples. The percentage of Ki­67 positive nuclei reduced dose-dependently following irradiation. In contrast no significant difference in LDH release was observed between irradiated groups and controls. There is clear inter- and intra-patient variability in response to irradiation when measuring a variety of parameters, which offers the potential for the approach to provide clinically valuable information.

Regional Association Analysis of MetaQTLs Delineates Candidate Grain Size Genes in Rice.

Molecular mapping studies which aim to identify genetic basis of diverse agronomic traits are vital for marker-assisted crop improvement. Numerous Quantitative Trait Loci (QTLs) mapped in rice span long genomic intervals with hundreds to thousands of genes, which limits their utilization for marker-assisted genetic enhancement of rice. Although potent, fine mapping of QTLs is challenging task as it requires screening of large number of segregants to identify suitable recombination events. Association mapping offers much higher resolution as compared to QTL mapping, but detects considerable number of spurious QTLs. Therefore, combined use of QTL and association mapping strategies can provide advantages associated with both these methods. In the current study, we utilized meta-analysis approach to identify metaQTLs associated with grain size/weight in diverse Indian indica and aromatic rice accessions. Subsequently, attempt has been made to narrow-down identified grain size/weight metaQTLs through individual SNP- as well as haplotype-based regional association analysis. The study identified six different metaQTL regions, three of which were successfully revalidated, and substantially scaled-down along with GS3 QTL interval (positive control) by regional association analysis. Consequently, two potential candidate genes within two reduced metaQTLs were identified based on their differential expression profiles in different tissues/stages of rice accessions during seed development. The developed strategy has broader practical utility for rapid delineation of candidate genes and natural alleles underlying QTLs associated with complex agronomic traits in rice as well as major crop plants enriched with useful genetic and genomic information.

A Multiple QTL-Seq Strategy Delineates Potential Genomic Loci Governing Flowering Time in Chickpea.

Identification of functionally relevant potential genomic loci using an economical, simpler and user-friendly genomics-assisted breeding strategy is vital for rapid genetic dissection of complex flowering time quantitative trait in chickpea. A high-throughput multiple QTL-seq strategy was employed in two inter (Cicer arietinum desi accession ICC 4958 × C reticulatum wild accession ICC 17160)- and intra (ICC 4958 × C. arietinum kabuli accession ICC 8261)-specific RIL mapping populations to identify the major QTL genomic regions governing flowering time in chickpea. The whole genome resequencing discovered 1635117 and 592486 SNPs exhibiting differentiation between early- and late-flowering mapping parents and bulks, constituted by pooling the homozygous individuals of extreme flowering time phenotypic trait from each of two aforesaid RIL populations. The multiple QTL-seq analysis using these mined SNPs in two RIL mapping populations narrowed-down two longer (907.1 kb and 1.99 Mb) major flowering time QTL genomic regions into the high-resolution shorter (757.7 kb and 1.39 Mb) QTL intervals on chickpea chromosome 4. This essentially identified regulatory as well as coding (non-synonymous/synonymous) novel SNP allelic variants from two efl1 (early flowering 1) and GI (GIGANTEA) genes regulating flowering time in chickpea. Interestingly, strong natural allelic diversity reduction (88-91%) of two known flowering genes especially mapped at major QTL intervals as compared to that of background genomic regions (where no flowering time QTLs were mapped; 61.8%) in cultivated vis-à-vis wild Cicer gene pools was evident inferring the significant impact of evolutionary bottlenecks on these loci during chickpea domestication. Higher association potential of coding non-synonymous and regulatory SNP alleles mined from efl1 (36-49%) and GI (33-42%) flowering genes for early and late flowering time differentiation among chickpea accessions was evident. The robustness and validity of two functional allelic variants-containing genes localized at major flowering time QTLs was apparent by their identification from multiple intra-/inter-specific mapping populations of chickpea. The functionally relevant molecular tags delineated can be of immense use for deciphering the natural allelic diversity-based domestication pattern of flowering time and expediting genomics-aided crop improvement to develop early flowering cultivars of chickpea.

Draft genome sequence of Cicer reticulatum L., the wild progenitor of chickpea provides a resource for agronomic trait improvement.

Cicer reticulatum L. is the wild progenitor of the fourth most important legume crop chickpea (C. arietinum L.). We assembled short-read sequences into 416 Mb draft genome of C. reticulatum and anchored 78% (327 Mb) of this assembly to eight linkage groups. Genome annotation predicted 25,680 protein-coding genes covering more than 90% of predicted gene space. The genome assembly shared a substantial synteny and conservation of gene orders with the genome of the model legume Medicago truncatula. Resistance gene homologs of wild and domesticated chickpeas showed high sequence homology and conserved synteny. Comparison of gene sequences and nucleotide diversity using 66 wild and domesticated chickpea accessions suggested that the desi type chickpea was genetically closer to the wild species than the kabuli type. Comparative analyses predicted gene flow between the wild and the cultivated species during domestication. Molecular diversity and population genetic structure determination using 15,096 genome-wide single nucleotide polymorphisms revealed an admixed domestication pattern among cultivated (desi and kabuli) and wild chickpea accessions belonging to three population groups reflecting significant influence of parentage or geographical origin for their cultivar-specific population classification. The assembly and the polymorphic sequence resources presented here would facilitate the study of chickpea domestication and targeted use of wild Cicer germplasms for agronomic trait improvement in chickpea.

A High-Resolution InDel (Insertion-Deletion) Markers-Anchored Consensus Genetic Map Identifies Major QTLs Governing Pod Number and Seed Yield in Chickpea.

Development and large-scale genotyping of user-friendly informative genome/gene-derived InDel markers in natural and mapping populations is vital for accelerating genomics-assisted breeding applications of chickpea with minimal resource expenses. The present investigation employed a high-throughput whole genome next-generation resequencing strategy in low and high pod number parental accessions and homozygous individuals constituting the bulks from each of two inter-specific mapping populations [(Pusa 1103 × ILWC 46) and (Pusa 256 × ILWC 46)] to develop non-erroneous InDel markers at a genome-wide scale. Comparing these high-quality genomic sequences, 82,360 InDel markers with reference to kabuli genome and 13,891 InDel markers exhibiting differentiation between low and high pod number parental accessions and bulks of aforementioned mapping populations were developed. These informative markers were structurally and functionally annotated in diverse coding and non-coding sequence components of genome/genes of kabuli chickpea. The functional significance of regulatory and coding (frameshift and large-effect mutations) InDel markers for establishing marker-trait linkages through association/genetic mapping was apparent. The markers detected a greater amplification (97%) and intra-specific polymorphic potential (58-87%) among a diverse panel of cultivated desi, kabuli, and wild accessions even by using a simpler cost-efficient agarose gel-based assay implicating their utility in large-scale genetic analysis especially in domesticated chickpea with narrow genetic base. Two high-density inter-specific genetic linkage maps generated using aforesaid mapping populations were integrated to construct a consensus 1479 InDel markers-anchored high-resolution (inter-marker distance: 0.66 cM) genetic map for efficient molecular mapping of major QTLs governing pod number and seed yield per plant in chickpea. Utilizing these high-density genetic maps as anchors, three major genomic regions harboring each of pod number and seed yield robust QTLs (15-28% phenotypic variation explained) were identified on chromosomes 2, 4, and 6. The integration of genetic and physical maps at these QTLs mapped on chromosomes scaled-down the long major QTL intervals into high-resolution short pod number and seed yield robust QTL physical intervals (0.89-2.94 Mb) which were essentially got validated in multiple genetic backgrounds of two chickpea mapping populations. The genome-wide InDel markers including natural allelic variants and genomic loci/genes delineated at major six especially in one colocalized novel congruent robust pod number and seed yield robust QTLs mapped on a high-density consensus genetic map were found most promising in chickpea. These functionally relevant molecular tags can drive marker-assisted genetic enhancement to develop high-yielding cultivars with increased seed/pod number and yield in chickpea.