Development and large-scale production of human milk fat analog by fermentation of microalgae.

Background: Human milk contains a complex mixture of triacylglycerols (TAG), making it challenging to recreate using common ingredients. Objective: The study aimed to develop an innovative fermentation technique to produce essential human milk TAG, effectively tackling a significant hurdle in infant nutrition. Method: An in-depth analysis of the literature has been conducted to identify the specific TAG to be targeted. We used a microalgal oil production platform and a two-step procedure to modify its fatty acid and TAG composition. The palmitic acid (16:0) content has been increased by classical strain improvement techniques, followed by a step involving the expression of a lysophosphatidic acid acyltransferase (LPAAT) sequence capable of esterifying 16:0 specifically at the internal position (sn-2 palmitate) of TAG. Once the strain was stabilized, the fermentation was scaled up in a 50-L reactor to yield several kilograms of biomass. Subsequently, the oil was extracted and refined using standard oil processing conditions. Liquid chromatography-mass spectrometry was employed to monitor the TAG profile and the region specificity of 16:0 at the internal position (sn-2 palmitate) of TAG. Results: The initial strain had a 16:0 level of 25% of total fatty acids, which was increased to 30% by classical strain improvement. Simultaneously, the oleic acid level decreased from 61% to 57% of total fatty acids. Upon expression of an exogenous LPAAT gene, the level of the 16:0 esterified in the internal position of the TAG (sn-2 palmitate) increased by a factor of 10, to reach 73% of total palmitic acid. Consequently, the concentration of oleic acid in the internal position decreased from 81% to 22% of total fatty acids, with TAG analysis confirming that the primary TAG species in the oil was 1,3-dioleoyl-2-palmitoyl-glycerol (OPO). The 50-L-scale fermentation trial confirmed the strain's ability to produce oil with a yield of >150 g of oil per liter of fermentation broth in a timeframe of 5 days, rendering the process scalable for larger-scale industrialization. Conclusion: We have demonstrated the feasibility of producing a suitable TAG composition that can be effectively integrated into the formulations of infant nutrition in combination with other fats and oils to meet the infant feeding requirements.

Genome sequence for the blue-flowered Andean shrub Iochroma cyaneum reveals extensive discordance across the berry clade of Solanaceae.

The tomato (Solanum lycopersicum L.) family, Solanaceae, is a model clade for a wide range of applied and basic research questions. Currently, reference-quality genomes are available for over 30 species from seven genera, and these include numerous crops as well as wild species [e.g., Jaltomata sinuosa (Miers) Mione and Nicotiana attenuata Torr. ex S. Watson]. Here we present the genome of the showy-flowered Andean shrub Iochroma cyaneum (Lindl.) M. L. Green, a woody lineage from the tomatillo (Physalis philadelphica Lam.) subfamily Physalideae. The assembled size of the genome (2.7 Gb) is more similar in size to pepper (Capsicum annuum L.) (2.6 Gb) than to other sequenced diploid members of the berry clade of Solanaceae [e.g., potato (Solanum tuberosum L.), tomato, and Jaltomata]. Our assembly recovers 92% of the conserved orthologous set, suggesting a nearly complete genome for this species. Most of the genomic content is repetitive (69%), with Gypsy elements alone accounting for 52% of the genome. Despite the large amount of repetitive content, most of the 12 I. cyaneum chromosomes are highly syntenic with tomato. Bayesian concordance analysis provides strong support for the berry clade, including I. cyaneum, but reveals extensive discordance along the backbone, with placement of chili pepper and Jaltomata being highly variable across gene trees. The I. cyaneum genome contributes to a growing wealth of genomic resources in Solanaceae and underscores the need for expanded sampling of diverse berry genomes to dissect major morphological transitions.

An adaptive teosinte mexicana introgression modulates phosphatidylcholine levels and is associated with maize flowering time.

Native Americans domesticated maize (Zea mays ssp. mays) from lowland teosinte parviglumis (Zea mays ssp. parviglumis) in the warm Mexican southwest and brought it to the highlands of Mexico and South America where it was exposed to lower temperatures that imposed strong selection on flowering time. Phospholipids are important metabolites in plant responses to low-temperature and phosphorus availability and have been suggested to influence flowering time. Here, we combined linkage mapping with genome scans to identify High PhosphatidylCholine 1 (HPC1), a gene that encodes a phospholipase A1 enzyme, as a major driver of phospholipid variation in highland maize. Common garden experiments demonstrated strong genotype-by-environment interactions associated with variation at HPC1, with the highland HPC1 allele leading to higher fitness in highlands, possibly by hastening flowering. The highland maize HPC1 variant resulted in impaired function of the encoded protein due to a polymorphism in a highly conserved sequence. A meta-analysis across HPC1 orthologs indicated a strong association between the identity of the amino acid at this position and optimal growth in prokaryotes. Mutagenesis of HPC1 via genome editing validated its role in regulating phospholipid metabolism. Finally, we showed that the highland HPC1 allele entered cultivated maize by introgression from the wild highland teosinte Zea mays ssp. mexicana and has been maintained in maize breeding lines from the Northern United States, Canada, and Europe. Thus, HPC1 introgressed from teosinte mexicana underlies a large metabolic QTL that modulates phosphatidylcholine levels and has an adaptive effect at least in part via induction of early flowering time.

Selective sorting of ancestral introgression in maize and teosinte along an elevational cline.

While often deleterious, hybridization can also be a key source of genetic variation and pre-adapted haplotypes, enabling rapid evolution and niche expansion. Here we evaluate these opposing selection forces on introgressed ancestry between maize (Zea mays ssp. mays) and its wild teosinte relative, mexicana (Zea mays ssp. mexicana). Introgression from ecologically diverse teosinte may have facilitated maize's global range expansion, in particular to challenging high elevation regions (> 1500 m). We generated low-coverage genome sequencing data for 348 maize and mexicana individuals to evaluate patterns of introgression in 14 sympatric population pairs, spanning the elevational range of mexicana, a teosinte endemic to the mountains of Mexico. While recent hybrids are commonly observed in sympatric populations and mexicana demonstrates fine-scale local adaptation, we find that the majority of mexicana ancestry tracts introgressed into maize over 1000 generations ago. This mexicana ancestry seems to have maintained much of its diversity and likely came from a common ancestral source, rather than contemporary sympatric populations, resulting in relatively low FST between mexicana ancestry tracts sampled from geographically distant maize populations. Introgressed mexicana ancestry in maize is reduced in lower-recombination rate quintiles of the genome and around domestication genes, consistent with pervasive selection against introgression. However, we also find mexicana ancestry increases across the sampled elevational gradient and that high introgression peaks are most commonly shared among high-elevation maize populations, consistent with introgression from mexicana facilitating adaptation to the highland environment. In the other direction, we find patterns consistent with adaptive and clinal introgression of maize ancestry into sympatric mexicana at many loci across the genome, suggesting that maize also contributes to adaptation in mexicana, especially at the lower end of its elevational range. In sympatric maize, in addition to high introgression regions we find many genomic regions where selection for local adaptation maintains steep gradients in introgressed mexicana ancestry across elevation, including at least two inversions: the well-characterized 14 Mb Inv4m on chromosome 4 and a novel 3 Mb inversion Inv9f surrounding the macrohairless1 locus on chromosome 9. Most outlier loci with high mexicana introgression show no signals of sweeps or local sourcing from sympatric populations and so likely represent ancestral introgression sorted by selection, resulting in correlated but distinct outcomes of introgression in different contemporary maize populations.

Conserved noncoding sequences provide insights into regulatory sequence and loss of gene expression in maize.

Thousands of species will be sequenced in the next few years; however, understanding how their genomes work, without an unlimited budget, requires both molecular and novel evolutionary approaches. We developed a sensitive sequence alignment pipeline to identify conserved noncoding sequences (CNSs) in the Andropogoneae tribe (multiple crop species descended from a common ancestor ∼18 million years ago). The Andropogoneae share similar physiology while being tremendously genomically diverse, harboring a broad range of ploidy levels, structural variation, and transposons. These contribute to the potential of Andropogoneae as a powerful system for studying CNSs and are factors we leverage to understand the function of maize CNSs. We found that 86% of CNSs were comprised of annotated features, including introns, UTRs, putative cis-regulatory elements, chromatin loop anchors, noncoding RNA (ncRNA) genes, and several transposable element superfamilies. CNSs were enriched in active regions of DNA replication in the early S phase of the mitotic cell cycle and showed different DNA methylation ratios compared to the genome-wide background. More than half of putative cis-regulatory sequences (identified via other methods) overlapped with CNSs detected in this study. Variants in CNSs were associated with gene expression levels, and CNS absence contributed to loss of gene expression. Furthermore, the evolution of CNSs was associated with the functional diversification of duplicated genes in the context of maize subgenomes. Our results provide a quantitative understanding of the molecular processes governing the evolution of CNSs in maize.

Gene regulatory effects of a large chromosomal inversion in highland maize.

Chromosomal inversions play an important role in local adaptation. Inversions can capture multiple locally adaptive functional variants in a linked block by repressing recombination. However, this recombination suppression makes it difficult to identify the genetic mechanisms underlying an inversion's role in adaptation. In this study, we used large-scale transcriptomic data to dissect the functional importance of a 13 Mb inversion locus (Inv4m) found almost exclusively in highland populations of maize (Zea mays ssp. mays). Inv4m was introgressed into highland maize from the wild relative Zea mays ssp. mexicana, also present in the highlands of Mexico, and is thought to be important for the adaptation of these populations to cultivation in highland environments. However, the specific genetic variants and traits that underlie this adaptation are not known. We created two families segregating for the standard and inverted haplotypes of Inv4m in a common genetic background and measured gene expression effects associated with the inversion across 9 tissues in two experimental conditions. With these data, we quantified both the global transcriptomic effects of the highland Inv4m haplotype, and the local cis-regulatory variation present within the locus. We found diverse physiological effects of Inv4m across the 9 tissues, including a strong effect on the expression of genes involved in photosynthesis and chloroplast physiology. Although we could not confidently identify the causal alleles within Inv4m, this research accelerates progress towards understanding this inversion and will guide future research on these important genomic features.

Filtering of target sequence capture individuals facilitates species tree construction in the plant subtribe Iochrominae (Solanaceae).

Advances in sequencing technology have made it possible to produce large multi-locus datasets required for species tree analyses. One challenge with constructing high throughput sequencing datasets, however, is that missing information is propagated at different steps in the sequence preparation process. To date, species tree studies have focused on filtering and removing errors that occur at particular loci. Given the way that high throughput sequencing datasets are constructed, however, large amounts of error or ambiguity may also manifest across individuals. Here we use a novel tree-based multivariate clustering method to identify and remove individuals with low phylogenetic signal in a nuclear sequence capture dataset for the Iochrominae clade (Solanaceae). Our results suggest that the low quality tips are the result of the library preparation process (e.g. unequal pooling) rather than poor capture due to phylogenetic distance from the reference species. After implementing the clustering approach and removing low quality tips, we construct an Iochrominae species tree that resolves a number of unknown relationships. We propose this pipeline as a valuable tool for species tree reconstruction with phylogenomic datasets containing variable levels of missing data.

A novel R3 MYB transcriptional repressor associated with the loss of floral pigmentation in Iochroma.

Losses of floral pigmentation represent one of the most common evolutionary transitions in flower color, yet the genetic basis for these changes has been elucidated in only a handful of cases. Here we used crossing studies, bulk-segregant RNA sequencing, phylogenetic analyses and functional tests to identify the gene(s) responsible for the transition to white flowers in Iochroma loxense. Crosses between I. loxense and its blue-flowered sister species, I. cyaneum, suggested that a single locus controls the flower color difference and that the white allele causes a nearly complete loss of pigmentation. Examining sequence variation across phenotypic pools from the crosses, we found that alleles at a novel R3 MYB transcription factor were tightly associated with flower color variation. This gene, which we term MYBL1, falls into a class of MYB transcriptional repressors and, accordingly, higher expression of this gene is associated with downregulation of multiple anthocyanin pigment pathway genes. We confirmed the repressive function of MYBL1 through stable transformation of Nicotiana. The mechanism underlying the evolution of white flowers in I. loxense differs from that uncovered in previous studies, pointing to multiple mechanisms for achieving fixed transitions in flower color intensity.

How to make a domesticate.

The Neolithic Revolution brought about the transition from hunting and gathering to sedentary societies, laying the foundation for the development of modern civilizations. The primary innovation that facilitated these changes was the domestication of plants and animals. In the case of plants, this involved the cultivation and selection of individuals with larger edible parts, easier harvesting, and decreased defenses, traits that allowed for the production of a food surplus and occupational specialization. Plant domestication is a process which started approximately 10,000 years ago and has thereafter been repeated independently in many locales around the world. Here, we offer a perspective that seeks to predict what factors influence the success of domestication, how many genes contributed to the process, where these genes originated and the implications for de novo domestication.

Diversification of R2R3-MYB Transcription Factors in the Tomato Family Solanaceae.

MYB transcription factors play an important role in regulating key plant developmental processes involving defense, cell shape, pigmentation, and root formation. Within this gene family, sequences containing an R2R3 MYB domain are the most abundant type and exhibit a wide diversity of functions. In this study, we identify 559 R2R3 MYB genes using whole genome data from four species of Solanaceae and reconstruct their evolutionary relationships. We compare the Solanaceae R2R3 MYBs to the well-characterized Arabidopsis thaliana sequences to estimate functional diversity and to identify gains and losses of MYB clades in the Solanaceae. We identify numerous R2R3 MYBs that do not appear closely related to Arabidopsis MYBs, and thus may represent clades of genes that have been lost along the Arabidopsis lineage or gained after the divergence of Rosid and Asterid lineages. Despite differences in the distribution of R2R3 MYBs across functional subgroups and species, the overall size of the R2R3 subfamily has changed relatively little over the roughly 50 million-year history of Solanaceae. We added our information regarding R2R3 MYBs in Solanaceae to other data and performed a meta-analysis to trace the evolution of subfamily size across land plants. The results reveal many shifts in the number of R2R3 genes, including a 54 % increase along the angiosperm stem lineage. The variation in R2R3 subfamily size across land plants is weakly positively correlated with genome size and strongly positively correlated with total number of genes. The retention of such a large number of R2R3 copies over long evolutionary time periods suggests that they have acquired new functions and been maintained by selection. Discovering the nature of this functional diversity will require integrating forward and reverse genetic approaches on an -omics scale.

A prospective study of the incidence of myocarditis/pericarditis and new onset cardiac symptoms following smallpox and influenza vaccination.

BACKGROUND: Although myocarditis/pericarditis (MP) has been identified as an adverse event following smallpox vaccine (SPX), the prospective incidence of this reaction and new onset cardiac symptoms, including possible subclinical injury, has not been prospectively defined. PURPOSE: The study's primary objective was to determine the prospective incidence of new onset cardiac symptoms, clinical and possible subclinical MP in temporal association with immunization. METHODS: New onset cardiac symptoms, clinical MP and cardiac specific troponin T (cTnT) elevations following SPX (above individual baseline values) were measured in a multi-center prospective, active surveillance cohort study of healthy subjects receiving either smallpox vaccine or trivalent influenza vaccine (TIV). RESULTS: New onset chest pain, dyspnea, and/or palpitations occurred in 10.6% of SPX-vaccinees and 2.6% of TIV-vaccinees within 30 days of immunization (relative risk (RR) 4.0, 95% CI: 1.7-9.3). Among the 1081 SPX-vaccinees with complete follow-up, 4 Caucasian males were diagnosed with probable myocarditis and 1 female with suspected pericarditis. This indicates a post-SPX incidence rate more than 200-times higher than the pre-SPX background population surveillance rate of myocarditis/pericarditis (RR 214, 95% CI 65-558). Additionally, 31 SPX-vaccinees without specific cardiac symptoms were found to have over 2-fold increases in cTnT (>99th percentile) from baseline (pre-SPX) during the window of risk for clinical myocarditis/pericarditis and meeting a proposed case definition for possible subclinical myocarditis. This rate is 60-times higher than the incidence rate of overt clinical cases. No clinical or possible subclinical myocarditis cases were identified in the TIV-vaccinated group. CONCLUSIONS: Passive surveillance significantly underestimates the true incidence of myocarditis/pericarditis after smallpox immunization. Evidence of subclinical transient cardiac muscle injury post-vaccinia immunization is a finding that requires further study to include long-term outcomes surveillance. Active safety surveillance is needed to identify adverse events that are not well understood or previously recognized.

Flowering asynchrony and mating system effects on reproductive assurance and mutualism persistence in fragmented fig-fig wasp populations.

PREMISE OF THE STUDY: Plants and animals may experience reproductive Allee effects in fragmented populations, and obligate pollination mutualisms may be especially sensitive to extinction risk via this density-dependent process. In this study we examine how a shift from within-crown reproductive synchrony to asynchrony influences reproductive assurance through contributions to selfing and outcrossing in small, spatially isolated populations of Ficus. METHODS: The research focuses on the monoecious fig F. petiolaris and consists of phenological analyses and genetic assessments of selfing and outcrossing for populations located in Baja California's Sonoran Desert. KEY RESULTS: Phenological censuses of eight populations revealed within-crown asynchrony in 44% of reproducing trees, with 16% having sufficient overlap of male and female flowering phases to permit selfing via the cycling of pollinating fig wasps within natal trees. In mating system analyses of two of these populations, however, multilocus outcrossing rates (t(m)) were indistinguishable from 1. This result, combined with low levels of inbreeding, indicates selfing to be absent or at best a minor contributor to reproductive assurance. CONCLUSIONS: The results indicate that the fitness benefits of within-crown asynchrony lie not with selfing, as commonly asserted, but with increased opportunities for outcross pollen transmission and receipt, changing our understanding of the mechanisms by which reproduction is facilitated and extinction risk minimized in naturally fragmented Ficus populations. Given the role of fig fruit as a keystone food resource in many tropical environments, trait variation leading to reproductive assurance in figs, such as within-crown asynchrony, has broader ecosystem-level implications.

The relationship between subclinical atherosclerosis and electrocardiographic abnormalities as biomarkers of cardiovascular risk.

Electrocardiographic findings indicating myocardial disease, such as left ventricular hypertrophy or ST-T wave abnormalities, or the presence of coronary artery calcium, indicating atherosclerotic coronary artery disease, are both biomarkers of future cardiovascular (CV) risk. Although the risk factors for myocardial and coronary artery disease are similar, their concomitant expression has implications for CV disease screening and prevention programmes. The relationship between the resting 12-lead ECG and subclinical atherosclerosis measured as coronary artery calcium (CAC) with electron beam tomography was examined in 937 healthy participants (aged 40-50 years) enrolled in a CV risk screening study. Electrocardiograms and CAC were interpreted in blinded fashion, using standard criteria. An abnormal ECG was coded in 268 (28.6%) participants, most commonly left ventricular hypertrophy (3.1%), delayed precordial R wave transition (5.7%), T-wave abnormalities (10.0%) and intraventricular conduction delay (10.4%). Although abnormal ECG findings were associated with CV risk variables, the prevalence of any CAC was similar in subjects with any ECG finding (43 of 268, 16.0%) compared with those with normal ECGs (125 of 669, 18.7%, p =NS). In a logistic model controlling for CV risk factors including systolic blood pressure, low-density lipoprotein cholesterol (LDL-C), body mass index (BMI), glycosylated haemoglobin, race, age and gender, significant associations with CAC were found for LDL-C, race and BMI. There was no significant relationship between CAC and ECG abnormalities (odds ratio 0.80, 95% confidence interval 0.54-1.20). In conclusion, electrocardiographic abnormalities and subclinical calcified atherosclerosis were not significantly associated with each other in this middle-aged screening population. This suggests these two biomarkers may be complementary towards broader detection of latent CV risk.