The limits of selection during maize domestication

Domestication and the evolution of crops: variable syndromes, complex genetic architectures, and ecological entanglements

Domestication can be considered a specialized mutualism in which a domesticator exerts control over the reproduction or propagation (fitness) of a domesticated species to gain resources or services. The evolution of crops by human-associated selection provides a powerful set of models to study recent evolutionary adaptations and their genetic bases. Moreover, the domestication and dispersal of crops such as rice, maize, and wheat during the Holocene transformed human social and political organization by serving as the key mechanism by which human societies fed themselves. Here we review major themes and identify emerging questions in three fundamental areas of crop domestication research: domestication phenotypes and syndromes, genetic architecture underlying crop evolution, and the ecology of domestication. Current insights on the domestication syndrome in crops largely come from research on cereal crops such as rice and maize, and recent work indicates distinct domestication phenotypes can arise from different domestication histories. While early studies on the genetics of domestication often identified single large-effect loci underlying major domestication traits, emerging evidence supports polygenic bases for many canonical traits such as shattering and plant architecture. Adaptation in human-constructed environments also influenced ecological traits in domesticates such as resource acquisition rates and interactions with other organisms such as root mycorrhizal fungi and pollinators. Understanding the ecological context of domestication will be key to developing resource-efficient crops and implementing more sustainable land management and cultivation practices.

Exploring the patterns of evolution: Core thoughts and focus on the saltational model

The Modern Synthesis, a pillar in biological thought, united Darwin's species origin concepts with Mendel's laws of character heredity, providing a comprehensive understanding of evolution within species. Highlighting phenotypic variation and natural selection, it elucidated the environment's role as a selective force, shaping populations over time. This framework integrated additional mechanisms, including genetic drift, random mutations, and gene flow, predicting their cumulative effects on microevolution and the emergence of new species. Beyond the Modern Synthesis, the Extended Evolutionary Synthesis expands perspectives by recognizing the role of developmental plasticity, non-genetic inheritance, and epigenetics. We suggest that these aspects coexist in the plant evolutionary process; in this context, we focus on the saltational model, emphasizing how saltation events, such as dichotomous saltation, chromosomal mutations, epigenetic phenomena, and polyploidy, contribute to rapid evolutionary changes. The saltational model proposes that certain evolutionary changes, such as the rise of new species, may result suddenly from single macromutations rather than from gradual changes in DNA sequences and allele frequencies within a species over time. These events, observed in domesticated and wild higher plants, provide well-defined mechanistic bases, revealing their profound impact on plant diversity and rapid evolutionary events. Notably, next-generation sequencing exposes the likely crucial role of allopolyploidy and autopolyploidy (saltational events) in generating new plant species, each characterized by distinct chromosomal complements. In conclusion, through this review, we offer a thorough exploration of the ongoing dissertation on the saltational model, elucidating its implications for our understanding of plant evolutionary processes and paving the way for continued research in this intriguing field.

Genetic basis controlling rice plant architecture and its modification for breeding

The shoot and root system architectures are fundamental for crop productivity. During the history of artificial selection of domestication and post-domestication breeding, the architecture of rice has significantly changed from its wild ancestor to fulfil requirements in agriculture. We review the recent studies on developmental biology in rice by focusing on components determining rice plant architecture; shoot meristems, leaves, tillers, stems, inflorescences and roots. We also highlight natural variations that affected these structures and were utilized in cultivars. Importantly, many core regulators identified from developmental mutants have been utilized in breeding as weak alleles moderately affecting these architectures. Given a surge of functional genomics and genome editing, the genetic mechanisms underlying the rice plant architecture discussed here will provide a theoretical basis to push breeding further forward not only in rice but also in other crops and their wild relatives.

Dominant complementation of biological pathways in maize hybrid lines is associated with heterosis

Allele-specific expressed genes (ASEGs) are widespread in maize hybrid lines and play important roles of complementation of biological pathways in heterosis. Heterosis (hybrid vigor) is an important phenomenon with both theoretical and practical value. However, our understanding of the genetic and molecular mechanisms behind heterosis is still limited. Here, we analyzed a comprehensive dataset of maize (Zea mays L.), including RNA-seq data from three hybrid-parent triplets (HPTs) and acetylated protein data from one HPT. The gene expression patterns exhibited extensive variation between the hybrids and their parents, and a substantial number of allele-specific expressed genes (ASEGs) were identified in the hybrids. Notably, ASEGs from different HPTs were significantly enriched in various conserved pathways. The parental alleles of ASEGs with fewer deleterious single-nucleotide polymorphisms were more likely to be expressed in hybrid lines than other parental alleles. ASEGs were mainly enriched in the functional gene ontology terms protein biosynthesis, photosynthesis, and metabolism. In addition, the ASEGs across the three HPTs were involved in key photosynthetic pathways and might enhance the photosynthetic efficiency of the hybrids. These findings suggest that ASEGs involved in complementary biological pathways in maize hybrids contribute to heterosis, shedding new light on the molecular mechanism of heterosis.

Wild potato ancestors as potential sources of resistance to the aphid Myzus persicae

BACKGROUND

Plant resistance to insects can be reduced by crop domestication which means their wild ancestors could provide novel sources of resistance. Thus, crossing wild ancestors with domesticated crops can potentially enhance their resistance against insects. However, a prerequisite for this is identification of sources of resistance. Here, we investigated the response of three wild potato (Solanum stoloniferum Schltdl.) accessions and cultivated potato (Solanum tuberosum) to aphid (Myzus persicae Sulzer) herbivory.

RESULTS

Results revealed that there was a significant reduction in aphid survival and reproduction on wild potato accessions (CGN18333, CGN22718, CGN23072) compared to cultivated (Desiree) potato plants. A similar trend was observed in olfactometer bioassay; the wild accessions had a repellent effect on adult aphids. In contrast, among the tested wild potato accessions, the parasitoid Diaeretiella rapae (M'Intosh) was significantly attracted to volatiles from CGN18333. Volatile analysis showed that wild accessions emitted significantly more volatiles compared to cultivated potato. Principal component analysis (PCA) of volatile data revealed that the volatile profiles of wild and cultivated potato are dissimilar. β-Bisabolene, (E)-β-farnesene, trans-α-bergamotene, d-limonene, (E,E)-4,8,12-trimethyl-1,3,7,11-tridecatetraene (TMTT), and p-cymen-7-ol were the main volatiles contributing to the emitted blends, suggesting possible involvement in the behavioural response of both M. persicae and D. rapae.

CONCLUSION

Our findings show that the tested wild accessions have the potential to be used to breed aphid-resistant potatoes. This opens new opportunities to reduce the aphid damage and to enhance the recruitment of natural enemies. © 2022 The Authors. Pest Management Science published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

An Efficient and Universal Protoplast Isolation Protocol Suitable for Transient Gene Expression Analysis and Single-Cell RNA Sequencing

The recent advent of single-cell RNA sequencing (scRNA-seq) has enabled access to the developmental landscape of a complex organ by monitoring the differentiation trajectory of every specialized cell type at the single-cell level. A main challenge in this endeavor is dissociating plant cells from the rigid cell walls and some species are recalcitrant to such cellular isolation. Here, we describe the establishment of a simple and efficient protocol for protoplast preparation in Chirita pumila, which includes two consecutive digestion processes with different enzymatic buffers. Using this protocol, we generated viable cell suspensions suitable for an array of expression analyses, including scRNA-seq. The universal application of this protocol was further tested by successfully isolating high-quality protoplasts from multiple organs (petals, fruits, tuberous roots, and gynophores) from representative species on the key branches of the angiosperm lineage. This work provides a robust method in plant science, overcoming barriers to isolating protoplasts in diverse plant species and opens a new avenue to study cell type specification, tissue function, and organ diversification in plants.

Fine mapping and candidate gene analysis of the up locus determining fruit orientation in pepper (Capsicum spp.)

The up locus determining fruit orientation was fine-mapped into a region with a physical length of ~169.51 kb on chromosome P12 in pepper. Capana12g000958, encoding a developmentally regulated G protein 2, was proposed as the strongest candidate via sequence comparison and expression analysis. Fruit orientation is an important horticultural and domesticated trait, which is controlled by a single semi-dominant gene (up) in pepper. However, the gene underlying up locus has not yet been identified. In this study, the previously detected major QTL UP12.1 was firstly verified using a backcross population (n = 225) stem from the cross of BB3 (C. annuum) and its wild relative Chiltepin (C. annuum var. glabriusculum) using BB3 as the recurrent parent. Then, a large BC₁F₂ population (n = 1827) was used for recombinant screening to delimit the up locus into an interval with ~ 169.51 kb in length. Sequence comparison and expression analysis suggested that Capana12g000958, encoding a developmentally regulated G protein 2, was the most likely candidate gene for the up locus. There is no difference within the coding sequences of Capana12g000958 between BB3 and Chiltepin, while a SNP in the upstream of Capana12g000958 showed a complete correlation with the fruit orientation among a panel of 40 diverse pepper inbred lines. These findings will form a basis for gene isolation and reveal of genetic mechanism underlying the fruit orientation domestication in pepper.

Recurrent Loss of abaA, a Master Regulator of Asexual Development in Filamentous Fungi, Correlates with Changes in Genomic and Morphological Traits

Gene regulatory networks (GRNs) drive developmental and cellular differentiation, and variation in their architectures gives rise to morphological diversity. Pioneering studies in Aspergillus fungi, coupled with subsequent work in other filamentous fungi, have shown that the GRN governed by the BrlA, AbaA, and WetA proteins controls the development of the asexual fruiting body or conidiophore. A specific aspect of conidiophore development is the production of phialides, conidiophore structures that are under the developmental control of AbaA and function to repetitively generate spores. Fungal genome sequencing has revealed that some filamentous fungi lack abaA, and also produce asexual structures that lack phialides, raising the hypothesis that abaA loss is functionally linked to diversity in asexual fruiting body morphology. To examine this hypothesis, we carried out an extensive search for the abaA gene across 241 genomes of species from the fungal subphylum Pezizomycotina. We found that abaA was independently lost in four lineages of Eurotiomycetes, including from all sequenced species within the order Onygenales, and that all four lineages that have lost abaA also lack the ability to form phialides. Genetic restoration of abaA from Aspergillus nidulans into Histoplasma capsulatum, a pathogenic species from the order Onygenales that lacks an endogenous copy of abaA, did not alter Histoplasma conidiation morphology but resulted in a marked increase in spore viability. We also discovered that species lacking abaA contain fewer AbaA binding motifs in the regulatory regions of orthologs of some AbaA target genes, suggesting that the asexual fruiting body GRN of organisms that have lost abaA has likely been rewired. Our results provide an illustration of how repeated losses of a key regulatory transcription factor have contributed to the diversity of an iconic fungal morphological trait.

Local adaptation contributes to gene expression divergence in maize

Gene expression links genotypes to phenotypes, so identifying genes whose expression is shaped by selection will be important for understanding the traits and processes underlying local adaptation. However, detecting local adaptation for gene expression will require distinguishing between divergence due to selection and divergence due to genetic drift. Here, we adapt a QST−FST framework to detect local adaptation for transcriptome-wide gene expression levels in a population of diverse maize genotypes. We compare the number and types of selected genes across a wide range of maize populations and tissues, as well as selection on cold-response genes, drought-response genes, and coexpression clusters. We identify a number of genes whose expression levels are consistent with local adaptation and show that genes involved in stress response show enrichment for selection. Due to its history of intense selective breeding and domestication, maize evolution has long been of interest to researchers, and our study provides insight into the genes and processes important for in local adaptation of maize.

Cis- and Trans-Regulatory Variations in the Domestication of the Chili Pepper Fruit

The process of domestication requires the rapid transformation of the wild morphology into the cultivated forms that humans select for. This process often takes place through changes in the regulation of genes, yet, there is no definite pattern on the role of cis- and trans-acting regulatory variations in the domestication of the fruit among crops. Using allele-specific expression and network analyses, we characterized the regulatory patterns and the inheritance of gene expression in wild and cultivated accessions of chili pepper, a crop with remarkable fruit morphological variation. We propose that gene expression differences associated to the cultivated form are best explained by cis-regulatory hubs acting through trans-regulatory cascades. We show that in cultivated chili, the expression of genes associated with fruit morphology is partially recessive with respect to those in the wild relative, consistent with the hybrid fruit phenotype. Decreased expression of fruit maturation and growth genes in cultivated chili suggest that selection for loss-of-function took place in its domestication. Trans-regulatory changes underlie the majority of the genes showing regulatory divergence and had larger effect sizes on gene expression than cis-regulatory variants. Network analysis of selected cis-regulated genes, including ARP9 and MED25, indicated their interaction with many transcription factors involved in organ growth and fruit ripening. Differentially expressed genes linked to cis-regulatory variants and their interactions with downstream trans-acting genes have the potential to drive the morphological differences observed between wild and cultivated fruits and provide an attractive mechanism of morphological transformation during the domestication of the chili pepper.

Designing future crops: challenges and strategies for sustainable agriculture

Crop production is facing unprecedented challenges. Despite the fact that the food supply has significantly increased over the past half-century, ~8.9 and 14.3% people are still suffering from hunger and malnutrition, respectively. Agricultural environments are continuously threatened by a booming world population, a shortage of arable land, and rapid changes in climate. To ensure food and ecosystem security, there is a need to design future crops for sustainable agriculture development by maximizing net production and minimalizing undesirable effects on the environment. The future crops design projects, recently launched by the National Natural Science Foundation of China and Chinese Academy of Sciences (CAS), aim to develop a roadmap for rapid design of customized future crops using cutting-edge technologies in the Breeding 4.0 era. In this perspective, we first introduce the background and missions of these projects. We then outline strategies to design future crops, such as improvement of current well-cultivated crops, de novo domestication of wild species and redomestication of current cultivated crops. We further discuss how these ambitious goals can be achieved by the recent development of new integrative omics tools, advanced genome-editing tools and synthetic biology approaches. Finally, we summarize related opportunities and challenges in these projects.

Tunability enhancement of gene regulatory motifs through competition for regulatory protein resources

Gene regulatory networks (GRNs) orchestrate the spatiotemporal levels of gene expression, thereby regulating various cellular functions ranging from embryonic development to tissue homeostasis. Some patterns called "motifs" recurrently appear in the GRNs. Owing to the prevalence of these motifs they have been subjected to much investigation, both in the context of understanding cellular decision making and engineering synthetic circuits. Mounting experimental evidence suggests that (1) the copy number of genes associated with these motifs varies, and (2) proteins produced from these genes bind to decoy binding sites on the genome as well as promoters driving the expression of other genes. Together, these two processes engender competition for protein resources within a cell. To unravel how competition for protein resources affects the dynamical properties of regulatory motifs, we propose a simple kinetic model that explicitly incorporates copy number variation (CNV) of genes and decoy binding of proteins. Using quasi-steady-state approximations, we theoretically investigate the transient and steady-state properties of three of the commonly found motifs: Autoregulation, toggle switch, and repressilator. While protein resource competition alters the timescales to reach the steady state for all these motifs, the dynamical properties of the toggle switch and repressilator are affected in multiple ways. For toggle switch, the basins of attraction of the known attractors are dramatically altered if one set of proteins binds to decoys more frequently than the other, an effect which gets suppressed as the copy number of the toggle switch is enhanced. For repressilators, protein sharing leads to an emergence of oscillation in regions of parameter space that were previously nonoscillatory. Intriguingly, both the amplitude and frequency of oscillation are altered in a nonlinear manner through the interplay of CNV and decoy binding. Overall, competition for protein resources within a cell provides an additional layer of regulation of gene regulatory motifs.

Identification of Ear Morphology Genes in Maize (Zea mays L.) Using Selective Sweeps and Association Mapping

The performance of maize hybrids largely depend on two parental inbred lines. Improving inbred lines using artificial selection is a key task in breeding programs. However, it is important to elucidate the effects of this selection on inbred lines. Altogether, 208 inbred lines from two maize heterosis groups, named Shaan A and Shaan B, were sequenced by the genotype-by-sequencing to detect genomic changes under selection pressures. In addition, we completed genome-wide association analysis in 121 inbred lines to identify candidate genes for ear morphology related traits. In a genome-wide selection scan, the inbred lines from Shaan A and Shaan B groups showed obvious population divergences and different selective signals distributed in 337 regions harboring 772 genes. Meanwhile, functional enrichment analysis showed those selected genes are mainly involved in regulating cell development. Interestingly, some ear morphology related traits showed significant differentiation between the inbred lines from the two heterosis groups. The genome-wide association analysis of ear morphology related traits showed that four associated genes were co-localized in the selected regions with high linkage disequilibrium. Our spatiotemporal pattern and gene interaction network results for the four genes further contribute to our understanding of the mechanisms behind ear and fruit length development. This study provides a novel insight into digging a candidate gene for complex traits using breeding materials. Our findings in relation to ear morphology will help accelerate future maize improvement.

Long-read bitter gourd (Momordica charantia) genome and the genomic architecture of nonclassic domestication

The genetic architecture of quantitative traits is determined by both Mendelian and polygenic factors, yet classic examples of plant domestication focused on selective sweep of newly mutated Mendelian genes. Here we report the chromosome-level genome assembly and the genomic investigation of a nonclassic domestication example, bitter gourd (Momordica charantia), an important Asian vegetable and medicinal plant of the family Cucurbitaceae. Population resequencing revealed the divergence between wild and South Asian cultivars about 6,000 y ago, followed by the separation of the Southeast Asian cultivars about 800 y ago, with the latter exhibiting more extreme trait divergence from wild progenitors and stronger signs of selection on fruit traits. Unlike some crops where the largest phenotypic changes and traces of selection happened between wild and cultivar groups, in bitter gourd large differences exist between two regional cultivar groups, likely reflecting the distinct consumer preferences in different countries. Despite breeding efforts toward increasing female flower proportion, a gynoecy locus exhibits complex patterns of balanced polymorphism among haplogroups, with potential signs of selective sweep within haplogroups likely reflecting artificial selection and introgression from cultivars back to wild accessions. Our study highlights the importance to investigate such nonclassic example of domestication showing signs of balancing selection and polygenic trait architecture in addition to classic selective sweep in Mendelian factors.

Comparative transcriptome analyses of fruit development among pears, peaches, and strawberries provide new insights into single sigmoid patterns

BACKGROUND

Pear fruit exhibit a single sigmoid pattern during development, while peach and strawberry fruits exhibit a double sigmoid pattern. However, little is known about the differences between these two patterns.

RESULTS

In this study, fruit weights were measured and paraffin sections were made from fruitlet to maturated pear, peach, and strawberry samples. Results revealed that both single and double sigmoid patterns resulted from cell expansion, but not cell division. Comparative transcriptome analyses were conducted among pear, peach, and strawberry fruits at five fruit enlargement stages. Comparing the genes involved in these intervals among peaches and strawberries, 836 genes were found to be associated with all three fruit enlargement stages in pears (Model I). Of these genes, 25 were located within the quantitative trait locus (QTL) regions related to fruit weight and 90 were involved in cell development. Moreover, 649 genes were associated with the middle enlargement stage, but not early or late enlargement in pears (Model II). Additionally, 22 genes were located within the QTL regions related to fruit weight and 63 were involved in cell development. Lastly, dual-luciferase assays revealed that the screened bHLH transcription factors induced the expression of cell expansion-related genes, suggesting that the two models explain the single sigmoid pattern.

CONCLUSIONS

Single sigmoid patterns are coordinately mediated by Models I and II, thus, a potential gene regulation network for the single sigmoid pattern was proposed. These results enhance our understanding of the molecular regulation of fruit size in Rosaceae.

Principles, Applications, and Biosafety of Plant Genome Editing Using CRISPR-Cas9

The terms genome engineering, genome editing, and gene editing, refer to modifications (insertions, deletions, substitutions) in the genome of a living organism. The most widely used approach to genome editing nowadays is based on Clustered Regularly Interspaced Short Palindromic Repeats and associated protein 9 (CRISPR-Cas9). In prokaryotes, CRISPR-Cas9 is an adaptive immune system that naturally protects cells from DNA virus infections. CRISPR-Cas9 has been modified to create a versatile genome editing technology that has a wide diversity of applications in medicine, agriculture, and basic studies of gene functions. CRISPR-Cas9 has been used in a growing number of monocot and dicot plant species to enhance yield, quality, and nutritional value, to introduce or enhance tolerance to biotic and abiotic stresses, among other applications. Although biosafety concerns remain, genome editing is a promising technology with potential to contribute to food production for the benefit of the growing human population. Here, we review the principles, current advances and applications of CRISPR-Cas9-based gene editing in crop improvement. We also address biosafety concerns and show that humans have been exposed to Cas9 protein homologues long before the use of CRISPR-Cas9 in genome editing.

Interaction Between Induced and Natural Variation at oil yellow1 Delays Reproductive Maturity in Maize

We previously demonstrated that maize (Zea mays) locus very oil yellow1 (vey1) encodes a putative cis-regulatory expression polymorphism at the magnesium chelatase subunit I gene (aka oil yellow1) that strongly modifies the chlorophyll content of the semi-dominant Oy1-N1989 mutants. The vey1 allele of Mo17 inbred line reduces chlorophyll content in the mutants leading to reduced photosynthetic output. Oy1-N1989 mutants in B73 reached reproductive maturity four days later than wild-type siblings. Enhancement of Oy1-N1989 by the Mo17 allele at the vey1 QTL delayed maturity further, resulting in detection of a flowering time QTL in two bi-parental mapping populations crossed to Oy1-N1989. The near isogenic lines of B73 harboring the vey1 allele from Mo17 delayed flowering of Oy1-N1989 mutants by twelve days. Just as previously observed for chlorophyll content, vey1 had no effect on reproductive maturity in the absence of the Oy1-N1989 allele. Loss of chlorophyll biosynthesis in Oy1-N1989 mutants and enhancement by vey1 reduced CO₂ assimilation. We attempted to separate the effects of photosynthesis on the induction of flowering from a possible impact of chlorophyll metabolites and retrograde signaling by manually reducing leaf area. Removal of leaves, independent of the Oy1-N1989 mutant, delayed flowering but surprisingly reduced chlorophyll contents of emerging leaves. Thus, defoliation did not completely separate the identity of the signal(s) that regulates flowering time from changes in chlorophyll content in the foliage. These findings illustrate the necessity to explore the linkage between metabolism and the mechanisms that connect it to flowering time regulation.

High-Resolution Mapping in Two RIL Populations Refines Major "QTL Hotspot" Regions for Seed Size and Shape in Soybean (Glycine max L.)

Seed size and shape are important traits determining yield and quality in soybean. However, the genetic mechanism and genes underlying these traits remain largely unexplored. In this regard, this study used two related recombinant inbred line (RIL) populations (ZY and K3N) evaluated in multiple environments to identify main and epistatic-effect quantitative trait loci (QTLs) for six seed size and shape traits in soybean. A total of 88 and 48 QTLs were detected through composite interval mapping (CIM) and mixed-model-based composite interval mapping (MCIM), respectively, and 15 QTLs were common among both methods; two of them were major (R2 > 10%) and novel QTLs (viz., qSW-1-1ZN and qSLT-20-1K3N). Additionally, 51 and 27 QTLs were identified for the first time through CIM and MCIM methods, respectively. Colocalization of QTLs occurred in four major QTL hotspots/clusters, viz., "QTL Hotspot A", "QTL Hotspot B", "QTL Hotspot C", and "QTL Hotspot D" located on Chr06, Chr10, Chr13, and Chr20, respectively. Based on gene annotation, gene ontology (GO) enrichment, and RNA-Seq analysis, 23 genes within four "QTL Hotspots" were predicted as possible candidates, regulating soybean seed size and shape. Network analyses demonstrated that 15 QTLs showed significant additive x environment (AE) effects, and 16 pairs of QTLs showing epistatic effects were also detected. However, except three epistatic QTLs, viz., qSL-13-3ZY, qSL-13-4ZY, and qSW-13-4ZY, all the remaining QTLs depicted no main effects. Hence, the present study is a detailed and comprehensive investigation uncovering the genetic basis of seed size and shape in soybeans. The use of a high-density map identified new genomic regions providing valuable information and could be the primary target for further fine mapping, candidate gene identification, and marker-assisted breeding (MAB).

Dynamic patterns of circular and linear RNAs in maize hybrid and parental lines

Hybrid vigor, also known as heterosis, has been widely utilized in agronomic production of maize (Zea mays L.) and other crops. However, the molecular mechanisms underlying heterosis are still not fully understood. To provide a more complete understanding of the transcriptomic dynamics associated with heterosis, we collected a comprehensive set of sequence data on linear mRNA transcripts and circular RNAs (circRNAs) from seedling leaves of two widely used maize inbred lines and their F₁ hybrid at the V4 growth stage. We detected over 25,000 expressed genes with more than 1200 circRNAs that showed dramatic and distinct variations in expression level across the three genotypes. Although most linear and circular transcripts exhibited additive expression in the hybrid, the expression of circRNAs was more likely to be nonadditive. Interestingly, the levels of linear transcripts and their corresponding circRNAs from the same loci showed a significant relationship and coordinated expression mode across all three genotypes. Notably, in the hybrid, allele-specific expression of linear transcripts was significantly associated with the expression of circRNAs from the same locus, suggesting potential regulatory cross talk between linear and circular transcripts. Our study provides a deeper understanding of dynamic variations for both the linear and circular transcriptome in a classical hybrid triplet of maize.

A population genomics appraisal suggests independent dispersals for bitter and sweet manioc in Brazilian Amazonia

Amazonia is a major world centre of plant domestication, but the genetics of domestication remains unclear for most Amazonian crops. Manioc (Manihot esculenta) is the most important staple food crop that originated in this region. Although manioc is relatively well-studied, little is known about the diversification of bitter and sweet landraces and how they were dispersed across Amazonia. We evaluated single nucleotide polymorphisms (SNPs) in wild and cultivated manioc to identify outlier SNPs putatively under selection and to assess the neutral genetic structure of landraces to make inferences about the evolution of the crop in Amazonia. Some outlier SNPs were in putative manioc genes possibly related to plant architecture, transcriptional regulation and responses to stress. The neutral SNPs revealed contrasting genetic structuring for bitter and sweet landraces. The outlier SNPs may be signatures of the genomic changes resulting from domestication, while the neutral genetic structure suggests independent dispersals for sweet and bitter manioc, possibly related to the earlier domestication and diversification of the former. Our results highlight the role of ancient peoples and current smallholders in the management and conservation of manioc genetic diversity, including putative genes and specific genetic resources with adaptive potential in the context of climate change.

Evolution of Plant Architecture in Oryza Driven by the PROG1 Locus

The genetic control of plant architecture in crops is critical for agriculture and understanding morphological evolution. This study showed that an open reading frame (ORF) of the rice domestication gene PROG1 appeared 3.4-3.9 million years ago (Mya). Subsequently, it acquired a novel protein-coding gene function in the genome of O. rufipogon (~0.3-0.4 Mya). This extremely young gene and its paralogous C2H2 genes located nearby define the prostrate architecture of O. rufipogon and, thus, are of adaptive significance for wild rice in swamp and water areas. However, selection for dense planting and high yield during rice domestication silenced the PROG1 gene and caused the loss of the RPAD locus containing functional C2H2 paralogs; hence, domesticated lines exhibit an erect plant architecture. Analysis of the stepwise origination process of PROG1 and its evolutionary genetics revealed that this zinc-finger coding gene may have rapidly evolved under positive selection and promoted the transition from non- or semi-prostrate growth to prostrate growth. A transgenic assay showed that PROG1 from O. rufipogon exerts a stronger function compared with PROG1 sequences from other Oryza species. However, the analysis of the expression levels of PROG1 in different Oryza species suggests that the transcriptional regulation of PROG1 has played an important role in its evolution. This study provides the first strong case showing how a fundamental morphological trait evolved in Oryza species driven by a gene locus.

Identifying loci with breeding potential across temperate and tropical adaptation via EigenGWAS and EnvGWAS

Understanding the genomic basis of adaptation in maize is important for gene discovery and the improvement of breeding germplasm, but much remains a mystery in spite of significant population genetics and archaeological research. Identifying the signals underpinning adaptation are challenging as adaptation often coincided with genetic drift, and the base genomic diversity of the species in massive. In this study, tGBS technology was used to genotype 1,143 diverse maize accessions including landraces collected from 20 countries and elite breeding lines of tropical lowland, highland, subtropical/midaltitude and temperate ecological zones. Based on 355,442 high‐quality single nucleotide polymorphisms, 13 genomic regions were detected as being under selection using the bottom‐up searching strategy, EigenGWAS. Of the 13 selection regions, 10 were first reported, two were associated with environmental parameters via EnvGWAS, and 146 genes were enriched. Combining large‐scale genomic and ecological data in this diverse maize panel, our study supports a polygenic adaptation model of maize and offers a framework to enhance our understanding of both the mechanistic basis and the evolutionary consequences of maize domestication and adaptation. The regions identified here are promising candidates for further, targeted exploration to identify beneficial alleles and haplotypes for deployment in maize breeding.

Posttranscriptional Regulation of RhBRC1 (Rosa hybrida BRANCHED1) in Response to Sugars is Mediated via its Own 3' Untranslated Region, with a Potential Role of RhPUF4 (Pumilio RNA-Binding Protein Family)

The shoot branching pattern is a determining phenotypic trait throughout plant development. During shoot branching, BRANCHED1 (BRC1) plays a master regulator role in bud outgrowth, and its transcript levels are regulated by various exogenous and endogenous factors. RhBRC1 (the homologous gene of BRC1 in Rosa hybrida) is a main branching regulator whose posttranscriptional regulation in response to sugar was investigated through its 3'UTR. Transformed Rosa calluses containing a construction composed of the CaMV35S promoter, the green fluorescent protein (GFP) reporter gene, and the 3'UTR of RhBRC1 (P35S:GFP::3'UTR_RhBRC1) were obtained and treated with various combinations of sugars and with sugar metabolism effectors. The results showed a major role of the 3'UTR of RhBRC1 in response to sugars, involving glycolysis/the tricarboxylic acid cycle (TCA) and the oxidative pentose phosphate pathway (OPPP). In Rosa vegetative buds, sequence analysis of the RhBRC1 3'UTR identified six binding motifs specific to the Pumilio/FBF RNA-binding protein family (PUF) and probably involved in posttranscriptional regulation. RhPUF4 was highly expressed in the buds of decapitated plants and in response to sugar availability in in-vitro-cultured buds. RhPUF4 was found to be close to AtPUM2, which encodes an Arabidopsis PUF protein. In addition, sugar-dependent upregulation of RhPUF4 was also found in Rosa calluses. RhPUF4 expression was especially dependent on the OPPP, supporting its role in OPPP-dependent posttranscriptional regulation of RhBRC1. These findings indicate that the 3'UTR sequence could be an important target in the molecular regulatory network of RhBRC1 and pave the way for investigating new aspects of RhBRC1 regulation.

Population genomic analysis of mango (Mangifera indica) suggests a complex history of domestication

Humans have domesticated diverse species from across the plant kingdom, yet much of our foundational knowledge of domestication has come from studies investigating relatively few of the most important annual food crops. Here, we examine the impacts of domestication on genetic diversity in a tropical perennial fruit species, mango (Mangifera indica). We used restriction site associated DNA sequencing to generate genomic single nucleotide polymorphism (SNP) data from 106 mango cultivars from seven geographical regions along with 52 samples of closely related species and unidentified cultivars to identify centers of mango genetic diversity and examine how post-domestication dispersal shaped the geographical distribution of diversity. We identify two gene pools of cultivated mango, representing Indian and Southeast Asian germplasm. We found no significant genetic bottleneck associated with the introduction of mango into new regions of the world. By contrast, we show that mango populations in introduced regions have elevated levels of diversity. Our results suggest that mango has a more complex history of domestication than previously supposed, perhaps including multiple domestication events, hybridization and regional selection. Our work has direct implications for mango breeding and genebank management, and also builds on recent efforts to understand how woody perennial crops respond to domestication.

Evaluating Maize Genotype Performance under Low Nitrogen Conditions Using RGB UAV Phenotyping Techniques

Maize is the most cultivated cereal in Africa in terms of land area and production, but low soil nitrogen availability often constrains yields. Developing new maize varieties with high and reliable yields using traditional crop breeding techniques in field conditions can be slow and costly. Remote sensing has become an important tool in the modernization of field-based high-throughput plant phenotyping (HTPP), providing faster gains towards the improvement of yield potential and adaptation to abiotic and biotic limiting conditions. We evaluated the performance of a set of remote sensing indices derived from red-green-blue (RGB) images along with field-based multispectral normalized difference vegetation index (NDVI) and leaf chlorophyll content (SPAD values) as phenotypic traits for assessing maize performance under managed low-nitrogen conditions. HTPP measurements were conducted from the ground and from an unmanned aerial vehicle (UAV). For the ground-level RGB indices, the strongest correlations to yield were observed with hue, greener green area (GGA), and a newly developed RGB HTPP index, NDLab (normalized difference Commission Internationale de I´Edairage (CIE)Lab index), while GGA and crop senescence index (CSI) correlated better with grain yield from the UAV. Regarding ground sensors, SPAD exhibited the closest correlation with grain yield, notably increasing in its correlation when measured in the vegetative stage. Additionally, we evaluated how different HTPP indices contributed to the explanation of yield in combination with agronomic data, such as anthesis silking interval (ASI), anthesis date (AD), and plant height (PH). Multivariate regression models, including RGB indices (R2 > 0.60), outperformed other models using only agronomic parameters or field sensors (R2 > 0.50), reinforcing RGB HTPP's potential to improve yield assessments. Finally, we compared the low-N results to the same panel of 64 maize genotypes grown under optimal conditions, noting that only 11% of the total genotypes appeared in the highest yield producing quartile for both trials. Furthermore, we calculated the grain yield loss index (GYLI) for each genotype, which showed a large range of variability, suggesting that low-N performance is not necessarily exclusive of high productivity in optimal conditions.

The Neutral Theory in Light of Natural Selection

In this perspective, we evaluate the explanatory power of the neutral theory of molecular evolution, 50 years after its introduction by Kimura. We argue that the neutral theory was supported by unreliable theoretical and empirical evidence from the beginning, and that in light of modern, genome-scale data, we can firmly reject its universality. The ubiquity of adaptive variation both within and between species means that a more comprehensive theory of molecular evolution must be sought.

Dynamic Patterns of Gene Expression Additivity and Regulatory Variation throughout Maize Development

Gene expression variation is a key component underlying phenotypic variation and heterosis. Transcriptome profiling was performed on 23 different tissues or developmental stages of two maize inbreds, B73 and Mo17, as well as their F₁ hybrid. The obtained large-scale datasets provided opportunities to monitor the developmental dynamics of differential expression, additivity for gene expression, and regulatory variation. The transcriptome can be divided into ∼30 000 genes that are expressed in at least one tissue of one inbred and an additional ∼10 000 ″silent" genes that are not expressed in any tissue of any genotype, 90% of which are non-syntenic relative to other grasses. Many (∼74%) of the expressed genes exhibit differential expression in at least one tissue. However, the majority of genes with differential expression do not exhibit consistent differential expression in different tissues. These genes often exhibit tissue-specific differential expression with equivalent expression in other tissues, and in many cases they switch the directionality of differential expression in different tissues. This suggests widespread variation for tissue-specific regulation of gene expression between the two maize inbreds B73 and Mo17. Nearly 5000 genes are expressed in only one parent in at least one tissue (single parent expression) and 97% of these genes are expressed at mid-parent levels or higher in the hybrid, providing extensive opportunities for hybrid complementation in heterosis. In general, additive expression patterns are much more common than non-additive patterns, and this trend is more pronounced for genes with strong differential expression or single parent expression. There is relatively little evidence for non-additive expression patterns that are maintained in multiple tissues. The analysis of allele-specific expression allowed classification of cis- and trans-regulatory variation. Genes with cis-regulatory variation often exhibit additive expression and tend to have more consistent regulatory variation throughout development. In contrast, genes with trans-regulatory variation are enriched for non-additive patterns and often show tissue-specific differential expression. Taken together, this study provides a deeper understanding of regulatory variation and the degree of additive gene expression throughout maize development. The dynamic nature of differential expression, additivity, and regulatory variation imply abundant variability for tissue-specific regulatory mechanisms and suggest that connections between transcriptome and phenome will require expression data from multiple tissues.

Regulatory Diversification of INDEHISCENT in the Capsella Genus Directs Variation in Fruit Morphology

Evolution of gene-regulatory sequences is considered the primary driver of morphological variation [1-3]. In animals, the diversity of body plans between distantly related phyla is due to the differential expression patterns of conserved "toolkit" genes [4]. In plants, variation in expression domains similarly underlie most of the reported diversity of organ shape both in natural evolution and in the domestication of crops [5-9]. The heart-shaped fruit from members of the Capsella genus is a morphological novelty that has evolved after Capsella diverged from Arabidopsis ∼8 mya [10]. Comparative studies of fruit growth in Capsella and Arabidopsis revealed that the difference in shape is caused by local control of anisotropic growth [11]. Here, we show that sequence variation in regulatory domains of the fruit-tissue identity gene, INDEHISCENT (IND), is responsible for expansion of its expression domain in the heart-shaped fruits from Capsella rubella. We demonstrate that expression of this CrIND gene in the apical part of the valves in Capsella contributes to the heart-shaped appearance. While studies on morphological diversity have revealed the importance of cis-regulatory sequence evolution, few examples exist where the downstream effects of such variation have been characterized in detail. We describe here how CrIND exerts its function on Capsella fruit shape by binding sequence elements of auxin biosynthesis genes to activate their expression and ensure auxin accumulation into highly localized maxima in the fruit valves. Thus, our data provide a direct link between changes in expression pattern and altered hormone homeostasis in the evolution of morphological novelty.

BRANCHED1: A Key Hub of Shoot Branching

Shoot branching is a key process for plant growth and fitness. Newly produced axes result from axillary bud outgrowth, which is at least partly mediated through the regulation of BRANCHED1 gene expression (BRC1/TB1/FC1). BRC1 encodes a pivotal bud-outgrowth-inhibiting transcription factor belonging to the TCP family. As the regulation of BRC1 expression is a hub for many shoot-branching-related mechanisms, it is influenced by endogenous (phytohormones and nutrients) and exogenous (light) inputs, which involve so-far only partly identified molecular networks. This review highlights the central role of BRC1 in shoot branching and its responsiveness to different stimuli, and emphasizes the different knowledge gaps that should be addressed in the near future.

Chinese lantern in Physalis is an advantageous morphological novelty and improves plant fitness

The origin of morphological novelties is an important but neglected issue of evolutionary biology. The fruit of the genus Physalis, a berry, is encapsulated by a novel morphological feature of the post-floral, accrescent calyx that is referred to as a Chinese lantern. The evolutionary developmental genetics of the Chinese lantern have been investigated in the last decade; however, the selective values of the morphological novelty remain elusive. Here, we measured the photosynthetic parameters of the fruiting calyces, monitored microclimatic variation within the Chinese lanterns during fruit development, performed floral-calyx-removal experiments, and recorded the fitness-related traits in Physalis floridana. Ultimately, we show that the green-fruiting calyx of Physalis has photosynthetic capabilities, thus serving as an energy source for fruit development. Moreover, the developing Chinese lantern provides a microclimate that benefits the development and maturation of berry and seed, and it improves plant fitness in terms of fruit/seed weight and number, and fruit maturation under low-temperature environments. Furthermore, the lantern structure facilitates the dispersal of fruits and seeds by water and wind. Our results suggest that the Chinese lantern morphology of Physalis is an evolutionary adaptive trait and improves plant fitness, thus providing new insight into the origin of morphological novelties.

Molecular analysis of mutant granule-bound starch synthase-I (waxy1) gene in diverse waxy maize inbreds

Waxy corn is popular beacuse of its high amylopectin due to mutation in granule-bound starch synthase-I or Waxy1 (Wx1) gene. Here, we characterized the wx1 allele among 24 diverse waxy inbreds using gene-based markers. A total of 29 alleles with average of 1.81 alleles/locus were observed. Major allele frequency varied from 0.42 to 1.00, with mean of 0.74. The polymorphism information content ranged from 0.00 to 0.56 (average 0.24). Three simple sequence repeat markers, viz., phi027, phi022 and phi061 were more polymorphic in the study. The mean heterozygosity was 0.04, which indicated attainment of higher levels of homozygosity. Dissimilarity coefficient varied from 0.00 to 0.90 with average of 0.51. Seventeen diverse haplotypes of wx1 allele were observed that was consistent with the pedigree. Cluster analyses grouped 24 genotypes into two main clusters each having sub-clusters. The information generated here possesses great potential for improvement of high amylopectin in maize through marker-assisted selection. This is the first report of molecular dissection of wx1 gene among the novel waxy inbreds developed in India.

Transcriptomics profiling in response to cold stress in cultivated rice and weedy rice

Weedy rice is an important germplasm resource for rice improvement because it has useful genes for many abiotic stresses including cold tolerance. We identified the cold tolerance and cold sensitivity of two weedy rice lines (WR 03-35 and WR 03-26) and two cultivated rice lines (Kongyu 131 and 9311). During the seedling stage of these lines, we used RNA-seq to measure changes in weedy rice and cultivated rice whole-genome transcriptome before and after cold treatment. We identified 14,213 and 14,730 differentially expressed genes (DEGs) in cold-tolerant genotypes (WR 03-35, Kongyu 131), and 9219 and 720 DEGs were observed in two cold-sensitive genotypes (WR 03-26, 9311). Many common and special DEGs were analyzed in cold-tolerant and cold-sensitive genotypes, respectively. Some typical genes related to cold stress such as the basic helix-loop-helix (bHLH) gene and leucine-rich repeat (LRR) domain gene etc. The number of these DEGs in cold-tolerant genotypes is more than those found in cold-sensitive genotypes. The gene ontology (GO) enrichment analyses showed significantly enriched terms for biological processes, cellular components and molecular functions. In addition, some genes related to several plant hormones such as abscisic acid (ABA), gibberellic acid (GA), auxin and ethylene were identified. To confirm the RNA-seq data, semi-quantitative RT-PCR and qRT-PCR were performed on 12 randomly selected DEGs. The expression patterns of RNA-seq on these genes corresponded with the semi-quantitative RT-PCR and qRT-PCR method. This study suggests the gene resources related to cold stress from weedy rice could be valuable for understanding the mechanisms involved in cold stress and rice breeding for improving cold tolerance.

Molecular cloning and sequence variance analysis of the TEOSINTE BRANCHED1 (TB1) gene in bermudagrass [Cynodon dactylon (L.) Pers]

TEOSINTE BRANCHED1 (TB1) encodes a TEOSINTE BRANCHED1, CYCLOIDEA, and PROLIFERATING CELL FACTOR (TCP) transcription factor that represses axillary bud outgrowth and lateral branch formation in plants. Previous studies have elucidated the essential tillering regulatory roles of TB1 in many grasses, including maize and rice; however, the functions of TB1 in turf grasses remain unclear. In this study, we cloned the CdTB1 gene from bermudagrass, an important turfgrass species, and characterized the transactivation function of the CdTB1 protein. Sequencing the CdTB1 gene locus in a mini-core germplasm collection of Chinese bermudagrasses led to the successful identification of 66 SNP and 2 indel mutations in the protein-coding region as well as 28 SNP and 11 indel mutations in the promoter region. Interestingly, mutations in the C-terminal transactivation domain of the CdTB1 protein had no significant influence on the transactivation activity, whereas a novel 335-bp insertion mutation located in the promoter region could significantly increase the expression of the CdTB1 gene. Furthermore, wild accessions of bermudagrass harboring the novel insertion mutation were found to have significantly reduced tillers compared with other accessions, suggesting a negative correlation between the mutation and tillering. The results of this study not only expanded our knowledge of TB1 gene expression regulation but also provided possible molecular markers to breed cultivars of turf and forage grasses with specific architectural features.

Genome-Wide Identification and Expression Profiling of the TCP Family Genes in Spike and Grain Development of Wheat (Triticum aestivum L.)

The TCP family genes are plant-specific transcription factors and play important roles in plant development. TCPs have been evolutionarily and functionally studied in several plants. Although common wheat (Triticum aestivum L.) is a major staple crop worldwide, no systematic analysis of TCPs in this important crop has been conducted. Here, we performed a genome-wide survey in wheat and found 66 TCP genes that belonged to 22 homoeologous groups. We then mapped these genes on wheat chromosomes and found that several TCP genes were duplicated in wheat including the ortholog of the maize TEOSINTE BRANCHED 1. Expression study using both RT-PCR and in situ hybridization assay showed that most wheat TCP genes were expressed throughout development of young spike and immature seed. Cis-acting element survey along promoter regions suggests that subfunctionalization may have occurred for homoeologous genes. Moreover, protein-protein interaction experiments of three TCP proteins showed that they can form either homodimers or heterodimers. Finally, we characterized two TaTCP9 mutants from tetraploid wheat. Each of these two mutant lines contained a premature stop codon in the A subgenome homoeolog that was dominantly expressed over the B subgenome homoeolog. We observed that mutation caused increased spike and grain lengths. Together, our analysis of the wheat TCP gene family provides a start point for further functional study of these important transcription factors in wheat.

Differential transcriptome patterns associated with early seedling development in a wild and a domesticated common bean (Phaseolus vulgaris L.) accession

Genes that control "Domestication Syndrome" traits were direct targets of selection, like those controlling increased seed size in the common bean. However, selection for this trait brought about unintentional selection on genes controlling seedling growth. We hypothesized that wild and domesticated plants have different early seedling growth patterns as an indirect consequence of selection for a larger seed size during domestication, and those differences resulted from changes in gene expression patterns of the wild ancestor. Large seeds pose a challenge to reserve remobilization during early heterotrophic growth, particularly during a transition towards more fertile alluvial soils. To address our hypothesis, we characterized the patterns of gene expression of cotyledon, root, and leaf tissues of 7-day old seedlings of a wild and a landrace accession of the common bean. Differential expression analyses detected genes with contrasting patterns of expression between the two genotypes in all three tissues. Some of the differentially expressed genes with contrasting genotypic patterns are known to have domestication-related signatures of selection. Among these genes were some transcription factors associated with key roles in development. These genes may represent targets of indirect selection and ultimately explain the growth phenotypic differences between wild and domesticated seedlings.

Specific LTR-Retrotransposons Show Copy Number Variations between Wild and Cultivated Sunflowers

The relationship between variation of the repetitive component of the genome and domestication in plant species is not fully understood. In previous work, variations in the abundance and proximity to genes of long terminal repeats (LTR)-retrotransposons of sunflower (Helianthus annuus L.) were investigated by Illumina DNA sequencingtocompare cultivars and wild accessions. In this study, we annotated and characterized 22 specific retrotransposon families whose abundance varies between domesticated and wild genotypes. These families mostly belonged to the Chromovirus lineage of the Gypsy superfamily and were distributed overall chromosomes. They were also analyzed in respect to their proximity to genes. Genes close to retrotransposon were classified according to biochemical pathways, and differences between domesticated and wild genotypes are shown. These data suggest that structural variations related to retrotransposons might have occurred to produce phenotypic variation between wild and domesticated genotypes, possibly by affecting the expression of genes that lie close to inserted or deleted retrotransposons and belong to specific biochemical pathways as those involved in plant stress responses.

Genetic diversity and evolution of reduced sulfur storage during domestication of maize

The domestication of maize has spanned a period of over 9000 years, during which time its wild relative teosinte underwent natural and artificial selection. We hypothesize that environmental conditions could have played a major role in this process. One factor of environmental variation is soil composition, which includes sulfur availability. Sulfur is reduced during photosynthesis and is used to synthesize cysteine and methionine, which drive the accumulation of δ10 (Zm00001d045937), δ18 (Zm00001d037436), β15 (Zm00001d035760), γ16 (Zm00001d005793), γ27 (Zm00001d020592), and γ50 (Zm00001d020591) zeins, representing the zein2 fraction (z2) of storage proteins in maize seeds. In this study, polymorphisms and haplotypes were detected based on six z2 genes in 60 maize and teosintes lines. Haplotypes were unevenly distributed, and abundant genetic diversity was found in teosintes. Polymorphism was highest in z2δ18, whereas for z2β15 single nucleotide polymorphism (SNP) density and insertion/deletion (indel) abundance were the lowest, indicating differential roles in seed evolution. Indels showed a clustered distribution, and most of these derived from teosintes. The indels not only led to tandem repeat polymorphisms, but also to frameshift mutations, which could also be used as null variants. In addition, neutral evolutionary tests, phylogenetic analyses, and population structures indicated that z2δ10 and z2γ50 had undergone natural selection. Indeed, a natural selection imprint could also be found with z2γ27 and z2γ16, whereas z2δ18 and z2β15 tended to be under neutral evolution. These results suggested that genetic diversity and evolution of a subset of sulfur-rich zeins could be under environmental adaptation during maize domestication.

The origin and early evolution of vascular plant shoots and leaves

The morphology of plant fossils from the Rhynie chert has generated longstanding questions about vascular plant shoot and leaf evolution, for instance, which morphologies were ancestral within land plants, when did vascular plants first arise and did leaves have multiple evolutionary origins? Recent advances combining insights from molecular phylogeny, palaeobotany and evo-devo research address these questions and suggest the sequence of morphological innovation during vascular plant shoot and leaf evolution. The evidence pinpoints testable developmental and genetic hypotheses relating to the origin of branching and indeterminate shoot architectures prior to the evolution of leaves, and demonstrates underestimation of polyphyly in the evolution of leaves from branching forms in 'telome theory' hypotheses of leaf evolution. This review discusses fossil, developmental and genetic evidence relating to the evolution of vascular plant shoots and leaves in a phylogenetic framework.This article is part of a discussion meeting issue 'The Rhynie cherts: our earliest terrestrial ecosystem revisited'.

Circular RNAs mediated by transposons are associated with transcriptomic and phenotypic variation in maize

Circular RNAs (circRNAs) are covalently closed RNA molecules. Recent studies have shown that circRNAs can arise from the transcripts of transposons. Given the prevalence of transposons in the maize genome and dramatic genomic variation driven by transposons, we hypothesize that transposons in maize may be involved in the formation of circRNAs and further modulate phenotypic variation. We performed circRNA-Seq on B73 seedling leaves and uncovered 2804 high-confidence maize circRNAs, which show distinct genomic features. Comprehensive analyses demonstrated that sequences related to LINE1-like elements (LLEs) and their Reverse Complementary Pairs (LLERCPs) are significantly enriched in the flanking regions of circRNAs. Interestingly, as the number of LLERCPs increase, the accumulation of circRNAs varies, whereas that of linear transcripts decreases. Furthermore, genes with LLERCP-mediated circRNAs are enriched among loci that are associated with phenotypic variation. These results suggest that circRNAs are likely to be involved in the modulation of phenotypic variation by LLERCPs. Further, we showed that the presence/absence variation of LLERCPs was associated with expression variation of circRNA-circ1690 and was related to ear height, potentially through the interplay between circRNAs and functional linear transcripts. Our first study of maize circRNAs uncovers a potential new way for transposons to modulate transcriptomic and phenotypic variations.

Evolutionary history of the NAM-B1 gene in wild and domesticated tetraploid wheat

BACKGROUND

The NAM-B1 gene in wheat has for almost three decades been extensively studied and utilized in breeding programs because of its significant impact on grain protein and mineral content and pleiotropic effects on senescence rate and grain size. First detected in wild emmer wheat, the wild-type allele of the gene has been introgressed into durum and bread wheat. Later studies have, however, also found the presence of the wild-type allele in some domesticated subspecies. In this study we trace the evolutionary history of the NAM-B1 in tetraploid wheat species and evaluate it as a putative domestication gene.

RESULTS

Genotyping of wild and landrace tetraploid accessions showed presence of only null alleles in durum. Domesticated emmer wheats contained both null alleles and the wild-type allele while wild emmers, with one exception, only carried the wild-type allele. One of the null alleles consists of a deletion that covers several 100 kb. The other null-allele, a one-basepair frame-shift insertion, likely arose among wild emmer. This allele was the target of a selective sweep, extending over several 100 kb.

CONCLUSIONS

The NAM-B1 gene fulfils some criteria for being a domestication gene by encoding a trait of domestication relevance (seed size) and is here shown to have been under positive selection. The presence of both wild-type and null alleles in domesticated emmer does, however, suggest the gene to be a diversification gene in this species. Further studies of genotype-environment interactions are needed to find out under what conditions selection on different NAM-B1 alleles have been beneficial.

Signatures of soft sweeps across the Dt1 locus underlying determinate growth habit in soya bean [Glycine max (L.) Merr.]

Determinate growth habit is an agronomically important trait associated with domestication in soya bean. Previous studies have demonstrated that the emergence of determinacy is correlated with artificial selection on four nonsynonymous mutations in the Dt1 gene. To better understand the signatures of the soft sweeps across the Dt1 locus and track the origins of the determinate alleles, we examined patterns of nucleotide variation in Dt1 and the surrounding genomic region of approximately 800 kb. Four local, asymmetrical hard sweeps on four determinate alleles, sized approximately 660, 120, 220 and 150 kb, were identified, which constitute the soft sweeps for the adaptation. These variable-sized sweeps substantially reflected the strength and timing of selection and indicated that the selection on the alleles had been completed rapidly within half a century. Statistics of EHH, iHS, H12 and H2/H1 based on haplotype data had the power to detect the soft sweeps, revealing distinct signatures of extensive long-range LD and haplotype homozygosity, and multiple frequent adaptive haplotypes. A haplotype network constructed for Dt1 and a phylogenetic tree based on its extended haplotype block implied independent sources of the adaptive alleles through de novo mutations or rare standing variation in quick succession during the selective phase, strongly supporting multiple origins of the determinacy. We propose that the adaptation of soya bean determinacy is guided by a model of soft sweeps and that this model might be indispensable during crop domestication or evolution.

Understanding and engineering plant form

A plant's form is an important determinant of its fitness and economic value. Here, we review strategies for producing plants with altered forms. Historically, the process of changing a plant's form has been slow in agriculture, requiring iterative rounds of growth and selection. We discuss modern techniques for identifying genes involved in the development of plant form and tools that will be needed to effectively design and engineer plants with altered forms. Synthetic genetic circuits are highlighted for their potential to generate novel plant forms. We emphasize understanding development as a prerequisite to engineering and discuss the potential role of computer models in translating knowledge about single genes or pathways into a more comprehensive understanding of development.

Domestication rewired gene expression and nucleotide diversity patterns in tomato

Plant domestication has led to considerable phenotypic modifications from wild species to modern varieties. However, although changes in key traits have been well documented, less is known about the underlying molecular mechanisms, such as the reduction of molecular diversity or global gene co-expression patterns. In this study, we used a combination of gene expression and population genetics in wild and crop tomato to decipher the footprints of domestication. We found a set of 1729 differentially expressed genes (DEG) between the two genetic groups, belonging to 17 clusters of co-expressed DEG, suggesting that domestication affected not only individual genes but also regulatory networks. Five co-expression clusters were enriched in functional terms involving carbohydrate metabolism or epigenetic regulation of gene expression. We detected differences in nucleotide diversity between the crop and wild groups specific to DEG. Our study provides an extensive profiling of the rewiring of gene co-expression induced by the domestication syndrome in one of the main crop species.

QTL mapping of domestication and diversifying selection related traits in round-fruited semi-wild Xishuangbanna cucumber (Cucumis sativus L. var. xishuangbannanesis)

QTL analysis revealed 11 QTL underlying flowering time and fruit size variation in the semi-wild Xishuangbanna cucumber, of which, FT6.2 and FS5.2 played the most important roles in determining photoperiod-dependent flowering time and round-fruit shape, respectively. Flowering time and fruit size are two important traits in domestication and diversifying selection in cucumber, but their genetic basis is not well understood. Here we reported QTL mapping results on flowering time and fruit size with F₂ and F_2:3 segregating populations derived from the cross between WI7200, a small fruited, early flowering primitive cultivated cucumber and WI7167, a round-fruited, later flowering semi-wild Xishuangbanna (XIS) cucumber. A linkage map with 267 microsatellite marker loci was developed with 138 F₂ plants. Phenotypic data of male and female flowering time, fruit length and diameter and three other traits (mature fruit weight and number, and seedling hypocotyl length) were collected in multiple environments. Three flowering time QTL, FT1.1, FT5.1 and FT6.2 were identified, in which FT6.2 played the most important role in conferring less photoperiod sensitive early flowering during domestication whereas FT1.1 seemed more influential in regulating flowering time within the cultivated cucumber. Eight consensus fruit size QTL distributed in 7 chromosomes were detected, each of which contributed to both longitudinal and radial growth in cucumber fruit development. Among them, FS5.2 on chromosome 5 exhibited the largest effect on the determination of round fruit shape that was characteristic of the WI7167 XIS cucumber. Possible roles of these flowering time and fruit size QTL in domestication of cucumber and crop evolution of the semi-wild XIS cucumber, as well as the genetic basis of round fruit shape in cucumber are discussed.

Evolution of reduced co-activator dependence led to target expansion of a starvation response pathway

Although combinatorial regulation is a common feature in gene regulatory networks, how it evolves and affects network structure and function is not well understood. In S. cerevisiae, the phosphate starvation (PHO) responsive transcription factors Pho4 and Pho2 are required for gene induction and survival during phosphate starvation. In the related human commensal C. glabrata, Pho4 is required but Pho2 is dispensable for survival in phosphate starvation and is only partially required for inducing PHO genes. Phylogenetic survey suggests that reduced dependence on Pho2 evolved in C. glabrata and closely related species. In S. cerevisiae, less Pho2-dependent Pho4 orthologs induce more genes. In C. glabrata, its Pho4 binds to more locations and induces three times as many genes as Pho4 in S. cerevisiae does. Our work shows how evolution of combinatorial regulation allows for rapid expansion of a gene regulatory network’s targets, possibly extending its physiological functions.

The diversity of life on Earth has intrigued generations of scientists and nature lovers alike. Research over recent decades has revealed that much of the diversity we can see did not require the invention of new genes. Instead, living forms diversified mostly by using old genes in new ways – for example, by changing when or where an existing gene became active. This kind of change is referred to as “regulatory evolution”.

A class of proteins called transcription factors are hot spots in regulatory evolution. These proteins recognize specific sequences of DNA to control the activity of other genes, and so represent the “readers” of the genetic information. Small changes to how a transcription factor is regulated, or the genes it targets, can lead to dramatic changes in an organism. Before we can understand how life on Earth evolved to be so diverse, scientists must first answer how transcription factors evolve and what consequences this has on their target genes.

So far, most studies of regulatory evolution have focused on networks of transcription factors and genes that control how an organism develops. He et al. have now studied a regulatory network that is behind a different process, namely how an organism responds to stress or starvation. These two types of regulatory networks are structured differently and work in different ways. These differences made He et al. wonder if the networks evolved differently too.

The chemical phosphate is an essential nutrient for all living things, and He et al. compared how two different species of yeast responded to a lack of phosphate. The key difference was how much a major transcription factor known as Pho4 depended on a so-called co-activator protein named Pho2 to carry out its role. Baker’s yeast (Saccharomyces cerevisiae), which is commonly used in laboratory experiments, requires both Pho4 and Pho2 to activate about 20 genes when inorganic phosphate is not available in its environment. However, in a related yeast species called Candida glabrata, Pho4 has evolved to depend less on Pho2. He et al. went on to show that, as well as being less dependent on Pho2, Pho4 in C. glabrata activates more than three times as many genes as Pho4 in S. cerevisiae does in the absence of phosphate. These additional gene targets for Pho4 in C. glabrata are predicted to extend the network’s activities, and allow it to regulate new process including the yeast’s responses to other types of stress and the building of the yeast’s cell wall.

Together these findings show a new way that regulatory networks can evolve, that is, by reducing its dependence on the co-activator, a transcription factor can expand the number of genes it targets. This has not been seen for regulatory networks related to development, suggesting that different networks can indeed evolve in different ways. Lastly, because disease-causing microbes are often stressed inside their hosts and C. glabrata sometimes infects humans, understanding how this yeast’s response to stress has evolved may lead to new ways to prevent and treat this infection.

Effect of transcription factor resource sharing on gene expression noise

Gene expression is intrinsically a stochastic (noisy) process with important implications for cellular functions. Deciphering the underlying mechanisms of gene expression noise remains one of the key challenges of regulatory biology. Theoretical models of transcription often incorporate the kinetics of how transcription factors (TFs) interact with a single promoter to impact gene expression noise. However, inside single cells multiple identical gene copies as well as additional binding sites can compete for a limiting pool of TFs. Here we develop a simple kinetic model of transcription, which explicitly incorporates this interplay between TF copy number and its binding sites. We show that TF sharing enhances noise in mRNA distribution across an isogenic population of cells. Moreover, when a single gene copy shares it’s TFs with multiple competitor sites, the mRNA variance as a function of the mean remains unaltered by their presence. Hence, all the data for variance as a function of mean expression collapse onto a single master curve independent of the strength and number of competitor sites. However, this result does not hold true when the competition stems from multiple copies of the same gene. Therefore, although previous studies showed that the mean expression follows a universal master curve, our findings suggest that different scenarios of competition bear distinct signatures at the level of variance. Intriguingly, the introduction of competitor sites can transform a unimodal mRNA distribution into a multimodal distribution. These results demonstrate the impact of limited availability of TF resource on the regulation of noise in gene expression.

Genetically identical cells, even when they are exposed to the same environmental conditions, display incredible diversity. Gene expression noise is attributed to be a key source of this phenotypic diversity. Transcriptional dynamics is a dominant source of expression noise. Although scores of theoretical and experimental studies have explored how noise is regulated at the level of transcription, most of them focus on the gene specific, cis regulatory elements, such as the number of transcription factor (TF) binding sites, their binding strength, etc. However, how the global properties of transcription, such as the limited availability of TFs impact noise in gene expression remains rather elusive. Here we build a theoretical model that incorporates the effect of limiting TF pool on gene expression noise. We find that competition between genes for TFs leads to enhanced variability in mRNA copy number across an isogenic population. Moreover, for gene copies sharing TFs with other competitor sites, mRNA variance as a function of the mean shows distinct imprints for one gene copy and multiple gene copies respectively. This stands in sharp contrast to the universal behavior found in mean expression irrespective of the different scenarios of competition. An interesting feature of competition is that introduction of competitor sites can transform a unimodal mRNA distribution into a multimodal distribution, which could lead to phenotypic variability.

The origin and evolution of fibromelanosis in domesticated chickens: Genomic comparison of Indonesian Cemani and Chinese Silkie breeds

Like Chinese Silkie, Indonesian Ayam Cemani exhibits fibromelanosis or dermal hyperpigmentation and possesses complex segmental duplications on chromosome 20 that involve the endothelin 3 gene, EDN3. A genomic region, DR1 of 127 kb, together with another region, DR2 of 171 kb, was duplicated by unequal crossing over, accompanied by inversion of one DR2. Quantitative PCR and copy number variation analyses on the Cemani genome sequence confirmed the duplication of EDN3. These genetic arrangements are identical in Cemani and Silkie, indicating a single origin of the genetic cause of Fm. The two DR1s harbor two distinct EDN3 haplotypes in a form of permanent heterozygosity, although they remain allelic in the ancestral Red Jungle Fowl population and some domesticated chicken breeds, with their allelic divergence time being as recent as 0.3 million years ago. In Cemani and Silkie breeds, artificial selection favoring the Fm phenotype has left an unambiguous record for selective sweep that extends in both directions from tandemly duplicated EDN3 loci. This highly homozygous tract is different in length between Cemani and Silkie, reflecting their distinct breeding histories. It is estimated that the Fm phenotype came into existence at least 6600-9100 years ago, prior to domestication of Cemani and Silkie, and that throughout domestication there has been intense artificial selection with strength s > 50% in each breed.

Gene co-expression network connectivity is an important determinant of selective constraint

While several studies have investigated general properties of the genetic architecture of natural variation in gene expression, few of these have considered natural, outbreeding populations. In parallel, systems biology has established that a general feature of biological networks is that they are scale-free, rendering them buffered against random mutations. To date, few studies have attempted to examine the relationship between the selective processes acting to maintain natural variation of gene expression and the associated co-expression network structure. Here we utilised RNA-Sequencing to assay gene expression in winter buds undergoing bud flush in a natural population of Populus tremula, an outbreeding forest tree species. We performed expression Quantitative Trait Locus (eQTL) mapping and identified 164,290 significant eQTLs associating 6,241 unique genes (eGenes) with 147,419 unique SNPs (eSNPs). We found approximately four times as many local as distant eQTLs, with local eQTLs having significantly higher effect sizes. eQTLs were primarily located in regulatory regions of genes (UTRs or flanking regions), regardless of whether they were local or distant. We used the gene expression data to infer a co-expression network and investigated the relationship between network topology, the genetic architecture of gene expression and signatures of selection. Within the co-expression network, eGenes were underrepresented in network module cores (hubs) and overrepresented in the periphery of the network, with a negative correlation between eQTL effect size and network connectivity. We additionally found that module core genes have experienced stronger selective constraint on coding and non-coding sequence, with connectivity associated with signatures of selection. Our integrated genetics and genomics results suggest that purifying selection is the primary mechanism underlying the genetic architecture of natural variation in gene expression assayed in flushing leaf buds of P. tremula and that connectivity within the co-expression network is linked to the strength of purifying selection.