Comparative transcriptome database for Camellia sinensis reveals genes related to the cold sensitivity and albino mechanism of 'Anji Baicha'.

Tea, a globally popular beverage, contains various beneficial secondary metabolites. Tea plants (Camellia sinensis) exhibit diverse genetic traits across cultivars, impacting yield, adaptability, morphology, and secondary metabolite composition. Many tea cultivars have been the subject of much research interest, which have led to the accumulation of publicly available RNA-seq data. As such, it has become possible to systematically summarize the characteristics of different cultivars at the transcriptomic level, identify functional genes, and infer gene functions through co-expression analysis. Here, the transcriptomes of 9 tea cultivars were assembled, and comparative analysis was conducted on the coding sequences of 13 cultivars. To give access to this data, we present TeaNekT (https://teanekt.sbs.ntu.edu.sg/), a web resource that facilitates the prediction of gene functions of various tea cultivars. We used TeaNekT to perform a cross-cultivar comparison of co-expressed gene clusters and tissue-specific gene expression. We observed that 'Anji Baicha' possesses the highest number of cultivar-specific genes and the second-highest number of expanded genes. These genes in 'Anji Baicha' tend to be enriched in functions associated with cold stress response, chloroplast thylakoid structure, and nitrogen metabolism. Notably, we identified three significantly expanded homologous genes in 'Anji Baicha' encoding the ICE1, SIZ1, and MAPKK2, which are closely associated with the cold sensitivity of 'Anji Baicha'. Additionally, one significantly expanded homologous gene in 'Anji Baicha' encoding regulatory factor RIQ may play a crucial role in the abnormal chloroplast structure and absence of thylakoid membranes in 'Anji Baicha'.

Analysis of Rice Blast Fungus Genetic Diversity and Identification of a Novel Blast Resistance OsDRq12 Gene.

The rice blast fungus Magnaporthe oryzae poses a significant challenge to maintaining rice production. Developing rice varieties with resistance to this disease is crucial for its effective control. To understand the genetic variability of blast isolates collected between 2015 and 2017, the 27 monogenic rice lines that carry specific resistance genes were used to evaluate blast disease reactions. Based on criteria such as viability, virulence, and reactions to resistance genes, 20 blast isolates were selected as representative strains. To identify novel resistance genes, a quantitative trait locus analysis was carried out utilizing a mixture of the 20 representative rice blast isolates and a rice population derived from crossing the blast-resistant cultivar 'Cheongcheong' with the blast-susceptible cultivar 'Nagdong'. This analysis revealed a significant locus, RM1227-RM1261 on chromosome 12, that is associated with rice blast resistance. Within this locus, 12 disease resistance-associated protein genes were identified. Among them, OsDRq12, a member of the nucleotide-binding, leucine-rich repeat disease resistance family, was chosen as the target gene for additional computational investigation. The findings of this study have significant implications for enhancing rice production and ensuring food security by controlling rice blast and developing resistant rice cultivars.

Phenotypic and transcriptomics characterization uncovers genes underlying tuber yield traits and gene expression marker development in potato under aeroponics.

MAIN CONCLUSION: Transcriptome analysis in potato varieties revealed genes associated with tuber yield-related traits and developed gene expression markers. This study aimed to identify genes involved in high tuber yield and its component traits in test potato varieties (Kufri Frysona, Kufri Khyati, and Kufri Mohan) compared to control (Kufri Sutlej). The aeroponic evaluation showed significant differences in yield-related traits in the varieties. Total RNA sequencing was performed using tuber and leaf tissues on the Illumina platform. The high-quality reads (QV > 25) mapping with the reference potato genomes revealed statistically significant (P < 0.05) differentially expressed genes (DEGs) into two categories: up-regulated (> 2 Log2 fold change) and down-regulated (< -2 Log2 fold change). DEGs were characterized by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. Collectively, we identified genes participating in sugar metabolism, stress response, transcription factors, phytohormones, kinase proteins, and other genes greatly affecting tuber yield and its related traits. A few selected genes were UDP-glucose glucosyltransferase, glutathion S-transferase, GDSL esterase/lipase, transcription factors (MYB, WRKY, bHLH63, and BURP), phytohormones (auxin-induced protein X10A, and GA20 oxidase), kinase proteins (Kunitz-type tuber invertase inhibitor, BRASSINOSTEROID INSENSITIVE 1-associated receptor kinase 1) and laccase. Based on the selected 17 peptide sequences representing 13 genes, a phylogeny tree and motifs were analyzed. Real time-quantitative polymerase chain reaction (RT-qPCR) analysis was used to validate the RNA-seq results. RT-qPCR based gene expression markers were developed for the genes such as 101 kDa heat shock protein, catechol oxidase B chloroplastic, cysteine protease inhibitor 1, Kunitz-type tuber invertase inhibitor, and laccase to identify high yielding potato genotypes. Thus, our study paved the path for potential genes associated with tuber yield traits in potato under aeroponics.

Defining genomic landscape for identification of potential candidate resistance genes associated with major rice diseases through MetaQTL analysis.

Rice production is severely affected by various diseases such as bacterial leaf blight (BLB), brown spot (BS), false smut (FS), foot rot (FR), rice blast (RB), and sheath blight (SB). In recent years, several quantitative trait loci (QTLs) studies involving different populations have been carried out, resulting in the identification of hundreds of resistance QTLs for each disease. These QTLs can be integrated and analyzed using meta-QTL (MQTL) analysis for better understanding of the genetic architecture underlying multiple disease resistance (MDR). This study involved an MQTL analysis on 661 QTLs (378, 161, 21, 41, 44, and 16 QTLs for SB, RB, BLB, BS, FS, and FR, respectively) retrieved from 50 individual studies published from 1995 to 2021. Of these, 503 QTLs were projected finally onto the consensus map saturated with 6,275 markers, resulting in 73 MQTLs, including 27 MDR-MQTLs conferring resistance to three or more diseases. Forty-seven MQTLs were validated using marker-trait associations identified in published genome-wide association studies. A total of 3,310 genes, including both R and defense genes, were also identified within some selected high-confidence MQTL regions that were investigated further for the syntenic relationship with barley, wheat, and maize genomes. Thirty-nine high-confidence candidate genes were selected based on their expression patterns and recommended for future studies involving functional validation, genetic engineering, and gene editing. Nineteen MQTLs were co-localized with 39 known R genes for BLB and RB diseases. These results could pave the way to utilize candidate genes in a marker-assisted breeding program for MDR in rice.

Mapping of Aegilops speltoides derived leaf rust and stripe rust resistance genes using 35K SNP array.

Wheat is an essential food commodity cultivated throughout the world. However, this crop faces continuous threats from fungal pathogens, leaf rust (LR) and stripe rust (YR). To continue feeding the growing population, these major destructors of wheat must be effectively countered by enhancing the genetic diversity of cultivated germplasm. In this study, an introgression line with hexaploid background (ILsp3603) carrying resistance against Pt pathotypes 77-5 (121R63-1), 77-9 (121R60-1) and Pst pathotypes 46S119 (46E159), 110S119 (110E159), 238S119 (238E159) was developed from donor wheat wild progenitor, Aegilops speltoides acc pau 3603. To understand the genetic basis of resistance and map these genes (named Lrsp3603 and Yrsp3603), inheritance studies were carried out in F6 and F7 mapping population, developed by crossing ILsp3603 with LR and YR susceptible cultivar WL711, which revealed a monogenic (single gene) inheritance pattern for each of these traits. Bulk segregant analysis combined with 35 K Axiom SNP array genotyping mapped both genes as separate entities on the short arm of chromosome 6B. A genetic linkage map, comprising five markers, 1 SNP, 1 PLUG and three gene based SSRs, covered a genetic distance of 12.65 cM. Lrsp3603 was flanked by markers Tag-SSR14 (located proximally at 2.42 cM) and SNP AX-94542331 (at 3.28 cM) while Yrsp3603 was mapped at one end closest to AX-94542331 at 6.62 cM distance. Functional annotation of Lrsp3603 target region (∼ 1 Mbp) revealed 10 gene IDs associated with disease resistance mechanisms including three encoding typical R gene domains.

Application of orthology and network biology to infer gene functions in non-model plants.

Approximately 60% of the genes and gene products in the model species Arabidopsis thaliana have been functionally characterized. In non-model plant species, the functional annotation of the gene space is largely based on homology, with the assumption that genes with shared common ancestry have conserved functions. However, the wide variety in possible morphological, physiological, and ecological differences between plant species gives rise to many species- and clade-specific genes, for which this transfer of knowledge is not possible. Other complications, such as difficulties with genetic transformation, the absence of large-scale mutagenesis methods, and long generation times, further lead to the slow characterization of genes in non-model species. Here, we discuss different resources that integrate plant gene function information. Different approaches that support the functional annotation of gene products, based on orthology or network biology, are described. While sequence-based tools to characterize the functional landscape in non-model species are maturing and becoming more readily available, easy-to-use network-based methods inferring plant gene functions are not as prevalent and have limited functionality.

Combined Genome-Wide Association Study and Transcriptome Analysis Reveal Candidate Genes for Resistance to Rust (Puccinia graminis) in Dactylis glomerata.

Rust disease is a common plant disease that can cause wilting, slow growth of plant leaves, and even affect the growth and development of plants. Orchardgrass (Dactylis glomerata L.) is native to temperate regions of Europe, which has been introduced as a superior forage grass in temperate regions worldwide. Orchardgrass has rich genetic diversity and is widely distributed in the world, which may contain rust resistance genes not found in other crops. Therefore, we collected a total of 333 orchardgrass accessions from different regions around the world. Through a genome-wide association study (GWAS) analysis conducted in four different environments, 91 genes that overlap or are adjacent to significant single nucleotide polymorphisms (SNPs) were identified as potential rust disease resistance genes. Combining transcriptome data from susceptible (PI292589) and resistant (PI251814) accessions, the GWAS candidate gene DG5C04160.1 encoding glutathione S-transferase (GST) was found to be important for orchardgrass rust (Puccinia graminis) resistance. Interestingly, by comparing the number of GST gene family members in seven species, it was found that orchardgrass has the most GST gene family members, containing 119 GST genes. Among them, 23 GST genes showed significant differential expression after inoculation with the rust pathogen in resistant and susceptible accessions; 82% of the genes still showed significantly increased expression 14 days after inoculation in resistant accessions, while the expression level significantly decreased in susceptible accessions. These results indicate that GST genes play an important role in orchardgrass resistance to rust (P. graminis) stress by encoding GST to reduce its oxidative stress response.

Identification and analysis of nine new flo2 allelic mutants in rice.

FLO2 is involved in grain development and storage substance synthesis in rice, and therefore can regulate grain size and quality. In this study, we identified 4 new flo2 allelic mutants with nonsense and frameshift mutation in the exon of 6, 10, 11 and 21 and 5 new flo2 allelic mutants with alternative splicing and frameshift mutation at the splicing site of intron 13, 14, 16 and 17. Compared with wild-type rice, the outer endosperm of flo2 mutants was transparent, and the inner endosperm was floury. Different mutation sites and types of FLO2 significantly decreased kernel width, thickness and weight to some extent. The contents of storage protein, starch, amylose and amylopectin showed significant decrease at different levels among 9 flo2 mutants. The expressions of most storage protein synthesis genes and starch synthesis-related genes were significantly down-regulated, and exhibited different ranges of variation among different flo2 mutants. This study could add helpful information for the roles of flo2 alleles in rice quality regulation and provide abundant germplasm resources for rice quality breeding.

Evolutionary and immune-activating character analyses of NLR genes in algae suggest the ancient origin of plant intracellular immune receptors.

Nucleotide-binding leucine-rich repeat (NLR) proteins are crucial intracellular immune receptors in plants, responsible for detecting invading pathogens and initiating defense responses. While previous studies on the evolution and function of NLR genes were mainly limited to land plants, the evolutionary trajectory and immune-activating character of NLR genes in algae remain less explored. In this study, genome-wide NLR gene analysis was conducted on 44 chlorophyte species across seven classes and seven charophyte species across five classes. A few but variable number of NLR genes, ranging from one to 20, were identified in five chlorophytes and three charophytes, whereas no NLR gene was identified from the remaining algal genomes. Compared with land plants, algal genomes possess fewer or usually no NLR genes, implying that the expansion of NLR genes in land plants can be attributed to their adaptation to the more complex terrestrial pathogen environments. Through phylogenetic analysis, domain composition analysis, and conserved motifs profiling of the NBS domain, we detected shared and lineage-specific features between NLR genes in algae and land plants, supporting the common origin and continuous evolution of green plant NLR genes. Immune-activation assays revealed that both TNL and RNL proteins from green algae can elicit hypersensitive responses in Nicotiana benthamiana, indicating the molecular basis for immune activation has emerged in the early evolutionary stage of different types of NLR proteins. In summary, the results from this study suggest that NLR proteins may have taken a role as intracellular immune receptors in the common ancestor of green plants.

Genome-wide SNPs and candidate genes underlying the genetic variations for protein and amino acids in pearl millet (Pennisetum glaucum) germplasm.

MAIN CONCLUSION: A total of 544 significant marker-trait associations and 286 candidate genes associated with total protein and 18 amino acids were identified. Thirty-three candidate genes were found near the strong marker trait associations (- log10P ≥ 5.5). Pearl millet (Pennisetum glaucum) is largely grown as a subsistence crop in South Asia and sub-Saharan Africa. It serves as a major source of daily protein intake in these regions. Despite its importance, no systematic effort has been made to study the genetic variations of protein and amino acid content in pearl millet germplasm. The present study was undertaken to dissect the global genetic variations of total protein and 18 essential and non-essential amino acids in pearl millet, using a set of 435 K Single Nucleotide Polymorphisms (SNPs) and 161 genotypes of the Pearl Millet Inbred Germplasm Association Panel (PMiGAP). A total of 544 significant marker-trait associations (at P < 0.0001; - log10P ≥ 4) were detected and 23 strong marker-trait associations were identified using Bonferroni's correction method. Forty-eight pleiotropic loci were found in the genome for the studied traits. In total, 286 candidate genes associated with total protein and 18 amino acids were identified. Thirty-three candidate genes were found near strongly associated SNPs. The associated markers and the candidate genes provide an insight into the genetic architecture of the traits studied and are going to be useful in breeding improved pearl millet varieties in the future. Availabilities of improved pearl millet varieties possessing higher protein and amino acid compositions will help combat the rising malnutrition problem via diet.

Genome-wide association studies meta-analysis uncovers NOJO and SGS3 novel genes involved in Arabidopsis thaliana primary root development and plasticity.

BACKGROUND: Arabidopsis thaliana primary root growth has become a model for evo-devo studies due to its simplicity and facility to record cell proliferation and differentiation. To identify new genetic components relevant to primary root growth, we used a Genome-Wide Association Studies (GWAS) meta-analysis approach using data published in the last decade. In this work, we performed intra and inter-studies analyses to discover new genetic components that could participate in primary root growth. METHODS AND RESULTS: We used 639 accessions from nine different studies under control conditions and performed different GWAS tests. We found that primary root growth changes were associated with 41 genes, of which six (14.6%) have been previously described as inhibitors or promoters of primary root growth. The knockdown lines of two genes, Suppressor of Gene Silencing (SGS3), involved in tasiRNA processing, and a gene with a Sterile Alpha Motif (SAM) motif named NOJOCH MOOTS (NOJO), confirmed their role as repressors of primary root growth, none has been shown to participate in this developmental process before. CONCLUSIONS: In summary, our GWAS analysis of different available studies identified new genes that participate in primary root growth; two of them were identified as repressors of primary root growth.

Leaf morphology related genes revealed by integrating Pan-transcriptome, GWAS and eQTL analyses in a Liriodendron population.

Leaf plays an indispensable role in plant development and growth. Although many known genes related to leaf morphology development have been identified, elucidating the complex genetic basis of leaf morphological traits remains a challenge. Liriodendron plants are common ornamental trees due to their unique leaf shapes, while the molecular mechanism underlying Liriodendron leaf morphogenesis has remained unknown. Herein, we firstly constructed a population-level pan-transcriptome of Liriodendron from 81 accessions to explore the expression presence or absence variations (ePAVs), global expression differences at the population level, as well as differentially expressed genes (DEGs) between the Liriodendron chinense and Liriodendron tulipifera accessions. Subsequently, we integrated a genome-wide association study (GWAS), expression quantitative trait loci (eQTL), and transcriptome-wide association study (TWAS) to identify candidate genes related to leaf morphology. Through GWAS analysis, we identified 18 and 17 significant allelic loci in the leaf size and leaf shape modules, respectively. In addition, we discerned 16 candidate genes in relation to leaf morphological traits via TWAS. Further, integrating the co-localization results of GWAS and eQTL, we determined two regulatory hotspot regions, hot88 and hot758, related to leaf size and leaf shape, respectively. Finally, co-expression analysis, eQTL, and linkage mapping together demonstrated that Lchi_4g10795 regulate their own expression levels through cis-eQTL to affect the expression of downstream genes and cooperatively participate in the development of Liriodendron leaf morphology. These findings will improve our understanding of the molecular regulatory mechanism of Liriodendron leaf morphogenesis and will also accelerate molecular breeding of Liriodendron.

PABLOG: a Primer Analysis tool using a Bee-Like approach on Orthologous Genes.

RNA-seq data is currently generated in numerous non-model organisms that lack a reference genome. Nevertheless, the confirmation of gene expression levels using RT-qPCR remains necessary, and the existing techniques do not seamlessly interface with the omics pipeline workflow. Developing primers for many targets by utilising orthologous genes can be a laborious, imprecise, and subjective process, particularly for plant species that are not commonly studied and do not have a known genome. We have developed a primer design tool, named PABLOG, that analyses the alignments generated from long or short RNA-seq reads and a reference orthologous gene. PABLOG scans, much like a bee searching several flowers for pollen, and presents a sorted list of potential exon-exon junction locations, ranked according to their reliability. Through computational analysis across the whole genomes of several non-model species, we demonstrate that PABLOG performs more effectively than other methods in identifying exon-exon junctions since it generates significantly fewer false-positive results. Examination of candidate regions at the gene level, in conjunction with laboratory studies, shows that the suggested primers successfully amplified particular targets in non-model plants without any presence of genomic contamination. Our tool includes a consensus sequence feature that enables the complete process of primer design, from aligning with the target gene to determining amplification parameters. The utility can be accessed via the GitHub repository located at: https://github.com/tools4plant-omics/PABLOG.

Validation and identification of promising gene specific markers governing foliar disease resistance in groundnut (Arachis hypogaea L.).

BACKGROUND: Groundnut is vulnerable to the major foliar fungal disease viz., late leaf spot (LLS) and rust in kharif season, which results in severe yield losses. Until now, LLS and rust resistance linked markers were developed based on GPBD 4 as a major donor source and were validated in its derivatives only, which restricted their use in marker assisted selection (MAS) involving other donors. METHODS AND RESULTS: The current study focused to validate LLS and rust resistance linked markers employing advanced breeding lines of F6 generation, derived from nine different crosses involving nine diverse parents, to identify potential markers for marker-assisted breeding of LLS and rust resistance in groundnut. Out of 28-trait linked markers used for validation, 8 were polymorphic (28.57%). Marker-trait association (MTA) and Single Marker Analysis (SMA) revealed that the SSR marker pPGPseq5D05 is significantly associated with both LLS (15.8% PVE) and rust (17.5% PVE) resistance, whereas, the marker IPAHM103 is tightly linked with rust resistance (26.8% PVE) alone. In silico analysis revealed that the marker gene for IPAHM103 is a zinc finger protein and the marker gene for pPGPseq5D05 is an ADP-ribosylation factor GTPase-activating protein. Both these protein products impart resistance or tolerance to biotic stress in crop plants. Two other markers namely, GMLQ975 and pPGPseq13A10 were also found to be associated with LLS resistance explaining MTA up to 60%. CONCLUSION: These gene specific markers will enable us to screen more number of germplasm lines or newly developed lines in MAS schemes for LLS and rust resistance using a wide range of resistant sources.

Adventitious embryonic causal gene FhRWP regulates multiple developmental phenotypes in citrus reproduction.

Citrus is a model plant for studying adventitious embryos, a form of asexual reproduction controlled by a single dominant gene, RWP. This gene has been identified as the causal gene for nucellar embryogenesis, but its function has not yet been fully understood. In this study, we used the fast-growing Fortunella hindsii as a system to explore chromatin accessibility during the nucellar embryony initiation, emphasizing elevated chromatin accessibility in polyembryonic (PO) genotypes compared to monoembryonic ones (MO). Notably, a higher level of accessible chromatin was observed in one allele of the promoter region of FhRWP, consistent with increased expression of the allele carrying the causal structural variant. By independently performing RNAi and gene editing experiments on PO genotypes, we found the downregulation of FhRWP expression could reduce the number of nucellar embryos, while its knockout resulted in abnormal axillary bud development. In overexpression experiments, FhRWP was identified as having the unique capability of inducing the embryogenic callus formation in MO stem segments, possibly through the regulation of the WUS-CLV signaling network and the ABA and cytokinin pathway, marking the inaugural demonstration of FhRWP's potential to reignite somatic cells' embryogenic fate. This study reveals the pleiotropic function of RWP in citrus and constructs a regulatory network during adventitious embryo formation, providing a new tool for bioengineering applications in plant regeneration.

Research on the Effects of the Relationship between Agronomic Traits and Dwarfing Genes on Yield in Colored Wheat.

This research focuses on 72 approved varieties of colored wheat from different provinces in China. Utilizing coefficients of variation, structural equation models, and correlation analyses, six agronomic traits of colored wheat were comprehensively evaluated, followed by further research on different dwarfing genes in colored wheat. Using the entropy method revealed that among the 72 colored wheat varieties, 10 were suitable for cultivation. Variety 70 was the top-performing variety, with a comprehensive index of 87.15%. In the final established structural equation model, each agronomic trait exhibited a positive direct effect on yield. Notably, plant height, spike length, and flag leaf width had significant impacts on yield, with path coefficients of 0.55, 0.40, and 0.27. Transcriptome analysis and real-time fluorescence quantitative polymerase chain reaction (RT-qPCR) validation were used to identify three dwarfing genes controlling plant height: Rht1, Rht-D1, and Rht8. Subsequent RT-qPCR validation clustering heatmap results indicated that Rht-D1 gene expression increased with the growth of per-acre yield. Rht8 belongs to the semi-dwarf gene category and has a significant positive effect on grain yield. However, the impact of Rht1, as a dwarfing gene, on agronomic traits varies. These research findings provide crucial references for the breeding of new varieties.

Deficiency of multiple RNA silencing-associated genes may contribute to the increased susceptibility of Nicotiana benthamiana to viruses.

KEY MESSAGE: Recently published high-quality reference genome assemblies indicate that, in addition to RDR1-deficiency, the loss of several key RNA silencing-associated genes may contribute to the hypersusceptibility of Nicotiana benthamiana to viruses.

Evolution, diversity, and function of the disease susceptibility gene Snn1 in wheat.

Septoria nodorum blotch (SNB), caused by Parastagonospora nodorum, is a disease of durum and common wheat initiated by the recognition of pathogen-produced necrotrophic effectors (NEs) by specific wheat genes. The wheat gene Snn1 was previously cloned, and it encodes a wall-associated kinase that directly interacts with the NE SnTox1 leading to programmed cell death and ultimately the development of SNB. Here, sequence analysis of Snn1 from 114 accessions including diploid, tetraploid, and hexaploid wheat species revealed that some wheat lines possess two copies of Snn1 (designated Snn1-B1 and Snn1-B2) approximately 120 kb apart. Snn1-B2 evolved relatively recently as a paralog of Snn1-B1, and both genes have undergone diversifying selection. Three point mutations associated with the formation of the first SnTox1-sensitive Snn1-B1 allele from a primitive wild wheat were identified. Four subsequent and independent SNPs, three in Snn1-B1 and one in Snn1-B2, converted the sensitive alleles to insensitive forms. Protein modeling indicated these four mutations could abolish Snn1-SnTox1 compatibility either through destabilization of the Snn1 protein or direct disruption of the protein-protein interaction. A high-throughput marker was developed for the absent allele of Snn1, and it was 100% accurate at predicting SnTox1-insensitive lines in both durum and spring wheat. Results of this study increase our understanding of the evolution, diversity, and function of Snn1-B1 and Snn1-B2 genes and will be useful for marker-assisted elimination of these genes for better host resistance.

Identification of anther thermotolerance genes by the integration of linkage and association analysis in maize.

Maize is one of the world's most important staple crops, yet its production is increasingly threatened by the rising frequency of high-temperature stress (HTS). To investigate the genetic basis of anther thermotolerance under field conditions, we performed linkage and association analysis to identify HTS response quantitative trait loci (QTL) using three recombinant inbred line (RIL) populations and an association panel containing 375 diverse maize inbred lines. These analyses resulted in the identification of 16 co-located large QTL intervals. Among the 37 candidate genes identified in these QTL intervals, five have rice or Arabidopsis homologs known to influence pollen and filament development. Notably, one of the candidate genes, ZmDUP707, has been subject to selection pressure during breeding. Its expression is suppressed by HTS, leading to pollen abortion and barren seeds. We also identified several additional candidate genes potentially underly QTL previously reported by other researchers. Taken together, our results provide a pool of valuable candidate genes that could be employed by future breeding programs aiming at enhancing maize HTS tolerance.

Population-level gene expression can repeatedly link genes to functions in maize.

Transcriptome-wide association studies (TWAS) can provide single gene resolution for candidate genes in plants, complementing genome-wide association studies (GWAS) but efforts in plants have been met with, at best, mixed success. We generated expression data from 693 maize genotypes, measured in a common field experiment, sampled over a 2-h period to minimize diurnal and environmental effects, using full-length RNA-seq to maximize the accurate estimation of transcript abundance. TWAS could identify roughly 10 times as many genes likely to play a role in flowering time regulation as GWAS conducted data from the same experiment. TWAS using mature leaf tissue identified known true-positive flowering time genes known to act in the shoot apical meristem, and trait data from a new environment enabled the identification of additional flowering time genes without the need for new expression data. eQTL analysis of TWAS-tagged genes identified at least one additional known maize flowering time gene through trans-eQTL interactions. Collectively these results suggest the gene expression resource described here can link genes to functions across different plant phenotypes expressed in a range of tissues and scored in different experiments.