Author Correction: Increasing CO2 threatens human nutrition.

An Amendment to this paper has been published and can be accessed via a link at the top of the paper.

Impact of multiple selective breeding programs on genetic diversity in soybean germplasm.

KEY MESSAGE: Independent soybean breeding programs shape genetic diversity from unimproved germplasm to modern cultivars in similar ways, but distinct breeding populations retain unique genetic variation, preserving additional diversity. From the domestication of wild soybean (Glycine soja Sieb. & Zucc.), over 3,000 years ago, to the modern soybean (Glycine max L. Merr) cultivars that provide much of the world's oil and protein, soybean populations have undergone fundamental changes. We evaluated the molecular impact of breeding and selection using 391 soybean accessions including US cultivars and their progenitors from the USDA Soybean Germplasm Collection (CGP), plus two new populations specifically developed to increase genetic diversity and high yield in two alternative gene pools: one derived from exotic G. max germplasm (AGP) and one derived from G. soja (SGP). Reduction in nucleotide genetic diversity (π) was observed with selection within gene pools, but artificial selection in the AGP maintained more diversity than in the CGP. The highest FST levels were seen between ancestral and elite lines in all gene pools, but specific nucleotide-level patterns varied between gene pools. Population structure analyses support that independent selection resulted in high-yielding elite lines with similar allelic compositions in the AGP and CGP. SGP, however, produced elite progeny that were well differentiated from, but lower yielding than, CGP elites. Both the AGP and SGP retained a significant number of private alleles that are absent in CGP. We conclude that the genomic diversity shaped by multiple selective breeding programs can result in gene pools of highly productive elite lines with similar allelic compositions in a genome-wide perspective. Breeding programs with different ancestral lines, however, can retain private alleles representing unique genetic diversity.

Soybean BARCSoySNP6K: An assay for soybean genetics and breeding research.

The limited number of recombinant events in recombinant inbred lines suggests that for a biparental population with a limited number of recombinant inbred lines, it is unnecessary to genotype the lines with many markers. For genomic prediction and selection, previous studies have demonstrated that only 1000-2000 genome-wide common markers across all lines/accessions are needed to reach maximum efficiency of genomic prediction in populations. Evaluation of too many markers will not only increase the cost but also generate redundant information. We developed a soybean (Glycine max) assay, BARCSoySNP6K, containing 6000 markers, which were carefully chosen from the SoySNP50K assay based on their position in the soybean genome and haplotype block, polymorphism among accessions and genotyping quality. The assay includes 5000 single nucleotide polymorphisms (SNPs) from euchromatic and 1000 from heterochromatic regions. The percentage of SNPs with minor allele frequency >0.10 was 95% and 91% in the euchromatic and heterochromatic regions, respectively. Analysis of progeny from two large families genotyped with SoySNP50K versus BARCSoySNP6K showed that the position of the common markers and number of unique bins along linkage maps were consistent based on the SNPs genotyped with the two assays; however, the rate of redundant markers was dramatically reduced with the BARCSoySNP6K. The BARCSoySNP6K assay is proven as an excellent tool for detecting quantitative trait loci, genomic selection and assessing genetic relationships. The assay is commercialized by Illumina Inc. and being used by soybean breeders and geneticists and the list of SNPs in the assay is an ideal resource for SNP genotyping by targeted amplicon sequencing.

Identification of SNP markers associated with soybean fatty acids contents by genome-wide association analyses.

Composition of fatty acids (FAs) in soybean seed is important for the quality and uses of soybean oil. Using gas chromatography, we have measured soybean FAs profiles of 621 soybean accessions (maturity groups I through IV) grown in five different environments; Columbus, OH (2015), Wooster, OH (2014 and 2015), Plymouth, NC (2015), and Urbana, IL (2015). Using publicly available SoySNP50K genotypic data and the FA profiles from this study, a genome-wide association analysis was completed with a compressed mixed linear model to identify 43 genomic regions significantly associated with a fatty acid at a genome wide significance threshold of 5%. Among these regions, one and three novel genomic regions associated with palmitic acid and stearic acid, respectively, were identified across all five environments. Additionally, nine novel environment-specific FA-related genomic regions were discovered providing new insights into the genetics of soybean FAs. Previously reported FA-related loci, such as FATB1a, SACPD-C, and KASIII, were also confirmed in this study. Our results will be useful for future functional studies and marker-assisted breeding for soybean FAs. Supplementary Information: The online version contains supplementary material available at 10.1007/s11032-021-01216-1.

Genomic prediction using training population design in interspecific soybean populations.

Agronomically important traits generally have complex genetic architecture, where many genes have a small and largely additive effect. Genomic prediction has been demonstrated to increase genetic gain and efficiency in plant breeding programs beyond marker-assisted selection and phenotypic selection. The objective of this study was to evaluate the impact of allelic origin, marker density, training population size, and cross-validation schemes on the accuracy of genomic prediction models in an interspecific soybean nested association mapping (NAM) panel. Three cross-validation schemes were used: (a) Within-Family (WF): training population and predictions are made exclusively within each family; (b) Across All families (AF): all the individuals from the three families were randomly assigned to either the training or validation set; (c) Leave one Family out (LFO): each family is predicted using a training set that contains the other two families. Predictive abilities increased with training population size up to 350 individuals, but no significant gains were noted beyond 250 individuals in the training population. The number of markers had a limited impact on the observed predictive ability across traits; increasing markers used in the model above 1000 revealed no significant increases in prediction accuracy. Predictive abilities for AF were not significantly different from the WF method, and predictive abilities across populations for the WF method had a range of 0.58 to 0.70 for maturity, protein, meal, and oil. Our results also showed encouraging prediction accuracies for grain yield (0.58-0.69) using the WF method. Partitioning genomic prediction between G. max and G. soja alleles revealed useful information to select material with a larger allele contribution from both parents and could accelerate allele introgression from exotic germplasm into the elite soybean gene pool. Supplementary Information: The online version contains supplementary material available at 10.1007/s11032-021-01203-6.

Nested association mapping of important agronomic traits in three interspecific soybean populations.

KEY MESSAGE: Glycine soja germplasm can be used to successfully introduce new alleles with the potential to add valuable new genetic diversity to the current elite soybean gene pool. Given the demonstrated narrow genetic base of the US soybean production, it is essential to identify beneficial alleles from exotic germplasm, such as wild soybean, to enhance genetic gain for favorable traits. Nested association mapping (NAM) is an approach to population development that permits the comparison of allelic effects of the same QTL in multiple parents. Seed yield, plant maturity, plant height and plant lodging were evaluated in a NAM panel consisting of 392 recombinant inbred lines derived from three biparental interspecific soybean populations in eight environments during 2016 and 2017. Nested association mapping, combined with linkage mapping, identified three major QTL for plant maturity in chromosomes 6, 11 and 12 associated with alleles from wild soybean resulting in significant increases in days to maturity. A significant QTL for plant height was identified on chromosome 13 with the allele increasing plant height derived from wild soybean. A significant grain yield QTL was detected on chromosome 17, and the allele from Glycine soja had a positive effect of 166 kg ha-1; RIL's with the wild soybean allele yielded on average 6% more than the lines carrying the Glycine max allele. These findings demonstrate the usefulness and potential of alleles from wild soybean germplasm to enhance important agronomic traits in a soybean breeding program.

Increasing CO2 threatens human nutrition.

Dietary deficiencies of zinc and iron are a substantial global public health problem. An estimated two billion people suffer these deficiencies, causing a loss of 63 million life-years annually. Most of these people depend on C3 grains and legumes as their primary dietary source of zinc and iron. Here we report that C3 grains and legumes have lower concentrations of zinc and iron when grown under field conditions at the elevated atmospheric CO2 concentration predicted for the middle of this century. C3 crops other than legumes also have lower concentrations of protein, whereas C4 crops seem to be less affected. Differences between cultivars of a single crop suggest that breeding for decreased sensitivity to atmospheric CO2 concentration could partly address these new challenges to global health.

Genome-wide association study of seed protein, oil and amino acid contents in soybean from maturity groups I to IV.

KEY MESSAGE: Genomic regions associated with seed protein, oil and amino acid contents were identified by genome-wide association analyses. Geographic distributions of haplotypes indicate scope of improvement of these traits. Soybean [Glycine max (L.) Merr.] protein and oil are used worldwide in feed, food and industrial materials. Increasing seed protein and oil contents is important; however, protein content is generally negatively correlated with oil content. We conducted a genome-wide association study using phenotypic data collected from five environments for 621 accessions in maturity groups I-IV and 34,014 markers to identify quantitative trait loci (QTL) for seed content of protein, oil and several essential amino acids. Three and five genomic regions were associated with seed protein and oil contents, respectively. One, three, one and four genomic regions were associated with cysteine, methionine, lysine and threonine content (g kg-1 crude protein), respectively. As previously shown, QTL on chromosomes 15 and 20 were associated with seed protein and oil contents, with both exhibiting opposite effects on the two traits, and the chromosome 20 QTL having the most significant effect. A multi-trait mixed model identified trait-specific QTL. A QTL on chromosome 5 increased oil with no effect on protein content, and a QTL on chromosome 10 increased protein content with little effect on oil content. The chromosome 10 QTL co-localized with maturity gene E2/GmGIa. Identification of trait-specific QTL indicates feasibility to reduce the negative correlation between protein and oil contents. Haplotype blocks were defined at the QTL identified on chromosomes 5, 10, 15 and 20. Frequencies of positive effect haplotypes varied across maturity groups and geographic regions, providing guidance on which alleles have potential to contribute to soybean improvement for specific regions.

Genetic dissection of domestication-related traits in soybean through genotyping-by-sequencing of two interspecific mapping populations.

KEY MESSAGE: A total of 132 domestication-related QTLs, of which 41 were novel, were identified through genotyping-by-sequencing of two Glycine max × Glycine soja populations. Soybean [Glycine max (L.) Merr.] was domesticated in East Asia from the wild progenitor Glycine soja. The domestication process led to many distinct morphological changes that adapt it to cultivation. These include larger seeds, erect growth, larger stem diameter, reduced pod shattering, and altered growth habit. The objective of this study was to identify QTLs controlling key domestication-related traits (DRTs) using interspecific mapping populations. A total of 151 RILs from Williams 82 × PI 468916 and 510 RILs from Williams 82 × PI 479752 were utilized for QTL mapping. These lines were genotyped using a genotyping-by-sequencing protocol which resulted in approximately 5000 polymorphic SNP markers. The number of QTLs detected for each of the eleven DRTs ranged between 0-4 QTLs in the smaller Williams 82 × PI 468916 population and 3-16 QTLs in the larger Williams 82 × PI 479752 population. A total of 132 QTLs were detected, of which 51 are associated with selective sweeps previously related to soybean domestication. These QTLs were detected across all 20 chromosomes within 42 genomic regions. This study identifies 41 novel QTLs not detected in previous studies using smaller populations while also confirming the quantitative nature for several of the important DRTs in soybeans. These results would enable more effective use of the wild germplasm for soybean improvement.

Ensuring due process in the IACUC and animal welfare setting: considerations in developing noncompliance policies and procedures for institutional animal care and use committees and institutional officials.

Every institution that is involved in research with animals is expected to have in place policies and procedures for the management of allegations of noncompliance with the Animal Welfare Act and the U.S. Public Health Service Policy on the Humane Care and Use of Laboratory Animals. We present here a model set of recommendations for institutional animal care and use committees and institutional officials to ensure appropriate consideration of allegations of noncompliance with federal Animal Welfare Act regulations that carry a significant risk or specific threat to animal welfare. This guidance has 3 overarching aims: 1) protecting the welfare of research animals; 2) according fair treatment and due process to an individual accused of noncompliance; and 3) ensuring compliance with federal regulations. Through this guidance, the present work seeks to advance the cause of scientific integrity, animal welfare, and the public trust while recognizing and supporting the critical importance of animal research for the betterment of the health of both humans and animals.-Hansen, B. C., Gografe, S., Pritt, S., Jen, K.-L. C., McWhirter, C. A., Barman, S. M., Comuzzie, A., Greene, M., McNulty, J. A., Michele, D. E., Moaddab, N., Nelson, R. J., Norris, K., Uray, K. D., Banks, R., Westlund, K. N., Yates, B. J., Silverman, J., Hansen, K. D., Redman, B. Ensuring due process in the IACUC and animal welfare setting: considerations in developing noncompliance policies and procedures for institutional animal care and use committees and institutional officials.

Rabies in a Dog Imported from Egypt - Connecticut, 2017.

In 2007, the United States successfully eliminated canine rabies virus variant. Globally, however, dogs remain the principal source of human rabies infections. Since 2007, three cases of canine rabies virus variant were reported in dogs imported into the United States, one each from India (2007), Iraq (2008), and Egypt (2015) (1-3). On December 20, 2017, a dog imported into the United States from Egypt was identified with rabies, representing the second case from Egypt in 3 years. An Egyptian-based animal rescue organization delivered four dogs from Cairo, Egypt, to a flight parent (a person solicited through social media, often not affiliated with the rescue organization, and usually compensated with an airline ticket), who transported the dogs to the United States. The flight parent arrived at John F. Kennedy International Airport (JFK) in New York City and, via transporters (persons who shuttle dogs from one state to another), transferred the dogs to foster families; the dogs ultimately were adopted in three states. The Connecticut Department of Public Health Laboratory (CDPHL) confirmed the presence of a canine rabies virus variant in one of the dogs, a male aged 6 months that was adopted by a Connecticut family. An investigation revealed the possibility of falsified rabies vaccination documentation presented on entry at JFK, allowing the unvaccinated dog entry to the United States. This report highlights the continuing risk posed by the importation of dogs inadequately vaccinated against rabies from high-risk countries and the difficulties in verifying any imported dog's health status and rabies vaccination history.

Dt2 is a gain-of-function MADS-domain factor gene that specifies semideterminacy in soybean.

Similar to Arabidopsis thaliana, the wild soybeans (Glycine soja) and many cultivars exhibit indeterminate stem growth specified by the shoot identity gene Dt1, the functional counterpart of Arabidopsis TERMINAL FLOWER1 (TFL1). Mutations in TFL1 and Dt1 both result in the shoot apical meristem (SAM) switching from vegetative to reproductive state to initiate terminal flowering and thus produce determinate stems. A second soybean gene (Dt2) regulating stem growth was identified, which, in the presence of Dt1, produces semideterminate plants with terminal racemes similar to those observed in determinate plants. Here, we report positional cloning and characterization of Dt2, a dominant MADS domain factor gene classified into the APETALA1/SQUAMOSA (AP1/SQUA) subfamily that includes floral meristem (FM) identity genes AP1, FUL, and CAL in Arabidopsis. Unlike AP1, whose expression is limited to FMs in which the expression of TFL1 is repressed, Dt2 appears to repress the expression of Dt1 in the SAMs to promote early conversion of the SAMs into reproductive inflorescences. Given that Dt2 is not the gene most closely related to AP1 and that semideterminacy is rarely seen in wild soybeans, Dt2 appears to be a recent gain-of-function mutation, which has modified the genetic pathways determining the stem growth habit in soybean.

Disialylated apolipoprotein C-III proteoform is associated with improved lipids in prediabetes and type 2 diabetes.

The apoC-III proteoform containing two sialic acid residues (apoC-III2) has different in vitro effects on lipid metabolism compared with asialylated (apoC-III0) or the most abundant monosialylated (apoC-III1) proteoforms. Cross-sectional and longitudinal associations between plasma apoC-III proteoforms (by mass spectrometric immunoassay) and plasma lipids were tested in two randomized clinical trials: ACT NOW, a study of pioglitazone in subjects with impaired glucose tolerance (n = 531), and RACED (n = 296), a study of intensive glycemic control and atherosclerosis in type 2 diabetes patients. At baseline, higher relative apoC-III2 and apoC-III2/apoC-III1 ratios were associated with lower triglycerides and total cholesterol in both cohorts, and with lower small dense LDL in the RACED. Longitudinally, changes in apoC-III2/apoC-III1 were inversely associated with changes in triglycerides in both cohorts, and with total and small dense LDL in the RACED. apoC-III2/apoC-III1 was also higher in patients treated with PPAR-γ agonists and was associated with reduced cardiovascular events in the RACED control group. Ex vivo studies of apoC-III complexes with higher apoC-III2/apoC-III1 showed attenuated inhibition of VLDL uptake by HepG2 cells and LPL-mediated lipolysis, providing possible functional explanations for the inverse association between a higher apoC-III2/apoC-III1 and hypertriglyceridemia, proatherogenic plasma lipid profiles, and cardiovascular risk.

Toward Analysis of Proteins in Single Cells: A Quantitative Approach Employing Isobaric Tags with MALDI Mass Spectrometry Realized with a Microfluidic Platform.

Protein identification and quantification in individual cells is essential to understand biological processes such as those involved in cell apoptosis, cancer, biomarker discovery, disease diagnostics, pathology, or therapy. Compared with present single cell genome analysis, probing the protein content of single cells has been hampered by the lack of a protein amplification technique. Here, we report the development of a quantitative mass spectrometric approach combined with microfluidic technology reaching the detection sensitivity of high abundant proteins in single cells. A microfluidic platform with a series of chambers and valves, ensuring a set of defined wells for absolute quantification of targeted proteins, was developed and combined with isotopic labeling strategies employing isobaric tags for relative and absolute quantitation (iTRAQ)-labels. To this aim, we adapted iTRAQ labeling to an on-chip protocol. Simultaneous protein digestion and labeling performed on the microfluidic platform rendered the labeling strategy compatible with all necessary manipulation steps on-chip, including the matrix delivery for MALDI-TOF analysis. We demonstrate this approach with the apoptosis related protein Bcl-2 and quantitatively assess the number of Bcl-2 molecules detected. We anticipate that this approach will eventually allow quantification of protein expression on the single cell level.

The association of plasma cystatin C proteoforms with diabetic chronic kidney disease.

BACKGROUND: Cystatin C (CysC) is an endogenous cysteine protease inhibitor that can be used to assess the progression of kidney function. Recent studies demonstrate that CysC is a more specific indicator of glomerular filtration rate (GFR) than creatinine. CysC in plasma exists in multiple proteoforms. The goal of this study was to clarify the association of native CysC, CysC missing N-terminal Serine (CysC des-S), and CysC without three N-terminal residues (CysC des-SSP) with diabetic chronic kidney disease (CKD). RESULTS: Using mass spectrometric immunoassay, the plasma concentrations of native CysC and the two CysC truncation proteoforms were examined in 111 individuals from three groups: 33 non-diabetic controls, 34 participants with type 2 diabetes (DM) and without CKD and 44 participants with diabetic CKD. Native CysC concentrations were 1.4 fold greater in CKD compared to DM group (p = 0.02) and 1.5 fold greater in CKD compared to the control group (p = 0.001). CysC des-S concentrations were 1.55 fold greater in CKD compared to the DM group (p = 0.002) and 1.9 fold greater in CKD compared to the control group (p = 0.0002). CysC des-SSP concentrations were 1.8 fold greater in CKD compared to the DM group (p = 0.008) and 1.52 fold greater in CKD compared to the control group (p = 0.002). In addition, the concentrations of CysC proteoforms were greater in the setting of albuminuria. The truncated CysC proteoform concentrations were associated with estimated GFR independent of native CysC concentrations. CONCLUSION: Our findings demonstrate a greater amount of CysC proteoforms in diabetic CKD. We therefore suggest assessing the role of cystatin C proteoforms in the progression of CKD.

Development of quantitative mass spectrometric immunoassay for serum amyloid A.

OBJECTIVE: Proteins can exist as multiple proteoforms in vivo that can have important roles in physiological and pathological states. METHODS: We present the development and characterization of mass spectrometric immunoassay (MSIA) for quantitative determination of serum amyloid A (SAA) proteoforms. RESULTS: Intra- and inter-day precision revealed CVs <10%. Against existing SAA ELISA, the developed MSIA showed good correlation according to the Altman-Bland plot. Individual concentrations of the SAA proteoforms across a cohort of 170 samples revealed 7 diverse SAA polymorphic types and 12 different proteoforms. CONCLUSION: The new SAA MSIA enables parallel analysis of SAA polymorphisms and quantification of all expressed SAA proteoforms, in a high-throughput and time-efficient manner.

Development of multiplex mass spectrometric immunoassay for detection and quantification of apolipoproteins C-I, C-II, C-III and their proteoforms.

The impetus for discovery and evaluation of protein biomarkers has been accelerated by recent development of advanced technologies for rapid and broad proteome analyses. Mass spectrometry (MS)-based protein assays hold great potential for in vitro biomarker studies. Described here is the development of a multiplex mass spectrometric immunoassay (MSIA) for quantification of apolipoprotein C-I (apoC-I), apolipoprotein C-II (apoC-II), apolipoprotein C-III (apoC-III) and their proteoforms. The multiplex MSIA assay was fast (∼ 40 min) and high-throughput (96 samples at a time). The assay was applied to a small cohort of human plasma samples, revealing the existence of multiple proteoforms for each apolipoprotein C. The quantitative aspect of the assay enabled determination of the concentration for each proteoform individually. Low-abundance proteoforms, such as fucosylated apoC-III, were detected in less than 20% of the samples. The distribution of apoC-III proteoforms varied among samples with similar total apoC-III concentrations. The multiplex analysis of the three apolipoproteins C and their proteoforms using quantitative MSIA represents a significant step forward toward better understanding of their physiological roles in health and disease.

Anthropometric measurement standardization in the US-affiliated pacific: Report from the Children's Healthy Living Program.

OBJECTIVES: Anthropometric standardization is essential to obtain reliable and comparable data from different geographical regions. The purpose of this study is to describe anthropometric standardization procedures and findings from the Children's Healthy Living (CHL) Program, a study on childhood obesity in 11 jurisdictions in the US-Affiliated Pacific Region, including Alaska and Hawai'i. METHODS: Zerfas criteria were used to compare the measurement components (height, waist, and weight) between each trainee and a single expert anthropometrist. In addition, intra- and inter-rater technical error of measurement (TEM), coefficient of reliability, and average bias relative to the expert were computed. RESULTS: From September 2012 to December 2014, 79 trainees participated in at least 1 of 29 standardization sessions. A total of 49 trainees passed either standard or alternate Zerfas criteria and were qualified to assess all three measurements in the field. Standard Zerfas criteria were difficult to achieve: only 2 of 79 trainees passed at their first training session. Intra-rater TEM estimates for the 49 trainees compared well with the expert anthropometrist. Average biases were within acceptable limits of deviation from the expert. Coefficient of reliability was above 99% for all three anthropometric components. CONCLUSIONS: Standardization based on comparison with a single expert ensured the comparability of measurements from the 49 trainees who passed the criteria. The anthropometric standardization process and protocols followed by CHL resulted in 49 standardized field anthropometrists and have helped build capacity in the health workforce in the Pacific Region. Am. J. Hum. Biol. 28:364-371, 2016. © 2015 Wiley Periodicals, Inc.

Genome-wide association study for flowering time, maturity dates and plant height in early maturing soybean (Glycine max) germplasm.

BACKGROUND: Soybean (Glycine max) is a photoperiod-sensitive and self-pollinated species. Days to flowering (DTF) and maturity (DTM), duration of flowering-to-maturity (DFTM) and plant height (PH) are crucial for soybean adaptability and yield. To dissect the genetic architecture of these agronomically important traits, a population consisting of 309 early maturity soybean germplasm accessions was genotyped with the Illumina Infinium SoySNP50K BeadChip and phenotyped in multiple environments. A genome-wide association study (GWAS) was conducted using a mixed linear model that involves both relative kinship and population structure. RESULTS: The linkage disequilibrium (LD) decayed slowly in soybean, and a substantial difference in LD pattern was observed between euchromatic and heterochromatic regions. A total of 27, 6, 18 and 27 loci for DTF, DTM, DFTM and PH were detected via GWAS, respectively. The Dt1 gene was identified in the locus strongly associated with both DTM and PH. Ten candidate genes homologous to Arabidopsis flowering genes were identified near the peak single nucleotide polymorphisms (SNPs) associated with DTF. Four of them encode MADS-domain containing proteins. Additionally, a pectin lyase-like gene was also identified in a major-effect locus for PH where LD decayed rapidly. CONCLUSIONS: This study identified multiple new loci and refined chromosomal regions of known loci associated with DTF, DTM, DFTM and/or PH in soybean. It demonstrates that GWAS is powerful in dissecting complex traits and identifying candidate genes although LD decayed slowly in soybean. The loci and trait-associated SNPs identified in this study can be used for soybean genetic improvement, especially the major-effect loci associated with PH could be used to improve soybean yield potential. The candidate genes may serve as promising targets for studies of molecular mechanisms underlying the related traits in soybean.

Mass Spectrometric Immunoassay for the qualitative and quantitative analysis of the cytokine Macrophage Migration Inhibitory Factor (MIF).

BACKGROUND: The cytokine MIF (Macrophage Migration Inhibitory Factor) has diverse physiological roles and is present at elevated concentrations in numerous disease states. However, its molecular heterogeneity has not been previously investigated in biological samples. Mass Spectrometric Immunoassay (MSIA) may help elucidate MIF post-translational modifications existing in vivo and provide additional clarity regarding its relationship to diverse pathologies. RESULTS: In this work, we have developed and validated a fully quantitative MSIA assay for MIF, and used it in the discovery and quantification of different proteoforms of MIF in serum samples, including cysteinylated and glycated MIF. The MSIA assay had a linear range of 1.56-50 ng/mL, and exhibited good precision, linearity, and recovery characteristics. The new assay was applied to a small cohort of human serum samples, and benchmarked against an MIF ELISA assay. CONCLUSIONS: The quantitative MIF MSIA assay provides a sensitive, precise and high throughput method to delineate and quantify MIF proteoforms in biological samples.