Complementing model species with model clades.

Model species continue to underpin groundbreaking plant science research. At the same time, the phylogenetic resolution of the land plant Tree of Life continues to improve. The intersection of these two research paths creates a unique opportunity to further extend the usefulness of model species across larger taxonomic groups. Here we promote the utility of the Arabidopsis thaliana model species, especially the ability to connect its genetic and functional resources, to species across the entire Brassicales order. We focus on the utility of using genomics and phylogenomics to bridge the evolution and diversification of several traits across the Brassicales to the resources in Arabidopsis, thereby extending scope from a model species by establishing a "model clade". These Brassicales-wide traits are discussed in the context of both the model species Arabidopsis thaliana and the family Brassicaceae. We promote the utility of such a "model clade" and make suggestions for building global networks to support future studies in the model order Brassicales.

Genome balance and dosage effect drive allopolyploid formation in Brassica.

Polyploidy is a major evolutionary force that has shaped plant diversity. However, the various pathways toward polyploid formation and interploidy gene flow remain poorly understood. Here, we demonstrated that the immediate progeny of allotriploid AAC Brassica (obtained by crossing allotetraploid Brassica napus and diploid Brassica rapa) was predominantly aneuploids with ploidal levels ranging from near-triploidy to near-hexaploidy, and their chromosome numbers deviated from the theoretical distribution toward increasing chromosome numbers, suggesting that they underwent selection. Karyotype and phenotype analyses showed that aneuploid individuals containing fewer imbalanced chromosomes had higher viability and fertility. Within three generations of self-fertilization, allotriploids mainly developed into near or complete allotetraploids similar to B. napus via gradually increasing chromosome numbers and fertility, suggesting that allotriploids could act as a bridge in polyploid formation, with aneuploids as intermediates. Self-fertilized interploidy hybrids ultimately generated new allopolyploids carrying different chromosome combinations, which may create a reproductive barrier preventing allotetraploidy back to diploidy and promote gene flow from diploids to allotetraploids. These results suggest that the maintenance of a proper genome balance and dosage drove the recurrent conversion of allotriploids to allotetraploids, which may contribute to the formation and evolution of polyploids.

The final piece of the Triangle of U: Evolution of the tetraploid Brassica carinata genome.

Ethiopian mustard (Brassica carinata) is an ancient crop with remarkable stress resilience and a desirable seed fatty acid profile for biofuel uses. Brassica carinata is one of six Brassica species that share three major genomes from three diploid species (AA, BB, and CC) that spontaneously hybridized in a pairwise manner to form three allotetraploid species (AABB, AACC, and BBCC). Of the genomes of these species, that of B. carinata is the least understood. Here, we report a chromosome scale 1.31-Gbp genome assembly with 156.9-fold sequencing coverage for B. carinata, completing the reference genomes comprising the classic Triangle of U, a classical theory of the evolutionary relationships among these six species. Our assembly provides insights into the hybridization event that led to the current B. carinata genome and the genomic features that gave rise to the superior agronomic traits of B. carinata. Notably, we identified an expansion of transcription factor networks and agronomically important gene families. Completion of the Triangle of U comparative genomics platform has allowed us to examine the dynamics of polyploid evolution and the role of subgenome dominance in the domestication and continuing agronomic improvement of B. carinata and other Brassica species.

Gene expression bias between the subgenomes of allopolyploid hybrids is an emergent property of the kinetics of expression.

Hybridization coupled to polyploidy, or allopolyploidy, has dramatically shaped the evolution of flowering plants, teleost fishes, and other lineages. Studies of recently formed allopolyploid plants have shown that the two subgenomes that merged to form that new allopolyploid do not generally express their genes equally. Instead, one of the two subgenomes expresses its paralogs more highly on average. Meanwhile, older allopolyploidy events tend to show biases in duplicate losses, with one of the two subgenomes retaining more genes than the other. Since reduced expression is a pathway to duplicate loss, understanding the origins of expression biases may help explain the origins of biased losses. Because we expect gene expression levels to experience stabilizing selection, our conceptual frameworks for how allopolyploid organisms form tend to assume that the new allopolyploid will show balanced expression between its subgenomes. It is then necessary to invoke phenomena such as differences in the suppression of repetitive elements to explain the observed expression imbalances. Here we show that, even for phenotypically identical diploid progenitors, the inherent kinetics of gene expression give rise to biases between the expression levels of the progenitor genes in the hybrid. Some of these biases are expected to be gene-specific and not give rise to global differences in progenitor gene expression. However, particularly in the case of allopolyploids formed from progenitors with different genome sizes, global expression biases favoring one subgenome are expected immediately on formation. Hence, expression biases are arguably the expectation upon allopolyploid formation rather than a phenomenon needing explanation. In the future, a deeper understanding of the kinetics of allopolyploidy may allow us to better understand both biases in duplicate losses and hybrid vigor.

Feralization: Confronting the Complexity of Domestication and Evolution.

Feral populations, those which successfully persist outside of cultivation or husbandry, provide unique opportunities to study the genomic impacts of domestication and local adaptation. We argue that by leveraging genomic resources designed for domestic counterparts, powerful phylogenetic and population genomic data collection and analyses can be designed to disentangle complex demographic processes.

The contributions from the progenitor genomes of the mesopolyploid Brassiceae are evolutionarily distinct but functionally compatible.

The members of the tribe Brassiceae share a whole-genome triplication (WGT), and one proposed model for its formation is a two-step pair of hybridizations producing hexaploid descendants. However, evidence for this model is incomplete, and the evolutionary and functional constraints that drove evolution after the hexaploidy are even less understood. Here, we report a new genome sequence of Crambe hispanica, a species sister to most sequenced Brassiceae. Using this new genome and three others that share the hexaploidy, we traced the history of gene loss after the WGT using the Polyploidy Orthology Inference Tool (POInT). We confirm the two-step formation model and infer that there was a significant temporal gap between those two allopolyploidizations, with about a third of the gene losses from the first two subgenomes occurring before the arrival of the third. We also, for the 90,000 individual genes in our study, make parental subgenome assignments, inferring, with measured uncertainty, from which of the progenitor genomes of the allohexaploidy each gene derives. We further show that each subgenome has a statistically distinguishable rate of homoeolog losses. There is little indication of functional distinction between the three subgenomes: the individual subgenomes show no patterns of functional enrichment, no excess of shared protein-protein or metabolic interactions between their members, and no biases in their likelihood of having experienced a recent selective sweep. We propose a "mix and match" model of allopolyploidy, in which subgenome origin drives homoeolog loss propensities but where genes from different subgenomes function together without difficulty.

Robust Yet Fragile: Expression Noise, Protein Misfolding, and Gene Dosage in the Evolution of Genomes.

The complex manner in which organisms respond to changes in their gene dosage has long fascinated geneticists. Oddly, although the existence of dominance implies that dosage reductions often have mild phenotypes, extra copies of whole chromosomes (aneuploidy) are generally strongly deleterious. Even more paradoxically, an extra copy of the genome is better tolerated than is aneuploidy. We review the resolution of this paradox, highlighting the roles of biochemistry, protein aggregation, and disruption of cellular microstructure in that explanation. Returning to life's curious combination of robustness and sensitivity to dosage changes, we argue that understanding how biological robustness evolved makes these observations less inexplicable. We propose that noise in gene expression and evolutionary strategies for its suppression play a role in generating dosage phenotypes. Finally, we outline an unappreciated mechanism for the preservation of duplicate genes, namely preservation to limit expression noise, arguing that it is particularly relevant in polyploid organisms.

Halophytes and heavy metals: A multi-omics approach to understand the role of gene and genome duplication in the abiotic stress tolerance of Cakile maritima.

PREMISE: The origin of diversity is a fundamental biological question. Gene duplications are one mechanism that provides raw material for the emergence of novel traits, but evolutionary outcomes depend on which genes are retained and how they become functionalized. Yet, following different duplication types (polyploidy and tandem duplication), the events driving gene retention and functionalization remain poorly understood. Here we used Cakile maritima, a species that is tolerant to salt and heavy metals and shares an ancient whole-genome triplication with closely related salt-sensitive mustard crops (Brassica), as a model to explore the evolution of abiotic stress tolerance following polyploidy. METHODS: Using a combination of ionomics, free amino acid profiling, and comparative genomics, we characterize aspects of salt stress response in C. maritima and identify retained duplicate genes that have likely enabled adaptation to salt and mild levels of cadmium. RESULTS: Cakile maritima is tolerant to both cadmium and salt treatments through uptake of cadmium in the roots. Proline constitutes greater than 30% of the free amino acid pool in C. maritima and likely contributes to abiotic stress tolerance. We find duplicated gene families are enriched in metabolic and transport processes and identify key transport genes that may be involved in C. maritima abiotic stress tolerance. CONCLUSIONS: These findings identify pathways and genes that could be used to enhance plant resilience and provide a putative understanding of the roles of duplication types and retention on the evolution of abiotic stress response.

Responses of selected plasma metabolites to a 2-day nutritional challenge of goats divergently selected for functional longevity.

Understanding the extent to which genetics × environment plays a role in shaping individual strategies to environmental challenges is of considerable interest for future selection of more resilient animals. The objective of this study was to evaluate the metabolic responses of goats, selected for functional longevity, to a nutritional challenge. This evaluation was based on plasma metabolites and was repeated across 2 experimental farms over 2 years. We carried out 6 different experimental trials from 2018 to 2022 (4 trials on site Bourges, 2018-2021 and 2 trials on site Grignon, 2021-2022) in which 267 first kidding goats, daughters of Alpine bucks divergently selected for functional longevity, longevity plus (n = 137), and longevity minus (n = 130), were exposed to a 2-d nutritional challenge in early lactation. The experiments consisted of a 5 or 7-d control period (prechallenge) on a standard lactation diet followed by a 2-d nutritional challenge with straw-only feeding, and then a 7- or 10-d recovery period on a standard lactation diet, for site Bourges and Grignon, respectively. During the challenge plasma metabolite composition was measured. Linear mixed models were used to analyze all traits, considering the individual as a random effect and the 2 × 2 treatments (i.e., genetic line and year nested in site) and litter size as fixed effects. The linear mixed-effects model using a piecewise arrangement was used to analyze the response/recovery profiles to the nutritional challenge. The estimates of the random effect for each individual, using the mixed models without the fixed effects of genetic line, were used in a sparse partial least square discriminant analysis (sPLS-DA) to compare the goat metabolism response to the challenge on a multivariate scale. The plasma glucose, BHB, nonesterified fatty acids (NEFA), and urea concentrations significantly responded to the 2-d nutritional challenge. Selection for functional longevity did not affect plasma glucose, NEFA, BHB, and urea response/recoveries to a 2-d nutritional challenge. However, site, trial, and litter size significantly affected these responses. Moreover, the plasma metabolites seem not to fully recover to prechallenge levels after the recovery phase. The sPLS-DA analysis did not discriminate between the 2 longevity lines. We observed meaningful between-individuals variability in plasma BHB, especially on the prechallenge and rate of response and rate of recovery from the 2-d nutritional challenge (CV = 26.2%, 36.1%, and 41.2%, repeatability = 0.749, 0.322, and 0.741, respectively). Plasma NEFA recovery from challenge also demonstrated high between-individuals variability (CV = 16.4%, repeatability = 0.323). Selection for functional longevity did not affect plasma metabolites responses to a 2-d nutritional challenge in dairy goats. Plasma NEFA and BHB response and recovery presented high between-individuals variability, indicating that individual adaptative characteristics to nutritional challenges are not related to the environmental conditions but to inherent individual characteristics.

MiRNome variations in milk fractions during feed restrictions of different intensities in dairy cows.

BACKGROUND: In dairy cows, diet is one factor that can affect their milk production and composition. However, the effect of feed restriction on milk miRNome has not yet been described. Indeed, milk is the body fluid with the highest RNA concentration, which includes numerous microRNA. Its presence in the four different milk fractions, whole milk, fat globules, mammary epithelial cells and extracellular vesicles, is still poorly documented. This study aimed to describe the effects of different feed restrictions on the miRNome composition of different milk fractions. RESULTS: Two feed restrictions were applied to lactating dairy cows, one of high intensity and one of moderate intensity. 2,896 mature microRNA were identified in the different milk fractions studied, including 1,493 that were already known in the bovine species. Among the 1,096 microRNA that were sufficiently abundant to be informative, the abundance of 1,027 of them varied between fractions: 36 of those were exclusive to one milk fraction. Feed restriction affected the abundance of 155 microRNA, with whole milk and milk extracellular vesicles being the most affected, whereas milk fat globules and exfoliated mammary epithelial cells were little or not affected at all. The high intensity feed restriction led to more microRNA variations in milk than moderate restriction. The target prediction of known microRNA that varied under feed restriction suggested the modification of some key pathways for lactation related to milk fat and protein metabolisms, cell cycle, and stress responses. CONCLUSIONS: This study highlighted that the miRNome of each milk fraction is specific, with mostly the same microRNA composition but with variations in abundance between fractions. These specific miRNomes were affected differently by feed restrictions, the intensity of which appeared to be a major factor modulating milk miRNomes. These findings offer opportunities for future research on the use of milk miRNA as biomarkers of energy status in dairy cows, which is affected by feed restrictions.

Patterns, mechanisms, and consequences of homoeologous exchange in allopolyploid angiosperms: a genomic and epigenomic perspective.

Allopolyploids result from hybridization between different evolutionary lineages coupled with genome doubling. Homoeologous chromosomes (chromosomes with common shared ancestry) may undergo recombination immediately after allopolyploid formation and continue over successive generations. The outcome of this meiotic pairing behavior is dynamic and complex. Homoeologous exchanges (HEs) may lead to the formation of unbalanced gametes, reduced fertility, and selective disadvantage. By contrast, HEs could act as sources of novel evolutionary substrates, shifting the relative dosage of parental gene copies, generating novel phenotypic diversity, and helping the establishment of neo-allopolyploids. However, HE patterns vary among lineages, across generations, and even within individual genomes and chromosomes. The causes and consequences of this variation are not fully understood, though interest in this evolutionary phenomenon has increased in the last decade. Recent technological advances show promise in uncovering the mechanistic basis of HEs. Here, we describe recent observations of the common patterns among allopolyploid angiosperm lineages, underlying genomic and epigenomic features, and consequences of HEs. We identify critical research gaps and discuss future directions with far-reaching implications in understanding allopolyploid evolution and applying them to the development of important phenotypic traits of polyploid crops.

Mangroves and multiplications: Influence of genome duplications on salt tolerance.

Whole-genome duplication (WGD) is an important force in plant diversification and novel environment adaptation. Various hypotheses have been proposed on the mechanism of how WGD influences this evolutionary process from chromosome recombination to genetic diversity (Qi et al., 2021). In this issue of Molecular Ecology, Xu et al. (2022) conduct a comparative transcriptome study across species to investigate the effects of WGD on plant adaptation and the differentiation of paralogues within two distantly related mangrove species (Xu et al., 2021). They identify retained WGD-derived duplicate pairs, assess their contributions to the salt tolerance phenotype, and propose a model for adaptation expansion in mangroves through the neofunctionalization of WGD-derived duplicate genes. Their results show that neofunctionalized copies of known salt tolerance-related genes show upregulated expression in fluctuating salinity (i.e., fresh and hypersaline) compared to the optimal medium salinity environment. This suggests that after WGD, retained duplicates increase genomic plasticity allowing for adaptation to new and unstable environments.

Multivariate analysis of milk metabolite measures shows potential for deriving new resilience phenotypes.

In a context of growing interest in breeding more resilient animals, a noninvasive indicator of resilience would be very valuable. We hypothesized that the time-course of concentrations of several milk metabolites through a short-term underfeeding challenge could reflect the variation of resilience mechanisms to such a challenge. We submitted 138 one-year-old primiparous goats, selected for extreme functional longevity (i.e., productive longevity corrected for milk yield [60 low longevity line goats and 78 high longevity line goats]), to a 2-d underfeeding challenge during early lactation. We measured the concentration of 13 milk metabolites and the activity of 1 enzyme during prechallenge, challenge, and recovery periods. Functional principal component analysis summarized the trends of milk metabolite concentration over time efficiently without preliminary assumptions concerning the shapes of the curves. We first ran a supervised prediction of the longevity line of the goats based on the milk metabolite curves. The partial least square analysis could not predict the longevity line accurately. We thus decided to explore the large overall variability of milk metabolite curves with an unsupervised clustering. The large year × facility effect on the metabolite concentrations was precorrected for. This resulted in 3 clusters of goats defined by different metabolic responses to underfeeding. The cluster that showed higher ß-hydroxybutyrate, cholesterol, and triacylglycerols increase during the underfeeding challenge was associated with poorer survival compared with the other 2 clusters. These results suggest that multivariate analysis of noninvasive milk measures show potential for deriving new resilience phenotypes.

Brassica rapa Domestication: Untangling Wild and Feral Forms and Convergence of Crop Morphotypes.

The study of domestication contributes to our knowledge of evolution and crop genetic resources. Human selection has shaped wild Brassica rapa into diverse turnip, leafy, and oilseed crops. Despite its worldwide economic importance and potential as a model for understanding diversification under domestication, insights into the number of domestication events and initial crop(s) domesticated in B. rapa have been limited due to a lack of clarity about the wild or feral status of conspecific noncrop relatives. To address this gap and reconstruct the domestication history of B. rapa, we analyzed 68,468 genotyping-by-sequencing-derived single nucleotide polymorphisms for 416 samples in the largest diversity panel of domesticated and weedy B. rapa to date. To further understand the center of origin, we modeled the potential range of wild B. rapa during the mid-Holocene. Our analyses of genetic diversity across B. rapa morphotypes suggest that noncrop samples from the Caucasus, Siberia, and Italy may be truly wild, whereas those occurring in the Americas and much of Europe are feral. Clustering, tree-based analyses, and parameterized demographic inference further indicate that turnips were likely the first crop type domesticated, from which leafy types in East Asia and Europe were selected from distinct lineages. These findings clarify the domestication history and nature of wild crop genetic resources for B. rapa, which provides the first step toward investigating cases of possible parallel selection, the domestication and feralization syndrome, and novel germplasm for Brassica crop improvement.

The Evolutionary History of Wild, Domesticated, and Feral Brassica oleracea (Brassicaceae).

Understanding the evolutionary history of crops, including identifying wild relatives, helps to provide insight for conservation and crop breeding efforts. Cultivated Brassica oleracea has intrigued researchers for centuries due to its wide diversity in forms, which include cabbage, broccoli, cauliflower, kale, kohlrabi, and Brussels sprouts. Yet, the evolutionary history of this species remains understudied. With such different vegetables produced from a single species, B. oleracea is a model organism for understanding the power of artificial selection. Persistent challenges in the study of B. oleracea include conflicting hypotheses regarding domestication and the identity of the closest living wild relative. Using newly generated RNA-seq data for a diversity panel of 224 accessions, which represents 14 different B. oleracea crop types and nine potential wild progenitor species, we integrate phylogenetic and population genetic techniques with ecological niche modeling, archaeological, and literary evidence to examine relationships among cultivars and wild relatives to clarify the origin of this horticulturally important species. Our analyses point to the Aegean endemic B. cretica as the closest living relative of cultivated B. oleracea, supporting an origin of cultivation in the Eastern Mediterranean region. Additionally, we identify several feral lineages, suggesting that cultivated plants of this species can revert to a wild-like state with relative ease. By expanding our understanding of the evolutionary history in B. oleracea, these results contribute to a growing body of knowledge on crop domestication that will facilitate continued breeding efforts including adaptation to changing environmental conditions.

Prospects of Feral Crop De Novo Redomestication.

Modern agriculture depends on a narrow variety of crop species, leaving global food and nutritional security highly vulnerable to the adverse effects of climate change and population expansion. Crop improvement using conventional and molecular breeding approaches leveraging plant genetic diversity using crop wild relatives (CWRs) has been one approach to address these issues. However, the rapid pace of the global change requires additional innovative solutions to adapt agriculture to meet global needs. Neodomestication-the rapid and targeted introduction of domestication traits using introgression or genome editing of CWRs-is being explored as a supplementary approach. These methods show promise; however, they have so far been limited in efficiency and applicability. We propose expanding the scope of neodomestication beyond truly wild CWRs to include feral crops as a source of genetic diversity for novel crop development, in this case 'redomestication'. Feral crops are plants that have escaped cultivation and evolved independently, typically adapting to their local environments. Thus, feral crops potentially contain valuable adaptive features while retaining some domestication traits. Due to their genetic proximity to crop species, feral crops may be easier targets for de novo domestication (i.e. neodomestication via genome editing techniques). In this review, we explore the potential of de novo redomestication as an application for novel crop development by genome editing of feral crops. This approach to efficiently exploit plant genetic diversity would access an underutilized reservoir of genetic diversity that could prove important in support of global food insecurity in the face of the climate change.

Nonadditive gene expression in polyploids.

Allopolyploidy involves hybridization and duplication of divergent parental genomes and provides new avenues for gene expression. The expression levels of duplicated genes in polyploids can show deviation from parental additivity (the arithmetic average of the parental expression levels). Nonadditive expression has been widely observed in diverse polyploids and comprises at least three possible scenarios: (a) The total gene expression level in a polyploid is similar to that of one of its parents (expression-level dominance); (b) total gene expression is lower or higher than in both parents (transgressive expression); and (c) the relative contribution of the parental copies (homeologs) to the total gene expression is unequal (homeolog expression bias). Several factors may result in expression nonadditivity in polyploids, including maternal-paternal influence, gene dosage balance, cis- and/or trans-regulatory networks, and epigenetic regulation. As our understanding of nonadditive gene expression in polyploids remains limited, a new generation of investigators should explore additional phenomena (i.e., alternative splicing) and use other high-throughput "omics" technologies to measure the impact of nonadditive expression on phenotype, proteome, and metabolome.

Anomalous Transport Signatures in Weyl Semimetals with Point Defects.

We present the first theoretical study of transport properties of Weyl semimetals with point defects. Focusing on a class of time-reversal symmetric Weyl lattice models, we show that dilute lattice vacancies induce a finite density of quasilocalized states at and near the nodal energy, causing strong modifications to the low-energy spectrum. This generates novel transport effects, namely, (i) an oscillatory behavior of the dc conductivity with the charge carrier density in the absence of magnetic fields, and (ii) a plateau-shaped dissipative optical response for photon frequencies below the interband threshold, E_{F}≲ℏω≲2E_{F}. Our results provide a path to engineer unconventional quantum transport effects in Weyl semimetals by means of common point defects.

Milk metabolites and fatty acids as noninvasive biomarkers of metabolic status and energy balance in early-lactation cows.

The objective was to study the effects of week of lactation (WOL) and experimental nutrient restriction on concentrations of selected milk metabolites and fatty acids (FA), and assess their potential as biomarkers of energy status in early-lactation cows. To study WOL effects, 17 multiparous Holstein cows were phenotyped from calving until 7 WOL while allowed ad libitum intake of a lactation diet. Further, to study the effects of nutrient restriction, 8 of these cows received a diet containing 48% straw (high-straw) for 4 d starting at 24 ± 3 days in milk (mean ± SD), and 8 cows maintained on the lactation diet were sampled to serve as controls. Blood and milk samples were collected weekly for the WOL data set, and daily from d -1 to 3 of nutrient restriction (or control) for the nutritional challenge data set. Milk ß-hydroxybutyrate (BHB), isocitrate, glucose, glucose-6-phosphate (glucose-6P), galactose, glutamate, creatinine, uric acid, and N-acetyl-ß-d-glucosaminidase activity (NAGase) were analyzed in p.m. and a.m. samples, and milk FA were analyzed in pooled p.m. and a.m. samples. Average energy balance (EB) per day ranged from -27 MJ/d to neutral when cows received the lactation total mixed ration, and from -109 to -87 ± 7 MJ/d for high-straw (least squares means ± standard error of the mean). Plasma nonesterified FA concentration was 1.67 ± 0.13 mM and BHB was 2.96 ± 0.39 mM on the d 3 of high-straw (least squares means ± standard error of the mean). Milk concentrations of BHB, glucose, glucose-6P, glutamate, and uric acid differed significantly between p.m. and a.m. milkings. Milk isocitrate, glucose-6P, creatinine, and NAGase decreased, whereas milk glucose and galactose increased with WOL. Changes in milk BHB, isocitrate, glucose, glucose-6P, and creatinine were concordant during early lactation and in response to nutrient restriction. Milk galactose and NAGase were modulated by WOL only, whereas glutamate and uric acid concentrations responded to nutrient restriction only. The high-straw increased milk concentrations of FA potentially mobilized from adipose tissue (e.g., C18:0 and cis-9 C18:1 and sum of odd- and branched-chain FA (OBCFA) with carbon chain greater than 16; ∑ OBCFA >C16), and decreased concentrations of FA synthesized de novo by the mammary gland (e.g., sum of FA with 6 to 15 carbons; ∑ C6:0 to C15:0). Similar observations were made during early lactation. Plasma nonesterified FA concentrations had the best single linear regression with EB (R2 = 0.62). Milk isocitrate, Σ C6:0 to C15:0. and cis-9 C18:1 had the best single linear regressions with EB (R2 ≥ 0.44). Milk BHB, isocitrate, galactose, glutamate, and creatinine explained up to 64% of the EB variation observed in the current study using multiple linear regression. Milk concentrations of ∑ C6:0 to C15:0, C18:0, cis-9 C18:1, and ∑ OBCFA >C16 presented some of the best correlations and regressions with other indicators of metabolic status, lipomobilization, and EB, and their responses were concordant during early lactation and during experimental nutrient restriction. Metabolites and FA secreted in milk may serve as noninvasive indicators of metabolic status and EB of early-lactation cows.

Modelling of gene loss propensity in the pangenomes of three Brassica species suggests different mechanisms between polyploids and diploids.

Plant genomes demonstrate significant presence/absence variation (PAV) within a species; however, the factors that lead to this variation have not been studied systematically in Brassica across diploids and polyploids. Here, we developed pangenomes of polyploid Brassica napus and its two diploid progenitor genomes B. rapa and B. oleracea to infer how PAV may differ between diploids and polyploids. Modelling of gene loss suggests that loss propensity is primarily associated with transposable elements in the diploids while in B. napus, gene loss propensity is associated with homoeologous recombination. We use these results to gain insights into the different causes of gene loss, both in diploids and following polyploidization, and pave the way for the application of machine learning methods to understanding the underlying biological and physical causes of gene presence/absence.