The B73 Maize Genome: Complexity, Diversity, and Dynamics

Loss of Lateral suppressor gene is associated with evolution of root nodule symbiosis in Leguminosae

BACKGROUND

Root nodule symbiosis (RNS) is a fascinating evolutionary event. Given that limited genes conferring the evolution of RNS in Leguminosae have been functionally validated, the genetic basis of the evolution of RNS remains largely unknown. Identifying the genes involved in the evolution of RNS will help to reveal the mystery.

RESULTS

Here, we investigate the gene loss event during the evolution of RNS in Leguminosae through phylogenomic and synteny analyses in 48 species including 16 Leguminosae species. We reveal that loss of the Lateral suppressor gene, a member of the GRAS-domain protein family, is associated with the evolution of RNS in Leguminosae. Ectopic expression of the Lateral suppressor (Ls) gene from tomato and its homolog MONOCULM 1 (MOC1) and Os7 from rice in soybean and Medicago truncatula result in almost completely lost nodulation capability. Further investigation shows that Lateral suppressor protein, Ls, MOC1, and Os7 might function through an interaction with NODULATION SIGNALING PATHWAY 2 (NSP2) and CYCLOPS to repress the transcription of NODULE INCEPTION (NIN) to inhibit the nodulation in Leguminosae. Additionally, we find that the cathepsin H (CTSH), a conserved protein, could interact with Lateral suppressor protein, Ls, MOC1, and Os7 and affect the nodulation.

CONCLUSIONS

This study sheds light on uncovering the genetic basis of the evolution of RNS in Leguminosae and suggests that gene loss plays an essential role.

Biochemical and Epigenetic Regulation of Glutamate Metabolism in Maize (Zea mays L.) Leaves under Salt Stress

The effect of salt stress (150 mM NaCl) on the expression of genes, methylation of their promoters, and enzymatic activity of glutamate dehydrogenase (GDH), glutamate decarboxylase (GAD), and the 2-oxoglutarate (2-OG)-dehydrogenase (2-OGDH) complex was studied in maize (Zea mays L.). GDH activity increased continuously under salt stress, being 3-fold higher after 24 h. This was accompanied by the appearance of a second isoform with lower electrophoretic mobility. The expression of the Gdh1 gene strongly increased after 6-12 h of incubation, which corresponded to the demethylation of its promoter, while Gdh2 gene expression slightly increased after 2-6 h and then decreased. GAD activity gradually increased in the first 12 h, and then returned to the control level. This corresponded to the increase of Gad expression and its demethylation. Salt stress led to a 2-fold increase in the activity of 2-OGDH during the first 6 h of NaCl treatment, then the activity returned to the control level. Expression of the genes Ogdh1 and Ogdh3 peaked after 1-2 h of incubation. After 6-8 h with NaCl, the expression of these genes declined below the control levels, which correlated with the higher methylation of their promoters. We conclude that salt stress causes a redirection of the 2-OG flux to the γ-aminobutyric acid shunt via its amination to glutamate, by altering the expression of the Gdh1 and Gdh2 genes, which likely promotes the assembly of the native GDH molecule having a different subunit composition and greater affinity for 2-OG.

First insight into the genomes of the Pulmonaria officinalis group (Boraginaceae) provided by repeatome analysis and comparative karyotyping

BACKGROUND

The genus Pulmonaria (Boraginaceae) represents a taxonomically complex group of species in which morphological similarity contrasts with striking karyological variation. The presence of different numbers of chromosomes in the diploid state suggests multiple hybridization/polyploidization events followed by chromosome rearrangements (dysploidy). Unfortunately, the phylogenetic relationships and evolution of the genome, have not yet been elucidated. Our study focused on the P. officinalis group, the most widespread species complex, which includes two morphologically similar species that differ in chromosome number, i.e. P. obscura (2n = 14) and P. officinalis (2n = 16). Ornamental cultivars, morphologically similar to P. officinalis (garden escapes), whose origin is unclear, were also studied. Here, we present a pilot study on genome size and repeatome dynamics of these closely related species in order to gain new information on their genome and chromosome structure.

RESULTS

Flow cytometry confirmed a significant difference in genome size between P. obscura and P. officinalis, corresponding to the number of chromosomes. Genome-wide repeatome analysis performed on genome skimming data showed that retrotransposons were the most abundant repeat type, with a higher proportion of Ty3/Gypsy elements, mainly represented by the Tekay lineage. Comparative analysis revealed no species-specific retrotransposons or striking differences in their copy number between the species. A new set of chromosome-specific cytogenetic markers, represented by satellite DNAs, showed that the chromosome structure in P. officinalis was more variable compared to that of P. obscura. Comparative karyotyping supported the hybrid origin of putative hybrids with 2n = 15 collected from a mixed population of both species and outlined the origin of ornamental garden escapes, presumably derived from the P. officinalis complex.

CONCLUSIONS

Large-scale genome size analysis and repeatome characterization of the two morphologically similar species of the P. officinalis group improved our knowledge of the genome dynamics and differences in the karyotype structure. A new set of chromosome-specific cytogenetic landmarks was identified and used to reveal the origin of putative hybrids and ornamental cultivars morphologically similar to P. officinalis.

Genetic analysis of QTLs for lysine content in four maize DH populations

BACKGROUND

Low levels of the essential amino acid lysine in maize endosperm is considered to be a major problem regarding the nutritional quality of food and feed. Increasing the lysine content of maize is important to improve the quality of food and feed nutrition. Although the genetic basis of quality protein maize (QPM) has been studied, the further exploration of the quantitative trait loci (QTL) underlying lysine content variation still needs more attention.

RESULTS

Eight maize inbred lines with increased lysine content were used to construct four double haploid (DH) populations for identification of QTLs related to lysine content. The lysine content in the four DH populations exhibited continuous and normal distribution. A total of 12 QTLs were identified in a range of 4.42-12.66% in term of individual phenotypic variation explained (PVE) which suggested the quantitative control of lysine content in maize. Five main genes involved in maize lysine biosynthesis pathways in the QTL regions were identified in this study.

CONCLUSIONS

The information presented will allow the exploration of candidate genes regulating lysine biosynthesis pathways and be useful for marker-assisted selection and gene pyramiding in high-lysine maize breeding programs.

Trait association and prediction through integrative k-mer analysis

Genome-wide association study (GWAS) with single nucleotide polymorphisms (SNPs) has been widely used to explore genetic controls of phenotypic traits. Alternatively, GWAS can use counts of substrings of length k from longer sequencing reads, k-mers, as genotyping data. Using maize cob and kernel color traits, we demonstrated that k-mer GWAS can effectively identify associated k-mers. Co-expression analysis of kernel color k-mers and genes directly found k-mers from known causal genes. Analyzing complex traits of kernel oil and leaf angle resulted in k-mers from both known and candidate genes. A gene encoding a MADS transcription factor was functionally validated by showing that ectopic expression of the gene led to less upright leaves. Evolution analysis revealed most k-mers positively correlated with kernel oil were strongly selected against in maize populations, while most k-mers for upright leaf angle were positively selected. In addition, genomic prediction of kernel oil, leaf angle, and flowering time using k-mer data resulted in a similarly high prediction accuracy to the standard SNP-based method. Collectively, we showed k-mer GWAS is a powerful approach for identifying trait-associated genetic elements. Further, our results demonstrated the bridging role of k-mers for data integration and functional gene discovery.

The Zea mays PeptideAtlas: A New Maize Community Resource

This study presents the Maize PeptideAtlas resource (www.peptideatlas.org/builds/maize) to help solve questions about the maize proteome. Publicly available raw tandem mass spectrometry (MS/MS) data for maize collected from ProteomeXchange were reanalyzed through a uniform processing and metadata annotation pipeline. These data are from a wide range of genetic backgrounds and many sample types and experimental conditions. The protein search space included different maize genome annotations for the B73 inbred line from MaizeGDB, UniProtKB, NCBI RefSeq, and for the W22 inbred line. 445 million MS/MS spectra were searched, of which 120 million were matched to 0.37 million distinct peptides. Peptides were matched to 66.2% of proteins in the most recent B73 nuclear genome annotation. Furthermore, most conserved plastid- and mitochondrial-encoded proteins (NCBI RefSeq annotations) were identified. Peptides and proteins identified in the other B73 genome annotations will improve maize genome annotation. We also illustrate the high-confidence detection of unique W22 proteins. N-terminal acetylation, phosphorylation, ubiquitination, and three lysine acylations (K-acetyl, K-malonyl, and K-hydroxyisobutyryl) were identified and can be inspected through a PTM viewer in PeptideAtlas. All matched MS/MS-derived peptide data are linked to spectral, technical, and biological metadata. This new PeptideAtlas is integrated in MaizeGDB with a peptide track in JBrowse.

Insertions and Deletions: Computational Methods, Evolutionary Dynamics, and Biological Applications

Insertions and deletions constitute the second most important source of natural genomic variation. Insertions and deletions make up to 25% of genomic variants in humans and are involved in complex evolutionary processes including genomic rearrangements, adaptation, and speciation. Recent advances in long-read sequencing technologies allow detailed inference of insertions and deletion variation in species and populations. Yet, despite their importance, evolutionary studies have traditionally ignored or mishandled insertions and deletions due to a lack of comprehensive methodologies and statistical models of insertions and deletion dynamics. Here, we discuss methods for describing insertions and deletion variation and modeling insertions and deletions over evolutionary time. We provide practical advice for tackling insertions and deletions in genomic sequences and illustrate our discussion with examples of insertions and deletion-induced effects in human and other natural populations and their contribution to evolutionary processes. We outline promising directions for future developments in statistical methodologies that would allow researchers to analyze insertions and deletion variation and their effects in large genomic data sets and to incorporate insertions and deletions in evolutionary inference.

The biosynthesis of storage reserves and auxin is coordinated by a hierarchical regulatory network in maize endosperm

Grain filling in maize (Zea mays) is intricately linked to cell development, involving the regulation of genes responsible for the biosynthesis of storage reserves (starch, proteins, and lipids) and phytohormones. However, the regulatory network coordinating these biological functions remains unclear. In this study, we identified 1744 high-confidence target genes co-regulated by the transcription factors (TFs) ZmNAC128 and ZmNAC130 (ZmNAC128/130) through chromatin immunoprecipitation sequencing coupled with RNA-seq analysis in the zmnac128/130 loss-of-function mutants. We further constructed a hierarchical regulatory network using DNA affinity purification sequencing analysis of downstream TFs regulated by ZmNAC128/130. In addition to target genes involved in the biosynthesis of starch and zeins, we discovered novel target genes of ZmNAC128/130 involved in the biosynthesis of lipids and indole-3-acetic acid (IAA). Consistently, the number of oil bodies, as well as the contents of triacylglycerol, and IAA were significantly reduced in zmnac128/130. The hierarchical regulatory network centered by ZmNAC128/130 revealed a significant overlap between the direct target genes of ZmNAC128/130 and their downstream TFs, particularly in regulating the biosynthesis of storage reserves and IAA. Our results indicated that the biosynthesis of storage reserves and IAA is coordinated by a multi-TFs hierarchical regulatory network in maize endosperm.

Identification of the Coexisting Virus-Derived siRNA in Maize and Rice Infected by Rice Black-Streaked Dwarf Virus

Rice black-streaked dwarf virus is transmitted by small brown planthoppers, which causes maize rough dwarf disease and rice black-streaked dwarf disease. This virus leads to slow growth or death of the host plants. During the coevolutionary arms race between viruses and plants, virus-derived small interfering RNAs (vsiRNAs) challenge the plant's defense response and inhibit host immunity through the RNA silencing system. However, it is currently unknown if rice black-streaked dwarf virus can produce the same siRNAs to mediate the RNA silencing in different infected species. In this study, four small RNA libraries and four degradome libraries were constructed by extracting total RNAs from the leaves of the maize (Zea mays) inbred line B73 and japonica rice (Oryza sativa) variety Nipponbare exposed to feeding by viruliferous and nonviruliferous small brown planthoppers. We analyzed the characteristics of small RNAs and explored virus-derived siRNAs in small RNA libraries through high-throughput sequencing. On analyzing the characteristics of small RNA, we noted that the size distributions of small RNAs were mainly 24 nt (19.74 to 62.00%), whereas those of vsiRNAs were mostly 21 nt (41.06 to 41.87%) and 22 nt (39.72 to 42.26%). The 5'-terminal nucleotides of vsiRNAs tended to be adenine or uracil. Exploring the distribution of vsiRNA hot spots on the viral genome segments revealed that the frequency of hotspots in B73 was higher than those in Nipponbare. Meanwhile, hotspots in the S9 and S10 virus genome segments were distributed similarly in both hosts. In addition, the target genes of small RNA were explored by degradome sequencing. Analyses of the regulatory pathway of these target genes unveiled that viral infection affected the ribosome-related target genes in maize and the target genes in the metabolism and biosynthesis pathways in rice. Here, 562 and 703 vsiRNAs were separately obtained in maize and rice and 73 vsiRNAs named as coexisting vsiRNAs (co-vsiRNAs) were detected in both hosts. Stem-loop PCR and real-time quantitative PCR confirmed that co-vsiRNA 3.1 and co-vsiRNA 3.5, derived from genome segment S3, simultaneously play a role in maize and rice and inhibited host gene expression. The study revealed that rice black-streaked dwarf virus can produce the same siRNAs in different species and provides a new direction for developing new antiviral strategies.

Genome-wide characterization and expression profiling of FARL (FHY3/FAR1) family genes in Zea mays

A significant role of the plant is played by the transcription factor FARL, which is light signal transduction as well as plant growth and development. Despite being transposases, FARL has developed a variety of dominant biological actions in evolution and speciation. On the other hand, little is known about the Zea mays FARL protein family. This study identifies and characterizes fifteen ZmFARL genes genome-wide, and RNA sequencing data was used to profile their expression. 105 FARL proteins from five plant species were classified into five groups based on sequence alignment and phylogeny. The ZmFARL genes' exon-intron and motif distribution were conserved based on their evolutionary group. The fifteen ZmFARL genes were distributed over seven of the ten Z. mays chromosomes, although no duplication was discovered. Cis-element analysis reveals that ZmFARL genes play a variety of activities, including tissue-specific, stress- and hormone-responsive expressions. Furthermore, the results of the RNA sequencing used to profile expression showed that the genes ZmFARL2 and ZmFARL5 were much more expressed than other genes in various tissues, particularly in leaf characteristics. The identification of likely genes involved in cellular activity in Z. mays and related species will be aided by the characterization of the FARL genes.

Maize mutant screens: from classical methods to new CRISPR-based approaches

Mutations play a pivotal role in shaping the trajectory and outcomes of a species evolution and domestication. Maize (Zea mays) has been a major staple crop and model for genetic research for more than 100 yr. With the arrival of site-directed mutagenesis and genome editing (GE) driven by the Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR), maize mutational research is once again in the spotlight. If we combine the powerful physiological and genetic characteristics of maize with the already available and ever increasing toolbox of CRISPR-Cas, prospects for its future trait engineering are very promising. This review aimed to give an overview of the progression and learnings of maize screening studies analyzing forward genetics, natural variation and reverse genetics to focus on recent GE approaches. We will highlight how each strategy and resource has contributed to our understanding of maize natural and induced trait variability and how this information could be used to design the next generation of mutational screenings.

Genome-wide identification and characterization of SLEEPER, a transposon-derived gene family and their expression pattern in Brassica napus L

BACKGROUND

The transposons of the hAT superfamily are the most widespread transposons ever known. SLEEPER genes encode domesticated transposases from the hAT superfamily, which may have lost their transposable functions during long-term evolution and transformed into host proteins that regulate plant growth and development.

RESULTS

This study identified 162 members of the SLEEPER gene family from Brassica napus. These members are widely distributed on 19 chromosomes, mainly in the Cn subgenome, and have promoters with various cis-acting elements related to hormone regulation, abiotic stress, and growth and development regulation. Most of the genes in this family contain similar conserved domains and motifs, and the closer the genes are distributed on evolutionary branches, the more similar their structures are. Transcriptome sequencing performed on tissues at different growth stages from B. napus line 3529 indicated that these genes had different expression patterns, and nearly half of the genes were not detectably expressed in all samples.

CONCLUSIONS

This study investigated the gene structure, expression patterns, evolutionary features, and gene localization of the SLEEPER family members to confirm the significance of these genes in the growth of B. napus, providing a reference for the study of transposon domestication and outstanding genetic resources for the genetic improvement of B. napus.

Triumphs of genomic-assisted breeding in crop improvement

Conventional breeding approaches have played a significant role in meeting the food demand remarkably well until now. However, the increasing population, yield plateaus in certain crops, and limited recombination necessitate using genomic resources for genomics-assisted crop improvement programs. As a result of advancements in the next-generation sequence technology, GABs have developed dramatically to characterize allelic variants and facilitate their rapid and efficient incorporation in crop improvement programs. Genomics-assisted breeding (GAB) has played an important role in harnessing the potential of modern genomic tools, exploiting allelic variation from genetic resources and developing cultivars over the past decade. The availability of pangenomes for major crops has been a significant development, albeit with varying degrees of completeness. Even though adopting these technologies is essentially determined on economic grounds and cost-effective assays, which create a wealth of information that can be successfully used to exploit the latent potential of crops. GAB has been instrumental in harnessing the potential of modern genomic resources and exploiting allelic variation for genetic enhancement and cultivar development. GAB strategies will be indispensable for designing future crops and are expected to play a crucial role in breeding climate-smart crop cultivars with higher nutritional value.

Evolutionary History of the DD41D Family of Tc1/Mariner Transposons in Two Mayetiola Species

Tc1/mariner elements are ubiquitous in eukaryotic genomes including insects. They are diverse and divided into families and sub-families. The DD34D family including mauritiana and irritans subfamilies have already been identified in two closely related species of Cecidomyiids M. destructor and M. hordei. In the current study the de novo and similarity-based methods allowed the identification for the first time of seven consensuses in M. destructor and two consensuses in M. hordei belonging to DD41D family whereas the in vitro method allowed the amplification of two and three elements in these two species respectively. Most of identified elements accumulated different mutations and long deletions spanning the N-terminal region of the transposase. Phylogenetic analyses showed that the DD41D elements were clustered in two groups belonging to rosa and Long-TIR subfamilies. The age estimation of the last transposition events of the identified Tc1/mariner elements in M. destructor showed different evolutionary histories. Indeed, irritans elements have oscillated between periods of silencing and reappearance while rosa and mauritiana elements have shown regular activity with large recent bursts. The study of insertion sites showed that they are mostly intronic and that some recently transposed elements occurred in genes linked to putative DNA-binding domains and enzymes involved in metabolic chains. Thus, this study gave evidence of the existence of DD41D family in two Mayetiola species and an insight on their evolutionary history.

Evolution of Einkorn wheat centromeres is driven by the mutualistic interplay of two LTR retrotransposons

BACKGROUND

Centromere function is highly conserved across eukaryotes, but the underlying centromeric DNA sequences vary dramatically between species. Centromeres often contain a high proportion of repetitive DNA, such as tandem repeats and/or transposable elements (TEs). Einkorn wheat centromeres lack tandem repeat arrays and are instead composed mostly of the two long terminal repeat (LTR) retrotransposon families RLG_Cereba and RLG_Quinta which specifically insert in centromeres. However, it is poorly understood how these two TE families relate to each other and if and how they contribute to centromere function and evolution.

RESULTS

Based on conservation of diagnostic motifs (LTRs, integrase and primer binding site and polypurine-tract), we propose that RLG_Cereba and RLG_Quinta are a pair of autonomous and non-autonomous partners, in which the autonomous RLG_Cereba contributes all the proteins required for transposition, while the non-autonomous RLG_Quinta contributes GAG protein. Phylogenetic analysis of predicted GAG proteins showed that the RLG_Cereba lineage was present for at least 100 million years in monocotyledon plants. In contrast, RLG_Quinta evolved from RLG_Cereba between 28 and 35 million years ago in the common ancestor of oat and wheat. Interestingly, the integrase of RLG_Cereba is fused to a so-called CR-domain, which is hypothesized to guide the integrase to the functional centromere. Indeed, ChIP-seq data and TE population analysis show only the youngest subfamilies of RLG_Cereba and RLG_Quinta are found in the active centromeres. Importantly, the LTRs of RLG_Quinta and RLG_Cereba are strongly associated with the presence of the centromere-specific CENH3 histone variant. We hypothesize that the LTRs of RLG_Cereba and RLG_Quinta contribute to wheat centromere integrity by phasing and/or placing CENH3 nucleosomes, thus favoring their persistence in the competitive centromere-niche.

CONCLUSION

Our data show that RLG_Cereba cross-mobilizes the non-autonomous RLG_Quinta retrotransposons. New copies of both families are specifically integrated into functional centromeres presumably through direct binding of the integrase CR domain to CENH3 histone variants. The LTRs of newly inserted RLG_Cereba and RLG_Quinta elements, in turn, recruit and/or phase new CENH3 deposition. This mutualistic interplay between the two TE families and the plant host dynamically maintains wheat centromeres.

Genomic and Transcriptomic Insights into the Evolution of C4 Photosynthesis in Grasses

C4 photosynthesis has independently evolved over 62 times within 19 angiosperm families. The recurrent evolution of C4 photosynthesis appears to contradict the complex anatomical and biochemical modifications required for the transition from C3 to C4 photosynthesis. In this study, we conducted an integrated analysis of genomics and transcriptomics to elucidate the molecular underpinnings of convergent C4 evolution in the grass family. Our genome-wide exploration of C4-related gene families suggests that the expansion of these gene families may have played an important role in facilitating C4 evolution in the grass family. A phylogenomic synteny network analysis uncovered the emergence of C4 genes in various C4 grass lineages from a common ancestral gene pool. Moreover, through a comparison between non-C4 and C4  PEPCs, we pinpointed 14 amino acid sites exhibiting parallel adaptations. These adaptations, occurring post the BEP-PACMAD divergence, shed light on why all C4 origins in grasses are confined to the PACMAD clade. Furthermore, our study revealed that the ancestor of Chloridoideae grasses possessed a more favorable molecular preadaptation for C4 functions compared to the ancestor of Panicoideae grasses. This molecular preadaptation potentially explains why C4 photosynthesis evolved earlier in Chloridoideae than in Panicoideae and why the C3-to-C4 transition occurred once in Chloridoideae but multiple times in Panicoideae. Additionally, we found that C4 genes share similar cis-elements across independent C4 lineages. Notably, NAD-ME subtype grasses may have retained the ancestral regulatory machinery of the C4  NADP-ME gene, while NADP-ME subtype grasses might have undergone unique cis-element modifications.

Genetic analyses of embryo homology and ontogeny in the model grass Zea mays subsp. mays

The homology of the single cotyledon of grasses and the ontogeny of the scutellum and coleoptile as the initial, highly modified structures of the grass embryo are investigated using leaf developmental genetics and targeted transcript analyses in the model grass Zea mays subsp. mays. Transcripts of leaf developmental genes are identified in both the initiating scutellum and the coleoptile, while mutations disrupting mediolateral leaf development also disrupt scutellum and coleoptile morphology, suggesting that these grass-specific organs are modified leaves. Higher-order mutations in WUSCHEL-LIKE HOMEOBOX3 (WOX3) genes, involved in mediolateral patterning of plant lateral organs, inform a model for the fusion of coleoptilar margins during maize embryo development. Genetic, RNA-targeting, and morphological evidence supports models for cotyledon evolution where the scutellum and coleoptile, respectively, comprise the distal and proximal domains of the highly modified, single grass cotyledon.

Engineering the future cereal crops with big biological data: toward intelligence-driven breeding by design

How to feed 10 billion human populations is one of the challenges that need to be addressed in the following decades, especially under an unpredicted climate change. Crop breeding, initiating from the phenotype-based selection by local farmers and developing into current biotechnology-based breeding, has played a critical role in securing the global food supply. However, regarding the changing environment and ever-increasing human population, can we breed outstanding crop varieties fast enough to achieve high productivity, good quality, and widespread adaptability? This review outlines the recent achievements in understanding cereal crop breeding, including the current knowledge about crop agronomic traits, newly developed techniques, crop big biological data research, and the possibility of integrating them for intelligence-driven breeding by design, which ushers in a new era of crop breeding practice and shapes the novel architecture of future crops. This review focuses on the major cereal crops, including rice, maize, and wheat, to explain how intelligence-driven breeding by design is becoming a reality.

Genome-wide identification and expression analysis of the U-box E3 ubiquitin ligase gene family related to bacterial wilt resistance in tobacco (Nicotiana tabacum L.) and eggplant (Solanum melongena L.)

The E3 enzyme in the UPS pathway is a crucial factor for inhibiting substrate specificity. In Solanaceae, the U-box E3 ubiquitin ligase has a complex relationship with plant growth and development, and plays a pivotal role in responding to various biotic and abiotic stresses. The analysis of the U-box gene family in Solanaceae and its expression profile under different stresses holds significant implications. A total of 116 tobacco NtU-boxs and 56 eggplant SmU-boxs were identified based on their respective genome sequences. Phylogenetic analysis of U-box genes in tobacco, eggplant, tomato, Arabidopsis, pepper, and potato revealed five distinct subgroups (I-V). Gene structure and protein motifs analysis found a high degree of conservation in both exon/intron organization and protein motifs among tobacco and eggplant U-box genes especially the members within the same subfamily. A total of 15 pairs of segmental duplication and 1 gene pair of tandem duplication were identified in tobacco based on the analysis of gene duplication events, while 10 pairs of segmental duplication in eggplant. It is speculated that segmental duplication events are the primary driver for the expansion of the U-box gene family in both tobacco and eggplant. The promoters of NtU-box and SmU-box genes contained cis-regulatory elements associated with cellular development, phytohormones, environment stress, and photoresponsive elements. Transcriptomic data analysis shows that the expression levels of the tobacco and eggplant U-box genes in different tissues and various abiotic stress conditions. Using cultivar Hongda of tobacco and cultivar Yanzhi of eggplant as materials, qRT-PCR analysis has revealed that 15 selected NtU-box genes and 8 SmU-box may play important roles in response to pathogen Ras invasion both in tobacco and eggplant.

Reclassification of Botryococcus braunii chemical races into separate species based on a comparative genomics analysis

The colonial green microalga Botryococcus braunii is well known for producing liquid hydrocarbons that can be utilized as biofuel feedstocks. B. braunii is taxonomically classified as a single species made up of three chemical races, A, B, and L, that are mainly distinguished by the hydrocarbons produced. We previously reported a B race draft nuclear genome, and here we report the draft nuclear genomes for the A and L races. A comparative genomic study of the three B. braunii races and 14 other algal species within Chlorophyta revealed significant differences in the genomes of each race of B. braunii. Phylogenomically, there was a clear divergence of the three races with the A race diverging earlier than both the B and L races, and the B and L races diverging from a later common ancestor not shared by the A race. DNA repeat content analysis suggested the B race had more repeat content than the A or L races. Orthogroup analysis revealed the B. braunii races displayed more gene orthogroup diversity than three closely related Chlamydomonas species, with nearly 24-36% of all genes in each B. braunii race being specific to each race. This analysis suggests the three races are distinct species based on sufficient differences in their respective genomes. We propose reclassification of the three chemical races to the following species names: Botryococcus alkenealis (A race), Botryococcus braunii (B race), and Botryococcus lycopadienor (L race).

Prioritizing Maize Metabolic Gene Regulators through Multi-Omic Network Integration

Elucidating gene regulatory networks is a major area of study within plant systems biology. Phenotypic traits are intricately linked to specific gene expression profiles. These expression patterns arise primarily from regulatory connections between sets of transcription factors (TFs) and their target genes. Here, we integrated 46 co-expression networks, 283 protein-DNA interaction (PDI) assays, and 16 million SNPs used to identify expression quantitative trait loci (eQTL) to construct TF-target networks. In total, we analyzed ∼4.6M interactions to generate four distinct types of TF-target networks: co-expression, PDI, trans -eQTL, and cis -eQTL combined with PDIs. To functionally annotate TFs based on their target genes, we implemented three different network integration strategies. We evaluated the effectiveness of each strategy through TF loss-of function mutant inspection and random network analyses. The multi-network integration allowed us to identify transcriptional regulators of several biological processes. Using the topological properties of the fully integrated network, we identified potential functionally redundant TF paralogs. Our findings retrieved functions previously documented for numerous TFs and revealed novel functions that are crucial for informing the design of future experiments. The approach here-described lays the foundation for the integration of multi-omic datasets in maize and other plant systems.

GRAPHICAL ABSTRACT

Mining for QTL controlling maize low-phosphorus response genes combined with deep resequencing of RIL parental genomes and in silico GWAS analysis

KEY MESSAGE

Extensive and comprehensive phenotypic data from a maize RIL population under both low- and normal-Pi treatments were used to conduct QTL mapping. Additionally, we integrated parental resequencing data from the RIL population, GWAS results, and transcriptome data to identify candidate genes associated with low-Pi stress in maize. Phosphorus (Pi) is one of the essential nutrients that greatly affect the maize yield. However, the genes underlying the QTL controlling maize low-Pi response remain largely unknown. In this study, a total of 38 traits at both seedling and maturity stages were evaluated under low- and normal-Pi conditions using a RIL population constructed from X178 (tolerant) and 9782 (sensitive), and most traits varied significantly between low- and normal-Pi treatments. Twenty-nine QTLs specific to low-Pi conditions were identified after excluding those with common intervals under both low- and normal-Pi conditions. Furthermore, 45 additional QTLs were identified based on the index value ((Trait_under_LowPi-Trait_under_NormalPi)/Trait_under_NormalPi) of each trait. These 74 QTLs collectively were classified as Pi-dependent QTLs. Additionally, 39 Pi-dependent QTLs were clustered in nine HotspotQTLs. The Pi-dependent QTL interval contained 19,613 unique genes, 6,999 of which exhibited sequence differences with non-synonymous mutation sites between X178 and 9782. Combined with in silico GWAS results, 277 consistent candidate genes were identified, with 124 genes located within the HotspotQTL intervals. The transcriptome analysis revealed that 21 genes, including the Pi transporter ZmPT7 and the strigolactones pathway-related gene ZmPDR1, exhibited consistent low-Pi stress response patterns across various maize inbred lines or tissues. It is noteworthy that ZmPDR1 in maize roots can be sharply up-regulated by low-Pi stress, suggesting its potential importance as a candidate gene for responding to low-Pi stress through the strigolactones pathway.

Quantitative trait locus analysis of gray leaf spot resistance in the maize IBM Syn10 DH population

KEY MESSAGE

The exploration and dissection of a set of QTLs and candidate genes for gray leaf spot disease resistance using two fully assembled parental genomes may help expedite maize resistance breeding. The fungal disease of maize known as gray leaf spot (GLS), caused by Cercospora zeae-maydis and Cercospora zeina, is a significant concern in China, Southern Africa, and the USA. Resistance to GLS is governed by multiple genes with an additive effect and is influenced by both genotype and environment. The most effective way to reduce the cost of production is to develop resistant hybrids. In this study, we utilized the IBM Syn 10 Doubled Haploid (IBM Syn10 DH) population to identify quantitative trait loci (QTLs) associated with resistance to gray leaf spot (GLS) in multiple locations. Analysis of seven distinct environments revealed a total of 58 QTLs, 49 of which formed 12 discrete clusters distributed across chromosomes 1, 2, 3, 4, 8 and 10. By comparing these findings with published research, we identified colocalized QTLs or GWAS loci within eleven clustering intervals. By integrating transcriptome data with genomic structural variations between parental individuals, we identified a total of 110 genes that exhibit both robust disparities in gene expression and structural alterations. Further analysis revealed 19 potential candidate genes encoding conserved resistance gene domains, including putative leucine-rich repeat receptors, NLP transcription factors, fucosyltransferases, and putative xyloglucan galactosyltransferases. Our results provide a valuable resource and linked loci for GLS marker resistance selection breeding in maize.

Genome-Wide Investigation of the CRF Gene Family in Maize and Functional Analysis of ZmCRF9 in Response to Multiple Abiotic Stresses

The cytokinin response factors (CRFs) are pivotal players in regulating plant growth, development, and responses to diverse stresses. Despite their significance, comprehensive information on CRF genes in the primary food crop, maize, remains scarce. In this study, a genome-wide analysis of CRF genes in maize was conducted, resulting in the identification of 12 members. Subsequently, we assessed the chromosomal locations, gene duplication events, evolutionary relationships, conserved motifs, and gene structures of all ZmCRF members. Analysis of ZmCRF promoter regions indicated the presence of cis-regulatory elements associated with plant growth regulation, hormone response, and various abiotic stress responses. The expression patterns of maize CRF genes, presented in heatmaps, exhibited distinctive patterns of tissue specificity and responsiveness to multiple abiotic stresses. qRT-PCR experiments were conducted on six selected genes and confirmed the involvement of ZmCRF genes in the plant's adaptive responses to diverse environmental challenges. In addition, ZmCRF9 was demonstrated to positively regulate cold and salt tolerance. Ultimately, we explored the putative interaction partners of ZmCRF proteins. In summary, this systematic overview and deep investigation of ZmCRF9 provides a solid foundation for further exploration into how these genes contribute to the complex interplay of plant growth, development, and responses to stress.

Comprehensive analysis of the Xya riparia genome uncovers the dominance of DNA transposons, LTR/Gypsy elements, and their evolutionary dynamics

Transposable elements (TEs) are DNA sequences that can move or replicate within a genome, and their study has become increasingly important in understanding genome evolution and function. The Tridactylidae family, including Xya riparia (pygmy mole cricket), harbors a variety of transposable elements (TEs) that have been insufficiently investigated. Further research is required to fully understand their diversity and evolutionary characteristics. Hence, we conducted a comprehensive repeatome analysis of X. riparia species using the chromosome-level assembled genome. The study aimed to comprehensively analyze the abundance, distribution, and age of transposable elements (TEs) in the genome. The results indicated that the genome was 1.67 Gb, with 731.63 Mb of repetitive sequences, comprising 27% of Class II (443.25 Mb), 16% of Class I (268.45 Mb), and 1% of unknown TEs (19.92 Mb). The study found that DNA transposons dominate the genome, accounting for approximately 60% of the total repeat size, with retrotransposons and unknown elements accounting for 37% and 3% of the genome, respectively. The members of the Gypsy superfamily were the most abundant amongst retrotransposons, accounting for 63% of them. The transposable superfamilies (LTR/Gypsy, DNA/nMITE, DNA/hAT, and DNA/Helitron) collectively constituted almost 70% of the total repeat size of all six chromosomes. The study further unveiled a significant linear correlation (Pearson correlation: r = 0.99, p-value = 0.00003) between the size of the chromosomes and the repetitive sequences. The average age of DNA transposon and retrotransposon insertions ranges from 25 My (million years) to 5 My. The satellitome analysis discovered 13 satellite DNA families that comprise about 0.15% of the entire genome. In addition, the transcriptional analysis of TEs found that DNA transposons were more transcriptionally active than retrotransposons. Overall, the study suggests that the genome of X. riparia is complex, characterized by a substantial portion of repetitive elements. These findings not only enhance our understanding of TE evolution within the Tridactylidae family but also provide a foundation for future investigations into the genomic intricacies of related species.

Identification of Highly Repetitive Enhancers with Long-range Regulation Potential in Barley via STARR-seq

Enhancers are DNA sequences that can strengthen transcription initiation. However, the global identification of plant enhancers is complicated due to uncertainty in the distance and orientation of enhancers, especially in species with large genomes. In this study, we performed self-transcribing active regulatory region sequencing (STARR-seq) for the first time to identify enhancers across the barley genome. A total of 7323 enhancers were successfully identified, and among 45 randomly selected enhancers, over 75% were effective as validated by a dual-luciferase reporter assay system in the lower epidermis of tobacco leaves. Interestingly, up to 53.5% of the barley enhancers were repetitive sequences, especially transposable elements (TEs), thus reinforcing the vital role of repetitive enhancers in gene expression. Both the common active mark H3K4me3 and repressive mark H3K27me3 were abundant among the barley STARR-seq enhancers. In addition, the functional range of barley STARR-seq enhancers seemed much broader than that of rice or maize and extended to ±100 kb of the gene body, and this finding was consistent with the high expression levels of genes in the genome. This study specifically depicts the unique features of barley enhancers and provides available barley enhancers for further utilization.

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.

Genetic response of a perennial grass to warm and wet environments interacts and is associated with trait means as well as plasticity

The potential for rapid evolution is an important mechanism allowing species to adapt to changing climatic conditions. Although such potential has been largely studied in various short-lived organisms, to what extent we can observe similar patterns in long-lived plant species, which often dominate natural systems, is largely unexplored. We explored the potential for rapid evolution in Festuca rubra, a long-lived grass with extensive clonal growth dominating in alpine grasslands. We used a field sowing experiment simulating expected climate change in our model region. Specifically, we exposed seeds from five independent seed sources to novel climatic conditions by shifting them along a natural climatic grid and explored the genetic profiles of established seedlings after 3 years. Data on genetic profiles of plants selected under different novel conditions indicate that different climate shifts select significantly different pools of genotypes from common seed pools. Increasing soil moisture was more important than increasing temperature or the interaction of the two climatic factors in selecting pressure. This can indicate negative genetic interaction in response to the combined effects or that the effects of different climates are interactive rather than additive. The selected alleles were found in genomic regions, likely affecting the function of specific genes or their expression. Many of these were also linked to morphological traits (mainly to trait plasticity), suggesting these changes may have a consequence on plant performance. Overall, these data indicate that even long-lived plant species may experience strong selection by climate, and their populations thus have the potential to rapidly adapt to these novel conditions.

Natural polymorphisms in ZMET2 encoding a DNA methyltransferase modulate the number of husk layers in maize

DNA methylation affects agronomic traits and the environmental adaptability of crops, but the natural polymorphisms in DNA methylation-related genes and their contributions to phenotypic variation in maize (Zea mays) remain elusive. Here, we show that a polymorphic 10-bp insertion/deletion variant in the 3'UTR of Zea methyltransferase2 (ZMET2) alters its transcript level and accounts for variation in the number of maize husk layers. ZMET2 encodes a chromomethylase and is required for maintaining genome-wide DNA methylation in the CHG sequence context. Disruption of ZMET2 increased the number of husk layers and resulted in thousands of differentially methylated regions, a proportion of which were also distinguishable in natural ZMET2 alleles. Population genetic analyses indicated that ZMET2 was a target of selection and might play a role in the spread of maize from tropical to temperate regions. Our results provide important insights into the natural variation of ZMET2 that confers both global and locus-specific effects on DNA methylation, which contribute to phenotypic diversity in maize.

Genome-wide association analysis of four yield-related traits using a maize (Zea mays L.) F1 population

Increasing the yield of maize F1 hybrid is one of the most important target for breeders. However, as a result of the genetic complexity and extremely low heritability, it is very difficult to directly dissect the genetic basis and molecular mechanisms of yield, and reports on genetic analysis of F1 hybrid yield are rare. Taking F1 hybrid as the research object and dividing the yield into different affect factors, this approach may be the best strategy for clarifying the genetic mechanism of yield. Therefore, in this study, a maize F1 population consisting of 300 hybrids with 17,652 single nucleotide polymorphisms (SNPs) markers was used for genome-wide association study (GWAS) to filtrate candidate genes associated with the four yield-related traits, i.e., kernel row number (KRN), kernel number per row (KNPR), ear tip-barrenness (ETB), and hundred kernel weight (HKW). Combined with the results of previous studies and functional annotation information of candidate genes, a total of six candidate genes were identified as being associated with the four traits, which were involved in plant growth and development, protein synthesis response, phytohormone biosynthesis and signal transduction. Our results improve the understanding of the genetic basis of the four yield-related traits and may be provide a new strategy for the genetic basis of maize yield.

Genome-wide analysis of AP2/ERF transcription factors that regulate fruit development of Chinese prickly ash

BACKGROUND

AP2/ERF is a large family of plant transcription factor proteins that play essential roles in signal transduction, plant growth and development, and responses to various stresses. The AP2/ERF family has been identified and verified by functional analysis in various plants, but so far there has been no comprehensive study of these factors in Chinese prickly ash. Phylogenetic, motif, and functional analyses combined with transcriptome analysis of Chinese prickly ash fruits at different developmental stages (30, 60, and 90 days after anthesis) were conducted in this study.

RESULTS

The analysis identified 146 ZbAP2/ERF genes that could be classified into 15 subgroups. The motif analysis revealed the presence of different motifs or elements in each group that may explain the functional differences between the groups. ZbERF13.2, ZbRAP2-12, and ZbERF2.1 showed high levels of expression in the early stages of fruit development. ZbRAP2-4, and ZbERF3.1 were significantly expressed at the fruit coloring stage (R2 and G2). ZbERF16 were significantly expressed at fruit ripening and expression level increased as the fruit continued to develop. Relative gene expression levels of 6 representative ZbAP2/ERFs assessed by RT-qPCR agreed with transcriptome analysis results.

CONCLUSIONS

These genes identified by screening can be used as candidate genes that affect fruit development. The results of the analysis can help guide future genetic improvement of Chinese prickly ash and enrich our understanding of AP2/ERF transcription factors and their regulatory functions in plants.

Metabolic and Functional Interactions of H2S and Sucrose in Maize Thermotolerance through Redox Homeodynamics

Hydrogen sulfide (H2S) is a novel gasotransmitter. Sucrose (SUC) is a source of cellular energy and a signaling molecule. Maize is the third most common food crop worldwide. However, the interaction of H2S and SUC in maize thermotolerance is not widely known. In this study, using maize seedlings as materials, the metabolic and functional interactions of H2S and SUC in maize thermotolerance were investigated. The data show that under heat stress, the survival rate and tissue viability were increased by exogenous SUC, while the malondialdehyde content and electrolyte leakage were reduced by SUC, indicating SUC could increase maize thermotolerance. Also, SUC-promoted thermotolerance was enhanced by H2S, while separately weakened by an inhibitor (propargylglycine) and a scavenger (hypotaurine) of H2S and a SUC-transport inhibitor (N-ethylmaleimide), suggesting the interaction of H2S and SUC in the development of maize thermotolerance. To establish the underlying mechanism of H2S-SUC interaction-promoted thermotolerance, redox parameters in mesocotyls of maize seedlings were measured before and after heat stress. The data indicate that the activity and gene expression of H2S-metabolizing enzymes were up-regulated by SUC, whereas H2S had no significant effect on the activity and gene expression of SUC-metabolizing enzymes. In addition, the activity and gene expression of catalase, glutathione reductase, ascorbate peroxidase, peroxidase, dehydroascorbate reductase, monodehydroascorbate reductase, and superoxide dismutase were reinforced by H2S, SUC, and their combination under non-heat and heat conditions to varying degrees. Similarly, the content of ascorbic acid, flavone, carotenoid, and polyphenol was increased by H2S, SUC, and their combination, whereas the production of superoxide radicals and the hydrogen peroxide level were impaired by these treatments to different extents. These results imply that the metabolic and functional interactions of H2S and sucrose signaling exist in the formation of maize thermotolerance through redox homeodynamics. This finding lays the theoretical basis for developing climate-resistant maize crops and improving food security.

Conserved functions of chromatin regulators in basal Archaeplastida

Chromatin is a dynamic network that regulates genome organization and gene expression. Different types of chromatin regulators are highly conserved among Archaeplastida, including unicellular algae, while some chromatin genes are only present in land plant genomes. Here, we review recent advances in understanding the function of conserved chromatin factors in basal land plants and algae. We focus on the role of Polycomb-group genes which mediate H3K27me3-based silencing and play a role in balancing gene dosage and regulating haploid-to-diploid transitions by tissue-specific repression of the transcription factors KNOX and BELL in many representatives of the green lineage. Moreover, H3K27me3 predominantly occupies repetitive elements which can lead to their silencing in a unicellular alga and basal land plants, while it covers mostly protein-coding genes in higher land plants. In addition, we discuss the role of nuclear matrix constituent proteins as putative functional lamin analogs that are highly conserved among land plants and might have an ancestral function in stress response regulation. In summary, our review highlights the importance of studying chromatin regulation in a wide range of organisms in the Archaeplastida.

Dissecting the genome sequence of a clinical isolated Cunninghamella bertholletiae Z2 strain with rich cytochrome P450 enzymes (Article)

Mucormycosis is receiving much more attention because of its high morbidity and extremely high mortality rate in immunosuppressed populations. In this study, we isolated a Cunnignhamella bertholletiae Z2 strain from a skin lesion of a 14 year, 9 months old girl with acute lymphoblastic leukemia who die of infection from the Z2 strain. Genome sequencing was performed after isolation and amplification of the Z2 strain to reveal potential virulence factors and pathogenic mechanisms. The results showed that the genome size of the Z2 strain is 30.9 Mb with 9213 genes. Mucoral specific virulence factor genes found are ARF, CalN, and CoTH, while no gliotoxin biosynthesis gene cluster was found, which is a known virulence factor in Aspergillus fumigatus adapted to the environment. The Z2 strain was found to have 69 cytochrome P450 enzymes, which are potential drug resistant targets. Sensitivity testing of Z2 showed it was only inhibited by amphotericin B and posaconazole. Detailed genomic information of the C. bertholletiae Z2 strain may provide useful data for treatment.

How diverse a monocentric chromosome can be? Repeatome and centromeric organization of Juncus effusus (Juncaceae)

Juncus is the largest genus of Juncaceae and was considered holocentric for a long time. Recent findings, however, indicated that 11 species from different clades of the genus have monocentric chromosomes. Thus, the Juncus centromere organization and evolution need to be reassessed. We aimed to investigate the major repetitive DNA sequences of two accessions of Juncus effusus and its centromeric structure by employing whole-genome analyses, fluorescent in situ hybridization, CENH3 immunodetection, and chromatin immunoprecipitation sequencing. We showed that the repetitive fraction of the small J. effusus genome (~270 Mbp/1C) is mainly composed of Class I and Class II transposable elements (TEs) and satellite DNAs. Three identified satellite DNA families were mainly (peri)centromeric, with two being associated with the centromeric protein CENH3, but not strictly centromeric. Two types of centromere organization were discerned in J. effusus: type 1 was characterized by a single CENH3 domain enriched with JefSAT1-155 or JefSAT2-180, whereas type 2 showed multiple CENH3 domains interrupted by other satellites, TEs or genes. Furthermore, while type 1 centromeres showed a higher degree of satellite identity along the array, type 2 centromeres had less homogenized arrays along the multiple CENH3 domains per chromosome. Although the analyses confirmed the monocentric organization of J. effusus chromosomes, our data indicate a more dynamic arrangement of J. effusus centromeres than observed for other plant species, suggesting it may constitute a transient state between mono- and holocentricity.

A chromosome-level genome of mango exclusively from long-read sequence data

Improvements in long-read sequencing techniques have greatly accelerated plant genome sequencing. Current de novo assemblies are routinely achieved by assembling long-read sequence data into contigs that are assembled to chromosome level by chromatin conformation capture. We report here a chromosome-level mango genome using only PacBio high-fidelity (HiFi) long reads. HiFi reads at high coverage (204x) resulted in the assembly of 17 chromosomes, each as a single contig with telomeres at both ends. The remaining three chromosomes were represented each by two contigs, with telomeres at one end and ribosomal repeats at the other end. Analyzing contig ends allowed them to be paired and linked to generate the remaining three complete chromosomes, telomere-to-telomere but with ribosomal repeats of uncertain length. The assembled genome was 365 Mb with 100% completeness as assessed by Benchmarking Universal Single-Copy Orthologs analysis. The haplotypes assembled demonstrated extensive structural differences. This approach using very high genome coverage may be useful for assembling high-quality genomes for many other plants.

Whole-genome sequencing of allotetraploid bermudagrass reveals the origin of Cynodon and candidate genes for salt tolerance

Bermudagrass (Cynodon dactylon) is a globally distributed, extensively used warm-season turf and forage grass with high tolerance to salinity and drought stress in alkaline environments. However, the origin of the species and genetic mechanisms for salinity tolerance in the species are basically unknown. Accordingly, we set out to study evolution divergence events in the Cynodon genome and to identify genes for salinity tolerance. We developed a 604.0 Mb chromosome-level polyploid genome sequence for bermudagrass 'A12359' (n = 18). The C. dactylon genome comprises 2 complete sets of homoeologous chromosomes, each with approximately 30 000 genes, and most genes are conserved as syntenic pairs. Phylogenetic study showed that the initial Cynodon species diverged from Oropetium thomaeum approximately 19.7-25.4 million years ago (Mya), the A and B subgenomes of C. dactylon diverged approximately 6.3-9.1 Mya, and the bermudagrass polyploidization event occurred 1.5 Mya on the African continent. Moreover, we identified 82 candidate genes associated with seven agronomic traits using a genome-wide association study, and three single-nucleotide polymorphisms were strongly associated with three salt resistance genes: RAP2-2, CNG channels, and F14D7.1. These genes may be associated with enhanced bermudagrass salt tolerance. These bermudagrass genomic resources, when integrated, may provide fundamental insights into evolution of diploid and tetraploid genomes and enhance the efficacy of comparative genomics in studying salt tolerance in Cynodon.

Transcriptome-wide expression landscape and starch synthesis pathway co-expression network in sorghum

The gene expression landscape across different tissues and developmental stages reflects their biological functions and evolutionary patterns. Integrative and comprehensive analyses of all transcriptomic data in an organism are instrumental to obtaining a comprehensive picture of gene expression landscape. Such studies are still very limited in sorghum, which limits the discovery of the genetic basis underlying complex agricultural traits in sorghum. We characterized the genome-wide expression landscape for sorghum using 873 RNA-sequencing (RNA-seq) datasets representing 19 tissues. Our integrative analysis of these RNA-seq data provides the most comprehensive transcriptomic atlas for sorghum, which will be valuable for the sorghum research community for functional characterizations of sorghum genes. Based on the transcriptome atlas, we identified 595 housekeeping genes (HKGs) and 2080 tissue-specific expression genes (TEGs) for the 19 tissues. We identified different gene features between HKGs and TEGs, and we found that HKGs have experienced stronger selective constraints than TEGs. Furthermore, we built a transcriptome-wide co-expression network (TW-CEN) comprising 35 modules with each module enriched in specific Gene Ontology terms. High-connectivity genes in TW-CEN tend to express at high levels while undergoing intensive selective pressure. We also built global and seed-preferential co-expression networks of starch synthesis pathways, which indicated that photosynthesis and microtubule-based movement play important roles in starch synthesis. The global transcriptome atlas of sorghum generated by this study provides an important functional genomics resource for trait discovery and insight into starch synthesis regulation in sorghum.

Establishing an optimized ATAC-seq protocol for the maize

The advent of next-generation sequencing in crop improvement offers unprecedented insights into the chromatin landscape closely linked to gene activity governing key traits in plant development and adaptation. Particularly in maize, its dynamic chromatin structure is found to collaborate with massive transcriptional variations across tissues and developmental stages, implying intricate regulatory mechanisms, which highlights the importance of integrating chromatin information into breeding strategies for precise gene controls. The depiction of maize chromatin architecture using Assay for Transposase Accessible Chromatin with high-throughput sequencing (ATAC-seq) provides great opportunities to investigate cis-regulatory elements, which is crucial for crop improvement. In this context, we developed an easy-to-implement ATAC-seq protocol for maize with fewer nuclei and simple equipment. We demonstrate a streamlined ATAC-seq protocol with four key steps for maize in which nuclei purification can be achieved without cell sorting and using only a standard bench-top centrifuge. Our protocol, coupled with the bioinformatic analysis, including validation by read length periodicity, key metrics, and correlation with transcript abundance, provides a precise and efficient assessment of the maize chromatin landscape. Beyond its application to maize, our testing design holds the potential to be applied to other crops or other tissues, especially for those with limited size and amount, establishing a robust foundation for chromatin structure studies in diverse crop species.

The quantification of southern corn leaf blight disease using deep UV fluorescence spectroscopy and autoencoder anomaly detection techniques

Southern leaf blight (SLB) is a foliar disease caused by the fungus Cochliobolus heterostrophus infecting maize plants in humid, warm weather conditions. SLB causes production losses to corn producers in different regions of the world such as Latin America, Europe, India, and Africa. In this paper, we demonstrate a non-destructive method to quantify the signs of fungal infection in SLB-infected corn plants using a deep UV (DUV) fluorescence spectrometer, with a 248.6 nm excitation wavelength, to acquire the emission spectra of healthy and SLB-infected corn leaves. Fluorescence emission spectra of healthy and diseased leaves were used to train an Autoencoder (AE) anomaly detection algorithm-an unsupervised machine learning model-to quantify the phenotype associated with SLB-infected leaves. For all samples, the signature of corn leaves consisted of two prominent peaks around 450 nm and 325 nm. However, SLB-infected leaves showed a higher response at 325 nm compared to healthy leaves, which was correlated to the presence of C. heterostrophus based on disease severity ratings from Visual Scores (VS). Specifically, we observed a linear inverse relationship between the AE error and the VS (R2 = 0.94 and RMSE = 0.935). With improved hardware, this method may enable improved quantification of SLB infection versus visual scoring based on e.g., fungal spore concentration per unit area and spatial localization.

Auxin efflux carrier PsPIN4 identified through genome-wide analysis as vital factor of petal abscission

PIN-FORMED (PIN) proteins, which function as efflux transporters, play many crucial roles in the polar transportation of auxin within plants. In this study, the exogenous applications of auxin IAA and TIBA were found to significantly prolong and shorten the florescence of tree peony (Paeonia suffruticosa Andr.) flowers. This finding suggests that auxin has some regulatory influence in petal senescence and abscission. Further analysis revealed a total of 8 PsPINs distributed across three chromosomes, which could be categorized into two classes based on phylogenetic and structural analysis. PsPIN1, PsPIN2a-b, and PsPIN4 were separated into the "long" PIN category, while PsPIN5, PsPIN6a-b, and PsPIN8 belonged to the "short" one. Additionally, the cis-regulatory elements of PsPIN promoters were associated with plant development, phytohormones, and environmental stress. These genes displayed tissue-specific expression, and phosphorylation sites were abundant throughout the protein family. Notably, PsPIN4 displayed distinct and elevated expression levels in roots, leaves, and flower organs. Expression patterns among the abscission zone (AZ) and adjacent areas during various flowering stages and IAA treatment indicate that PsPIN4 likely influences the initiation of peony petal abscission. The PsPIN4 protein was observed to be co-localized on both the plasma membrane and the cell nucleus. The ectopic expression of PsPIN4 reversed the premature flower organs abscission in the Atpin4 and significantly protracted florescence when introduced to Col Arabidopsis. Our findings established a strong basis for further investigation of PIN gene biological functions, particularly concerning intrinsic relationship between PIN-mediated auxin polar.

Enhanced pan-genomic resources at the maize genetics and genomics database

Pan-genomes, encompassing the entirety of genetic sequences found in a collection of genomes within a clade, are more useful than single reference genomes for studying species diversity. This is especially true for a species like Zea mays, which has a particularly diverse and complex genome. Presenting pan-genome data, analyses, and visualization is challenging, especially for a diverse species, but more so when pan-genomic data is linked to extensive gene model and gene data, including classical gene information, markers, insertions, expression and proteomic data, and protein structures as is the case at MaizeGDB. Here, we describe MaizeGDB's expansion to include the genic subset of the Zea pan-genome in a pan-gene data center featuring the maize genomes hosted at MaizeGDB, and the outgroup teosinte Zea genomes from the Pan-Andropoganeae project. The new data center offers a variety of browsing and visualization tools, including sequence alignment visualization, gene trees and other tools, to explore pan-genes in Zea that were calculated by the pipeline Pandagma. Combined, these data will help maize researchers study the complexity and diversity of Zea, and to use the comparative functions to validate pan-gene relationships for a selected gene model.

Are cereal grasses a single genetic system?

In 1993, a passionate and provocative call to arms urged cereal researchers to consider the taxon they study as a single genetic system and collaborate with each other. Since then, that group of scientists has seen their discipline blossom. In an attempt to understand what unity of genetic systems means and how the notion was borne out by later research, we survey the progress and prospects of cereal genomics: sequence assemblies, population-scale sequencing, resistance gene cloning and domestication genetics. Gene order may not be as extraordinarily well conserved in the grasses as once thought. Still, several recurring themes have emerged. The same ancestral molecular pathways defining plant architecture have been co-opted in the evolution of different cereal crops. Such genetic convergence as much as cross-fertilization of ideas between cereal geneticists has led to a rich harvest of genes that, it is hoped, will lead to improved varieties.

Helitrons: genomic parasites that generate developmental novelties

Helitrons, classified as DNA transposons, employ rolling-circle intermediates for transposition. Distinguishing themselves from other DNA transposons, they leave the original template element unaltered during transposition, which has led to their characterization as 'peel-and-paste elements'. Helitrons possess the ability to capture and mobilize host genome fragments, with enormous consequences for host genomes. This review discusses the current understanding of Helitrons, exploring their origins, transposition mechanism, and the extensive repercussions of their activity on genome structure and function. We also explore the evolutionary conflicts stemming from Helitron-transposed gene fragments and elucidate their domestication for regulating responses to environmental challenges. Looking ahead, further research in this evolving field promises to bring interesting discoveries on the role of Helitrons in shaping genomic landscapes.

Identification and characterization of ACR gene family in maize for salt stress tolerance

Background

Members of the ACR gene family are commonly involved in various physiological processes, including amino acid metabolism and stress responses. In recent decades, significant progress has been made in the study of ACR genes in plants. However, little is known about their characteristics and function in maize.

Methods

In this study, ACR genes were identified from the maize genome, and their molecular characteristics, gene structure, gene evolution, gene collinearity analysis, cis-acting elements were analyzed. qRT-PCR technology was used to verify the expression patterns of the ZmACR gene family in different tissues under salt stress. In addition, Ectopic expression technique of ZmACR5 in Arabidopsis thaliana was utilized to identify its role in response to salt stress.

Results

A total of 28 ZmACR genes were identified, and their molecular characteristics were extensively described. Two gene pairs arising from segmented replication events were detected in maize, and 18 collinear gene pairs were detected between maize and 3 other species. Through phylogenetic analysis, three subgroups were revealed, demonstrating distinct divergence between monocotyledonous and dicotyledonous plants. Analysis of ZmACR cis-acting elements revealed the optional involvement of ZmACR genes in light response, hormone response and stress resistance. Expression analysis of 8 ZmACR genes under salt treatment clearly revealed their role in the response to salt stress. Ectopic overexpression of ZmACR5 in Arabidopsis notably reduced salt tolerance compared to that of the wild type under salt treatment, suggesting that ZmACR5 has a negative role in the response to salt stress.

Conclusion

Taken together, these findings confirmed the involvement of ZmACR genes in regulating salt stress and contributed significantly to our understanding of the molecular function of ACR genes in maize, facilitating further research in this field.

In vitro functional analysis and in silico structural modelling of pathogen-secreted polyglycine hydrolases

Polyglycine hydrolases are fungal effectors composed of an N-domain with unique sequence and structure and a C-domain that resembles β-lactamases, with serine protease activity. These secreted fungal proteins cleave Gly-Gly bonds within a polyglycine sequence in corn ChitA chitinase. The polyglycine hydrolase N-domain (PND) function is unknown. In this manuscript we provide evidence that the PND does not directly participate in ChitA cleavage. In vitro analysis of site-directed mutants in conserved residues of the PND of polyglycine hydrolase Es-cmp did not specifically impair protease activity. Furthermore, in silico structural models of three ChitA-bound polyglycine hydrolases created by High Ambiguity Driven protein-protein DOCKing (HADDOCK) did not predict significant interactions between the PND and ChitA. Together these results suggest that the PND has another function. To determine what types of PND-containing proteins exist in nature we performed a computational analysis of Foldseek-identified PND-containing proteins. The analysis showed that proteins with PNDs are present throughout biology as either single domain proteins or fused to accessory domains that are diverse but are usually proteases or kinases.

Genomic analysis reveals phylogeny of Zygophyllales and mechanism for water retention of a succulent xerophyte

Revealing the genetic basis for stress-resistant traits in extremophile plants will yield important information for crop improvement. Zygophyllum xanthoxylum, an extant species of the ancient Mediterranean, is a succulent xerophyte that can maintain a favorable water status under desert habitats; however, the genetic basis of this adaptive trait is poorly understood. Furthermore, the phylogenetic position of Zygophyllales, to which Z. xanthoxylum belongs, remains controversial. In this study, we sequenced and assembled the chromosome-level genome of Z. xanthoxylum. Phylogenetic analysis showed that Zygophyllales and Myrtales form a separated taxon as a sister to the clade comprising fabids and malvids, clarifying the phylogenetic position of Zygophyllales at whole-genome scale. Analysis of genomic and transcriptomic data revealed multiple critical mechanisms underlying the efficient osmotic adjustment using Na+ and K+ as "cheap" osmolytes that Z. xanthoxylum has evolved through the expansion and synchronized expression of genes encoding key transporters/channels and their regulators involved in Na+/K+ uptake, transport, and compartmentation. It is worth noting that ZxCNGC1;1 (cyclic nucleotide-gated channels) and ZxCNGC1;2 constituted a previously undiscovered energy-saving pathway for Na+ uptake. Meanwhile, the core genes involved in biosynthesis of cuticular wax also featured an expansion and upregulated expression, contributing to the water retention capacity of Z. xanthoxylum under desert environments. Overall, these findings boost the understanding of evolutionary relationships of eudicots, illustrate the unique water retention mechanism in the succulent xerophyte that is distinct from glycophyte, and thus provide valuable genetic resources for the improvement of stress tolerance in crops and insights into the remediation of sodic lands.

Multiscale physiological responses to nitrogen supplementation of maize hybrids

Maize (Zea mays) production systems are heavily reliant on the provision of managed inputs such as fertilizers to maximize growth and yield. Hence, the effective use of nitrogen (N) fertilizer is crucial to minimize the associated financial and environmental costs, as well as maximize yield. However, how to effectively utilize N inputs for increased grain yields remains a substantial challenge for maize growers that requires a deeper understanding of the underlying physiological responses to N fertilizer application. We report a multiscale investigation of five field-grown maize hybrids under low or high N supplementation regimes that includes the quantification of phenolic and prenyl-lipid compounds, cellular ultrastructural features, and gene expression traits at three developmental stages of growth. Our results reveal that maize perceives the lack of supplemented N as a stress and, when provided with additional N, will prolong vegetative growth. However, the manifestation of the stress and responses to N supplementation are highly hybrid-specific. Eight genes were differentially expressed in leaves in response to N supplementation in all tested hybrids and at all developmental stages. These genes represent potential biomarkers of N status and include two isoforms of Thiamine Thiazole Synthase involved in vitamin B1 biosynthesis. Our results uncover a detailed view of the physiological responses of maize hybrids to N supplementation in field conditions that provides insight into the interactions between management practices and the genetic diversity within maize.

Chromosome-level genome assembly of the diploid oat species Avena longiglumis

Diploid wild oat Avena longiglumis has nutritional and adaptive traits which are valuable for common oat (A. sativa) breeding. The combination of Illumina, Nanopore and Hi-C data allowed us to assemble a high-quality chromosome-level genome of A. longiglumis (ALO), evidenced by contig N50 of 12.68 Mb with 99% BUSCO completeness for the assembly size of 3,960.97 Mb. A total of 40,845 protein-coding genes were annotated. The assembled genome was composed of 87.04% repetitive DNA sequences. Dotplots of the genome assembly (PI657387) with two published ALO genomes were compared to indicate the conservation of gene order and equal expansion of all syntenic blocks among three genome assemblies. Two recent whole-genome duplication events were characterized in genomes of diploid Avena species. These findings provide new knowledge for the genomic features of A. longiglumis, give information about the species diversity, and will accelerate the functional genomics and breeding studies in oat and related cereal crops.

Study of Dispersed Repeats in the Cyanidioschyzon merolae Genome

In this study, we applied the iterative procedure (IP) method to search for families of highly diverged dispersed repeats in the genome of Cyanidioschyzon merolae, which contains over 16 million bases. The algorithm included the construction of position weight matrices (PWMs) for repeat families and the identification of more dispersed repeats based on the PWMs using dynamic programming. The results showed that the C. merolae genome contained 20 repeat families comprising a total of 33,938 dispersed repeats, which is significantly more than has been previously found using other methods. The repeats varied in length from 108 to 600 bp (522.54 bp in average) and occupied more than 72% of the C. merolae genome, whereas previously identified repeats, including tandem repeats, have been shown to constitute only about 28%. The high genomic content of dispersed repeats and their location in the coding regions suggest a significant role in the regulation of the functional activity of the genome.