Stepping towards pollen DNA metabarcoding: A breakthrough in forensic sciences.

This review is engaged in determining the capability of plant pollen as a significant source of evidence for the linkage between suspects and crime location in forensic sciences. Research and review articles were collected from Google Scholar, the Web of Science, and PubMed. Articles were searched using specific keywords such as "Forensic Palynology," "Pollen metabarcoding," "Plant forensics," and "Pollen" AND "criminal investigation." Boolean logic was also utilized to narrow the articles to be included in this review article. Through the literature and exploratory research, it has been observed in the current study that with advancements in technology, forensic palynology has found its application in creating an association between the crime scene and suspected individuals to have a link to it, as pollen DNA is a long-lasting investigative tool that can effectively help forensic investigations. Moreover, the literature shows that the DNA of pollen and spores has helped forensic scientists link suspects to crime scenes, and the introduction of pollen DNA metabarcoding tools has eased the efforts of palynologists to analyze pollen DNA. The introduction of DNA metabarcoding techniques to analyze pollen from plants has helped identify the geological locations of the plants and ultimately identify the culprit.

Antimicrobial resistome and mobilome in the urban river affected by combined sewer overflows and wastewater treatment effluent.

The dissemination of antimicrobial resistance in the environment is an emerging global health problem. Wastewater treatment effluent and combined sewer overflows (CSOs) are major sources of antimicrobial resistance in urban rivers. This study aimed to clarify the effect of municipal wastewater treatment effluent and CSO on antimicrobial resistance genes (ARGs), mobile gene elements, and the microbial community in an urban river. The ARG abundance per 16S-based microbial population in the target river was 0.37-0.54 and 0.030-0.097 during the CSO event and dry weather, respectively. During the CSO event, the antimicrobial resistome in the river shifted toward a higher abundance of ARGs to clinically important drug classes, including macrolide, fluoroquinolone, and ß-lactam, whereas ARGs to sulfonamide and multidrug by efflux pump were relatively abundant in dry weather. The abundance of intI1 and tnpA genes were highly associated with the total ARG abundance, suggesting their potential application as an indicator for estimating resistome contamination. Increase of prophage during the CSO event suggested that impact of CSO has a greater potential for horizontal gene transfer (HGT) via transduction. Consequently, CSO not only increases the abundance of ARGs to clinically important antimicrobials but also possibly enhances potential of HGT in urban rivers.

Genetic Risk Score Constructed from Polymorphisms in the VEGFA, TBX5, and SMAD7 Genes Provides Novel Insights into the Molecular Mechanisms of the Tetralogy of Fallot and Ventricular Septal Defect (Case-Control Study from the Pakistani Population).

Congenital heart defects are common and complex birth-defect malformations in developed and developing countries. It is a multifactorial disease that involves the interaction of either gene-gene or gene-environment. This comparative study was the first report on the genotypic-phenotypic correlation in the Pakistani population. The single nucleotide polymorphisms (SNPs) were further tested for association with maternal diabetes mellitus or hypertension. In addition, the cumulative genetic risk score (GRS) for low to moderately-associated SNPs was calculated for each study subject, which can ultimately guide us for better therapeutic options and prevention strategies. According to the predefined selection criteria, 376 subjects were recruited. The multiplex mini-sequencing genotyping technique opted for the cost-effective genotyping of selected loci. The association of variants with the disease was examined using logistic regression analysis. The statistical and graphical analysis was conducted using SPSS, Haploview, SNPStats, and GraphPad Prism. The results for all SNPs analysis suggested a nonsignificant association with overall congenital heart defect risk except rs3809923. However, interestingly on stratified analysis variants, rs3809923 and rs3809922 showed an association only with tetralogy of Fallot. The remaining risk factor analysis for maternal hypertension and diabetes mellitus association with SNPs were nonsignificant. The GRS was the first time constructed for this low to moderately-associated variants. Interestingly, the cumulative GRS was significantly different from the control group revealing the cumulative effect of these polymorphisms panel in patients. In conclusion, the use of GRS in the clinical setting can predict better risk association and patient outcomes.

Impact of prostate radiotherapy on survival outcomes in clinically node-positive prostate cancer: A multicentre retrospective analysis.

PURPOSE: To evaluate clinical outcomes for cN1M0 prostate cancer treated with varied modalities. MATERIALS AND METHODS: Men with radiological stage cN1M0 prostate cancer on conventional imaging, treated from 2011-2019 with various modalities across four centres in the UK were included. Demographics, tumour grade and stage, and treatment details were collected. Biochemical and radiological progression-free survival (bPFS, rPFS) and overall survival (OS) were estimated using Kaplan Meier analyses. Potential factors impacting survival were tested with univariable log-rank test and multivariable Cox-proportional hazards model. RESULTS: Total 337 men with cN1M0 prostate cancer were included, 47% having Gleason grade group 5 disease. Treatment modalities included androgen deprivation therapy (ADT) in 98.9% men, either alone (19%) or in combinations including prostate radiotherapy (70%), pelvic nodal radiotherapy (38%), docetaxel (22%), or surgery (7%). At median follow up of 50 months, 5-year bPFS, rPFS, and OS were 62.7%, 71.0%, and 75.8% respectively. Prostate radiotherapy was associated with significantly higher bPFS (74.1% vs 34.2%), rPFS (80.7% vs 44.3%) and OS (86.7% vs 56.2%) at five years (log rank p < 0.001 each). On multivariable analysis including age, Gleason grade group, tumour stage, ADT duration, docetaxel, and nodal radiotherapy, benefit of prostate radiotherapy persisted for bPFS [HR 0.33 (95% CI 0.18-0.62)], rPFS [HR 0.25 (0.12-0.51)], and OS [HR 0.27 (0.13-0.58)] (p < 0.001 each). Impact of nodal radiotherapy or docetaxel was not established due to small subgroups. CONCLUSION: Addition of prostate radiotherapy to ADT in cN1M0 prostate cancer yielded improved disease control and overall survival independent of other tumour and treatment factors.

Corrigendum: Contaminants from a former Croatian coal sludge dictate the structure of microbiota in the estuarine (Rasa Bay) sediment and soil.

[This corrects the article DOI: 10.3389/fmicb.2023.1126612.].

Association of maternal hypertension and diabetes with variants of the NKX2-5, LEFTY1 and LEFTY2 genes in children with congenital heart defects: a case-control study from Pakistani Population.

BACKGROUND: Globally, congenital heart defect (CHD) is the most common congenital malformation, responsible for higher morbidity and mortality in the pediatric population. It is a complex multifactorial disease influenced by gene-environment and gene-gene interactions. The current study was the first attempt to study these polymorphisms in common clinical phenotypes of CHD in Pakistan and the association between maternal hypertension and diabetes with single nucleotide polymorphisms (SNPs) in children. METHODS: A total of 376 subjects were recruited in this current case-control study. Six variants from three genes were analyzed by cost-effective multiplex PCR and genotyped by minisequencing. Statistical analysis was done by GraphPad prism and Haploview. The association of SNPs and CHD was determined using logistic regression. RESULTS: The risk allele frequency was higher in cases as compared to healthy subjects, but the results were not significant for rs703752. However, stratification analysis suggested that rs703752 was significantly associated with the tetralogy of Fallot. The rs2295418 was significantly associated with maternal hypertension (OR = 16.41, p = 0.003), while a weak association was present between maternal diabetes and rs360057 (p = 0.08). CONCLUSION: In conclusion, variants in transcriptional and signaling genes were associated with Pakistani pediatric CHD patients that showed varied susceptibility between different clinical phenotypes of CHD. In addition, this study was the first report regarding the significant association between maternal hypertension and the LEFTY2 gene variant.

Contaminants from a former Croatian coal sludge dictate the structure of microbiota in the estuarine (Rasa Bay) sediment and soil.

Introduction: Croatian superhigh-organic-sulfur Rasa coal had been mined for nearly 400 years. The release of hazardous trace elements (HTEs) and toxic organic pollutants (TOPs) into the local environment by coal mining, preparation, and combustion activities has resulted in pollution. Methods: In this study, the diversity and composition of microbial communities in estuarine sediment and soil samples as well as community function responses to the pollutants were investigated. Results: The results showed that PAH degradation does occur following 60 years of natural attenuation, the location is still heavily polluted by polycyclic aromatic hydrocarbons (PAHs) and HTEs. Microbial analyses have shown that high concentrations of PAHs have reduced the diversity and abundance of microbial communities. The pollution exerted an adverse, long-term impact on the microbial community structure and function in the brackish aquatic ecosystem. Microorganisms associated with the degradation of PAHs and sulfur-containing compounds have been enriched although the diversity and abundance of the microbial community have reduced. Fungi which are believed to be the main PAH degrader may play an important role initially, but the activity remains lower thereafter. It is the high concentrations of coal-derived PAHs, rather than HTEs, that have reduced the diversity and abundance of microbial communities and shaped the structure of the local microbiota. Discussion: This study could provide a basis for the monitoring and restoration of ecosystems impacted by coal mining activities considering the expected decommission of a large number of coal plants on a global scale in the coming years due to growing global climate change concerns.

Role of Vascular Endothelial Growth Factor as a Potential Biomarker in Congenital Heart Defects: A Systematic Review.

Background: Several studies have investigated the role of vascular endothelial growth factor (VEGF) variants, serum levels, and correlations with other extrinsic factors in congenital heart defects (CHDs); however, the findings need confirmation. The present systematic review evaluates the association between CHDs and genetic polymorphisms and serum expressions. Methods: Relevant literature was searched through electronic databases using keywords and MeSH terms. VEGF activity was comparatively assessed between cyanotic and acyanotic CHDs, and the association between different polymorphisms and heart defects was evaluated. Results: We ultimately evaluated 12 studies regarding the association between VEGF serum patterns and found that serum VEGF levels were upregulated or downregulated in correlation with hypoxia and hemoglobin levels and were significantly associated with cyanotic CHDs compared with acyanotic CHDs. Our results also showed a significant role for different single-nucleotide polymorphisms, including rs699947, rs2010963, and rs3025039. Conclusion: The findings of the current study suggested a significant association between CHDs and VEGF genetic polymorphisms or varied serum levels. Nevertheless, more comprehensive studies may provide conclusive results and valuable insights into the pathogenesis of CHDs and relevant treatment strategies.

Single Nucleotide Polymorphisms Of Transforming Growth Factor-Β1 Gene As Potential Asthma Susceptible Variants In Punjabi Population Of Pakistan.

Background: Candidate gene approach based on case-control model is a valuable strategy to determined disease related genetic variants. Two single nucleotide polymorphisms rs1800469 and rs2241715 in TGF ß1gene have been reported to affect the asthmatic status in different populations. The main focus of this research was to find any relationship between these SNPs and asthma in Pakistani population. Methods: Using case-control model, a total of 108 individuals including 52 asthma patients and 56 healthy controls were screened to find asthma susceptibility of variants rs1800469 and rs2241715. These SNPs were genotyped using SNaPshot minisequencing assay followed by capillary electrophoresis using ABI 3130xl genetic analyzer platform. The statistical analysis of genetic data was performed by using SPSS 21, SHEsis online platform and SNPStats online web software. Results: No association with asthma was seen in allelic model for both SNPs but genotypes analyzed under codominant, dominant, over dominant and recessive models of inheritance revealed that SNP rs2241715 is strongly associated with asthma under genotypic model. Conclusion: rs2241715 was found to be a genetic risk factor for asthma in Pakistani population.

Identification of QTLs for rice grain size and weight by high-throughput SNP markers in the IR64 x Sadri population.

Grain appearance is one of the most important attributes of rice. It is determined by grain size, shape, and weight, which in turn influences the rice yield and market value. In this study, QTLs for grain length, grain width, grain length/width ratio, and grain weight were mapped using the high-throughput indica/indica SNP platforms. The population of the mega indica variety IR64 and the high-quality aromatic variety Sadri from Iran was phenotyped. Based on this phenotypic data, plants of 94 F2:3 families including both parents were selected. A linkage map analysis of 210 SNP markers identified 14 QTLs controlling the grain length, grain width, length/width ratio, and 1,000 grain weight. Among these 14, one important region containing the QTLs for all the four studies' traits was mapped on chromosome 8. It was derived from Sadri for the decreased length/width ratio and increased grain weight. This study demonstrated the speed and efficiency in using multiplex SNP genotyping for QTL analysis. Moreover, this study identified four novel QTLs (qGL8, qTGW8, qLWR8, and qGW8) sharing the same position on chromosome 8 which were linked with grain quality characteristics between one indica and one aromatic variety. It will enable more precise marker-assisted selection for grain weight, shape, and size. Further in-depth studies are required to dissect this region of interest and identify the related gene(s).

CrAssphage as an indicator of human-fecal contamination in water environment and virus reduction in wastewater treatment.

Viral indicators of human-fecal contamination in wastewaters and environmental waters have been getting much attention in the past decade. Cross-assembly phage (crAssphage) is the most abundant DNA virus in human feces. Recently, the usefulness of crAssphage as a microbial source tracking and water quality monitoring tool for human-fecal contamination has been highlighted. Here, we conducted a comprehensive review on crAssphage in water, focusing on detection methodology, concentration range in various waters and wastewaters, specificity to human-fecal contamination, and reduction in wastewater treatment systems. This review highlights that crAssphage is globally distributed in wastewaters and various fecal-contaminated water bodies at high concentrations without seasonal fluctuations. CrAssphage is highly specific to human-fecal contamination and is rarely found in animal feces. It also has a good potential as a performance indicator to ensure virus reduction in wastewater treatment systems. Accordingly, crAssphage could be an effective tool for monitoring of human-fecal contamination and potential presence of fecal pathogenic microbes in environmental waters. Bridging the research gaps highlighted in this review would make crAssphage a powerful tool to support the control of water-related health risks.

Integrated Genomic Analysis Identifies ANKRD36 Gene as a Novel and Common Biomarker of Disease Progression in Chronic Myeloid Leukemia.

Background: Chronic myeloid leukemia (CML) is initiated in bone marrow due to chromosomal translocation t(9;22) leading to fusion oncogene BCR-ABL. Targeting BCR-ABL by tyrosine kinase inhibitors (TKIs) has changed fatal CML into an almost curable disease. Despite that, TKIs lose their effectiveness due to disease progression. Unfortunately, the mechanism of CML progression is poorly understood and common biomarkers for CML progression are unavailable. This study was conducted to find novel biomarkers of CML progression by employing whole-exome sequencing (WES). Materials and Methods: WES of accelerated phase (AP) and blast crisis (BC) CML patients was carried out, with chronic-phase CML (CP-CML) patients as control. After DNA library preparation and exome enrichment, clustering and sequencing were carried out using Illumina platforms. Statistical analysis was carried out using SAS/STAT software version 9.4, and R package was employed to find mutations shared exclusively by all AP-/BC-CML patients. Confirmation of mutations was carried out using Sanger sequencing and protein structure modeling using I-TASSER followed by mutant generation and visualization using PyMOL. Results: Three novel genes (ANKRD36, ANKRD36B and PRSS3) were mutated exclusively in all AP-/BC-CML patients. Only ANKRD36 gene mutations (c.1183_1184 delGC and c.1187_1185 dupTT) were confirmed by Sanger sequencing. Protein modeling studies showed that mutations induce structural changes in ANKRD36 protein. Conclusions: Our studies show that ANKRD36 is a potential common biomarker and drug target of early CML progression. ANKRD36 is yet uncharacterized in humans. It has the highest expression in bone marrow, specifically myeloid cells. We recommend carrying out further studies to explore the role of ANKRD36 in the biology and progression of CML.

The role of NKX2-5 gene polymorphisms in congenital heart disease (CHD): a systematic review and meta-analysis.

BACKGROUND: The gene NKX2-5 is a key transcription factor that plays an essential role in normal cardiac development. Although some recent studies have studied the role of polymorphisms in the NKX2-5 gene in congenital heart diseases (CHDs), the results were not consistent and remained uncertain. Therefore, we conduct a review of literature and investigate the association of genetic polymorphisms with CHDs. RESULTS: We selected seventeen studies regarding the association of NKX2-5 gene rs2277923 polymorphism with CHDs. Overall, in all the tested genetic models, the 63A > G polymorphism was not significantly associated with increased congenital heart defects risk. We used pooled odds ratios (OR) to calculate the association of CHDs with rs2277923 including allelic model: OR 1.00, 95% CI 0.82-1.21; homozygote model: OR 0.95, 95%CI 0.68-1.33, recessive model: OR 0.89 CI 0.70-1.13, heterozygote model: OR: 1.09, 95%CI 0.87-1.37, dominant model: OR 1.08 CI 0.82-1.42 and overdominant model: OR 1.17 CI 1.01-1.35. In addition, our analysis suggests that no publication bias exists in this meta-analysis. CONCLUSIONS: Our findings suggested that 63A > G polymorphism in the NKX2-5 gene was not significantly associated with congenital heart defects. However, in the future, more studies with increased sample size are required that may provide us more definite conclusions.

Familial DNA analysis and criminal investigation: Usage, downsides and privacy concerns.

Since the discovery of Deoxyribonucleic acid (DNA) capability in forensic investigation, it has been an important part of the criminal justice system. In most criminal cases DNA profile originating from evidence sample collected from the crime scene is compared with the DNA profile from the reference sample. However, when a reference sample is not available for comparison, familial DNA analysis can provide important investigation leads in a criminal investigation process by identifying an individual. Moreover, this analysis is also proving effective in the identification of ethnicity and ancestry of an individual. A number of different methodologies and software are being used for familial DNA analysis. This review describes the importance of familial DNA analysis, methodologies used for familial DNA searching and identification, and its advantages in forensic. Moreover, ethical, legal and social issues associated with familial DNA analysis have also been discussed along with future directions for the proper implementation of this technology.

Prevention programmes and prenatal diagnosis for beta thalassemia in Pakistan: A narrative review.

Beta thalassemia in Pakistan is a serious health concern with an estimated 5-8% carrier frequency and birth of 5000 major children every year in the country. The treatment of beta thalassemia major patients poses a great economic burden; hence, the ideal approach towards this disease should encompass effective prevention services. At present only one government funded project "Punjab Thalassemia Prevention Programme" existed in Punjab province, and providing free of cost services for beta thalassemia screening and prenatal diagnosis. Complete blood count and haemoglobin electrophoresis remains the preliminary test for screening, while chorionic villi sampling and amplification refractory mutation system method have been most widely used for molecular diagnosis of beta thalassemia. Modern molecular techniques, non-invasive prenatal diagnosis, and pre-implantation diagnosis are in trial phases. In this review we have discussed the available diagnostic facilities and status of prevention programmes for beta thalassemia in Pakistan as well as future perspectives.

Major histocompatibility complex class II polymorphic variants are associated with asthma predisposition in the Punjabi population of Lahore, Pakistan.

INTRODUCTION: Various genome wide association studies have manifested that Major Histocompatibility Complex (MHC) region on chromosome 6p21 houses many potential candidate genes for asthma. OBJECTIVE: This Case-Control association study was planned to determine the association of 10 Single Nucleotide Polymorphisms (SNPs), residing within and around MHC genes' region on chromosome 6p21, with Asthma in Punjabi population of Lahore, Pakistan. METHODS: A total of 161 subjects, 61 physician-diagnosed asthma patients and 100 age-matched healthy controls, were recruited from Lahore, a city in Punjab. Ten single nucleotide polymorphisms (SNPs) (rs9378249, rs2070600, rs404860, rs6689, rs1049124, rs1063355, rs1049225, rs1049219, rs7773955 and rs928976) located within or near AGER, NOTCH and HLA genes in MHC region, were genotyped in both patients and controls using single base extension reaction and capillary electrophoresis-based genetic analyser. Statistical models were applied using SHEsis Plus. Results were adjusted for various cofactors (age, gender and environment) and by applying multiple corrections. Haplotype and linkage disequilibrium analyses were performed on Haploview software v4.1. RESULTS: Three of the studied SNPs rs1049124, rs1049219 and rs7773955 show independent significant association with asthma under allelic and genotypic models. Two of the haplotypes, H7 and H13, "CTAATTT" and "CCACTAT", respectively, for rs2070600, rs404860, rs6689, rs1049124, rs1063355, rs1049219 and rs7773955, are found to be significantly associated with the disease. CONCLUSION: This study reports association of SNP variants residing on HLA-DQB1 and HLA-DQA2 genes and haplotypes H7 and H13 on genomic region 6p21 with Asthma in the Punjabi population of Lahore, Pakistan.

Discovery and Protein Modeling Studies of Novel Compound Mutations Causing Resistance to Multiple Tyrosine Kinase Inhibitors in Chronic Myeloid Leukemia.

OBJECTIVE: BCR-ABL fusion oncogene is the hallmark of chronic myeloid leukemia (CML), causing genomic instability which leads to accumulation of mutations in BCR-ABL as well as other genes. BCR-ABL mutations are the cause of tyrosine kinase inhibitors (TKIs) resistance in CML. Recently, compound BCR-ABL mutations have been reported to resist all FDA approved TKIs. Therefore, finding novel compound BCR-ABL mutations can help and clinically manage CML. Therefore, our objective was to find out novel drug-resistant compound BCR-ABL mutations in CML and carry out their protein modelling studies. METHODOLOGY: Peripheral blood samples were collected from ten imatinib resistant CML patients receiving nilotinib treatment. BCR-ABL transcript mutations were investigated by employing capillary sequencing. Patient follow-up was carried out using European LeukemiaNet guidelines. Protein modeling  studies were carried out for new compound mutations using PyMol to see the effects of mutations at structural level. RESULTS: A novel compound mutation (K245N mutation along with G250W mutation) and previously known T351I utation was detected in two of the nilotinib resistance CML patients respectively while in the rest of 8 nilotinib responders, no resistant mutations were detected. Protein modelling studies indicated changes in BCR-ABL mutant protein which may have negatively impacted its binding with nilotinib leading to drug resistance. CONCLUSION: We report a novel nilotinib resistant BCR-ABL compound mutation (K245N along with G250W mutation) which impacts structural modification in BCR-ABL mutant protein leading to drug resistance. As compound mutations pose a new threat by causing resistance to all FDA approved tyrosine kinase inhibitors in BCR-ABL+ leukemias, our study opens a new direction for in vitro characterization of novel BCR-ABL compound mutations and their resistant to second  generation and third generation TKIs.

Grain Fe and Zn contents linked SSR markers based genetic diversity in rice.

Rice is critical for sustainable food and nutritional security; however, nominal micronutrient quantities in grains aggravate malnutrition in rice-eating poor populations. In this study, we evaluated genetic diversity in grain iron (Fe) and zinc (Zn) contents using trait-linked simple sequence repeat (SSR) markers in the representative subset of a large collection of local and exotic rice germplasm. Results demonstrated that aromatic fine grain accessions contained relatively higher Fe and Zn contents in brown rice (BR) than coarse grain accessions and a strong positive correlation between both mineral elements. Genotyping with 24 trait-linked SSR markers identified 21 polymorphic markers, among which 17 demonstrated higher gene diversity and polymorphism information content (PIC) values, strongly indicating that markers used in current research were moderate to highly informative for evaluating the genetic diversity. Population structure, principal coordinate and phylogenetic analyses classified studied rice accessions into two fine grain specific and one fine and coarse grain admixture subpopulations. Single marker analysis recognized four ZnBR and single FeBR significant marker-trait associations (MTAs) contributing 15.41-39.72% in total observed phenotypic variance. Furthermore, high grain Fe and Zn contents linked marker alleles from significant MTAs were also identified. Collectively, these results indicate a wide genetic diversity exist in grain Fe and Zn contents of studied rice accessions and reveal perspective for marker-assisted biofortification breeding.

Reproductive tissues-specific meta-QTLs and candidate genes for development of heat-tolerant rice cultivars.

KEY MESSAGE: By integrating genetics and genomics data, reproductive tissues-specific and heat stress responsive 35 meta-QTLs and 45 candidate genes were identified, which could be exploited through marker-assisted breeding for fast-track development of heat-tolerant rice cultivars. Rice holds the key to future food security. In rice-growing areas, temperature has already reached an optimum level for growth, hence, any further increase due to global climate change could significantly reduce rice yield. Several mapping studies have identified a plethora of reproductive tissue-specific and heat stress associated inconsistent quantitative trait loci (QTL), which could be exploited for improvement of heat tolerance. In this study, we performed a meta-analysis on previously reported QTLs and identified 35 most consistent meta-QTLs (MQTLs) across diverse genetic backgrounds and environments. Genetic and physical intervals of nearly 66% MQTLs were narrower than 5 cM and 2 Mb respectively, indicating hotspot genomic regions for heat tolerance. Comparative analyses of MQTLs underlying genes with microarray and RNA-seq based transcriptomic data sets revealed a core set of 45 heat-responsive genes, among which 24 were reproductive tissue-specific and have not been studied in detail before. Remarkably, all these genes corresponded to various stress associated functions, ranging from abiotic stress sensing to regulating plant stress responses, and included heat-shock genes (OsBiP2, OsMed37_1), transcription factors (OsNAS3, OsTEF1, OsWRKY10, OsWRKY21), transmembrane transporters (OsAAP7A, OsAMT2;1), sugar metabolizing (OsSUS4, α-Gal III) and abiotic stress (OsRCI2-7, SRWD1) genes. Functional data evidences from Arabidopsis heat-shock genes also suggest that OsBIP2 may be associated with thermotolerance of pollen tubes under heat stress conditions. Furthermore, promoters of identified genes were enriched with heat, dehydration, pollen and sugar responsive cis-acting regulatory elements, proposing a common regulatory mechanism might exist in rice for mitigating reproductive stage heat stress. These findings strongly support our results and provide new candidate genes for fast-track development of heat-tolerant rice cultivars.