Substantial contribution of genetic variation in the expression of transcription factors to phenotypic variation revealed by eRD-GWAS

Biotechnological studies towards improvement of finger millet using multi-omics approaches

A plethora of studies have uncovered numerous important genes with agricultural significance in staple crops. However, when it comes to orphan crops like minor millet, genomic research lags significantly behind that of major crops. This situation has promoted a focus on exploring research opportunities in minor millets, particularly in finger millet, using cutting-edge methods. Finger millet, a coarse cereal known for its exceptional nutritional content and ability to withstand environmental stresses represents a promising climate-smart and nutritional crop in the battle against escalating environmental challenges. The existing traditional improvement programs for finger millet are insufficient to address global hunger effectively. The lack of utilization of high-throughput platforms, genome editing, haplotype breeding, and advanced breeding approaches hinders the systematic multi-omics studies on finger millet, which are essential for pinpointing crucial genes related to agronomically important and various stress responses. The growing environmental uncertainties have widened the gap between the anticipated and real progress in crop improvement. To overcome these challenges a combination of cutting-edge multi-omics techniques such as high-throughput sequencing, speed breeding, mutational breeding, haplotype-based breeding, genomic selection, high-throughput phenotyping, pangenomics, genome editing, and more along with integration of deep learning and artificial intelligence technologies are essential to accelerate research efforts in finger millet. The scarcity of multi-omics approaches in finger millet leaves breeders with limited modern tools for crop enhancement. Therefore, leveraging datasets from previous studies could prove effective in implementing the necessary multi-omics interventions to enrich the genetic resource in finger millet.

TWAS facilitates gene-scale trait genetic dissection through gene expression, structural variations, and alternative splicing in soybean

A genome-wide association study (GWAS) identifies trait-associated loci, but identifying the causal genes can be a bottleneck, due in part to slow decay of linkage disequilibrium (LD). A transcriptome-wide association study (TWAS) addresses this issue by identifying gene expression-phenotype associations or integrating gene expression quantitative trait loci with GWAS results. Here, we used self-pollinated soybean (Glycine max [L.] Merr.) as a model to evaluate the application of TWAS to the genetic dissection of traits in plant species with slow LD decay. We generated RNA sequencing data for a soybean diversity panel and identified the genetic expression regulation of 29 286 soybean genes. Different TWAS solutions were less affected by LD and were robust to the source of expression, identifing known genes related to traits from different tissues and developmental stages. The novel pod-color gene L2 was identified via TWAS and functionally validated by genome editing. By introducing a new exon proportion feature, we significantly improved the detection of expression variations that resulted from structural variations and alternative splicing. As a result, the genes identified through our TWAS approach exhibited a diverse range of causal variations, including SNPs, insertions or deletions, gene fusion, copy number variations, and alternative splicing. Using this approach, we identified genes associated with flowering time, including both previously known genes and novel genes that had not previously been linked to this trait, providing insights complementary to those from GWAS. In summary, this study supports the application of TWAS for candidate gene identification in species with low rates of LD decay.

Transcriptome-wide association study coupled with eQTL analysis reveals the genetic connection between gene expression and flowering time in Arabidopsis

Genome-wide association study (GWAS) has improved our understanding of complex traits, but challenges exist in distinguishing causation versus association caused by linkage disequilibrium. Instead, transcriptome-wide association studies (TWAS) detect direct associations between expression levels and phenotypic variations, providing an opportunity to better prioritize candidate genes. To assess the feasibility of TWAS, we investigated the association between transcriptomes, genomes, and various traits in Arabidopsis, including flowering time. The associated genes formerly known to regulate growth allometry or metabolite production were first identified by TWAS. Next, for flowering time, six TWAS-newly identified genes were functionally validated. Analysis of the expression quantitative trait locus (eQTL) further revealed a trans-regulatory hotspot affecting the expression of several TWAS-identified genes. The hotspot covers the FRIGIDA (FRI) gene body, which possesses multiple haplotypes differentially affecting the expression of downstream genes, such as FLOWERING LOCUS C (FLC) and SUPPRESSOR OF OVEREXPRESSION OF CO 1 (SOC1). We also revealed multiple independent paths towards the loss of function of FRI in natural accessions. Altogether, this study demonstrates the potential of combining TWAS with eQTL analysis to identify important regulatory modules of FRI-FLC-SOC1 for quantitative traits in natural populations.

Integration of eQTL and machine learning to dissect causal genes with pleiotropic effects in genetic regulation networks of seed cotton yield

The dissection of a gene regulatory network (GRN) that complements the genome-wide association study (GWAS) locus and the crosstalk underlying multiple agronomical traits remains a major challenge. In this study, we generate 558 transcriptional profiles of lint-bearing ovules at one day post-anthesis from a selective core cotton germplasm, from which 12,207 expression quantitative trait loci (eQTLs) are identified. Sixty-six known phenotypic GWAS loci are colocalized with 1,090 eQTLs, forming 38 functional GRNs associated predominantly with seed yield. Of the eGenes, 34 exhibit pleiotropic effects. Combining the eQTLs within the seed yield GRNs significantly increases the portion of narrow-sense heritability. The extreme gradient boosting (XGBoost) machine learning approach is applied to predict seed cotton yield phenotypes on the basis of gene expression. Top-ranking eGenes (NF-YB3, FLA2, and GRDP1) derived with pleiotropic effects on yield traits are validated, along with their potential roles by correlation analysis, domestication selection analysis, and transgenic plants.

A complete telomere-to-telomere assembly of the maize genome

A complete telomere-to-telomere (T2T) finished genome has been the long pursuit of genomic research. Through generating deep coverage ultralong Oxford Nanopore Technology (ONT) and PacBio HiFi reads, we report here a complete genome assembly of maize with each chromosome entirely traversed in a single contig. The 2,178.6 Mb T2T Mo17 genome with a base accuracy of over 99.99% unveiled the structural features of all repetitive regions of the genome. There were several super-long simple-sequence-repeat arrays having consecutive thymine-adenine-guanine (TAG) tri-nucleotide repeats up to 235 kb. The assembly of the entire nucleolar organizer region of the 26.8 Mb array with 2,974 45S rDNA copies revealed the enormously complex patterns of rDNA duplications and transposon insertions. Additionally, complete assemblies of all ten centromeres enabled us to precisely dissect the repeat compositions of both CentC-rich and CentC-poor centromeres. The complete Mo17 genome represents a major step forward in understanding the complexity of the highly recalcitrant repetitive regions of higher plant genomes.

The elite variations in germplasms for soybean breeding

The genetic base of soybean cultivars (Glycine max (L.) Merr.) has been narrowed through selective domestication and specific breeding improvement, similar to other crops. This presents challenges in breeding new cultivars with improved yield and quality, reduced adaptability to climate change, and increased susceptibility to diseases. On the other hand, the vast collection of soybean germplasms offers a potential source of genetic variations to address those challenges, but it has yet to be fully leveraged. In recent decades, rapidly improved high-throughput genotyping technologies have accelerated the harness of elite variations in soybean germplasm and provided the important information for solving the problem of a narrowed genetic base in breeding. In this review, we will overview the situation of maintenance and utilization of soybean germplasms, various solutions provided for different needs in terms of the number of molecular markers, and the omics-based high-throughput strategies that have been used or can be used to identify elite alleles. We will also provide an overall genetic information generated from soybean germplasms in yield, quality traits, and pest resistance for molecular breeding.

The 2020 derecho revealed limited overlap between maize genes associated with root lodging and root system architecture

Roots anchor plants in soil, and the failure of anchorage (i.e., root lodging) is a major cause of crop yield loss. Anchorage is often assumed to be driven by root system architecture. We made use of a natural experiment to measure the overlap between the genetic regulation of root system architecture and anchorage. After one of the most devastating derechos ever recorded in August 2020, we phenotyped root lodging in a maize (Zea mays) diversity panel consisting of 369 genotypes grown in six environments affected by the derecho. Genome-wide association studies and transcriptome-wide association studies identified 118 candidate genes associated with root lodging. Thirty-four percent (40/118) of these were homologs of genes from Arabidopsis (Arabidopsis thaliana) that affect traits such as root morphology and lignin content, expected to affect root lodging. Finally, Gene Ontology enrichment analysis of the candidate genes and their predicted interaction partners at the transcriptional and translational levels revealed the complex regulatory networks of physiological and biochemical pathways underlying root lodging in maize. Limited overlap between genes associated with lodging resistance and root system architecture in this diversity panel suggests that anchorage depends in part on factors other than gross characteristics of root system architecture.

Methodology for Constructing a Knowledgebase for Plant Gene Regulation Information

The amount of biological data is growing at a rapid pace as many high-throughput omics technologies and data pipelines are developed. This is resulting in the growth of databases for DNA and protein sequences, gene expression, protein accumulation, structural, and localization information. The diversity and multi-omics nature of such bioinformatic data requires well-designed databases for flexible organization and presentation. Besides general-purpose online bioinformatic databases, users need narrowly focused online databases to quickly access a meaningful collection of related data for their research. Here, we describe the methodology used to implement a plant gene regulatory knowledgebase, with data, query, and tool features, as well as the ability to expand to accommodate future datasets. We exemplify this methodology for the GRASSIUS knowledgebase, but it is applicable to developing and updating similar plant gene regulatory knowledgebases. GRASSIUS organizes and presents gene regulatory data from grass species with a central focus on maize (Zea mays). The main class of data presented include not only the families of transcription factors (TFs) and co-regulators (CRs) but also protein-DNA interaction data, where available.

Association mapping across a multitude of traits collected in diverse environments in maize

Classical genetic studies have identified many cases of pleiotropy where mutations in individual genes alter many different phenotypes. Quantitative genetic studies of natural genetic variants frequently examine one or a few traits, limiting their potential to identify pleiotropic effects of natural genetic variants. Widely adopted community association panels have been employed by plant genetics communities to study the genetic basis of naturally occurring phenotypic variation in a wide range of traits. High-density genetic marker data-18M markers-from 2 partially overlapping maize association panels comprising 1,014 unique genotypes grown in field trials across at least 7 US states and scored for 162 distinct trait data sets enabled the identification of of 2,154 suggestive marker-trait associations and 697 confident associations in the maize genome using a resampling-based genome-wide association strategy. The precision of individual marker-trait associations was estimated to be 3 genes based on a reference set of genes with known phenotypes. Examples were observed of both genetic loci associated with variation in diverse traits (e.g., above-ground and below-ground traits), as well as individual loci associated with the same or similar traits across diverse environments. Many significant signals are located near genes whose functions were previously entirely unknown or estimated purely via functional data on homologs. This study demonstrates the potential of mining community association panel data using new higher-density genetic marker sets combined with resampling-based genome-wide association tests to develop testable hypotheses about gene functions, identify potential pleiotropic effects of natural genetic variants, and study genotype-by-environment interaction.

Combining GWAS and TWAS to identify candidate causal genes for tocochromanol levels in maize grain

Tocochromanols (tocopherols and tocotrienols, collectively vitamin E) are lipid-soluble antioxidants important for both plant fitness and human health. The main dietary sources of vitamin E are seed oils that often accumulate high levels of tocopherol isoforms with lower vitamin E activity. The tocochromanol biosynthetic pathway is conserved across plant species but an integrated view of the genes and mechanisms underlying natural variation of tocochromanol levels in seed of most cereal crops remains limited. To address this issue, we utilized the high mapping resolution of the maize Ames panel of ∼1,500 inbred lines scored with 12.2 million single-nucleotide polymorphisms to generate metabolomic (mature grain tocochromanols) and transcriptomic (developing grain) data sets for genetic mapping. By combining results from genome- and transcriptome-wide association studies, we identified a total of 13 candidate causal gene loci, including 5 that had not been previously associated with maize grain tocochromanols: 4 biosynthetic genes (arodeH2 paralog, dxs1, vte5, and vte7) and a plastid S-adenosyl methionine transporter (samt1). Expression quantitative trait locus (eQTL) mapping of these 13 gene loci revealed that they are predominantly regulated by cis-eQTL. Through a joint statistical analysis, we implicated cis-acting variants as responsible for colocalized eQTL and GWAS association signals. Our multiomics approach provided increased statistical power and mapping resolution to enable a detailed characterization of the genetic and regulatory architecture underlying tocochromanol accumulation in maize grain and provided insights for ongoing biofortification efforts to breed and/or engineer vitamin E and antioxidant levels in maize and other cereals.

Genome-wide mediation analysis: an empirical study to connect phenotype with genotype via intermediate transcriptomic data in maize

Mapping genotype to phenotype is an essential topic in genetics and genomics research. As the Omics data become increasingly available, 2-variable methods have been widely applied to associate genotype with the phenotype (genome-wide association study), gene expression with the phenotype (transcriptome-wide association study), and genotype with gene expression. However, signals detected by these 2-variable association methods suffer from low mapping resolution or inexplicit causality between genotype and phenotype, making it challenging to interpret and validate the molecular mechanisms of the underlying genomic variations and the candidate genes. Under the context of genetics research, we hypothesized a causal chain from genotype to phenotype partially mediated by intermediate molecular processes, i.e. gene expression. To test this hypothesis, we applied the high-dimensional mediation analysis, a class of causal inference method with an assumed causal chain from the exposure to the mediator to the outcome, and implemented it with a maize association panel (N = 280 lines). Using 40 publicly available agronomy traits, 66 newly generated metabolite traits, and published RNA-seq data from 7 different tissues, our empirical study detected 736 unique mediating genes. Noticeably, 83/736 (11%) genes were identified in mediating more than 1 trait, suggesting the prevalence of pleiotropic mediating effects. We demonstrated that several identified mediating genes are consistent with their known functions. In addition, our results provided explicit hypotheses for functional validation and suggested that the mediation analysis is a powerful tool to integrate Omics data to connect genotype to phenotype.

Quantitative Resistance Loci to Southern Rust Mapped in a Temperate Maize Diversity Panel

Southern rust is a severe foliar disease of maize (Zea mays) resulting from infection with the obligate biotrophic fungus Puccinia polysora. This disease reduces photosynthetic productivity, which in turn reduces yields, with the greatest yield losses (up to 50%) associated with earlier onset infections. P. polysora urediniospores overwinter only in tropical and subtropical regions but cause outbreaks when environmental conditions favor initial infection. Increased temperatures and humidity during the growing season combined with an increased frequency of moderate winters are likely to increase the frequency of severe southern rust outbreaks in the U.S. Corn Belt. In summer 2020, a severe outbreak of southern rust was observed in eastern Nebraska, United States. We scored a replicated maize association panel planted in Lincoln, NE for disease severity and found that disease incidence and severity showed significant variation among maize genotypes. Genome-wide association studies identified four loci associated with significant quantitative variation in disease severity. These loci were associated with candidate genes with plausible links to quantitative disease resistance. A transcriptome-wide association study identified additional genes associated with disease severity. Together, these results indicate that substantial diversity in resistance to southern rust exists among current temperate-adapted maize germplasm, including several candidate loci that may explain the observed variation in resistance to southern rust.[Formula: see text] Copyright © 2022 The Author(s). This is an open access article distributed under the CC BY 4.0 International license.

Genetic Variation in ZmPAT7 Contributes to Tassel Branch Number in Maize

Tassel branch number (TBN) is one of the important agronomic traits that contribute to the efficiency of seed production and has been selected strongly during the modern maize breeding process. However, the genetic mechanisms of TBN in maize are not entirely clear. In this study, we used a B73 × CML247 recombination inbred lines (RILs) population to detect quantitative trait loci (QTLs) for TBN. A total of four QTLs (qTBN2a, qTBN2b, qTBN4, and qTBN6) and six candidate genes were identified through expression analysis. Further, one of the candidates (GRMZM2G010011, ZmPAT7) encoding an S-acyltransferase was selected to validate its function by CRISPR-Cas9 technology, and its loss-of-function lines showed a significant increase in TBN. A key SNP(-101) variation in the promoter of ZmPAT7 was significantly associated with TBN. A total of 17 distant eQTLs associated with the expression of ZmPAT7 were identified in expression quantitative trait loci (eQTL) analysis, and ZmNAC3 may be a major factor involved in regulating ZmPAT7. These findings of the present study promote our understanding of the genetic basis of tassel architecture and provide new gene resources for maize breeding improvement.

Identification and Validation of Candidate Genes from Genome-Wide Association Studies

Exploiting the statistical associations coming out from a GWAS experiment to identify and validate candidate genes may be potentially difficult and time consuming. To fill the gap between the identification of candidate genes toward their functional validation onto the trait performance, the prioritization of variants underlying the GWAS-associated regions is necessary. In parallel, recent developments in genomics and statistical methods have been achieved notably in human genetic and they are accordingly being adopted in plant breeding toward the study of the genetic architecture of traits to sustain genetic gains. In this chapter, we aim at providing both theoretical and practical aspects underlying three main options including (1) the MetaGWAS analysis, (2) the statistical fine mapping and (3) the integration of functional data toward the identification and validation of candidate genes from a GWAS experiment.

Fast-forward breeding for a food-secure world

Crop production systems need to expand their outputs sustainably to feed a burgeoning human population. Advances in genome sequencing technologies combined with efficient trait mapping procedures accelerate the availability of beneficial alleles for breeding and research. Enhanced interoperability between different omics and phenotyping platforms, leveraged by evolving machine learning tools, will help provide mechanistic explanations for complex plant traits. Targeted and rapid assembly of beneficial alleles using optimized breeding strategies and precise genome editing techniques could deliver ideal crops for the future. Realizing desired productivity gains in the field is imperative for securing an adequate future food supply for 10 billion people.

TWAS results are complementary to and less affected by linkage disequilibrium than GWAS

A genome-wide association study (GWAS) is used to identify genetic markers associated with phenotypic variation. In contrast, a transcriptome-wide association study (TWAS) detects associations between gene expression levels and phenotypic variation. It has previously been shown that in the cross-pollinated species, maize (Zea mays), GWAS, and TWAS identify complementary sets of trait-associated genes, many of which exhibit characteristics of true positives. Here, we extend this conclusion to the self-pollinated species, Arabidopsis thaliana and soybean (Glycine max). Linkage disequilibrium (LD) can result in the identification, via GWAS, of false-positive associations. In all three analyzed plant species, most trait-associated genes identified via TWAS are well separated physically from other candidate genes. Hence, TWAS is less affected by LD than is GWAS, demonstrating that TWAS is particularly well suited for association studies in genomes with slow rates of LD decay, such as soybean. TWAS is reasonably robust to the plant organs/tissues used to determine expression levels. In summary, this study confirms that TWAS is a promising approach for accurate gene-level association mapping in plants that is complementary to GWAS, and established that TWAS can exhibit substantial advantages relative to GWAS in species with slow rates of LD decay.

Using high-throughput multiple optical phenotyping to decipher the genetic architecture of maize drought tolerance

BACKGROUND

Drought threatens the food supply of the world population. Dissecting the dynamic responses of plants to drought will be beneficial for breeding drought-tolerant crops, as the genetic controls of these responses remain largely unknown.

RESULTS

Here we develop a high-throughput multiple optical phenotyping system to noninvasively phenotype 368 maize genotypes with or without drought stress over a course of 98 days, and collected multiple optical images, including color camera scanning, hyperspectral imaging, and X-ray computed tomography images. We develop high-throughput analysis pipelines to extract image-based traits (i-traits). Of these i-traits, 10,080 were effective and heritable indicators of maize external and internal drought responses. An i-trait-based genome-wide association study reveals 4322 significant locus-trait associations, representing 1529 quantitative trait loci (QTLs) and 2318 candidate genes, many that co-localize with previously reported maize drought responsive QTLs. Expression QTL (eQTL) analysis uncovers many local and distant regulatory variants that control the expression of the candidate genes. We use genetic mutation analysis to validate two new genes, ZmcPGM2 and ZmFAB1A, which regulate i-traits and drought tolerance. Moreover, the value of the candidate genes as drought-tolerant genetic markers is revealed by genome selection analysis, and 15 i-traits are identified as potential markers for maize drought tolerance breeding.

CONCLUSION

Our study demonstrates that combining high-throughput multiple optical phenotyping and GWAS is a novel and effective approach to dissect the genetic architecture of complex traits and clone drought-tolerance associated genes.

An integrated framework reinstating the environmental dimension for GWAS and genomic selection in crops

Identifying mechanisms and pathways involved in gene-environment interplay and phenotypic plasticity is a long-standing challenge. It is highly desirable to establish an integrated framework with an environmental dimension for complex trait dissection and prediction. A critical step is to identify an environmental index that is both biologically relevant and estimable for new environments. With extensive field-observed complex traits, environmental profiles, and genome-wide single nucleotide polymorphisms for three major crops (maize, wheat, and oat), we demonstrated that identifying such an environmental index (i.e., a combination of environmental parameter and growth window) enables genome-wide association studies and genomic selection of complex traits to be conducted with an explicit environmental dimension. Interestingly, genes identified for two reaction-norm parameters (i.e., intercept and slope) derived from flowering time values along the environmental index were less colocalized for a diverse maize panel than for wheat and oat breeding panels, agreeing with the different diversity levels and genetic constitutions of the panels. In addition, we showcased the usefulness of this framework for systematically forecasting the performance of diverse germplasm panels in new environments. This general framework and the companion CERIS-JGRA analytical package should facilitate biologically informed dissection of complex traits, enhanced performance prediction in breeding for future climates, and coordinated efforts to enrich our understanding of mechanisms underlying phenotypic variation.

Expression genome-wide association study identifies that phosphatidylinositol-derived signalling regulates ALUMINIUM SENSITIVE3 expression under aluminium stress in the shoots of Arabidopsis thaliana

To identify unknown regulatory mechanisms leading to aluminium (Al)-induction of the Al tolerance gene ALS3, we conducted an expression genome-wide association study (eGWAS) for ALS3 in the shoots of 95 Arabidopsis thaliana accessions in the presence of Al. The eGWAS was conducted using a mixed linear model with 145,940 genome-wide single nucleotide polymorphisms (SNPs) and the association results were validated using reverse genetics. We found that many SNPs from the eGWAS were associated with genes related to phosphatidylinositol metabolism as well as stress signal transduction, including Ca²⁺signals, inter-connected in a co-expression network. Of these, PLC9, CDPK32, ANAC071, DIR1, and a hypothetical protein (AT4G10470) possessed amino acid sequence/ gene expression level polymorphisms that were significantly associated with ALS3 expression level variation. Furthermore, T-DNA insertion mutants of PLC9, CDPK32, and ANAC071 suppressed shoot ALS3 expression in the presence of Al. This study clarified the regulatory mechanisms of ALS3 expression in the shoot and provided genetic evidence of the involvement of phosphatidylinositol-derived signal transduction under Al stress.

Pituitary Tumor Transforming Gene 1 Orchestrates Gene Regulatory Variation in Mouse Ventral Midbrain During Aging

Dopaminergic neurons in the midbrain are of particular interest due to their role in diseases such as Parkinson’s disease and schizophrenia. Genetic variation between individuals can affect the integrity and function of dopaminergic neurons but the DNA variants and molecular cascades modulating dopaminergic neurons and other cells types of ventral midbrain remain poorly defined. Three genetically diverse inbred mouse strains – C57BL/6J, A/J, and DBA/2J – differ significantly in their genomes (∼7 million variants), motor and cognitive behavior, and susceptibility to neurotoxins. To further dissect the underlying molecular networks responsible for these variable phenotypes, we generated RNA-seq and ChIP-seq data from ventral midbrains of the 3 mouse strains. We defined 1000–1200 transcripts that are differentially expressed among them. These widespread differences may be due to altered activity or expression of upstream transcription factors. Interestingly, transcription factors were significantly underrepresented among the differentially expressed genes, and only one transcription factor, Pttg1, showed significant differences between all three strains. The changes in Pttg1 expression were accompanied by consistent alterations in histone H3 lysine 4 trimethylation at Pttg1 transcription start site. The ventral midbrain transcriptome of 3-month-old C57BL/6J congenic Pttg1^–/– mutants was only modestly altered, but shifted toward that of A/J and DBA/2J in 9-month-old mice. Principle component analysis (PCA) identified the genes underlying the transcriptome shift and deconvolution of these bulk RNA-seq changes using midbrain single cell RNA-seq data suggested that the changes were occurring in several different cell types, including neurons, oligodendrocytes, and astrocytes. Taken together, our results show that Pttg1 contributes to gene regulatory variation between mouse strains and influences mouse midbrain transcriptome during aging.

Heritability jointly explained by host genotype and microbiome: will improve traits prediction?

As we observe the $70$th anniversary of the publication by Robertson that formalized the notion of 'heritability', geneticists remain puzzled by the problem of missing/hidden heritability, where heritability estimates from genome-wide association studies (GWASs) fall short of that from twin-based studies. Many possible explanations have been offered for this discrepancy, including existence of genetic variants poorly captured by existing arrays, dominance, epistasis and unaccounted-for environmental factors; albeit these remain controversial. We believe a substantial part of this problem could be solved or better understood by incorporating the host's microbiota information in the GWAS model for heritability estimation and may also increase human traits prediction for clinical utility. This is because, despite empirical observations such as (i) the intimate role of the microbiome in many complex human phenotypes, (ii) the overlap between genetic variants associated with both microbiome attributes and complex diseases and (iii) the existence of heritable bacterial taxa, current GWAS models for heritability estimate do not take into account the contributory role of the microbiome. Furthermore, heritability estimate from twin-based studies does not discern microbiome component of the observed total phenotypic variance. Here, we summarize the concept of heritability in GWAS and microbiome-wide association studies, focusing on its estimation, from a statistical genetics perspective. We then discuss a possible statistical method to incorporate the microbiome in the estimation of heritability in host GWAS.

Maize adaptation across temperate climates was obtained via expression of two florigen genes

Expansion of the maize growing area was central for food security in temperate regions. In addition to the suppression of the short-day requirement for floral induction, it required breeding for a large range of flowering time that compensates the effect of South-North gradients of temperatures. Here we show the role of a novel florigen gene, ZCN12, in the latter adaptation in cooperation with ZCN8. Strong eQTLs of ZCN8 and ZCN12, measured in 327 maize lines, accounted for most of the genetic variance of flowering time in platform and field experiments. ZCN12 had a strong effect on flowering time of transgenic Arabidopsis thaliana plants; a path analysis showed that it directly affected maize flowering time together with ZCN8. The allelic composition at ZCN QTLs showed clear signs of selection by breeders. This suggests that florigens played a central role in ensuring a large range of flowering time, necessary for adaptation to temperate areas.

Genome-Phenome Wide Association in Maize and Arabidopsis Identifies a Common Molecular and Evolutionary Signature

Linking natural genetic variation to trait variation can help determine the functional roles ofdifferent genes. Variations of one or several traits are often assessed separately. High-throughput phenotyping and data mining can capture dozens or hundreds of traits from the same individuals. Here, we test the association between markers within a gene and many traits simultaneously. This genome-phenome wide association study (GPWAS) is both a multi-marker and multi-trait test. Genes identified using GPWAS with 260 phenotypic traits in maize were enriched for genes independently linked to phenotypic variation. Traits associated with classical mutants were consistent with reported phenotypes for mutant alleles. Genes linked to phenomic variation in maize using GPWAS shared molecular, population genetic, and evolutionary features with classical mutants in maize. Genes linked to phenomic variation in Arabidopsis using GPWAS are significantly enriched in genes with known loss-of-function phenotypes. GPWAS may be an effective strategy to identify genes in which loss-of-function alleles produce mutant phenotypes. The shared signatures present in classical mutants and genes identified using GPWAS may be markers for genes with a role in specifying plant phenotypes generally or pleiotropy specifically.

Meta Gene Regulatory Networks in Maize Highlight Functionally Relevant Regulatory Interactions

The regulation of gene expression is central to many biological processes. Gene regulatory networks (GRNs) link transcription factors (TFs) to their target genes and represent maps of potential transcriptional regulation. Here, we analyzed a large number of publically available maize (Zea mays) transcriptome data sets including >6000 RNA sequencing samples to generate 45 coexpression-based GRNs that represent potential regulatory relationships between TFs and other genes in different populations of samples (cross-tissue, cross-genotype, and tissue-and-genotype samples). While these networks are all enriched for biologically relevant interactions, different networks capture distinct TF-target associations and biological processes. By examining the power of our coexpression-based GRNs to accurately predict covarying TF-target relationships in natural variation data sets, we found that presence/absence changes rather than quantitative changes in TF gene expression are more likely associated with changes in target gene expression. Integrating information from our TF-target predictions and previous expression quantitative trait loci (eQTL) mapping results provided support for 68 TFs underlying 74 previously identified trans-eQTL hotspots spanning a variety of metabolic pathways. This study highlights the utility of developing multiple GRNs within a species to detect putative regulators of important plant pathways and provides potential targets for breeding or biotechnological applications.

Shared Genetic Control of Root System Architecture between Zea mays and Sorghum bicolor

Determining the genetic control of root system architecture (RSA) in plants via large-scale genome-wide association study (GWAS) requires high-throughput pipelines for root phenotyping. We developed Core Root Excavation using Compressed-air (CREAMD), a high-throughput pipeline for the cleaning of field-grown roots, and Core Root Feature Extraction (COFE), a semiautomated pipeline for the extraction of RSA traits from images. CREAMD-COFE was applied to diversity panels of maize (Zea mays) and sorghum (Sorghum bicolor), which consisted of 369 and 294 genotypes, respectively. Six RSA-traits were extracted from images collected from >3,300 maize roots and >1,470 sorghum roots. Single nucleotide polymorphism (SNP)-based GWAS identified 87 TAS (trait-associated SNPs) in maize, representing 77 genes and 115 TAS in sorghum. An additional 62 RSA-associated maize genes were identified via expression read depth GWAS. Among the 139 maize RSA-associated genes (or their homologs), 22 (16%) are known to affect RSA in maize or other species. In addition, 26 RSA-associated genes are coregulated with genes previously shown to affect RSA and 51 (37% of RSA-associated genes) are themselves transe-quantitative trait locus for another RSA-associated gene. Finally, the finding that RSA-associated genes from maize and sorghum included seven pairs of syntenic genes demonstrates the conservation of regulation of morphology across taxa.

Dynamic Patterns of Transcript Abundance of Transposable Element Families in Maize

Transposable Elements (TEs) are mobile elements that contribute the majority of DNA sequences in the maize genome. Due to their repetitive nature, genomic studies of TEs are complicated by the difficulty of properly attributing multi-mapped short reads to specific genomic loci. Here, we utilize a method to attribute RNA-seq reads to TE families rather than particular loci in order to characterize transcript abundance for TE families in the maize genome. We applied this method to assess per-family expression of transposable elements in >800 published RNA-seq libraries representing a range of maize development, genotypes, and hybrids. While a relatively small proportion of TE families are transcribed, expression is highly dynamic with most families exhibiting tissue-specific expression. A large number of TE families were specifically detected in pollen and endosperm, consistent with reproductive dynamics that maintain silencing of TEs in the germ line. We find that B73 transcript abundance is a poor predictor of TE expression in other genotypes and that transcript levels can differ even for shared TEs. Finally, by assessing recombinant inbred line and hybrid transcriptomes, complex patterns of TE transcript abundance across genotypes emerged. Taken together, this study reveals a dynamic contribution of TEs to maize transcriptomes.

Transcriptome-Wide Association Supplements Genome-Wide Association in Zea mays

Modern improvement of complex traits in agricultural species relies on successful associations of heritable molecular variation with observable phenotypes. Historically, this pursuit has primarily been based on easily measurable genetic markers. The recent advent of new technologies allows assaying and quantifying biological intermediates (hereafter endophenotypes) which are now readily measurable at a large scale across diverse individuals. The usefulness of endophenotypes for delineating the regulatory landscape of the genome and genetic dissection of complex trait variation remains underexplored in plants. The work presented here illustrated the utility of a large-scale (299-genotype and seven-tissue) gene expression resource to dissect traits across multiple levels of biological organization. Using single-tissue- and multi-tissue-based transcriptome-wide association studies (TWAS), we revealed that about half of the functional variation acts through altered transcript abundance for maize kernel traits, including 30 grain carotenoid abundance traits, 20 grain tocochromanol abundance traits, and 22 field-measured agronomic traits. Comparing the efficacy of TWAS with genome-wide association studies (GWAS) and an ensemble approach that combines both GWAS and TWAS, we demonstrated that results of TWAS in combination with GWAS increase the power to detect known genes and aid in prioritizing likely causal genes. Using a variance partitioning approach in the largely independent maize Nested Association Mapping (NAM) population, we also showed that the most strongly associated genes identified by combining GWAS and TWAS explain more heritable variance for a majority of traits than the heritability captured by the random genes and the genes identified by GWAS or TWAS alone. This not only improves the ability to link genes to phenotypes, but also highlights the phenotypic consequences of regulatory variation in plants.

Transcription factors involved in abiotic stress responses in Maize (Zea mays L.) and their roles in enhanced productivity in the post genomics era

BACKGROUND

Maize (Zea mays L.) is a principal cereal crop cultivated worldwide for human food, animal feed, and more recently as a source of biofuel. However, as a direct consequence of water insufficiency and climate change, frequent occurrences of both biotic and abiotic stresses have been reported in various regions around the world, and recently, this has become a constant threat in increasing global maize yields. Plants respond to abiotic stresses by utilizing the activities of transcription factors (TFs), which are families of genes coding for specific TF proteins. TF target genes form a regulon that is involved in the repression/activation of genes associated with abiotic stress responses. Therefore, it is of utmost importance to have a systematic study on each TF family, the downstream target genes they regulate, and the specific TF genes involved in multiple abiotic stress responses in maize and other staple crops.

METHOD

In this review, the main TF families, the specific TF genes and their regulons that are involved in abiotic stress regulation will be briefly discussed. Great emphasis will be given on maize abiotic stress improvement throughout this review, although other examples from different plants like rice, Arabidopsis, wheat, and barley will be used.

RESULTS

We have described in detail the main TF families in maize that take part in abiotic stress responses together with their regulons. Furthermore, we have also briefly described the utilization of high-efficiency technologies in the study and characterization of TFs involved in the abiotic stress regulatory networks in plants with an emphasis on increasing maize production. Examples of these technologies include next-generation sequencing, microarray analysis, machine learning, and RNA-Seq.

CONCLUSION

In conclusion, it is expected that all the information provided in this review will in time contribute to the use of TF genes in the research, breeding, and development of new abiotic stress tolerant maize cultivars.

ZmMADS69 functions as a flowering activator through the ZmRap2.7-ZCN8 regulatory module and contributes to maize flowering time adaptation

Flowering time is a major determinant of the local adaptation of plants. Although numerous loci affecting flowering time have been mapped in maize, their underlying molecular mechanisms and roles in adaptation remain largely unknown. Here, we report the identification and characterization of MADS-box transcription factor ZmMADS69 that functions as a flowering activator through the ZmRap2.7-ZCN8 regulatory module and contributes to adaptation. We show that ZmMADS69 underlies a quantitative trait locus controlling the difference in flowering time between maize and its wild ancestor, teosinte. Maize ZmMADS69 allele is expressed at a higher level at floral transition and confers earlier flowering than the teosinte allele under long days and short days. Overexpression of ZmMADS69 causes early flowering, while a transposon insertion mutant of ZmMADS69 exhibits delayed flowering. ZmMADS69 shows pleiotropic effects for multiple traits of agronomic importance. ZmMADS69 functions upstream of the flowering repressor ZmRap2.7 to downregulate its expression, thereby relieving the repression of the florigen gene ZCN8 and causing early flowering. Population genetic analyses showed that ZmMADS69 was a target of selection and may have played an important role as maize spread from the tropics to temperate zones. Our findings provide important insights into the regulation and adaptation of flowering time.

Evolutionary and ecological functional genomics, from lab to the wild

Plant phenotypes are the result of both genetic and environmental forces that act to modulate trait expression. Over the last few years, numerous approaches in functional genomics and systems biology have led to a greater understanding of plant phenotypic variation and plant responses to the environment. These approaches, and the questions that they can address, have been loosely termed evolutionary and ecological functional genomics (EEFG), and have been providing key insights on how plants adapt and evolve. In particular, by bringing these studies from the laboratory to the field, EEFG studies allow us to gain greater knowledge of how plants function in their natural contexts.

Insect Transcription Factors: A Landscape of Their Structures and Biological Functions in Drosophila and beyond

Transcription factors (TFs) play essential roles in the transcriptional regulation of functional genes, and are involved in diverse physiological processes in living organisms. The fruit fly Drosophila melanogaster, a simple and easily manipulated organismal model, has been extensively applied to study the biological functions of TFs and their related transcriptional regulation mechanisms. It is noteworthy that with the development of genetic tools such as CRISPR/Cas9 and the next-generation genome sequencing techniques in recent years, identification and dissection the complex genetic regulatory networks of TFs have also made great progress in other insects beyond Drosophila. However, unfortunately, there is no comprehensive review that systematically summarizes the structures and biological functions of TFs in both model and non-model insects. Here, we spend extensive effort in collecting vast related studies, and attempt to provide an impartial overview of the progress of the structure and biological functions of current documented TFs in insects, as well as the classical and emerging research methods for studying their regulatory functions. Consequently, considering the importance of versatile TFs in orchestrating diverse insect physiological processes, this review will assist a growing number of entomologists to interrogate this understudied field, and to propel the progress of their contributions to pest control and even human health.

Empirical Comparisons of Different Statistical Models To Identify and Validate Kernel Row Number-Associated Variants from Structured Multi-parent Mapping Populations of Maize

Advances in next generation sequencing technologies and statistical approaches enable genome-wide dissection of phenotypic traits via genome-wide association studies (GWAS). Although multiple statistical approaches for conducting GWAS are available, the power and cross-validation rates of many approaches have been mostly tested using simulated data. Empirical comparisons of single variant (SV) and multi-variant (MV) GWAS approaches have not been conducted to test if a single approach or a combination of SV and MV is effective, through identification and cross-validation of trait-associated loci. In this study, kernel row number (KRN) data were collected from a set of 6,230 entries derived from the Nested Association Mapping (NAM) population and related populations. Three different types of GWAS analyses were performed: 1) single-variant (SV), 2) stepwise regression (STR) and 3) a Bayesian-based multi-variant (BMV) model. Using SV, STR, and BMV models, 257, 300, and 442 KRN-associated variants (KAVs) were identified in the initial GWAS analyses. Of these, 231 KAVs were subjected to genetic validation using three unrelated populations that were not included in the initial GWAS. Genetic validation results suggest that the three GWAS approaches are complementary. Interestingly, KAVs in low recombination regions were more likely to exhibit associations in independent populations than KAVs in recombinationally active regions, probably as a consequence of linkage disequilibrium. The KAVs identified in this study have the potential to enhance our understanding of the genetic basis of ear development.

Intragenic Meiotic Crossovers Generate Novel Alleles with Transgressive Expression Levels

Meiotic recombination is an evolutionary force that generates new genetic diversity upon which selection can act. Whereas multiple studies have assessed genome-wide patterns of recombination and specific cases of intragenic recombination, few studies have assessed intragenic recombination genome-wide in higher eukaryotes. We identified recombination events within or near genes in a population of maize recombinant inbred lines (RILs) using RNA-sequencing data. Our results are consistent with case studies that have shown that intragenic crossovers cluster at the 5′ ends of some genes. Further, we identified cases of intragenic crossovers that generate transgressive transcript accumulation patterns, that is, recombinant alleles displayed higher or lower levels of expression than did nonrecombinant alleles in any of ∼100 RILs, implicating intragenic recombination in the generation of new variants upon which selection can act. Thousands of apparent gene conversion events were identified, allowing us to estimate the genome-wide rate of gene conversion at SNP sites (4.9 × 10⁻⁵). The density of syntenic genes (i.e., those conserved at the same genomic locations since the divergence of maize and sorghum) exhibits a substantial correlation with crossover frequency, whereas the density of nonsyntenic genes (i.e., those which have transposed or been lost subsequent to the divergence of maize and sorghum) shows little correlation, suggesting that crossovers occur at higher rates in syntenic genes than in nonsyntenic genes. Increased rates of crossovers in syntenic genes could be either a consequence of the evolutionary conservation of synteny or a biological process that helps to maintain synteny.