To investigate the natural variability of leaf metabolism and enzymatic activity in a maize inbred population, statistical and network analyses were employed on metabolite and enzyme profiles. The test of coefficient of variation showed that sugars and amino acids displayed opposite trends in their variance within the population, consistently with their related enzymes. The overall higher CV values for metabolites as compared to the tested enzymes are indicative for their greater phenotypic plasticity. H(2) tests revealed galactinol (1) and asparagine (0.91) as the highest scorers among metabolites and nitrate reductase (0.73), NAD-glutamate dehydrogenase (0.52), and phosphoglucomutase (0.51) among enzymes. The overall low H(2) scores for metabolites and enzymes are suggestive for a great environmental impact or gene-environment interaction. Correlation-based network generation followed by community detection analysis, partitioned the network into three main communities and one dyad, (i) reflecting the different levels of phenotypic plasticity of the two molecular classes as observed for the CV values and (ii) highlighting the concerted changes between classes of chemically related metabolites. Community 1 is composed mainly of enzymes and specialized metabolites, community 2' is enriched in N-containing compounds and phosphorylated-intermediates. The third community contains mainly organic acids and sugars. Cross-community linkages are supported by aspartate, by the photorespiration amino acids glycine and serine, by the metabolically related GABA and putrescine, and by citrate. The latter displayed the strongest node-betweenness value (185.25) of all nodes highlighting its fundamental structural role in the connectivity of the network by linking between different communities and to the also strongly connected enzyme aldolase.

Linking genotype with phenotype is a fundamental goal in biology and requires robust data for both. Recent advances in plant-genome sequencing have expedited comparisons among multiple-related individuals. The abundance of structural genomic within-species variation that has been discovered indicates that a single reference genome cannot represent the complete sequence diversity of a species, leading to the expansion of the pan-genome concept. For high-resolution forward genetics, this unprecedented access to genomic variation should be paralleled and integrated with phenotypic characterization of genetic diversity. We developed a multi-parental framework for trait dissection in melon (Cucumis melo), leveraging a novel pan-genome constructed for this highly variable cucurbit crop. A core subset of 25 diverse founders (MelonCore25), consisting of 24 accessions from the two widely cultivated subspecies of C. melo, encompassing 12 horticultural groups, and 1 feral accession was sequenced using a combination of short- and long-read technologies, and their genomes were assembled de novo. The construction of this melon pan-genome exposed substantial variation in genome size and structure, including detection of ~300 000 structural variants and ~9 million SNPs. A half-diallel derived set of 300 F2 populations, representing all possible MelonCore25 parental combinations, was constructed as a framework for trait dissection through integration with the pan-genome. We demonstrate the potential of this unified framework for genetic analysis of various melon traits, including rind color intensity and pattern, fruit sugar content, and resistance to fungal diseases. We anticipate that utilization of this integrated resource will enhance genetic dissection of important traits and accelerate melon breeding.

Molecular markers allowed breeders to mendelize quantitative trait loci (QTL) providing another demonstration that quantitative traits are governed by the same principles as single qualitative genes. This research extends the QTL analysis to two and three QTL and tests our ability to mendelize an oligogenic trait. In tomato, agricultural yield is determined by the weight of the fruits harvested per unit area and the total soluble solids (% Brix)-sugars and acids. The current study explores the segregation of multiple independent yield-related QTL that were identified and mapped using introgression lines (IL) of Solanum pennellii in cultivated processing tomato (S. lycopersicum). We screened 45 different double and triple IL-QTL combinations for agricultural yield, to identify QTL pyramids that behaved in an additive manner and were suitable substrate for mendelizing an oligogenic trait. A pyramid of three independent QTL that significantly improved Brix(∗)Yield (BXY - the soluble solids output per unit area) compared to M82 was selected. In the progenies of the tri-hybrid we bred using markers a nearly isogenic 'immortalized F2.' While the common mode of QTL-QTL interactions across the 45 IL-QTLs combinations was less than additive, the three QTLs in the selected triple-stack performed in an additive manner which made it an exceptional material for breeding. This study demonstrates that using the phenotypic effect of all 27 possible QTL-alleles combinations it is possible to make reliable predictions about the genotypes that will maximize the yield.

The broad variability of Cucumis melo (melon, Cucurbitaceae) presents a challenge to conventional classification and organization within the species. To shed further light on the infraspecific relationships within C. melo, we compared genotypic and metabolomic similarities among 44 accessions representative of most of the cultivar-groups. Genotyping-by-sequencing (GBS) provided over 20,000 single-nucleotide polymorphisms (SNPs). Metabolomics data of the mature fruit flesh and rind provided over 80,000 metabolomic and elemental features via an orchestra of six complementary metabolomic platforms. These technologies probed polar, semi-polar, and non-polar metabolite fractions as well as a set of mineral elements and included both flavor- and taste-relevant volatile and non-volatile metabolites. Together these results enabled an estimate of "metabolomic/elemental distance" and its correlation with the genetic GBS distance of melon accessions. This study indicates that extensive and non-targeted metabolomics/elemental characterization produced classifications that strongly, but not completely, reflect the current and extensive genetic classification. Certain melon Groups, such as Inodorous, clustered in parallel with the genetic classifications while other genome to metabolome/element associations proved less clear. We suggest that the combined genomic, metabolic, and element data reflect the extensive sexual compatibility among melon accessions and the breeding history that has, for example, targeted metabolic quality traits, such as taste and flavor.

Earliness and ripening behavior are important attributes of fruits on and off the vine, and affect quality and preference of both growers and consumers. Fruit ripening is a complex physiological process that involves metabolic shifts affecting fruit color, firmness, and aroma production. Melon is a promising model crop for the study of fruit ripening, as the full spectrum of climacteric behavior is represented across the natural variation. Using Recombinant Inbred Lines (RILs) population derived from the parental lines "Dulce" (reticulatus, climacteric) and "Tam Dew" (inodorus, non-climacteric) that vary in earliness and ripening traits, we mapped QTLs for ethylene emission, fruit firmness and days to flowering and maturity. To further annotate the main QTL intervals and identify candidate genes, we used Oxford Nanopore long-read sequencing in combination with Illumina short-read resequencing, to assemble the parental genomes de-novo. In addition to 2.5 million genome-wide SNPs and short InDels detected between the parents, we also highlight here the structural variation between these lines and the reference melon genome. Through systematic multi-layered prioritization process, we identified 18 potential polymorphisms in candidate genes within multi-trait QTLs. The associations of selected SNPs with earliness and ripening traits were further validated across a panel of 177 diverse melon accessions and across a diallel population of 190 F1 hybrids derived from a core subset of 20 diverse parents. The combination of advanced genomic tools with diverse germplasm and targeted mapping populations is demonstrated as a way to leverage forward genetics strategies to dissect complex horticulturally important traits.

Central carbon metabolism (CCM) is a fundamental component of life. The participating genes and enzymes are thought to be structurally and functionally conserved across and within species. Association mapping utilizes a rich history of mutation and recombination to achieve high resolution mapping. Therefore, applying association mapping in maize (Zea mays ssp. mays), the most diverse model crop species, to study the genetics of CCM is a particularly attractive system.

Cucumis melo is highly diverse for fruit traits providing wide breeding and genetic research opportunities, including genome-wide association (GWA) analysis. We used a collection of 177 accessions representing the two C. melo subspecies and 11 horticultural groups for detailed characterization of fruit traits variation and evaluation of the potential of GWA for trait mapping in melon. Through genotyping-by-sequencing, 23,931 informative SNPs were selected for genome-wide analyses. We found that linkage-disequilibrium decays at ~100 Kb in this collection and that population structure effect on association results varies between traits. We mapped several monogenic traits to narrow intervals overlapping with known causative genes, demonstrating the potential of diverse collections and GWA for mapping Mendelian traits to a candidate-gene level in melon. We further report on mapping of fruit shape quantitative trait loci (QTLs) and comparison with multiple previous QTL studies. Expansion of sample size and a more balanced representation of taxonomic groups might improve efficiency for simple traits dissection. But, as in other plant species, integrated linkage-association multi-allelic approaches are likely to produce better combination of statistical power, diversity capture and mapping resolution in melon. Our data can be utilized for selection of the most appropriate accessions for such approaches.

Color and pigment contents are important aspects of fruit quality and consumer acceptance of cucurbit crops. Here, we describe the independent mapping and cloning of a common causative APRR2 gene regulating pigment accumulation in melon and watermelon. We initially show that the APRR2 transcription factor is causative for the qualitative difference between dark and light green rind in both crops. Further analyses establish the link between sequence or expression level variations in the CmAPRR2 gene and pigment content in the rind and flesh of mature melon fruits. A genome-wide association study (GWAS) of young fruit rind color in a panel composed of 177 diverse melon accessions did not result in any significant association, leading to an earlier assumption that multiple genes are involved in shaping the overall phenotypic variation in this trait. Through resequencing of 25 representative accessions and allelism tests between light rind accessions, we show that multiple independent single nucleotide polymorphisms in the CmAPRR2 gene are causative of the light rind phenotype. The multi-haplotypic nature of this gene explains the lack of detection power obtained through genotyping by sequencing-based GWAS and confirms the pivotal role of this gene in shaping fruit color variation in melon. This study demonstrates the power of combining bi- and multi-allelic designs with deep sequencing, to resolve lack of power due to high haplotypic diversity and low allele frequencies. Due to its central role and broad effect on pigment accumulation in fruits, the APRR2 gene is an attractive target for carotenoid bio-fortification of cucurbit crops.

Carotenoids have various roles in plant physiology. Plant carotenoids are synthesized in plastids and are highly abundant in the chromoplasts of ripening fleshy fruits. Considerable research efforts have been devoted to elucidating mechanisms that regulate carotenoid biosynthesis, yet, little is known about the mechanism that triggers storage capacity, mainly through chromoplast differentiation. The Orange gene (OR) product stabilizes phytoene synthase protein (PSY) and triggers chromoplast differentiation. OR underlies carotenoid accumulation in orange cauliflower and melon. The OR's 'golden SNP', found in melon, alters the highly evolutionary conserved Arginine108 to Histidine and controls β-carotene accumulation in melon fruit, in a mechanism yet to be elucidated. We have recently shown that similar carotenogenic metabolic flux is active in non-orange and orange melon fruit. This flux probably leads to carotenoid turnover but known carotenoid turnover products are not detected in non-orange fruit. Arrest of this metabolic flux, using chemical inhibitors or mutations, induces carotenoid accumulation and biogenesis of chromoplasts, regardless of the allelic state of OR. We suggest that the 'golden SNP' induces β-carotene accumulation probably by negatively affecting the capacity to synthesize downstream compounds. The accumulation of carotenoids induces chromoplast biogenesis through a metabolite-induced mechanism. Carotenogenic turnover flux can occur in non-photosynthetic tissues, which do not accumulate carotenoids. Arrest of this flux by the 'golden SNP' or other flux-arrest mutations is a potential tool for the biofortification of agricultural products with carotenoids.

Years of selection for desirable fruit quality traits in dessert watermelon (Citrullus lanatus) has resulted in a narrow genetic base in modern cultivars. Development of novel genomic and genetic resources offers great potential to expand genetic diversity and improve important traits in watermelon. Here, we report a high-quality genome sequence of watermelon cultivar 'Charleston Gray', a principal American dessert watermelon, to complement the existing reference genome from '97103', an East Asian cultivar. Comparative analyses between genomes of 'Charleston Gray' and '97103' revealed genomic variants that may underlie phenotypic differences between the two cultivars. We then genotyped 1365 watermelon plant introduction (PI) lines maintained at the U.S. National Plant Germplasm System using genotyping-by-sequencing (GBS). These PI lines were collected throughout the world and belong to three Citrullus species, C. lanatus, C. mucosospermus and C. amarus. Approximately 25 000 high-quality single nucleotide polymorphisms (SNPs) were derived from the GBS data using the 'Charleston Gray' genome as the reference. Population genomic analyses using these SNPs discovered a close relationship between C. lanatus and C. mucosospermus and identified four major groups in these two species correlated to their geographic locations. Citrullus amarus was found to have a distinct genetic makeup compared to C. lanatus and C. mucosospermus. The SNPs also enabled identification of genomic regions associated with important fruit quality and disease resistance traits through genome-wide association studies. The high-quality 'Charleston Gray' genome and the genotyping data of this large collection of watermelon accessions provide valuable resources for facilitating watermelon research, breeding and improvement.

The ability to study human post-implantation development remains limited owing to ethical and technical challenges associated with intrauterine development after implantation1. Embryo-like models with spatially organized morphogenesis and structure of all defining embryonic and extra-embryonic tissues of the post-implantation human conceptus (that is, the embryonic disc, the bilaminar disc, the yolk sac, the chorionic sac and the surrounding trophoblast layer) remain lacking1,2. Mouse naive embryonic stem cells have recently been shown to give rise to embryonic and extra-embryonic stem cells capable of self-assembling into post-gastrulation structured stem-cell-based embryo models with spatially organized morphogenesis (called SEMs)3. Here we extend those findings to humans using only genetically unmodified human naive embryonic stem cells (cultured in human enhanced naive stem cell medium conditions)4. Such human fully integrated and complete SEMs recapitulate the organization of nearly all known lineages and compartments of post-implantation human embryos, including the epiblast, the hypoblast, the extra-embryonic mesoderm and the trophoblast layer surrounding the latter compartments. These human complete SEMs demonstrated developmental growth dynamics that resemble key hallmarks of post-implantation stage embryogenesis up to 13-14 days after fertilization (Carnegie stage 6a). These include embryonic disc and bilaminar disc formation, epiblast lumenogenesis, polarized amniogenesis, anterior-posterior symmetry breaking, primordial germ-cell specification, polarized yolk sac with visceral and parietal endoderm formation, extra-embryonic mesoderm expansion that defines a chorionic cavity and a connecting stalk, and a trophoblast-surrounding compartment demonstrating syncytium and lacunae formation. This SEM platform will probably enable the experimental investigation of previously inaccessible windows of human early post implantation up to peri-gastrulation development.

The voltage-gated Kv2.1 potassium channel encoded by KCNB1 produces the major delayed rectifier potassium current in pyramidal neurons. Recently, de novo heterozygous missense KCNB1 mutations have been identified in three patients with epileptic encephalopathy and a patient with neurodevelopmental disorder. However, the frequency of KCNB1 mutations in infantile epileptic patients and their effects on neuronal activity are yet unknown. We searched whole exome sequencing data of a total of 437 patients with infantile epilepsy, and found novel de novo heterozygous missense KCNB1 mutations in two patients showing psychomotor developmental delay and severe infantile generalized seizures with high-amplitude spike-and-wave electroencephalogram discharges. The mutation located in the channel voltage sensor (p.R306C) disrupted sensitivity and cooperativity of the sensor, while the mutation in the channel pore domain (p.G401R) selectively abolished endogenous Kv2 currents in transfected pyramidal neurons, indicating a dominant-negative effect. Both mutants inhibited repetitive neuronal firing through preventing production of deep interspike voltages. Thus KCNB1 mutations can be a rare genetic cause of infantile epilepsy, and insufficient firing of pyramidal neurons would disturb both development and stability of neuronal circuits, leading to the disease phenotypes.

Fungal pathogens are a significant threat to crops worldwide. The soil fungus, Magnaporthiopsis maydis, severely affects sensitive maize hybrids by causing the rapid wilting of plants at the maturity stage. Similarly, the soil fungus, Macrophomina phaseolina, develops in a variety of host plants, which leads to rot and plant mortality. The presence of both pathogens together in diseased cotton plants in Israel suggests possible interactions between them. Here, these relationships were tested in a series of experiments accompanied by real-time PCR tracking in maize and cotton. Despite the fact that neither of the pathogens was superior in a growth plate confrontation assay, their co-inoculum had a significant influence under field conditions. In maize sprouts and fully matured plants, infection by both pathogens (compared to inoculation with each of them alone) led to lesser amounts of M. maydis DNA but to increased amounts of M. phaseolina DNA levels. These results were obtained under a restricted water regime, while optimal water irrigation led to less pronounced differences. In water-stressed cotton sprouts, infection with both pathogens led to an increase in DNA amounts of each of the pathogens. Whereas the M. maydis DNA levels in the double infection remain high at the end of the season, a reduction in the amount of M. phaseolina DNA was observed. The double infection caused an increase in growth parameters in maize and cotton and decreased levels of dehydration in maize plants accompanied by an increase in yield production. Dehydration symptoms were minor in cotton under an optimal water supply. However, under a restricted water regime, the double infection abolished the harmful effect of M. phaseolina on the plants' development and yield. These findings are the first report of interactions between these two pathogens in maize and cotton, and they encourage expanding the study to additional plant hosts and examining the potential involvement of other pathogens.

The survival and function of tissues depend on appropriate vascularization. Blood vessels of the tissues supply oxygen, and nutrients and remove waste and byproducts. Incorporating blood vessels into engineered tissues is essential for overcoming diffusion limitations, improving tissue function, and thus facilitating the fabrication of thick tissues. Here, we present a modified ECM bioink, with enhanced mechanical properties and endothelial cell-specific adhesion motifs, to serve as a building material for 3D printing of a multiscale blood vessel network. The bioink is composed of natural ECM and alginate conjugated with a laminin adhesion molecule motif (YIGSR). The hybrid hydrogel was characterized for its mechanical properties, biochemical content, and ability to interact with endothelial cells. The pristine and modified hydrogels were mixed with induced pluripotent stem cells derived endothelial cells (iPSCs-ECs) and used to print large blood vessels with capillary beds in between.