<?xml version="1.0" encoding="UTF-8" standalone="no"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" article-type="research-article" dtd-version="2.3" xml:lang="EN">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Plant Sci.</journal-id>
<journal-title>Frontiers in Plant Science</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Plant Sci.</abbrev-journal-title>
<issn pub-type="epub">1664-462X</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.3389/fpls.2023.1218042</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Plant Science</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Genetic basis of maize maternal haploid induction beyond <italic>MATRILINEAL</italic> and <italic>ZmDMP</italic>
</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name>
<surname>Trentin</surname>
<given-names>Henrique Uliana</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<xref ref-type="author-notes" rid="fn003">
<sup>&#x2020;</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/2319997"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Krause</surname>
<given-names>Matheus Dalsente</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<xref ref-type="author-notes" rid="fn003">
<sup>&#x2020;</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/2294556"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Zunjare</surname>
<given-names>Rajkumar Uttamrao</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<xref ref-type="aff" rid="aff3">
<sup>3</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/402905"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Almeida</surname>
<given-names>Vin&#xed;cius Costa</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<xref ref-type="aff" rid="aff4">
<sup>4</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/2412986"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Peterlini</surname>
<given-names>Edicarlos</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<xref ref-type="aff" rid="aff5">
<sup>5</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/2366946"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Rotarenco</surname>
<given-names>Valeriy</given-names>
</name>
<xref ref-type="aff" rid="aff6">
<sup>6</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/2411891"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Frei</surname>
<given-names>Ursula Karoline</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/2269754"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Beavis</surname>
<given-names>William Dale</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/1224751"/>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>L&#xfc;bberstedt</surname>
<given-names>Thomas</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<xref ref-type="author-notes" rid="fn001">
<sup>*</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/416547"/>
</contrib>
</contrib-group>
<aff id="aff1">
<sup>1</sup>
<institution>Bayer Crop Science</institution>, <addr-line>Coxilha, RS</addr-line>, <country>Brazil</country>
</aff>
<aff id="aff2">
<sup>2</sup>
<institution>Department of Agronomy, Iowa State University</institution>, <addr-line>Ames, IA</addr-line>, <country>United States</country>
</aff>
<aff id="aff3">
<sup>3</sup>
<institution>Division of Genetics, ICAR-Indian Agricultural Research Institute</institution>, <addr-line>New Delhi</addr-line>, <country>India</country>
</aff>
<aff id="aff4">
<sup>4</sup>
<institution>Federal University of Vi&#xe7;osa</institution>, <addr-line>Vi&#xe7;osa, MG</addr-line>, <country>Brazil</country>
</aff>
<aff id="aff5">
<sup>5</sup>
<institution>Department of Agronomy, State University of Maring&#xe1;</institution>, <addr-line>Maring&#xe1;, PR</addr-line>, <country>Brazil</country>
</aff>
<aff id="aff6">
<sup>6</sup>
<institution>VAR BREEDING Ltd.</institution>, <addr-line>Chisinau</addr-line>, <country>Moldova</country>
</aff>
<author-notes>
<fn fn-type="edited-by">
<p>Edited by: Hakan Ozkan, &#xc7;ukurova University, T&#xfc;rkiye</p>
</fn>
<fn fn-type="edited-by">
<p>Reviewed by: Zhang Hongwei, Chinese Academy of Agricultural Sciences, China; Chen Chen, Beijing Academy of Agriculture and Forestry Sciences, China; Manje S. Gowda, The International Maize and Wheat Improvement Center (CIMMYT), Kenya; Haochuan Li, Henan Agricultural University, China</p>
</fn>
<fn fn-type="corresp" id="fn001">
<p>*Correspondence: Thomas L&#xfc;bberstedt, <email xlink:href="mailto:thomasl@iastate.edu">thomasl@iastate.edu</email>
</p>
</fn>
<fn fn-type="equal" id="fn003">
<p>&#x2020;These authors have contributed equally to this work</p>
</fn>
</author-notes>
<pub-date pub-type="epub">
<day>04</day>
<month>10</month>
<year>2023</year>
</pub-date>
<pub-date pub-type="collection">
<year>2023</year>
</pub-date>
<volume>14</volume>
<elocation-id>1218042</elocation-id>
<history>
<date date-type="received">
<day>06</day>
<month>05</month>
<year>2023</year>
</date>
<date date-type="accepted">
<day>05</day>
<month>09</month>
<year>2023</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2023 Trentin, Krause, Zunjare, Almeida, Peterlini, Rotarenco, Frei, Beavis and L&#xfc;bberstedt</copyright-statement>
<copyright-year>2023</copyright-year>
<copyright-holder>Trentin, Krause, Zunjare, Almeida, Peterlini, Rotarenco, Frei, Beavis and L&#xfc;bberstedt</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/">
<p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p>
</license>
</permissions>
<abstract>
<p>In maize, doubled haploid (DH) lines are created <italic>in vivo</italic> through crosses with maternal haploid inducers. Their induction ability, usually expressed as haploid induction rate (HIR), is known to be under polygenic control. Although two major genes (<italic>MTL</italic> and <italic>ZmDMP</italic>) affecting this trait were recently described, many others remain unknown. To identify them, we designed and performed a SNP based (~9007) genome-wide association study using a large and diverse panel of 159 maternal haploid inducers. Our analyses identified a major gene near <italic>MTL</italic>, which is present in all inducers and necessary to disrupt haploid induction. We also found a significant quantitative trait loci (QTL) on chromosome 10 using a case-control mapping approach, in which 793 noninducers were used as controls. This QTL harbors a kokopelli ortholog, whose role in maternal haploid induction was recently described in Arabidopsis. QTL with smaller effects were identified on six of the ten maize chromosomes, confirming the polygenic nature of this trait. These QTL could be incorporated into inducer breeding programs through marker-assisted selection approaches. Further improving HIR is important to reduce the cost of DH line production.</p>
</abstract>
<kwd-group>
<kwd>maize maternal haploid induction rate</kwd>
<kwd>maize haploid inducers</kwd>
<kwd>genome-wide association study</kwd>
<kwd>maize doubled haploid technique</kwd>
<kwd>single fertilization</kwd>
</kwd-group>
<counts>
<fig-count count="6"/>
<table-count count="3"/>
<equation-count count="6"/>
<ref-count count="90"/>
<page-count count="12"/>
<word-count count="6720"/>
</counts>
<custom-meta-wrap>
<custom-meta>
<meta-name>section-in-acceptance</meta-name>
<meta-value>Plant Breeding</meta-value>
</custom-meta>
</custom-meta-wrap>
</article-meta>
</front>
<body>
<sec id="s1" sec-type="intro">
<title>Introduction</title>
<p>Today&#x2019;s maize maternal haploid inducers have substantially higher induction rates than paternal haploid inducers, and are, therefore, the preferred tool for DH line production (<xref ref-type="bibr" rid="B29">Gilles et&#xa0;al., 2017a</xref>). The induction ability of an inducer is estimated by the rate of seeds with haploid embryo over the total number of seeds produced in cross-pollinations with a donor genotype. This rate is referred to as haploid induction rate (HIR) and is usually expressed as a percentage. In paternal inducers, where HIRs of up to 6% were observed (<xref ref-type="bibr" rid="B44">Kindiger and Hamann, 1993</xref>), haploid induction is attributed to a single gene, called <italic>indeterminate gametophyte1</italic> (<italic>ig1</italic>) (<xref ref-type="bibr" rid="B41">Kermicle, 1969</xref>). In maternal inducers, where breeding efforts raised HIRs from 3.2% to 14.5% (<xref ref-type="bibr" rid="B15">Coe, 1959</xref>; <xref ref-type="bibr" rid="B72">Rotarenco et&#xa0;al., 2010</xref>), HIR is under polygenic control (<xref ref-type="bibr" rid="B45">Lashermes and Beckert, 1988</xref>; <xref ref-type="bibr" rid="B19">Deimling et&#xa0;al., 1997</xref>; <xref ref-type="bibr" rid="B68">R&#xf6;ber, 1999</xref>; <xref ref-type="bibr" rid="B63">Prigge et&#xa0;al., 2012a</xref>; <xref ref-type="bibr" rid="B64">Prigge et&#xa0;al., 2012b</xref>; <xref ref-type="bibr" rid="B52">Liu et&#xa0;al., 2015</xref>). The first study designed to investigate the genetic nature of maternal haploid induction was performed by <xref ref-type="bibr" rid="B45">Lashermes and Beckert (1988)</xref>, who observed variation in the induction ability of different segregating populations derived from crosses between inducer Stock 6 (<xref ref-type="bibr" rid="B15">Coe, 1959</xref>) and various noninducer lines.</p>
<p>Two quantitative trait loci (QTL), located on chromosomes 1 and 2, were identified in a F<sub>3</sub> population derived from Stock 6, jointly explaining 17.9% of the phenotypic variance (<xref ref-type="bibr" rid="B19">Deimling et&#xa0;al., 1997</xref>). Studies conducted by <xref ref-type="bibr" rid="B68">R&#xf6;ber (1999)</xref> and <xref ref-type="bibr" rid="B5">Barret et&#xa0;al. (2008)</xref> confirmed the presence of a major QTL on chromosome 1 and reported strong segregation distortion against the inducer allele at this locus. A QTL mapping study conducted by <xref ref-type="bibr" rid="B64">Prigge et&#xa0;al. (2012b)</xref>, which used multiple filial generations derived from four different bi-parental crosses, located this major QTL (<italic>qhir1</italic>) to bin 1.04 and found additional QTL on chromosomes 3 (<italic>qhir2</italic> and <italic>qhir3</italic>), 4 (<italic>qhir4</italic>), 5 (<italic>qhir5</italic> and <italic>qhir6</italic>), 7 (<italic>qhir7</italic>) and 9 (<italic>qhir8</italic>). <italic>qhir8</italic> explained more than 20% of the genetic variance in all filial generations of the cross between inducers CAUHOI and UH400 (<xref ref-type="bibr" rid="B64">Prigge et&#xa0;al., 2012b</xref>). <xref ref-type="bibr" rid="B24">Dong et&#xa0;al. (2013)</xref> subsequently fine mapped <italic>qhir1</italic> to a 243 kb region and observed large variation in HIR among F<sub>2:3</sub> progeny homozygous for this locus. Their findings are in agreement with observations made by <xref ref-type="bibr" rid="B64">Prigge et&#xa0;al. (2012b)</xref>, who noted that while <italic>qhir1</italic> is required for the formation of seeds with haploid embryos, multiple segregating alleles affect HIR. A conditional haplotype extension (CHE) test designed to detect selective sweeps under the assumption of confounding trait expression and population structure, identified two distinct loci (<italic>qhir11</italic> and <italic>qhir12</italic>) within <italic>qhir1</italic> (<xref ref-type="bibr" rid="B31">Hu et&#xa0;al., 2016</xref>). Subsequently, <xref ref-type="bibr" rid="B59">Nair et&#xa0;al. (2017)</xref> demonstrated the importance of <italic>qhir11</italic> only in haploid induction by genetically delineating the <italic>qhir1</italic> locus.</p>
<p>Subsequently independent researchers discovered that a frameshift mutation in a gene underlying <italic>qhir11</italic>, and coding for a patatin-like phospholipase, triggers the formation of haploid seeds (<xref ref-type="bibr" rid="B28">Gilles et&#xa0;al., 2017b</xref>; <xref ref-type="bibr" rid="B40">Kelliher et&#xa0;al., 2017</xref>; <xref ref-type="bibr" rid="B51">Liu et&#xa0;al., 2017</xref>). This frameshift mutation is caused by a 4-bp insertion in the last exon of gene Zm00001d029412, which leads to an early stop codon and results in a truncated protein, losing its proper localization at the pollen endo-plasma membrane (<xref ref-type="bibr" rid="B27">Gilles et&#xa0;al., 2021</xref>). The almost simultaneous discovery of the function of this gene by different researchers resulted in different names: <italic>MATRILINEAL (MTL)</italic> (<xref ref-type="bibr" rid="B40">Kelliher et&#xa0;al., 2017</xref>), <italic>Zea mays Phospholipase A1 (ZmPLA1)</italic> (<xref ref-type="bibr" rid="B51">Liu et&#xa0;al., 2017</xref>) and <italic>NOT LIKE DAD (NLD)</italic> (<xref ref-type="bibr" rid="B28">Gilles et&#xa0;al., 2017b</xref>). By inducing mutations close to the 4-bp insertion site in a noninducer, <xref ref-type="bibr" rid="B40">Kelliher et&#xa0;al. (2017)</xref> observed HIRs ranging from 4.0% to 12.5% in T1 and derived plants. This suggest that mutations in <italic>MTL</italic> are capable to generate inducers with high HIRs. However, the same 4-bp insertion was found in all inducers sequenced by <xref ref-type="bibr" rid="B40">Kelliher et&#xa0;al. (2017)</xref>; <xref ref-type="bibr" rid="B51">Liu et&#xa0;al. (2017)</xref> and <xref ref-type="bibr" rid="B28">Gilles et&#xa0;al. (2017b)</xref>, which have HIRs ranging from 2.0% to 10.0%. This is consistent with the presumed polygenic control of HIR (<xref ref-type="bibr" rid="B45">Lashermes and Beckert, 1988</xref>; <xref ref-type="bibr" rid="B19">Deimling et&#xa0;al., 1997</xref>; <xref ref-type="bibr" rid="B68">R&#xf6;ber, 1999</xref>; <xref ref-type="bibr" rid="B64">Prigge et&#xa0;al., 2012b</xref>; <xref ref-type="bibr" rid="B52">Liu et&#xa0;al., 2015</xref>). <xref ref-type="bibr" rid="B90">Zhong et&#xa0;al. (2019)</xref> identified Zm00001d044822, which encodes for a DUF679 domain membrane protein, to be the underlying cause of <italic>qhir8</italic>. The authors named this gene <italic>ZmDMP</italic> and confirmed its ability to enhance the HIR of inducers possessing <italic>mtl/zmpla1/nld</italic>, as previously reported by <xref ref-type="bibr" rid="B52">Liu et&#xa0;al. (2015)</xref>. They also observed that in genotypes homozygous for the wild type <italic>MTL/ZmPLA1/NLD</italic> allele, knockout mutations in <italic>ZmDMP</italic> can disrupt haploid induction at very low frequencies (0.1-0.3%). This was an interesting observation because it indicated that mutations in different genes can promote and/or have synergetic effect on maternal haploid induction (<xref ref-type="bibr" rid="B37">Jacquier et&#xa0;al., 2020</xref> and <xref ref-type="bibr" rid="B36">Jacquier et&#xa0;al., 2021</xref>).</p>
<p>The savings from employing inducers with higher HIRs can be substantial and are important to breeding programs where haploid selection is visually performed through the <italic>R1-nj</italic> color marker (<xref ref-type="bibr" rid="B13">Chase and Nanda, 1965</xref>). This dominant allele leads to anthocyanin production in the scutellum and aleurone layer of seeds in which double fertilization and embryonic development occurred normally. However, since embryos of haploid seeds do not contain inducer chromosomes, they also lack anthocyanin pigmentation. This difference in embryo pigmentation allows differentiation of haploid and diploid seeds. The savings from employing more efficient inducers can be significant, since the number of donor seeds that will be planted, and whose plants would be detasseled, harvested, dried, and screened for haploid seeds directly depends on the inducers&#x2019; HIR and the number of DH lines that one wants to extract from that donor population. To gain better understanding of the genetic basis of HIR, we designed a genome-wide association study (GWAS) using a diverse set of haploid inducers adapted to Midwestern U.S. conditions. The objectives of this study were to (i) identify promising maternal inducers to be used as parents of new inducer breeding populations, (ii) locate and assess the effect of QTL affecting HIR, and (iii) gain better understanding of the biological processes underlying haploid induction by identifying candidate genes within the detected QTL.</p>
</sec>
<sec id="s2" sec-type="materials|methods">
<title>Materials and methods</title>
<sec id="s2_1">
<title>Germplasm</title>
<p>We used a diverse panel of North American and European haploid inducers to identify QTL affecting HIR. The North American panel contains inducers in the background of public and Ex-PVP lines, such as A637, B73, DK78004, LH82, Mo17, PHG50 and Va35 (<xref ref-type="supplementary-material" rid="SF3">
<bold>Table S1</bold>
</xref>). All lines from the North American panel have the cross between RWS (<xref ref-type="bibr" rid="B69">R&#xf6;ber et&#xa0;al., 2005</xref>) and RWK (<xref ref-type="bibr" rid="B69">R&#xf6;ber et&#xa0;al., 2005</xref>) as their source of <italic>mtl/zmpla1/nld</italic>. The European panel contains lines such as MHI (<xref ref-type="bibr" rid="B10">Chalyk, 1999</xref>), PHI-3 (<xref ref-type="bibr" rid="B72">Rotarenco et&#xa0;al., 2010</xref>), RWS, RWK and eight F<sub>7</sub> progeny derived from the cross of RWS and PHI-3, which were selected for high HIR. In total, 159 inducers were evaluated for HIR. The HIR of all inducers was evaluated in crosses with the commercial hybrid Viking 60-01N, from Albert Lea Seed Company (Minnesota, US). This hybrid was chosen as a donor because it possesses good inducibility and allows clear expression of the <italic>R1-nj</italic> marker. Inducibility is the ability of the donor parent to generate haploid seeds, and multiple studies indicated that the source germplasm impacts HIR (<xref ref-type="bibr" rid="B65">Randolph, 1940</xref>; <xref ref-type="bibr" rid="B12">Chase, 1952</xref>; <xref ref-type="bibr" rid="B45">Lashermes and Beckert, 1988</xref>; <xref ref-type="bibr" rid="B25">Eder and Chalyk, 2002</xref>; <xref ref-type="bibr" rid="B20">De La Fuente et&#xa0;al., 2018</xref>; <xref ref-type="bibr" rid="B58">Nair et&#xa0;al., 2020</xref>; <xref ref-type="bibr" rid="B67">Ren et&#xa0;al., 2022</xref>; <xref ref-type="bibr" rid="B81">Trentin et&#xa0;al., 2022</xref>).</p>
<p>A panel of 793 non-inducer genotypes was considered for a case-control GWAS (cc-GWAS) (<xref ref-type="supplementary-material" rid="SF3">
<bold>Table S1</bold>
</xref>). This panel was composed mainly of DH or highly homozygous lines derived from various genetic backgrounds, such as B73, BS39 (<xref ref-type="bibr" rid="B84">Verzegnazzi et&#xa0;al., 2021</xref>; <xref ref-type="bibr" rid="B73">Santos et&#xa0;al., 2022</xref>), several cycles of the Iowa Stiff Stalk Synthetic Population (BSSS) (Ledesma et&#xa0;al<italic>.</italic>, in preparation), BGEM (Iowa State University/USDA-ARS Germplasm Enhancement Maize project) and Ex-PVP lines. No haploid plants were observed when these genotypes were grown in the field, which is in accordance with <xref ref-type="bibr" rid="B11">Chase (1947)</xref> observation that the occurrence of haploid plants in elite inbred lines happens at very low frequencies (&lt;0.1%).</p>
</sec>
<sec id="s2_2">
<title>Field plot design</title>
<p>The inducers and donor used for this experiment were sown at the Iowa State University, Agronomy and Agricultural Engineering Farm, located in Boone (Iowa), during the summer of 2018. The trial was grown under rainfed conditions, following the recommended practices for maize production in Central Iowa. Pre- and post-emergence herbicides along with manual weeding were used to control invasive species. Urea ammonium nitrate was applied in the area before sowing. Two blocks of inducers and donor were sown at different planting dates to ensure that enough seeds would be generated to obtain reliable estimates of HIR. The first block of donor seeds was sown on May 8<sup>th</sup>, and the first block of inducers seeds was sown on May 21<sup>st</sup>. The second block of donor seeds was sown along with the first block of inducers seeds, while the second block of inducers seeds was sown on May 31<sup>st</sup>. Sowing of inducer blocks was delayed because most inducers have a significantly earlier maturity than the donor.</p>
<p>With the delayed planting of the haploid inducers, it was possible to pollinate the first donor block exclusively with the first inducer block, and the second donor block exclusively with the second inducer block. Each inducer and donor block was composed of subblocks containing 16 plots. Inducers were not randomized among subblocks because the great difference in vigor among them would adversely affect other traits for which data were collected for a companion study (<xref ref-type="bibr" rid="B82">Trentin et&#xa0;al., 2023</xref>). This rationale also justified sowing closely related inducers side-by-side within each subblock. For instance, if hybrid and inbred inducers were randomized, the shading caused by differences in plant height would be detrimental to the growth and development of inbred inducers. Plots were 5.5 meters long, 0.75 meters wide, and were sown with 25 seeds. Inducer and donor blocks were sown side-by-side, and pollen from inducers in a given subblock was carried to the adjacent donor subblock. Multicolored tags with easy-to-match codes were used to ensure that pollen from each inducer plot was placed in the corresponding donor plot. Bulk pollen was collected in tassel bags and used to pollinate at least 10 ears of the donor, which were covered before silk emergence using wax bags. Ears were harvested when seeds reached the black layer stage and were air-dried for one week.</p>
</sec>
<sec id="s2_3">
<title>Data collection</title>
<p>Visual haploid selection was performed using the <italic>R1-nj</italic> marker, and the number of putative haploid and diploid seeds of each ear was recorded. Embryo and endosperm abortion, which are correlated with HIR (<xref ref-type="bibr" rid="B64">Prigge et&#xa0;al., 2012b</xref>; <xref ref-type="bibr" rid="B88">Xu et&#xa0;al., 2013</xref>; <xref ref-type="bibr" rid="B59">Nair et&#xa0;al., 2017</xref>), occur at different stages of seed development (<xref ref-type="bibr" rid="B88">Xu et&#xa0;al., 2013</xref>). This makes the identification of embryo and endosperm aborted seeds quite subjective, and for this reason, these two classes of seeds were not considered as a category in the statistical analysis of the data.</p>
</sec>
<sec id="s2_4">
<title>Ploidy determination of the putative haploid seeds</title>
<p>Putative haploid seeds produced by each inducer at each planting date were bulked in a single envelope, whereas putative diploid seeds were discarded. The ploidy of the putative haploid seeds was verified by cutting seeds in half and observing the presence of anthocyanin pigmentation in the embryonic region. With the true number of haploid seeds, a weighted correction factor <inline-formula>
<mml:math display="inline" id="im1">
<mml:mrow>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mi>C</mml:mi>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> for the <inline-formula>
<mml:math display="inline" id="im2">
<mml:mrow>
<mml:msup>
<mml:mi>i</mml:mi>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mi>h</mml:mi>
</mml:mrow>
</mml:msup>
<mml:mo>&#xa0;</mml:mo>
</mml:mrow>
</mml:math>
</inline-formula> inducer was calculated as follows:</p>
<disp-formula>
<mml:math display="block" id="M1">
<mml:mrow>
<mml:mi>C</mml:mi>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>n</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mfrac>
<mml:mrow>
<mml:mi>T</mml:mi>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mi>P</mml:mi>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>+</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:msub>
<mml:mi>n</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mfrac>
<mml:mrow>
<mml:mi>T</mml:mi>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mi>P</mml:mi>
<mml:msub>
<mml:mi>H</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>n</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>+</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:msub>
<mml:mi>n</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
</disp-formula>
<p>Where the indices <inline-formula>
<mml:math display="inline" id="im3">
<mml:mn>1</mml:mn>
</mml:math>
</inline-formula> and <inline-formula>
<mml:math display="inline" id="im4">
<mml:mn>2</mml:mn>
</mml:math>
</inline-formula> stand for the first and second planting dates, respectively, <inline-formula>
<mml:math display="inline" id="im5">
<mml:mrow>
<mml:mi>T</mml:mi>
<mml:mi>H</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> is the true number of haploid seeds, <inline-formula>
<mml:math display="inline" id="im6">
<mml:mrow>
<mml:mi>P</mml:mi>
<mml:mi>H</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> is the putative number of haploid seeds, and <inline-formula>
<mml:math display="inline" id="im7">
<mml:mi>n</mml:mi>
</mml:math>
</inline-formula> is the number of harvested ears. The <inline-formula>
<mml:math display="inline" id="im8">
<mml:mrow>
<mml:mi>C</mml:mi>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> was then used to correct the number of putative haploid seeds for the <inline-formula>
<mml:math display="inline" id="im9">
<mml:mrow>
<mml:msup>
<mml:mi>i</mml:mi>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mi>h</mml:mi>
</mml:mrow>
</mml:msup>
<mml:mo>&#xa0;</mml:mo>
</mml:mrow>
</mml:math>
</inline-formula> inducer in each planting date, where the difference between the counts before and after the correction was re-classified as diploid seeds.</p>
</sec>
<sec id="s2_5">
<title>Statistical analysis of phenotypic data</title>
<p>Two modeling strategies were implemented to analyze and understand the variability of the phenotypic data by considering planting dates as blocks (<xref ref-type="bibr" rid="B77">Silva, 2017</xref>; <xref ref-type="bibr" rid="B16">Couto et&#xa0;al., 2019</xref>). In both cases, the corrected number of haploid seeds for the <inline-formula>
<mml:math display="inline" id="im10">
<mml:mrow>
<mml:msup>
<mml:mi>i</mml:mi>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mi>h</mml:mi>
</mml:mrow>
</mml:msup>
<mml:mo>&#xa0;</mml:mo>
</mml:mrow>
</mml:math>
</inline-formula> inducer, in the <inline-formula>
<mml:math display="inline" id="im11">
<mml:mrow>
<mml:msup>
<mml:mi>j</mml:mi>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mi>h</mml:mi>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:math>
</inline-formula> planting date, at the <inline-formula>
<mml:math display="inline" id="im12">
<mml:mrow>
<mml:msup>
<mml:mi>k</mml:mi>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mi>h</mml:mi>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:math>
</inline-formula> harvested ear, was modeled as <inline-formula>
<mml:math display="inline" id="im13">
<mml:mrow>
<mml:msub>
<mml:mi>Y</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x223c;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mi>B</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>n</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>m</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>l</mml:mi>
<mml:mo>&#xa0;</mml:mo>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>m</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:msub>
<mml:mi>&#x3c0;</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>, where <inline-formula>
<mml:math display="inline" id="im14">
<mml:mrow>
<mml:msub>
<mml:mi>Y</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>, <inline-formula>
<mml:math display="inline" id="im15">
<mml:mrow>
<mml:msub>
<mml:mi>m</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>, and <inline-formula>
<mml:math display="inline" id="im16">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c0;</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> are the number of haploid seeds, the total number of seeds, and the probability of successful production of haploid seeds, respectively.</p>
<p>The first strategy consisted of fitting a generalized linear model (GLM) with a <italic>logit</italic> link function, expressed:</p>
<disp-formula>
<label>(Model 1)</label>
<mml:math display="block" id="M2">
<mml:mrow>
<mml:mi>l</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>g</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>t</mml:mi>
<mml:mo>&#xa0;</mml:mo>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c0;</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mo>=</mml:mo>
<mml:mi>l</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>g</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c0;</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#x2212;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:msub>
<mml:mi>&#x3c0;</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mo>=</mml:mo>
<mml:mi>&#x3bc;</mml:mi>
<mml:mo>+</mml:mo>
<mml:msub>
<mml:mi>&#x3b2;</mml:mi>
<mml:mi>j</mml:mi>
</mml:msub>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>+</mml:mo>
<mml:msub>
<mml:mi>&#x3c4;</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</disp-formula>
<p>Where <inline-formula>
<mml:math display="inline" id="im17">
<mml:mi>&#x3bc;</mml:mi>
</mml:math>
</inline-formula> is the intercept, <inline-formula>
<mml:math display="inline" id="im18">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3b2;</mml:mi>
<mml:mi>j</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the fixed effect of the <inline-formula>
<mml:math display="inline" id="im19">
<mml:mrow>
<mml:msup>
<mml:mi>j</mml:mi>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mi>h</mml:mi>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:math>
</inline-formula> planting date, and <inline-formula>
<mml:math display="inline" id="im20">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c4;</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the fixed effect of the <inline-formula>
<mml:math display="inline" id="im21">
<mml:mrow>
<mml:msup>
<mml:mi>i</mml:mi>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mi>h</mml:mi>
</mml:mrow>
</mml:msup>
<mml:mo>&#xa0;</mml:mo>
</mml:mrow>
</mml:math>
</inline-formula> inducer. The dispersion parameter <inline-formula>
<mml:math display="inline" id="im22">
<mml:mrow>
<mml:mi>&#x3d5;</mml:mi>
<mml:mo>&#xa0;</mml:mo>
</mml:mrow>
</mml:math>
</inline-formula> from Model 1 was estimated by:</p>
<disp-formula>
<mml:math display="block" id="M3">
<mml:mrow>
<mml:mi>&#x3d5;</mml:mi>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>=</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msubsup>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>k</mml:mi>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>=</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi>n</mml:mi>
</mml:msubsup>
<mml:msubsup>
<mml:mi>r</mml:mi>
<mml:mi>k</mml:mi>
<mml:mn>2</mml:mn>
</mml:msubsup>
</mml:mrow>
<mml:mrow>
<mml:mi>n</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>p</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
</disp-formula>
<p>Where <inline-formula>
<mml:math display="inline" id="im23">
<mml:mrow>
<mml:msub>
<mml:mi>r</mml:mi>
<mml:mi>k</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the Pearson residual of the <inline-formula>
<mml:math display="inline" id="im24">
<mml:mrow>
<mml:msup>
<mml:mi>k</mml:mi>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mi>h</mml:mi>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:math>
</inline-formula> observation, and n-p is the degree of freedom for the deviance residuals. If Pearson&#x2019;s Chi-Square Statistic suggests a lack of fit from Model 1, then a Quasi-Likelihood inferential approach that accounts for overdispersion was assumed.</p>
<p>The second strategy consisted of fitting a generalized linear mixed model (GLMM) with a <italic>logit</italic> link function, expressed by the following linear predictor:</p>
<disp-formula>
<label>(Model 2)</label>
<mml:math display="block" id="M4">
<mml:mrow>
<mml:mi>l</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>g</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>t</mml:mi>
<mml:mo>&#xa0;</mml:mo>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c0;</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mo>=</mml:mo>
<mml:mi>l</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>g</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c0;</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#x2212;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:msub>
<mml:mi>&#x3c0;</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mo>=</mml:mo>
<mml:mi>&#x3bc;</mml:mi>
<mml:mo>+</mml:mo>
<mml:msub>
<mml:mi>&#x3b2;</mml:mi>
<mml:mi>j</mml:mi>
</mml:msub>
<mml:mo>+</mml:mo>
<mml:msub>
<mml:mi>&#x3c4;</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>+</mml:mo>
<mml:msub>
<mml:mi>&#x3f5;</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</disp-formula>
<p>Where <inline-formula>
<mml:math display="inline" id="im25">
<mml:mi>&#x3bc;</mml:mi>
</mml:math>
</inline-formula> is the intercept, <inline-formula>
<mml:math display="inline" id="im26">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3b2;</mml:mi>
<mml:mi>j</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the random effect of the <inline-formula>
<mml:math display="inline" id="im27">
<mml:mrow>
<mml:msup>
<mml:mi>j</mml:mi>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mi>h</mml:mi>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:math>
</inline-formula> planting date with <inline-formula>
<mml:math display="inline" id="im28">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3b2;</mml:mi>
<mml:mi>j</mml:mi>
</mml:msub>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#x223c;</mml:mo>
<mml:mi>N</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>0</mml:mn>
<mml:mo>,</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:msubsup>
<mml:mi>&#x3c3;</mml:mi>
<mml:mi>&#x3b2;</mml:mi>
<mml:mn>2</mml:mn>
</mml:msubsup>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>, <inline-formula>
<mml:math display="inline" id="im29">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c4;</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the fixed effect of the <inline-formula>
<mml:math display="inline" id="im30">
<mml:mrow>
<mml:msup>
<mml:mi>i</mml:mi>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mi>h</mml:mi>
</mml:mrow>
</mml:msup>
<mml:mo>&#xa0;</mml:mo>
</mml:mrow>
</mml:math>
</inline-formula> inducer, and <inline-formula>
<mml:math display="inline" id="im31">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3f5;</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the random effect of the plot with <inline-formula>
<mml:math display="inline" id="im32">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3f5;</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#x223c;</mml:mo>
<mml:mi>N</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>0</mml:mn>
<mml:mo>,</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:msubsup>
<mml:mi>&#x3c3;</mml:mi>
<mml:mi>&#x3f5;</mml:mi>
<mml:mn>2</mml:mn>
</mml:msubsup>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>. The following two models account for the random effects of individual observations (<inline-formula>
<mml:math display="inline" id="im33">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mrow>
<mml:mi>k</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>, i.e., harvested ears) and both individual observations and plot variation, respectively, as follows:</p>
<disp-formula>
<label>(Model 3)</label>
<mml:math display="block" id="M5">
<mml:mrow>
<mml:mi>l</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>g</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>t</mml:mi>
<mml:mo>&#xa0;</mml:mo>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c0;</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mo>=</mml:mo>
<mml:mi>l</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>g</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c0;</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#x2212;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:msub>
<mml:mi>&#x3c0;</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mo>=</mml:mo>
<mml:mi>&#x3bc;</mml:mi>
<mml:mo>+</mml:mo>
<mml:msub>
<mml:mi>&#x3b2;</mml:mi>
<mml:mi>j</mml:mi>
</mml:msub>
<mml:mo>+</mml:mo>
<mml:msub>
<mml:mi>&#x3c4;</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>+</mml:mo>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mrow>
<mml:mi>k</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</disp-formula>
<disp-formula>
<label>(Model 4)</label>
<mml:math display="block" id="M6">
<mml:mrow>
<mml:mi>l</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>g</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>t</mml:mi>
<mml:mo>&#xa0;</mml:mo>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c0;</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mo>=</mml:mo>
<mml:mi>l</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>g</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c0;</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#x2212;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:msub>
<mml:mi>&#x3c0;</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mo>=</mml:mo>
<mml:mi>&#x3bc;</mml:mi>
<mml:mo>+</mml:mo>
<mml:msub>
<mml:mi>&#x3b2;</mml:mi>
<mml:mi>j</mml:mi>
</mml:msub>
<mml:mo>+</mml:mo>
<mml:msub>
<mml:mi>&#x3c4;</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>+</mml:mo>
<mml:msub>
<mml:mi>&#x3f5;</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>+</mml:mo>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mrow>
<mml:mi>k</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</disp-formula>
<p>Where <inline-formula>
<mml:math display="inline" id="im34">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c9;</mml:mi>
<mml:mrow>
<mml:mi>k</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mo>&#xa0;</mml:mo>
</mml:mrow>
</mml:msub>
<mml:mo>&#x223c;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mi>N</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>0</mml:mn>
<mml:mo>,</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:msubsup>
<mml:mi>&#x3c3;</mml:mi>
<mml:mi>&#x3c9;</mml:mi>
<mml:mn>2</mml:mn>
</mml:msubsup>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>. All the other model terms were previously defined.</p>
<p>Three criteria were used to select the best-fit model: (a) the goodness-of-fit via the Akaike Information Criterion (AIC; <xref ref-type="bibr" rid="B1">Akaike, 1974</xref>), (b) the Bayesian Information Criterion (BIC; <xref ref-type="bibr" rid="B75">Schwarz, 1978</xref>); and (c) visual inspection using half-normal plots with a simulated envelope (<xref ref-type="bibr" rid="B22">Dem&#xe9;trio et&#xa0;al., 2014</xref>) using the R package <italic>hnp</italic> (<xref ref-type="bibr" rid="B55">Moral et&#xa0;al., 2017</xref>). Plots of the Pearson residuals were also compared between models. Maximum likelihood estimates for the fixed effects of inducers <inline-formula>
<mml:math display="inline" id="im35">
<mml:mrow>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mover accent="true">
<mml:mi>Y</mml:mi>
<mml:mo>^</mml:mo>
</mml:mover>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>.</mml:mo>
<mml:mo>.</mml:mo>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> from the selected model were further used for GWAS analysis. Asymptotic confidence intervals were obtained using the package <italic>emmeans</italic> (<xref ref-type="bibr" rid="B47">Lenth, 2020</xref>). All statistical analyzes were carried out using the <italic>lme4</italic> package (<xref ref-type="bibr" rid="B6">Bates et&#xa0;al., 2015</xref>) in the R environment (<xref ref-type="bibr" rid="B66">R Core Team, 2022</xref>).</p>
</sec>
<sec id="s2_6">
<title>Genotyping and quality control</title>
<p>Leaf samples were collected when seedlings were at the V3 stage. After lyophilization, they were sent to the International Maize and Wheat Improvement Center (CIMMYT), where they were genotyped with the Diversity Arrays Technology (DArTSeq) platform (<xref ref-type="bibr" rid="B42">Kilian et&#xa0;al., 2012</xref>). SNP calling was performed using the DArTsoft analytical pipepeline (<uri xlink:href="https://www.diversityarrays.com">https://www.diversityarrays.com</uri>) and the version 4 of B73&#x2019;s reference genome (AGPv4) as a reference. A total of 32,929 SNP markers were initially identified, and 9,007 remained after filtering for a call rate of at least 70% and minor allele frequency (MAF) of 1% (<uri xlink:href="https://dr.lib.iastate.edu/home">https://dr.lib.iastate.edu/home</uri>). The software Beagle 5.0 (<xref ref-type="bibr" rid="B8">Browning et&#xa0;al., 2018</xref>) was subsequently used for the imputation of missing data.</p>
</sec>
<sec id="s2_7">
<title>Genome-wide association analyses</title>
<p>We performed marker-trait association using quantitative and binary response variables. The quantitative GWAS was implemented in GAPIT (<xref ref-type="bibr" rid="B49">Lipka et&#xa0;al., 2012</xref>) with the multiple loci models (i) BLINK (<xref ref-type="bibr" rid="B34">Huang et&#xa0;al., 2019</xref>), (ii) FarmCPU (<xref ref-type="bibr" rid="B50">Liu et&#xa0;al., 2016</xref>), and (iii) MLMM (<xref ref-type="bibr" rid="B76">Segura et&#xa0;al., 2012</xref>). The estimated HIR from the best-fit phenotypic model was normalized (i.e., transformed) using the <italic>bestNormalize</italic> package in R (<xref ref-type="bibr" rid="B61">Peterson, 2021</xref>). The binary GWAS was a case-control and included the 159 inducers considered &#x201c;cases&#x201d;, and the panel of 793 noninducers considered &#x201c;controls&#x201d; (<xref ref-type="supplementary-material" rid="SF3">
<bold>Table S1</bold>
</xref>). In total, 952 genotypes were included. The phenotypes, a binary random variable, were coded as ones for the cases and zeros for the controls, and a mixed logistic regression model was implemented in the R package <italic>milorGWAS</italic> (<xref ref-type="bibr" rid="B54">Milet et&#xa0;al., 2020</xref>). The reader is referred to the cited references for details on these GWAS models.</p>
<p>For the quantitative GWAS, population/family structure was accounted by the additive genomic relationship G (<xref ref-type="bibr" rid="B83">VanRaden, 2008</xref>), and the first seven principal components (PCs) obtained from principal component analysis (PCA, built-in R function <italic>prcomp</italic>) of the numerical matrix of SNP markers, calculated via singular value decomposition (<xref ref-type="fig" rid="f1">
<bold>Figures&#xa0;1A, B</bold>
</xref>). For the cc-GWAS, the cases and controls grouping of genotypes were completely confounded by the first two PCs (<xref ref-type="fig" rid="f1">
<bold>Figure&#xa0;1C</bold>
</xref>), and hence population structure was only modeled by the matrix G. It is unlikely such correction can fully adjust for the population structure; hence, a permutation analysis was performed as follow: the phenotypes (zeros and ones) were shuffled under the null hypothesis of no association between markers and phenotypes, and the lowest <italic>p-value</italic> recorded. This process was repeated 10,000 times, and the 95% quantile of the distribution of the <italic>-log10(p-values)</italic> was used as a reference threshold. In addition to the permutation analysis for the cc-GWAS, SNPs were declared statistically significant for both GWAS analyses at an FDR-adjusted <italic>p-value</italic> of less than 0.05.</p>
<fig id="f1" position="float">
<label>Figure&#xa0;1</label>
<caption>
<p>Scatterplot of two first principal components (PCs) from principal component analysis (PCA) of SNP matrix for <bold>(A)</bold> quantitative and <bold>(C)</bold> binary GWAS. The percentage of explained variance for each PC in the quantitative GWAS is shown in <bold>(B)</bold>.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fpls-14-1218042-g001.tif"/>
</fig>
<p>Lastly, we calculated the genome-wide (i.e., across chromosomes) pairwise linkage disequilibrium (LD) using the function <italic>LD.decay</italic> from the <italic>sommer</italic> package (<xref ref-type="bibr" rid="B17">Covarrubias-Pazaran, 2016</xref>). The LD decay was determined by locally-weighted polynomial regression (LOESS, build-in R function) with the physical distance between markers as a function of the coefficient of determination <inline-formula>
<mml:math display="inline" id="im36">
<mml:mrow>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:msup>
<mml:mi>r</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>. The threshold <inline-formula>
<mml:math display="inline" id="im37">
<mml:mrow>
<mml:msup>
<mml:mi>r</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:math>
</inline-formula> value when estimating LD decay was calculated with the <inline-formula>
<mml:math display="inline" id="im38">
<mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mn>95</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mi>h</mml:mi>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:math>
</inline-formula> percentile of the LD distribution between unlinked markers [<italic>LD.decay(&#x2026;, unlinked = TRUE, gamma = 0.95)</italic>]. The search browser in MaizeGDB (<ext-link ext-link-type="uri" xlink:href="https://www.maizegdb.org/gbrowse">https://www.maizegdb.org/gbrowse</ext-link>) was used to identify the putative candidate genes in the interval of significant SNPs. Maize B73 RefGen_v4 version was used to locate the candidate genes.</p>
</sec>
</sec>
<sec id="s3" sec-type="results">
<title>Results</title>
<sec id="s3_1">
<title>Phenotypic analyses</title>
<p>Estimates of HIR with asymptotic confidence intervals were obtained with Model 4 (<xref ref-type="fig" rid="f2">
<bold>Figure&#xa0;2</bold>
</xref>; <xref ref-type="supplementary-material" rid="SF4">
<bold>Table S2</bold>
</xref>). For the selection of the best-fit phenotypic model, the AIC values ranged from 14702.24 (Model 4) to 16932.52 (Model 1), and for the BIC criterion, from 15668.51 (Model 4) to 17886.93 (Model 1). In addition, Model 1 was adjusted with the Quasi-likelihood inference due to overdispersion (<xref ref-type="table" rid="T1">
<bold>Table&#xa0;1</bold>
</xref>). Thus, both selection criteria suggested that Model 4, which included random effects for planting dates, plots, and individual observations (i.e., harvested ears), was the best-fit model for the observed data. The visual inspection of the half-normal and Pearson residuals plots also showed the superior fit of Model 4 (<xref ref-type="supplementary-material" rid="SF1">
<bold>Figures S1</bold>
</xref> and <xref ref-type="supplementary-material" rid="SF2">
<bold>S2</bold>
</xref>). Prior to GWAS analyses, the estimated HIR values were transformed to ensure normality (<xref ref-type="fig" rid="f3">
<bold>Figure&#xa0;3</bold>
</xref>).</p>
<fig id="f2" position="float">
<label>Figure&#xa0;2</label>
<caption>
<p>Estimated HIR with asymptotic confidence intervals estimated from Model 4.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fpls-14-1218042-g002.tif"/>
</fig>
<table-wrap id="T1" position="float">
<label>Table&#xa0;1</label>
<caption>
<p>Goodness-of-fit for Models 1 to 4. Model 4 is the best-fit model.</p>
</caption>
<table frame="hsides">
<thead>
<tr>
<th valign="top" rowspan="2" align="center">Model</th>
<th valign="top" rowspan="2" align="center">Log-likelihood</th>
<th valign="top" rowspan="2" align="center">AIC</th>
<th valign="top" rowspan="2" align="center">BIC</th>
<th valign="top" colspan="3" align="center">Variance components</th>
<th valign="top" rowspan="2" align="center">Dispersion<break/>parameter*</th>
</tr>
<tr>
<th valign="top" align="center">Block</th>
<th valign="top" align="center">Plot</th>
<th valign="top" align="center">Ears</th>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" align="left">1</td>
<td valign="top" align="center">-8305.26</td>
<td valign="top" align="center">16932.52</td>
<td valign="top" align="center">17886.93</td>
<td valign="top" align="center">&#x2013;</td>
<td valign="top" align="center">&#x2013;</td>
<td valign="top" align="center">&#x2013;</td>
<td valign="top" align="center">3.05</td>
</tr>
<tr>
<td valign="top" align="left">2</td>
<td valign="top" align="center">-7889.48</td>
<td valign="top" align="center">16102.95</td>
<td valign="top" align="center">17063.30</td>
<td valign="top" align="center">0.01</td>
<td valign="top" align="center">0.10</td>
<td valign="top" align="center">&#x2013;</td>
<td valign="top" align="center">&#x2013;</td>
</tr>
<tr>
<td valign="top" align="left">3</td>
<td valign="top" align="center">-7257.02</td>
<td valign="top" align="center">14838.03</td>
<td valign="top" align="center">15798.38</td>
<td valign="top" align="center">0.01</td>
<td valign="top" align="center">&#x2013;</td>
<td valign="top" align="center">0.26</td>
<td valign="top" align="center">&#x2013;</td>
</tr>
<tr>
<td valign="top" align="left">4</td>
<td valign="top" align="center">-7188.12</td>
<td valign="top" align="center">14702.24</td>
<td valign="top" align="center">15668.51</td>
<td valign="top" align="center">0.01</td>
<td valign="top" align="center">0.08</td>
<td valign="top" align="center">0.20</td>
<td valign="top" align="center">&#x2013;</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<fn>
<p>*Model fitted with Quasi-likelihood approach, "-" symbol indicates no data.</p>
</fn>
</table-wrap-foot>
</table-wrap>
<fig id="f3" position="float">
<label>Figure&#xa0;3</label>
<caption>
<p>Data transformation of the estimated HIR from Model 4 to perform the quantitative GWAS.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fpls-14-1218042-g003.tif"/>
</fig>
</sec>
<sec id="s3_2">
<title>Genome-wide association analyses</title>
<p>The average rate of LD decay at the <inline-formula>
<mml:math display="inline" id="im39">
<mml:mrow>
<mml:msup>
<mml:mi>r</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:math>
</inline-formula> estimated threshold of 0.096 was of 113 Kb (<xref ref-type="fig" rid="f4">
<bold>Figure&#xa0;4</bold>
</xref>). The reference threshold obtained from the permutation analysis for the cc-GWAS was of 4.89, whereas the <italic>-log10(p-value)</italic> of the statistically significant SNP marker on chromosome 10 was of 5.30, with an adjusted FDR value of 0.04 (<xref ref-type="fig" rid="f5">
<bold>Figure&#xa0;5D</bold>
</xref>). In the quantitative GWAS, significant SNPs were detected on chromosomes 1, 2, 3, 7 and 8 (<xref ref-type="fig" rid="f5">
<bold>Figures&#xa0;5A&#x2013;C</bold>
</xref>). The significant SNPs detected in chromosome 1 were located in positions 66.6 Mb (S1_66636144), 69.3 Mb (S1_69321282) and 76.1 Mb (S1_76160596), flaking the region of <italic>MTL/ZmPLA1/NLD</italic>, which is located in position 69.4 Mb. The significant SNP detected on chromosome 2 (S2_220376487) is located at position 220.3 Mb and does not overlap with the QTL detected by <xref ref-type="bibr" rid="B19">Deimling et&#xa0;al. (1997)</xref> on the same chromosome. The significant SNP detected on chromosome 3 (S3_13701123) is located at position 13.7 Mb and does not overlap with any of chromosome 3 QTL&#x2019;s previously identified by <xref ref-type="bibr" rid="B64">Prigge et&#xa0;al. (2012b)</xref>. The significant SNP detected on chromosome 7 (S7_131226237) is located at position 131.2 Mb, and also does not overlap with the chromosome 7 QTL previously identified by <xref ref-type="bibr" rid="B64">Prigge et&#xa0;al. (2012b)</xref>. The two significant SNP detected on chromosome 8 (S8_138779408 and S8_174792234) are located at positions 138.7 and 174.7 Mb, respectively. There are no reports of HIR QTL on this chromosome in the literature. The cc-GWAS detected a significant SNP (S10_141729953) at position 141.7 Mb on chromosome 10. To the best of our knowledge, there are no reports of HIR QTL on chromosome 10 in the literature. A summary overview of the statistically significant SNPs is given in <xref ref-type="table" rid="T2">
<bold>Table&#xa0;2</bold>
</xref>, and the allelic distribution as a function of the estimated HIR in <xref ref-type="fig" rid="f6">
<bold>Figure&#xa0;6</bold>
</xref>. For the binary GWAS, a classification table differentiating the number of cases and controls according to the genotyped allele is presented (<xref ref-type="table" rid="T3">
<bold>Table&#xa0;3</bold>
</xref>).</p>
<fig id="f4" position="float">
<label>Figure&#xa0;4</label>
<caption>
<p>LD decay with an estimated threshold value of 0.096. The fitted LOESS curve is shown in orange.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fpls-14-1218042-g004.tif"/>
</fig>
<fig id="f5" position="float">
<label>Figure&#xa0;5</label>
<caption>
<p>Manhattan and qq-plot of the quantitative GWAS using Blink <bold>(A)</bold>, FarmCPU <bold>(B)</bold>, MLMM <bold>(C)</bold>, and the case-control GWAS <bold>(D)</bold>.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fpls-14-1218042-g005.tif"/>
</fig>
<table-wrap id="T2" position="float">
<label>Table&#xa0;2</label>
<caption>
<p>Summary of the statistically significant SNPs containing name, chromosome, position in base pairs, alleles, minor allele frequency (MAF), and models&#x2019; name. C-C stands for cc-GWAS.</p>
</caption>
<table frame="hsides">
<thead>
<tr>
<th valign="top" align="center">SNP marker</th>
<th valign="top" align="center">Chromosome</th>
<th valign="top" align="center">Position</th>
<th valign="top" align="center">Alleles</th>
<th valign="top" align="center">MAF<xref ref-type="table-fn" rid="fnT2_1">
<sup>a</sup>
</xref>
</th>
<th valign="top" align="center">BLINK</th>
<th valign="top" align="center">FarmCPU</th>
<th valign="top" align="center">MLMM</th>
<th valign="top" align="center">C-C</th>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" align="left">S1_66636144</td>
<td valign="top" align="center">1</td>
<td valign="top" align="left">66636144</td>
<td valign="top" align="center">A/G</td>
<td valign="top" align="center">0.21</td>
<td valign="top" align="center">x</td>
<td valign="top" align="center"/>
<td valign="top" align="center">x</td>
<td valign="top" align="center"/>
</tr>
<tr>
<td valign="top" align="left">S1_69321282</td>
<td valign="top" align="center">1</td>
<td valign="top" align="left">69321282</td>
<td valign="top" align="center">C/G</td>
<td valign="top" align="center">0.13</td>
<td valign="top" align="center">x</td>
<td valign="top" align="center"/>
<td valign="top" align="center"/>
<td valign="top" align="center"/>
</tr>
<tr>
<td valign="top" align="left">S1_76160596</td>
<td valign="top" align="center">1</td>
<td valign="top" align="left">76160596</td>
<td valign="top" align="center">G/T</td>
<td valign="top" align="center">0.22</td>
<td valign="top" align="center"/>
<td valign="top" align="center">x</td>
<td valign="top" align="center"/>
<td valign="top" align="center"/>
</tr>
<tr>
<td valign="top" align="left">S2_220376487</td>
<td valign="top" align="center">2</td>
<td valign="top" align="left">220376487</td>
<td valign="top" align="center">A/G</td>
<td valign="top" align="center">0.18</td>
<td valign="top" align="center">x</td>
<td valign="top" align="center"/>
<td valign="top" align="center"/>
<td valign="top" align="center"/>
</tr>
<tr>
<td valign="top" align="left">S3_13701123</td>
<td valign="top" align="center">3</td>
<td valign="top" align="left">13701123</td>
<td valign="top" align="center">A/T</td>
<td valign="top" align="center">0.08</td>
<td valign="top" align="center"/>
<td valign="top" align="center">x</td>
<td valign="top" align="center"/>
<td valign="top" align="center"/>
</tr>
<tr>
<td valign="top" align="left">S7_131226237</td>
<td valign="top" align="center">7</td>
<td valign="top" align="left">31226237</td>
<td valign="top" align="center">C/T</td>
<td valign="top" align="center">0.03</td>
<td valign="top" align="center"/>
<td valign="top" align="center">x</td>
<td valign="top" align="center"/>
<td valign="top" align="center"/>
</tr>
<tr>
<td valign="top" align="left">S8_138779408</td>
<td valign="top" align="center">8</td>
<td valign="top" align="left">138779408</td>
<td valign="top" align="center">A/G</td>
<td valign="top" align="center">0.08</td>
<td valign="top" align="center"/>
<td valign="top" align="center">x</td>
<td valign="top" align="center"/>
<td valign="top" align="center"/>
</tr>
<tr>
<td valign="top" align="left">S8_174792234</td>
<td valign="top" align="center">8</td>
<td valign="top" align="left">74792234</td>
<td valign="top" align="center">A/C</td>
<td valign="top" align="center">0.06</td>
<td valign="top" align="center">x</td>
<td valign="top" align="center"/>
<td valign="top" align="center"/>
<td valign="top" align="center"/>
</tr>
<tr>
<td valign="top" align="left">S10_141729953</td>
<td valign="top" align="center">10</td>
<td valign="top" align="left">141729953</td>
<td valign="top" align="center">G/T</td>
<td valign="top" align="center">0.13</td>
<td valign="top" align="center"/>
<td valign="top" align="center"/>
<td valign="top" align="center"/>
<td valign="top" align="center">x</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<fn id="fnT2_1">
<label>a</label>
<p>MAF was calculated respective to the population used (i.e., 159 or 952 genotypes).</p>
</fn>
</table-wrap-foot>
</table-wrap>
<fig id="f6" position="float">
<label>Figure&#xa0;6</label>
<caption>
<p>Allelic distribution as a function of the estimated HIR at each statistically significant SNP.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fpls-14-1218042-g006.tif"/>
</fig>
<table-wrap id="T3" position="float">
<label>Table&#xa0;3</label>
<caption>
<p>Counts of number of genotypes for SNP marker S10_141729953 from cc-GWAS analysis.</p>
</caption>
<table frame="hsides">
<thead>
<tr>
<th valign="top" align="center">Allele</th>
<th valign="top" align="center">Cases</th>
<th valign="top" align="center">Controls</th>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" align="left">G</td>
<td valign="top" align="center">11</td>
<td valign="top" align="center">789</td>
</tr>
<tr>
<td valign="top" align="left">T</td>
<td valign="top" align="center">146</td>
<td valign="top" align="center">3</td>
</tr>
</tbody>
</table>
</table-wrap>
</sec>
<sec id="s3_3">
<title>Identification and annotation of candidate genes</title>
<p>The physical locations of the nine significant SNPs identified were recorded using the B73 RefGen_v4 based on the LD decay distance. The candidate genes were identified by the genome browser search available at Maizegdb (<ext-link ext-link-type="uri" xlink:href="https://www.maizegdb.org/gbrowse">https://www.maizegdb.org/gbrowse</ext-link>). A total of 41 candidate genes with descriptions were found, out of which, 37 were protein-coding, and the remaining 4 were not annotated. The chromosomal distribution of these candidate genes included 12 located on chromosome 1, 8 on chromosome 2, 5 on chromosome 3, 5 on chromosome 7, 3 on chromosome 8, and 8 on chromosome 10. The functional domain of each of the candidate genes is provided in <xref ref-type="supplementary-material" rid="SF5">
<bold>Table S3</bold>
</xref>.</p>
</sec>
</sec>
<sec id="s4" sec-type="discussion">
<title>Discussion</title>
<p>The widespread use of the DH technique in maize was influenced by the increase of HIRs to manageable levels and the addition of anthocyanin-based phenotypic markers in maternal inducers. These two fields have recently seen remarkable advancements. For instance, <xref ref-type="bibr" rid="B14">Chen et&#xa0;al. (2022)</xref> recently developed an efficient and accurate (99.1% accuracy) haploid selection marker through co-expression of two transcription factor genes (<italic>ZmC1</italic> and <italic>ZmR2</italic>) in novel haploid inducer line called &#x2018;Maize Anthocyanin Gene InduCer 1 (MAGIC1)&#x2019;. Shortly after, <xref ref-type="bibr" rid="B87">Wang et&#xa0;al. (2023)</xref> showed that expression of <italic>RUBY</italic> enabled 100% accuracy of haploid identification based on deep betalain pigmentation. The development of a group of inducers with numerous anthocyanin markers and HIRs above 12.0% allowed for effective haploid generation and screening (<xref ref-type="bibr" rid="B72">Rotarenco et&#xa0;al., 2010</xref>). The identification of the <italic>MTL/ZmPLA1/NLD</italic> (<xref ref-type="bibr" rid="B28">Gilles et&#xa0;al., 2017b</xref>; <xref ref-type="bibr" rid="B40">Kelliher et&#xa0;al., 2017</xref>; <xref ref-type="bibr" rid="B51">Liu et&#xa0;al., 2017</xref>) DUF679 (<xref ref-type="bibr" rid="B90">Zhong et&#xa0;al., 2019</xref>), <italic>ZmPLD3</italic> (<xref ref-type="bibr" rid="B48">Li et&#xa0;al., 2021</xref>) and <italic>ZmPOD65</italic> (<xref ref-type="bibr" rid="B38">Jiang et&#xa0;al., 2022</xref>) genes contributed to a better understanding of the molecular mechanisms of haploid induction (<xref ref-type="bibr" rid="B37">Jacquier et&#xa0;al., 2020</xref>). Previous research (<xref ref-type="bibr" rid="B45">Lashermes and Beckert, 1988</xref>; <xref ref-type="bibr" rid="B19">Deimling et&#xa0;al., 1997</xref>; <xref ref-type="bibr" rid="B68">R&#xf6;ber, 1999</xref>; <xref ref-type="bibr" rid="B64">Prigge et&#xa0;al., 2012b</xref>; <xref ref-type="bibr" rid="B52">Liu et&#xa0;al., 2015</xref>) revealed evidence for the polygenic regulation of HIR and identified various QTL controlling this trait. To increase the effectiveness of inducer breeding and, as a result, lower the cost of DH line production, it is critical to identify and validate the influence of QTL affecting HIR. To achieve this goal, we designed and performed a GWAS with a large and diverse panel of inducers. In a previous GWAS for HIR, the induction ability of 53 inducers was evaluated in distinct donors (<xref ref-type="bibr" rid="B31">Hu et&#xa0;al., 2016</xref>). Comparing the induction ability obtained through crosses with different donors is not ideal, since inducibility and HIR effects might be confounded (<xref ref-type="bibr" rid="B65">Randolph, 1940</xref>; <xref ref-type="bibr" rid="B12">Chase, 1952</xref>; <xref ref-type="bibr" rid="B45">Lashermes and Beckert, 1988</xref>; <xref ref-type="bibr" rid="B25">Eder and Chalyk, 2002</xref>; <xref ref-type="bibr" rid="B20">De La Fuente et&#xa0;al., 2018</xref>). In our study, 159 inducers were crossed to the same donor at two different planting dates.</p>
<sec id="s4_1">
<title>Genome-wide association analyses</title>
<p>To this date, just a few studies were performed to identify QTL affecting the HIR of maize maternal haploid inducers (<xref ref-type="bibr" rid="B19">Deimling et&#xa0;al., 1997</xref>; <xref ref-type="bibr" rid="B5">Barret et&#xa0;al., 2008</xref>; <xref ref-type="bibr" rid="B64">Prigge et&#xa0;al., 2012b</xref>; <xref ref-type="bibr" rid="B31">Hu et&#xa0;al., 2016</xref>). However, multiple studies have been performed to fine-map or validate the function of previously identified QTL (<xref ref-type="bibr" rid="B24">Dong et&#xa0;al., 2013</xref>; <xref ref-type="bibr" rid="B52">Liu et&#xa0;al., 2015</xref>; <xref ref-type="bibr" rid="B28">Gilles et&#xa0;al., 2017b</xref>; <xref ref-type="bibr" rid="B40">Kelliher et&#xa0;al., 2017</xref>; <xref ref-type="bibr" rid="B51">Liu et&#xa0;al., 2017</xref>; <xref ref-type="bibr" rid="B59">Nair et&#xa0;al., 2017</xref>; <xref ref-type="bibr" rid="B90">Zhong et&#xa0;al., 2019</xref>; <xref ref-type="bibr" rid="B48">Li et&#xa0;al., 2021</xref>).</p>
<p>In this study, we used four different GWAS methods, including BLINK, FarmCPU, MLMM, and cc-GWAS, to identify potential associations with single nucleotide polymorphisms (SNPs). We identified a total of nine SNPs that were significantly associated with the trait. They correspond to different chromosomal origins, validating the polygenic control of HIR in maize, as previously reported (<xref ref-type="bibr" rid="B45">Lashermes and Beckert, 1988</xref>; <xref ref-type="bibr" rid="B19">Deimling et&#xa0;al., 1997</xref>; <xref ref-type="bibr" rid="B68">R&#xf6;ber, 1999</xref>; <xref ref-type="bibr" rid="B64">Prigge et&#xa0;al., 2012b</xref>; <xref ref-type="bibr" rid="B52">Liu et&#xa0;al., 2015</xref>). For instance, two SNPs (S8_138779408 and S8_174792234) were identified on chromosome-8 and a SNP (S10_14172995) chromosome-10. No QTLs were previously reported on these chromosomes. However, some of the significant SNPs that we identified are closely linked to the previously discovered QTLs for HIR. The marker S1_69321282, at position 69.3 Mb on chromosome 1, is closely linked to the <italic>MTL/ZmPLA1/NLD</italic> gene at position 69.4 Mb (identified in the region <italic>qhir11</italic>). Similarly, marker S1_76160596 (position 76.1 Mb), is located inside the <italic>qhir12</italic> region which <xref ref-type="bibr" rid="B31">Hu et&#xa0;al. (2016)</xref> found a higher CHE score for <italic>qhir12</italic> than for <italic>qhir11</italic>. Interestingly, <xref ref-type="bibr" rid="B31">Hu et&#xa0;al. (2016)</xref> also reported that the <italic>qhir12</italic> haplotype was exclusively found in inducers, whereas the <italic>qhir11</italic> haplotype could be found in 2.7% of the non-inducers. <xref ref-type="bibr" rid="B59">Nair et&#xa0;al. (2017)</xref> showed the significance of <italic>qhir11</italic> region compared <italic>qhir12</italic> for inducing haploid embryo in maize. At the same time, it was discovered that <italic>qhir11</italic> harbours the <italic>MTL/ZmPLA1/NLD</italic> allele that possesses major effect on haploid induction in maize. The differences in germplasm, genotyping method, statistical models and sample size might explain the divergence between earlier studies and our results. The inbred panel used in this study essentially carried the mutant <italic>MTL/ZmPLA1/NLD</italic> allele and therefore it is obvious that the gene would not be detected in our quantitative GWAS.</p>
<p>We also found two significant markers, S2_220376487 at 220.3 Mb and S7_13226237 at 132.2 Mb on chromosome 2 and chromosome 7, respectively, in regions that does not overlap QTL identified on chromosome 2 (<xref ref-type="bibr" rid="B19">Deimling et&#xa0;al., 1997</xref>) and chromosome 7 (<xref ref-type="bibr" rid="B64">Prigge et&#xa0;al., 2012b</xref>). Although on the same chromosome, we report these new genomic regions that could influence haploid induction. PVE (phenotypic variance explained) is an important metric used in GWAS to quantify the contribution of SNP for trait value and indicator of the biological relevance of the SNP. PVE for the significant SNPs varied from 5.9% to 37% explaining substantial variance for HIR detected in our study. SNPs S1_66636144, located at position 66.6 Mb on chromosome 1 (linked to <italic>MTL/ZmPLA1/NLD</italic> gene), S8_174792234, located at position 174.7 Mb on chromosome 8 and S7_131226237, located at position 131.2 Mb on chromosome 7, with higher PVE values of 37%, 25% and 21%, respectively, are likely to be more important for the trait and are more likely to have a direct impact on the HIR. Notably, detection of S1_66636144 at position 66.6 Mb on chromosome 1 by BLINK and MLMM models makes it as potential region for investigation, in addition to SNP S8_174792234, at position 174.7 Mb on chromosome 8 which is reported first time for the HIR. A smaller p-value suggests that the association is less likely to be due to chance and more likely to be real. Apart from S1_66636144 at 66.6 MB on chromosome 1, the marker S8_174792234 at position 174.7 Mb on chromosome 8 possessed the lowest p-value followed by S1_69321282 on chromosome 1 at 69.3 Mb (tightly linked to <italic>MTL/ZmPLA1/NLD</italic> gene) compared to the two other markers (S1_66636144 and S1_76160596 at 66.6 and 76.1 Mb, respectively) detected on chromosome 1. Given their high significance, it is interesting that they have not been detected in other mapping studies (<xref ref-type="bibr" rid="B19">Deimling et&#xa0;al., 1997</xref>; <xref ref-type="bibr" rid="B5">Barret et&#xa0;al., 2008</xref>; <xref ref-type="bibr" rid="B64">Prigge et&#xa0;al., 2012b</xref>; <xref ref-type="bibr" rid="B31">Hu et&#xa0;al., 2016</xref>). This failure of detection might have been caused by the absence of genetic variation at these QTL in the germplasm used in other studies.</p>
<p>In a cc-GWAS, we compared inducers (cases) to noninducers (controls) to identify genetic variants that are associated with the haploid induction. As a result, we identified a significant novel genetic region, close to marker S10_141729953, which is located at position 141.7 Mb of chromosome 10. Therefore, the identification of S10_141729953 as associated with HIR can help improve our understanding of the genetic factors involved in this trait. The failure of detection of this region might have been due to the absence of genetic variation at these QTL in the germplasm used in other studies (<xref ref-type="bibr" rid="B45">Lashermes and Beckert, 1988</xref>; <xref ref-type="bibr" rid="B19">Deimling et&#xa0;al., 1997</xref>; <xref ref-type="bibr" rid="B68">R&#xf6;ber, 1999</xref>; <xref ref-type="bibr" rid="B64">Prigge et&#xa0;al., 2012b</xref>; <xref ref-type="bibr" rid="B52">Liu et&#xa0;al., 2015</xref>). For marker S10_141729953, this seems to be a plausible explanation since the G allele was almost exclusively found in MHI and in its progeny. Another explanation might be that their effect is dependent on the presence of the <italic>mtl/zmpla1/nld</italic> allele, similarly to what is described for the <italic>zmdmp</italic> allele. This also seems to be the case for inducers carrying the G allele at marker S10_141729953. From those, the only two which do not have the common inducer haplotype at the <italic>MTL/ZmPLA1/NLD</italic> region (HUT45 and HUT71), display an HIR&lt; 1.0% at all planting dates. The population structure for the cc-GWAS analysis was only controlled with the additive genomic relationship matrix due to the confounding between the cases and controls grouping of genotypes with the first two principal components. Thus, further research is needed to ensure that marker S10_141729953 was not a spurious association.</p>
</sec>
<sec id="s4_2">
<title>Candidate genes</title>
<p>We searched the putative candidate genes in the interval of 56kb on both sides of nine significant markers identified in the study. The aldose 1-epimerase enzyme-coding gene Zm00001d029341 has been shown to be expressed in meiotic tassel (<xref ref-type="bibr" rid="B78">Stelpflug et&#xa0;al., 2016</xref>). This could be one of the candidate genes for inducing haploid in maize. The gene Zm00001d029558, which encodes for a pathogenesis-related protein 1, contains the marker S1_76160596. Pollen allergen properties of pathogenesis-related proteins have been reported in a variety of species (<xref ref-type="bibr" rid="B32">Huang et&#xa0;al., 1997</xref> and <xref ref-type="bibr" rid="B33">Huang et&#xa0;al., 2016</xref>). However, there was no evidence of Zm00001d029558 expression in pollen (<xref ref-type="bibr" rid="B85">Walley et&#xa0;al., 2016</xref>). S1_76160596, the same marker, is found 19.9 kb downstream of the gene Zm00001d029559, which encodes an EID1-like F-box protein 2. EID1 (Empfindlicher im Dunkelroten Licht 1) is a F-box protein that belongs to the SKP1-Cullin1-F-box (SCF) proteomic complexes that target protein degradation (<xref ref-type="bibr" rid="B53">Marrocco et&#xa0;al., 2006</xref>). F-box-like protein 17 (FBL17) in Arabidopsis targets the degradation of two CDKA;1 inhibitors, KRP6 and KRP7 (<xref ref-type="bibr" rid="B43">Kim et&#xa0;al., 2008</xref>). Arabidopsis fbl17 mutants occasionally fail to divide and produce pollen grains with a single sperm cell, resulting in single fertilization events (<xref ref-type="bibr" rid="B43">Kim et&#xa0;al., 2008</xref>). A-type cyclin-dependent kinase 1 (CDKA;1) is a homolog of the human serine/threonine protein kinase cdc2. Arabidopsis CDKA;1 mutation delay pollen generative cell division (<xref ref-type="bibr" rid="B4">Aw et&#xa0;al., 2010</xref>), resulting in the delivery of a single sperm cell to the ovary. It is possible that gene Zm00001d029559 interferes with maize&#x2019;s CDKA;1 via the same pathway described in Arabidopsis, influencing HIR by increasing the frequency of single fertilization events, which is one of two processes thought to cause haploid induction.</p>
<p>The marker S10_141729953 is in the intragenic region of gene Zm00001d026242, which was reported to be either exclusively or highly expressed in pollen grains (<xref ref-type="bibr" rid="B18">Davidson et&#xa0;al., 2011</xref>; <xref ref-type="bibr" rid="B85">Walley et&#xa0;al., 2016</xref>). Zm00001d026242 currently does not have an annotation, but the National Center for Biotechnology Information (NCBI) nucleotide Basic Local Alignment Search Tool (BLAST) shows that its nucleotide sequence is 88.9% and 83.0% identical to the <italic>Sorghum bicolor</italic>&#x2019;s and maize&#x2019;s <italic>kokopelli</italic> (<italic>KPL</italic>) mRNAs, respectively. <italic>Kokopelli</italic> was first described in <italic>Arabidopsis thaliana</italic>, where it was named after the Native American deity of male fertility (<xref ref-type="bibr" rid="B71">Ron et&#xa0;al., 2010</xref>). In maize, <italic>kokopelli</italic> is located on chromosome 2 (NCBI&#x2019;s sequence ID: XP_008671112.1), whereas in sorghum it is located on chromosome 6 (NCBI&#x2019;s sequence ID: XP_021319809.1). <italic>Kokopelli</italic> encodes for a natural cis-antisense siRNA (cis-nat-siRNA), which is a class of small regulatory RNAs found in various eukaryotic species. Even though two apparently normal sperms are delivered to embryos sacs by the pollen tubes of the Arabidopsis <italic>kpl</italic> mutant, single fertilization was observed in approximately 40% of the ovules fertilized (<xref ref-type="bibr" rid="B71">Ron et&#xa0;al., 2010</xref>). In Arabidopsis, <italic>KPL</italic> forms a sperm-specific nat-siRNA pair with the inversely transcribed gene <italic>ARIADNE14</italic> (<italic>ARI14</italic>), which encodes for a putative ubiquitin E3 ligase, and whose overexpression decreases seed set (<xref ref-type="bibr" rid="B71">Ron et&#xa0;al., 2010</xref>). Protein degradation occurs during gametogenesis in plants and animals and is required for successful fertilization (<xref ref-type="bibr" rid="B23">Dickins et&#xa0;al., 2002</xref>; <xref ref-type="bibr" rid="B79">Suzumori et&#xa0;al., 2003</xref>; <xref ref-type="bibr" rid="B70">Roest et&#xa0;al., 2004</xref>; <xref ref-type="bibr" rid="B71">Ron et&#xa0;al., 2010</xref>; <xref ref-type="bibr" rid="B86">Wang et&#xa0;al., 2013</xref>). Altogether, these results suggest a possible involvement of Zm00001d026242 in fertilization and haploid induction. Recently, <xref ref-type="bibr" rid="B35">Jacquier et&#xa0;al. (2023)</xref> validated the role of <italic>Atkpl</italic> (KOKOPELLI, AT5G63720) mutant for triggering <italic>in planta</italic> maternal haploid induction.</p>
</sec>
<sec id="s4_3">
<title>Applicability of the findings to inducer breeding</title>
<p>The detection of candidate genes in six maize chromosomes is in agreement with other studies that suggested a polygenic control of HIR (<xref ref-type="bibr" rid="B45">Lashermes and Beckert, 1988</xref>; <xref ref-type="bibr" rid="B19">Deimling et&#xa0;al., 1997</xref>; <xref ref-type="bibr" rid="B68">R&#xf6;ber, 1999</xref>; <xref ref-type="bibr" rid="B64">Prigge et&#xa0;al., 2012b</xref>; <xref ref-type="bibr" rid="B52">Liu et&#xa0;al., 2015</xref>; <xref ref-type="bibr" rid="B59">Nair et&#xa0;al., 2017</xref>). The presence of genes with strong effects, like <italic>MTL/ZmPLA1/NLD</italic>, <italic>ZmDMP</italic>, <italic>ZmPLD3</italic>, <italic>ZmPOD65</italic> and the QTL identified in this study on chromosomes 1, 8 and 10, suggest that some form of MAS, such as F<sub>2</sub>-enrichment or marker-assisted backcrossing, might be efficient breeding strategies for inducer development. Employing genomic selection or predicted cross value (<xref ref-type="bibr" rid="B30">Han et&#xa0;al., 2017</xref>) during or after the fixation of these major genes, can help to capture the best combination of small effect QTL (<xref ref-type="bibr" rid="B3">Almeida et&#xa0;al., 2020</xref>). Phenotypic selection for HIR is very laborious and time-consuming, since multiple cross-pollinations need to be performed to obtain reliable estimates of HIR. Moreover, visual discrimination of haploid and diploid seeds is time-consuming and error-prone, since <italic>R1-nj</italic> expression is affected by environmental conditions, seed morphology and anthocyanin inhibitor genes present in the donor background (<xref ref-type="bibr" rid="B9">Burr et&#xa0;al., 1996</xref>; <xref ref-type="bibr" rid="B21">Della Vedova et&#xa0;al., 2005</xref>; <xref ref-type="bibr" rid="B69">R&#xf6;ber et&#xa0;al., 2005</xref>; <xref ref-type="bibr" rid="B39">Kebede et&#xa0;al., 2011</xref>; <xref ref-type="bibr" rid="B62">Prigge et&#xa0;al., 2011</xref>; <xref ref-type="bibr" rid="B63">Prigge et&#xa0;al., 2012a</xref>; <xref ref-type="bibr" rid="B64">Prigge et&#xa0;al., 2012b</xref>).</p>
<p>The development of high oil content (OC) inducers allows automated haploid selection, which greatly facilitates accurate estimation of HIRs. However, OC is a quantitatively inherited character (<xref ref-type="bibr" rid="B56">Moreno-Gonzalez et&#xa0;al., 1975</xref>; <xref ref-type="bibr" rid="B7">Berke and Rocheford, 1995</xref>; <xref ref-type="bibr" rid="B46">Laurie et&#xa0;al., 2004</xref>; <xref ref-type="bibr" rid="B57">M&#xf4;ro et&#xa0;al., 2012</xref>), as are other traits important to inducers, such as plant height, tassel size and disease tolerance (<xref ref-type="bibr" rid="B26">Gallavotti et&#xa0;al., 2004</xref>; <xref ref-type="bibr" rid="B74">Satoh-Nagasawa et&#xa0;al., 2006</xref>; <xref ref-type="bibr" rid="B60">Peiffer et&#xa0;al., 2014</xref>; <xref ref-type="bibr" rid="B80">Tembo et&#xa0;al., 2014</xref>; <xref ref-type="bibr" rid="B89">Zhang et&#xa0;al., 2016</xref>; <xref ref-type="bibr" rid="B3">Almeida et&#xa0;al., 2020</xref>). Employing genomic selection to simultaneously improve multiple polygenic traits using optimization criteria (<xref ref-type="bibr" rid="B2">Akdemir et&#xa0;al., 2019</xref>) dilutes the costs of genotyping and becomes more cost-effective as the number of traits that otherwise would have to be phenotyped increase. Nevertheless, visual discrimination of haploid seeds through the <italic>R1-nj</italic> anthocyanin marker is still performed in the majority of medium to small scale breeding programs employing the DH technique. Therefore, the economic benefits of using inducers with high HIR is still important.</p>
</sec>
</sec>
<sec id="s5" sec-type="data-availability">
<title>Data availability statement</title>
<p>The original contributions presented in the study are included in the article/<xref ref-type="supplementary-material" rid="SF1"><bold>Supplementary Material</bold></xref>. The genetic polymorphism data is available at <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.25380/iastate.24185784.v2">https://doi.org/10.25380/iastate.24185784.v2</ext-link>. Further inquiries can be directed to the corresponding author.</p>
</sec>
<sec id="s6" sec-type="author-contributions">
<title>Author contributions</title>
<p>HT, UF, and TL designed this project; HT, UF, VR, and TL developed the inducer panel used in this study; HT, VA, and EP performed most of the HIR phenotypic selection; MK and WB analyzed the phenotypic and genomic data; HT and RZ executed the candidate gene analysis; HT, MK, RZ and TL wrote this manuscript; TL supervised the research work; All authors contributed to the article and approved the submitted version.</p>
</sec>
</body>
<back>
<sec id="s7" sec-type="funding-information">
<title>Funding</title>
<p>This research was funded by USDA&#x2019;s National Institute of Food and Agriculture (NIFA) Project, No. IOW04314, IOW01018, and IOW05510; and NIFA award 2018-51181-28419. Funding for this work was also provided by the R.F. Baker Center for Plant Breeding, Plant Sciences Institute, and K.J. Frey Chair in Agronomy, at Iowa State University. The Brazilian National Council for Scientific and Technological Development (CNPq) also provided funding for this work.</p>
</sec>
<ack>
<title>Acknowledgments</title>
<p>RZ gratefully acknowledges the funds provided by SERB-International Research Experience Fellowship, Department of Science and Technology, India (File No. SIR/2022/000112). We appreciate Dr. Thomas Widiez for his constructive and insightful remarks on the manuscript.</p>
</ack>
<sec id="s8" sec-type="COI-statement">
<title>Conflict of interest</title>
<p>Author VR was employed by VAR BREEDING Ltd.</p>
<p>The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
<sec id="s9" sec-type="disclaimer">
<title>Publisher&#x2019;s note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec>
<sec id="s10" sec-type="supplementary-material">
<title>Supplementary material</title>
<p>The Supplementary Material for this article can be found online at: <ext-link ext-link-type="uri" xlink:href="https://www.frontiersin.org/articles/10.3389/fpls.2023.1218042/full#supplementary-material">https://www.frontiersin.org/articles/10.3389/fpls.2023.1218042/full#supplementary-material</ext-link>
</p>
<supplementary-material xlink:href="Image_1.jpeg" id="SF1" mimetype="image/jpeg">
<label>Supplementary Figure&#xa0;1</label>
<caption>
<p>Half-normal plots with a simulated envelope for <bold>(A)</bold> Model 1, <bold>(B)</bold> Model 1 with Quasi-likelihood inference, <bold>(C)</bold> Model 2, <bold>(D)</bold> Model 3, and <bold>(E)</bold> Model 4.</p>
</caption>
</supplementary-material>
<supplementary-material xlink:href="Image_2.jpeg" id="SF2" mimetype="image/jpeg">
<label>Supplementary Figure&#xa0;2</label>
<caption>
<p>Index Pearson residual plots for <bold>(A)</bold> Model 1, <bold>(B)</bold> Model 2, <bold>(C)</bold> Model 3, and <bold>(D)</bold> Model 4.</p>
</caption>
</supplementary-material>
<supplementary-material xlink:href="Table_1.xlsx" id="SF3" mimetype="application/vnd.openxmlformats-officedocument.spreadsheetml.sheet"/>
<supplementary-material xlink:href="Table_2.xlsx" id="SF4" mimetype="application/vnd.openxmlformats-officedocument.spreadsheetml.sheet"/>
<supplementary-material xlink:href="Table_3.docx" id="SF5" mimetype="application/vnd.openxmlformats-officedocument.wordprocessingml.document"/>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Akaike</surname> <given-names>H.</given-names>
</name>
</person-group> (<year>1974</year>). <article-title>A new look at the statistical model identification</article-title>. <source>IEEE Trans. Automatic Control</source> <volume>19</volume>, <fpage>716</fpage>&#x2013;<lpage>723</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1109/TAC.1974.1100705</pub-id>
</citation>
</ref>
<ref id="B2">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Akdemir</surname> <given-names>D.</given-names>
</name>
<name>
<surname>Beavis</surname> <given-names>W.</given-names>
</name>
<name>
<surname>Fritsche-Neto</surname> <given-names>R.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Multi-objective optimized genomic breeding strategies for sustainable food improvement</article-title>. <source>Heredity</source> <volume>122</volume>, <fpage>672</fpage>&#x2013;<lpage>683</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1038/s41437-018-0147-1</pub-id>
</citation>
</ref>
<ref id="B3">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Almeida</surname> <given-names>V. C.</given-names>
</name>
<name>
<surname>Trentin</surname> <given-names>H. U.</given-names>
</name>
<name>
<surname>Frei</surname> <given-names>U. K.</given-names>
</name>
<name>
<surname>L&#xfc;bberstedt</surname> <given-names>T.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Genomic prediction of maternal haploid induction rate in maize</article-title>. <source>Plant Genome</source> <volume>13</volume> (<issue>1</issue>), <elocation-id>e20014</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.1002/tpg2.20014</pub-id>
</citation>
</ref>
<ref id="B4">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Aw</surname> <given-names>S. J.</given-names>
</name>
<name>
<surname>Hamamura</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Chen</surname> <given-names>Z.</given-names>
</name>
<name>
<surname>Schnittger</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Bergeret</surname> <given-names>F.</given-names>
</name>
</person-group> (<year>2010</year>). <article-title>Sperm entry is sufficient to trigger division of the central cell but the paternal genome is required for endosperm development in Arabidopsis</article-title>. <source>Development</source> <volume>137</volume>, <fpage>2683</fpage>&#x2013;<lpage>2690</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1242/dev.052928</pub-id>
</citation>
</ref>
<ref id="B5">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Barret</surname> <given-names>P.</given-names>
</name>
<name>
<surname>Brinkmann</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Beckert</surname> <given-names>M.</given-names>
</name>
</person-group> (<year>2008</year>). <article-title>A major locus expressed in the male gametophyte with incomplete penetrance is responsible for in <italic>situ</italic> gynogenesis in maize</article-title>. <source>Theor. Appl. Genet.</source> <volume>117</volume>, <fpage>581</fpage>&#x2013;<lpage>594</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1007/s00122-008-0803-6</pub-id>
</citation>
</ref>
<ref id="B6">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bates</surname> <given-names>D.</given-names>
</name>
<name>
<surname>M&#xe4;chler</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Bolker</surname> <given-names>B.</given-names>
</name>
<name>
<surname>Walker</surname> <given-names>S.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Fitting linear mixed-effects models using lme4</article-title>. <source>J. Stat.</source> <volume>67</volume>, <fpage>1</fpage>&#x2013;<lpage>48</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.18637/jss.v067.i01</pub-id>
</citation>
</ref>
<ref id="B7">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Berke</surname> <given-names>T. G.</given-names>
</name>
<name>
<surname>Rocheford</surname> <given-names>T. R.</given-names>
</name>
</person-group> (<year>1995</year>). <article-title>Quantitative trait loci for flowering, plant and ear height, and kernel traits in maize</article-title>. <source>Crop Sci.</source> <volume>35</volume>, <fpage>1542</fpage>&#x2013;<lpage>1549</lpage>. doi: <pub-id pub-id-type="doi">10.2135/cropsci1995.0011183X003500060004x</pub-id>
</citation>
</ref>
<ref id="B8">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Browning</surname> <given-names>B. L.</given-names>
</name>
<name>
<surname>Zhou</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Browning</surname> <given-names>S. R.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>A one-penny imputed genome from next-generation reference panels</article-title>. <source>Am. J. Hum. Genet.</source> <volume>103</volume>, <fpage>338</fpage>&#x2013;<lpage>348</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.ajhg.2018.07.015</pub-id>
</citation>
</ref>
<ref id="B9">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Burr</surname> <given-names>F.</given-names>
</name>
<name>
<surname>Burr</surname> <given-names>B.</given-names>
</name>
<name>
<surname>Scheffler</surname> <given-names>B. E.</given-names>
</name>
<name>
<surname>Blewitt</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Wienand</surname> <given-names>U.</given-names>
</name>
<name>
<surname>Matzet</surname> <given-names>E. C.</given-names>
</name>
</person-group> (<year>1996</year>). <article-title>The maize repressor-like gene intensifier1 shares homology with the r1/b1 multigene family of transcription factors and exhibits missplicing</article-title>. <source>Plant Cell</source> <volume>8</volume>, <fpage>1249</fpage>&#x2013;<lpage>1259</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1105/tpc.8.8.1249</pub-id>
</citation>
</ref>
<ref id="B10">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chalyk</surname> <given-names>S. T.</given-names>
</name>
</person-group> (<year>1999</year>). <article-title>Creating new haploid-inducing lines of maize</article-title>. <source>Maize Genet. Coop. Newsl.</source> <volume>73</volume>, <fpage>53</fpage>&#x2013;<lpage>54</lpage>.</citation>
</ref>
<ref id="B11">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chase</surname> <given-names>S. S.</given-names>
</name>
</person-group> (<year>1947</year>). <article-title>Techniques for isolating haploid plants</article-title>. <source>Am. J. Bot.</source> <volume>34</volume>, <fpage>579</fpage>&#x2013;<lpage>609</lpage>.</citation>
</ref>
<ref id="B12">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chase</surname> <given-names>S. S.</given-names>
</name>
</person-group> (<year>1952</year>). <article-title>Production of homozygous diploids of maize from monoploids 1</article-title>. <source>Agron. J.</source> <volume>44</volume>, <fpage>263</fpage>&#x2013;<lpage>267</lpage>. doi: <pub-id pub-id-type="doi">10.2134/agronj1952.00021962004400050010x</pub-id>
</citation>
</ref>
<ref id="B13">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chase</surname> <given-names>S. S.</given-names>
</name>
<name>
<surname>Nanda</surname> <given-names>D. K.</given-names>
</name>
</person-group> (<year>1965</year>). <article-title>Screening for monoploids of maize by use of a purple embryo marker</article-title>. <source>Maize Genet. Coop. Newsl.</source> <volume>39</volume>, <fpage>59</fpage>&#x2013;<lpage>60</lpage>.</citation>
</ref>
<ref id="B14">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chen</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Liu</surname> <given-names>X.</given-names>
</name>
<name>
<surname>Li</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Liu</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Zhang</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Luo</surname> <given-names>L.</given-names>
</name>
<etal/>
</person-group>. (<year>2022</year>). <article-title>Co-expression of transcription factors ZmC1 and ZmR2 establishes an efficient and accurate haploid embryo identification system in maize</article-title>. <source>Plant J.</source> <volume>111</volume>, <fpage>1296</fpage>&#x2013;<lpage>1307</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1111/tpj.15888</pub-id>
</citation>
</ref>
<ref id="B15">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Coe</surname> <given-names>E. H.</given-names>
</name>
</person-group> (<year>1959</year>). <article-title>A line of maize with high haploid frequency</article-title>. <source>Am. Nat.</source> <volume>93</volume>, <fpage>381</fpage>&#x2013;<lpage>382</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1086/282187</pub-id>
</citation>
</ref>
<ref id="B16">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Couto</surname> <given-names>E. G. O.</given-names>
</name>
<name>
<surname>Cury</surname> <given-names>M. N.</given-names>
</name>
<name>
<surname>Bandeira, Souza</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Granato</surname> <given-names>I. S. C.</given-names>
</name>
<name>
<surname>Vidotti</surname> <given-names>M. S.</given-names>
</name>
<name>
<surname>Domingos, Garbuglio</surname> <given-names>D.</given-names>
</name>
<etal/>
</person-group>. (<year>2019</year>). <article-title>Effect of F1 and F2 generations on genetic variability and working steps of doubled haploid production in maize</article-title>. <source>PloS One</source> <volume>14</volume>, <elocation-id>e0224631</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.1371/journal.pone.0224631</pub-id>
</citation>
</ref>
<ref id="B17">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Covarrubias-Pazaran</surname> <given-names>G.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>Genome-assisted prediction of quantitative traits using the R package sommer</article-title>. <source>PloS One</source> <volume>11</volume>, <elocation-id>e0156744</elocation-id>. doi: <pub-id pub-id-type="doi">10.1371/journal.pone.0156744</pub-id>
</citation>
</ref>
<ref id="B18">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Davidson</surname> <given-names>R. M.</given-names>
</name>
<name>
<surname>Hansey</surname> <given-names>C. N.</given-names>
</name>
<name>
<surname>Gowda</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Childs</surname> <given-names>K. L.</given-names>
</name>
<name>
<surname>Lin</surname> <given-names>H.</given-names>
</name>
<name>
<surname>Vaillancourt</surname> <given-names>B.</given-names>
</name>
<etal/>
</person-group>. (<year>2011</year>). <article-title>Utility of RNA sequencing for analysis of maize reproductive transcriptomes</article-title>. <source>Plant Genome</source> <volume>4</volume>, <fpage>191</fpage>&#x2013;<lpage>203</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.3835/plantgenome2011.05.0015</pub-id>
</citation>
</ref>
<ref id="B19">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Deimling</surname> <given-names>S.</given-names>
</name>
<name>
<surname>R&#xf6;ber</surname> <given-names>F. K.</given-names>
</name>
<name>
<surname>Geiger</surname> <given-names>H. H.</given-names>
</name>
</person-group> (<year>1997</year>). <article-title>Methodology and genetics of in <italic>vivo</italic> haploid induction in maize</article-title>. <source>Vortr&#xe4;ge f&#xfc;r Pflanzenz&#xfc;chtung</source> <volume>38</volume>, <fpage>203</fpage>&#x2013;<lpage>224</lpage>.</citation>
</ref>
<ref id="B20">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>De La Fuente</surname> <given-names>G. N.</given-names>
</name>
<name>
<surname>Frei</surname> <given-names>U. K.</given-names>
</name>
<name>
<surname>Trampe</surname> <given-names>B.</given-names>
</name>
<name>
<surname>Nettleton</surname> <given-names>D.</given-names>
</name>
<name>
<surname>Zhang</surname> <given-names>W.</given-names>
</name>
<name>
<surname>L&#xfc;bberstedt</surname> <given-names>T.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>A diallel analysis of a maize donor population response to in <italic>vivo</italic> maternal haploid induction: I. inducibility</article-title>. <source>Crop Sci.</source> <volume>58</volume>, <fpage>1830</fpage>&#x2013;<lpage>1837</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.2135/cropsci2017.05.0285</pub-id>
</citation>
</ref>
<ref id="B21">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Della Vedova</surname> <given-names>C. B.</given-names>
</name>
<name>
<surname>Lorbiecke</surname> <given-names>R.</given-names>
</name>
<name>
<surname>Kirsch</surname> <given-names>H.</given-names>
</name>
<name>
<surname>Schulte</surname> <given-names>M. B.</given-names>
</name>
<name>
<surname>Scheets</surname> <given-names>K.</given-names>
</name>
<name>
<surname>Borchert</surname> <given-names>L. M.</given-names>
</name>
<etal/>
</person-group>. (<year>2005</year>). <article-title>The dominant inhibitory chalcone synthase allele C2-Idf (inhibitor diffuse) from <italic>Zea mays</italic> (L.) acts via an endogenous RNA silencing mechanism</article-title>. <source>Genetics</source> <volume>170</volume>, <fpage>1989</fpage>&#x2013;<lpage>2002</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1534/genetics.105.043406</pub-id>
</citation>
</ref>
<ref id="B22">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Dem&#xe9;trio</surname> <given-names>C. G. B.</given-names>
</name>
<name>
<surname>Hinde</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Moral</surname> <given-names>R. A.</given-names>
</name>
</person-group> (<year>2014</year>). <source>Models for overdispersed data in entomology bt - ecological modelling applied to entomology</source>. Eds. <person-group person-group-type="editor">
<name>
<surname>Ferreira</surname> <given-names>C. P.</given-names>
</name>
<name>
<surname>Godoy</surname> <given-names>W. A. C.</given-names>
</name>
</person-group> (<publisher-loc>Cham</publisher-loc>: <publisher-name>Springer International Publishing</publisher-name>), <fpage>219</fpage>&#x2013;<lpage>259</lpage>.</citation>
</ref>
<ref id="B23">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Dickins</surname> <given-names>R. A.</given-names>
</name>
<name>
<surname>Frew</surname> <given-names>I. J.</given-names>
</name>
<name>
<surname>House</surname> <given-names>C. M.</given-names>
</name>
<name>
<surname>O'Bryan</surname> <given-names>M. K.</given-names>
</name>
<name>
<surname>Holloway</surname> <given-names>A. J.</given-names>
</name>
<name>
<surname>Haviv</surname> <given-names>I.</given-names>
</name>
<etal/>
</person-group>. (<year>2002</year>). <article-title>The ubiquitin ligase component Siah1a is required for completion of meiosis I in male mice</article-title>. <source>Mol. Cell. Biol.</source> <volume>22</volume>, <fpage>2294</fpage>&#x2013;<lpage>2303</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1128/MCB.22.7.2294-2303.2002</pub-id>
</citation>
</ref>
<ref id="B24">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Dong</surname> <given-names>X.</given-names>
</name>
<name>
<surname>Xu</surname> <given-names>X.</given-names>
</name>
<name>
<surname>Miao</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Li</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Zhang</surname> <given-names>D.</given-names>
</name>
<name>
<surname>Mi</surname> <given-names>X.</given-names>
</name>
<etal/>
</person-group>. (<year>2013</year>). <article-title>Fine mapping of qhir1 influencing in <italic>vivo</italic> haploid induction in maize</article-title>. <source>Theor. Appl. Genet.</source> <volume>126</volume>, <fpage>1713</fpage>&#x2013;<lpage>1720</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1007/s00122-013-2086-9</pub-id>
</citation>
</ref>
<ref id="B25">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Eder</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Chalyk</surname> <given-names>S.</given-names>
</name>
</person-group> (<year>2002</year>). <article-title>
<italic>In vivo</italic> haploid induction in maize</article-title>. <source>Theor. Appl. Genet.</source> <volume>104</volume>, <fpage>703</fpage>&#x2013;<lpage>708</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1007/s00122-001-0773-4</pub-id>
</citation>
</ref>
<ref id="B26">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Gallavotti</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Zhao</surname> <given-names>Q.</given-names>
</name>
<name>
<surname>Kyozuka</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Meeley</surname> <given-names>R. B.</given-names>
</name>
</person-group> (<year>2004</year>). <article-title>The role of barren stalk1 in the architecture of maize</article-title>. <source>Nature</source> <volume>432</volume>, <fpage>1</fpage>&#x2013;<lpage>6</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1038/nature03122.1</pub-id>
</citation>
</ref>
<ref id="B27">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Gilles</surname> <given-names>L. M.</given-names>
</name>
<name>
<surname>Calhau</surname> <given-names>A. R. M.</given-names>
</name>
<name>
<surname>Padula</surname> <given-names>V. L.</given-names>
</name>
<name>
<surname>Jacquier</surname> <given-names>N. M.A.</given-names>
</name>
<name>
<surname>Lionnet</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Martinant</surname> <given-names>J.</given-names>
</name>
<etal/>
</person-group>. (<year>2021</year>). <article-title>Lipid anchoring and electrostatic interactions target NOT-LIKE-DAD to pollen endo-plasma membrane</article-title>. <source>J. Cell Biol.</source> <volume>220</volume> (<issue>10</issue>), <elocation-id>e202010077</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.1083/jcb.202010077</pub-id>
</citation>
</ref>
<ref id="B28">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Gilles</surname> <given-names>L. M.</given-names>
</name>
<name>
<surname>Khaled</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Laffaire</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Chaignon</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Gendrot</surname> <given-names>G.</given-names>
</name>
<name>
<surname>Laplaige</surname> <given-names>J.</given-names>
</name>
<etal/>
</person-group>. (<year>2017</year>b). <article-title>Loss of pollen-specific phospholipase NOT LIKE DAD triggers gynogenesis in maize</article-title>. <source>EMBO J.</source> <volume>36</volume>, <fpage>707</fpage>&#x2013;<lpage>717</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.15252/embj.201796603</pub-id>
</citation>
</ref>
<ref id="B29">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Gilles</surname> <given-names>L. M.</given-names>
</name>
<name>
<surname>Martinant</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Rogowsky</surname> <given-names>P. M.</given-names>
</name>
<name>
<surname>Widiez</surname> <given-names>T.</given-names>
</name>
</person-group> (<year>2017</year>a). <article-title>Haploid induction in plants</article-title>. <source>Curr. Biol.</source> <volume>27</volume> (<issue>20</issue>), <fpage>R1095</fpage>&#x2013;<lpage>R1097</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.cub.2017.07.055</pub-id>
</citation>
</ref>
<ref id="B30">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Han</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Cameron</surname> <given-names>J. N.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Beavis</surname> <given-names>W. D.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>The predicted cross value for genetic introgression of multiple alleles</article-title>. <source>Genetics</source> <volume>205</volume> (<issue>4</issue>), <fpage>1409</fpage>&#x2013;<lpage>1423</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1534/genetics.116.197095</pub-id>
</citation>
</ref>
<ref id="B31">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hu</surname> <given-names>H.</given-names>
</name>
<name>
<surname>Schrag</surname> <given-names>T. A.</given-names>
</name>
<name>
<surname>Peis</surname> <given-names>R.</given-names>
</name>
<name>
<surname>Unterseer</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Schipprack</surname> <given-names>W.</given-names>
</name>
<name>
<surname>Chen</surname> <given-names>S.</given-names>
</name>
<etal/>
</person-group>. (<year>2016</year>). <article-title>The genetic basis of haploid induction in maize identified with a novel genome-wide association method</article-title>. <source>Genetics</source> <volume>202</volume>, <fpage>1267</fpage>&#x2013;<lpage>1276</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1534/genetics.115.184234</pub-id>
</citation>
</ref>
<ref id="B32">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Huang</surname> <given-names>J. C.</given-names>
</name>
<name>
<surname>Chang</surname> <given-names>F. C.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>C. S.</given-names>
</name>
</person-group> (<year>1997</year>). <article-title>Characterization of a lily tapetal transcript that shares sequence similarity with a class of intracellular pathogenesis-related (IPR) proteins</article-title>. <source>Plant Mol. Biol.</source> <volume>34</volume>, <fpage>681</fpage>&#x2013;<lpage>686</lpage>. doi: <pub-id pub-id-type="doi">10.1023/A:1005824306560</pub-id>
</citation>
</ref>
<ref id="B33">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Huang</surname> <given-names>L. F.</given-names>
</name>
<name>
<surname>Lin</surname> <given-names>K. H.</given-names>
</name>
<name>
<surname>He</surname> <given-names>S. L.</given-names>
</name>
<name>
<surname>Chen</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Jiang</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Chen</surname> <given-names>B.</given-names>
</name>
<etal/>
</person-group>. (<year>2016</year>). <article-title>Multiple patterns of regulation and overexpression of a ribonuclease-like pathogenesis-related protein gene, OsPR10a, conferring disease resistance in rice and Arabidopsis</article-title>. <source>PloS One</source> <volume>11</volume>, <elocation-id>e0156414</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.1371/journal.pone.0156414</pub-id>
</citation>
</ref>
<ref id="B34">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Huang</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Liu</surname> <given-names>X.</given-names>
</name>
<name>
<surname>Zhou</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Summers</surname> <given-names>R. M.</given-names>
</name>
<name>
<surname>Zhang</surname> <given-names>Z.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>BLINK: a package for the next level of genome-wide association studies with both individuals and markers in the millions</article-title>. <source>Gigascience</source> <volume>8</volume>, <elocation-id>154</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.1093/gigascience/giy154</pub-id>
</citation>
</ref>
<ref id="B35">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Jacquier</surname> <given-names>N. M. A.</given-names>
</name>
<name>
<surname>Calhau</surname> <given-names>A. R. M.</given-names>
</name>
<name>
<surname>Fierlej</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Martinant</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Rogowsky</surname> <given-names>P. M.</given-names>
</name>
<name>
<surname>Gilles</surname> <given-names>L. M.</given-names>
</name>
<etal/>
</person-group>. (<year>2023</year>). <article-title>In planta haploid induction by kokopelli mutants</article-title>. <source>Plant Physiol.</source> <volume>193</volume> (<issue>1</issue>), <fpage>182</fpage>&#x2013;<lpage>185</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1093/plphys/kiad328</pub-id>
</citation>
</ref>
<ref id="B36">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Jacquier</surname> <given-names>N. M. A.</given-names>
</name>
<name>
<surname>Gilles</surname> <given-names>L. M.</given-names>
</name>
<name>
<surname>Martinant</surname> <given-names>J. P.</given-names>
</name>
<name>
<surname>Rogowsky</surname> <given-names>P. M.</given-names>
</name>
<name>
<surname>Widiez</surname> <given-names>T.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Maize In Planta haploid inducer lines, a cornerstone for doubled haploid technology</article-title>. <source>Methods Mol. Biol.</source> <volume>2288</volume>, <fpage>25</fpage>&#x2013;<lpage>48</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1007/978-1-0716-1335-1_2</pub-id>
</citation>
</ref>
<ref id="B37">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Jacquier</surname> <given-names>N. M.</given-names>
</name>
<name>
<surname>Gilles</surname> <given-names>L. M.</given-names>
</name>
<name>
<surname>Pyott</surname> <given-names>D. E.</given-names>
</name>
<name>
<surname>Martinant</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Rogowsky</surname> <given-names>P. M.</given-names>
</name>
<name>
<surname>Widiez</surname> <given-names>T.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Puzzling out plant reproduction by haploid induction for innovations in plant breeding</article-title>. <source>Nat. Plants</source> <volume>6</volume> (<issue>6</issue>), <fpage>610</fpage>&#x2013;<lpage>619</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1038/s41477-020-0664-9</pub-id>
</citation>
</ref>
<ref id="B38">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Jiang</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Sun</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Li</surname> <given-names>R.</given-names>
</name>
<name>
<surname>Yan</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Chen</surname> <given-names>W.</given-names>
</name>
<name>
<surname>Guo</surname> <given-names>L.</given-names>
</name>
<etal/>
</person-group>. (<year>2022</year>). <article-title>A reactive oxygen species burst causes haploid induction in maize</article-title>. <source>Mol. Plant</source> <volume>15</volume> (<issue>6</issue>), <fpage>943</fpage>&#x2013;<lpage>955</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.molp.2022.04.001</pub-id>
</citation>
</ref>
<ref id="B39">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kebede</surname> <given-names>A. Z.</given-names>
</name>
<name>
<surname>Dhillon</surname> <given-names>B. S.</given-names>
</name>
<name>
<surname>Schipprack</surname> <given-names>W.</given-names>
</name>
<name>
<surname>Araus</surname> <given-names>J. L.</given-names>
</name>
<name>
<surname>B&#xe4;nziger</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Semagn</surname> <given-names>K.</given-names>
</name>
<etal/>
</person-group>. (<year>2011</year>). <article-title>Effect of source germplasm and season on the in <italic>vivo</italic> haploid induction rate in tropical maize</article-title>. <source>Euphytica</source> <volume>180</volume>, <fpage>219</fpage>&#x2013;<lpage>226</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1007/s10681-011-0376-3</pub-id>
</citation>
</ref>
<ref id="B40">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kelliher</surname> <given-names>T.</given-names>
</name>
<name>
<surname>Starr</surname> <given-names>D.</given-names>
</name>
<name>
<surname>Richbourg</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Chintamanani</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Delzer</surname> <given-names>B.</given-names>
</name>
<name>
<surname>Nuccio</surname> <given-names>M. L.</given-names>
</name>
<etal/>
</person-group>. (<year>2017</year>). <article-title>MATRILINEAL, a sperm-specific phospholipase, triggers maize haploid induction</article-title>. <source>Nature</source> <volume>542</volume>, <fpage>105</fpage>&#x2013;<lpage>109</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1038/nature20827</pub-id>
</citation>
</ref>
<ref id="B41">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Kermicle</surname> <given-names>J. L.</given-names>
</name>
</person-group> (<year>1969</year>). <article-title>Androgenesis conditioned by a mutation in maize abstract</article-title>. <source>Science</source> <volume>166</volume> (<issue>3911</issue>), <fpage>1422</fpage>&#x2013;<lpage>1424</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1126/science.166.3911.1422</pub-id>
</citation>
</ref>
<ref id="B42">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kilian</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Wenzl</surname> <given-names>P.</given-names>
</name>
<name>
<surname>Huttner</surname> <given-names>E.</given-names>
</name>
<name>
<surname>Carling</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Xia</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Blois</surname> <given-names>H.</given-names>
</name>
<etal/>
</person-group>. (<year>2012</year>). <article-title>Diversity Arrays Technology: A generic genome profiling technology on open platforms</article-title>. <source>Methods Mol. Biol.</source> <volume>888</volume>, <fpage>67</fpage>&#x2013;<lpage>89</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1007/978-1-61779-870-2_5</pub-id>
</citation>
</ref>
<ref id="B43">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kim</surname> <given-names>H. J.</given-names>
</name>
<name>
<surname>Oh</surname> <given-names>S. A.</given-names>
</name>
<name>
<surname>Brownfield</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Brownfield</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Hong</surname> <given-names>S. H.</given-names>
</name>
<name>
<surname>Ryu</surname> <given-names>H.</given-names>
</name>
<etal/>
</person-group>. (<year>2008</year>). <article-title>Control of plant germline proliferation by SCFFBL17 degradation of cell cycle inhibitors</article-title>. <source>Nature</source> <volume>455</volume>, <fpage>1134</fpage>&#x2013;<lpage>1137</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1038/nature07289</pub-id>
</citation>
</ref>
<ref id="B44">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kindiger</surname> <given-names>B.</given-names>
</name>
<name>
<surname>Hamann</surname> <given-names>S.</given-names>
</name>
</person-group> (<year>1993</year>). <article-title>Generation of haploids in maize: A modification of the indeterminate gametophyte (ig) system</article-title>. <source>Crop Sci.</source> <volume>33</volume>, <fpage>342</fpage>&#x2013;<lpage>344</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.2135/cropsci1993.0011183X003300020027x</pub-id>
</citation>
</ref>
<ref id="B45">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lashermes</surname> <given-names>P.</given-names>
</name>
<name>
<surname>Beckert</surname> <given-names>M.</given-names>
</name>
</person-group> (<year>1988</year>). <article-title>Genetic control of maternal haploidy in maize (<italic>Zea mays</italic> L.) and selection of haploid inducing lines</article-title>. <source>Theor. Appl. Genet.</source> <volume>76</volume>, <fpage>405</fpage>&#x2013;<lpage>410</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1007/BF00265341</pub-id>
</citation>
</ref>
<ref id="B46">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Laurie</surname> <given-names>C. C.</given-names>
</name>
<name>
<surname>Chasalow</surname> <given-names>S. D.</given-names>
</name>
<name>
<surname>Ledeaux</surname> <given-names>J. R.</given-names>
</name>
<name>
<surname>McCarroll</surname> <given-names>R.</given-names>
</name>
<name>
<surname>Bush</surname> <given-names>D.</given-names>
</name>
<name>
<surname>Hauge</surname> <given-names>B.</given-names>
</name>
<etal/>
</person-group>. (<year>2004</year>). <article-title>The genetic architecture of response to long-term artificial selection for oil concentration in the maize kernel</article-title>. <source>Genetics</source> <volume>2155</volume>, <fpage>2141</fpage>&#x2013;<lpage>2155</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1534/genetics.104.029686</pub-id>
</citation>
</ref>
<ref id="B47">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Lenth</surname> <given-names>R.</given-names>
</name>
</person-group> (<year>2020</year>). <source>Emmeans: estimated marginal means, aka least-squares means</source>, <italic>R package</italic> version 1.5.1. Available at: <uri xlink:href="https://cran.r-project.org/src/contrib/Archive/emmeans/">https://cran.r-project.org/src/contrib/Archive/emmeans/</uri>.</citation>
</ref>
<ref id="B48">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Li</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Lin</surname> <given-names>Z.</given-names>
</name>
<name>
<surname>Yue</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Zhao</surname> <given-names>H.</given-names>
</name>
<name>
<surname>Fei</surname> <given-names>X.</given-names>
</name>
<name>
<surname>Lizhu</surname> <given-names>E.</given-names>
</name>
<etal/>
</person-group>. (<year>2021</year>). <article-title>Loss-of-function alleles of ZmPLD3 cause haploid induction in maize</article-title>. <source>Nat. Plants</source> <volume>7</volume>, <fpage>1579</fpage>&#x2013;<lpage>1588</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1038/s41477-021-01037-2</pub-id>
</citation>
</ref>
<ref id="B49">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lipka</surname> <given-names>A. E.</given-names>
</name>
<name>
<surname>Tian</surname> <given-names>F.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>Q.</given-names>
</name>
<name>
<surname>Peiffer</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Li</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Bradbury</surname> <given-names>P. J.</given-names>
</name>
<etal/>
</person-group>. (<year>2012</year>). <article-title>GAPIT: genome association and prediction integrated tool</article-title>. <source>Bioinformatics</source> <volume>28</volume>, <fpage>2397</fpage>&#x2013;<lpage>2399</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1093/bioinformatics/bts444</pub-id>
</citation>
</ref>
<ref id="B50">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Liu</surname> <given-names>X.</given-names>
</name>
<name>
<surname>Huang</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Fan</surname> <given-names>B.</given-names>
</name>
<name>
<surname>Buckler</surname> <given-names>E. S.</given-names>
</name>
<name>
<surname>Zhang</surname> <given-names>Z.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>Iterative usage of fixed and random effect models for powerful and efficient genome-wide association studies</article-title>. <source>PloS Genet.</source> <volume>12</volume>, <elocation-id>e1005767</elocation-id>. doi: <pub-id pub-id-type="doi">10.1371/journal.pgen.1005767</pub-id>
</citation>
</ref>
<ref id="B51">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Liu</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Li</surname> <given-names>X.</given-names>
</name>
<name>
<surname>Meng</surname> <given-names>D.</given-names>
</name>
<name>
<surname>Zhong</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Chen</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Dong</surname> <given-names>X.</given-names>
</name>
<etal/>
</person-group>. (<year>2017</year>). <article-title>A 4-bp insertion at ZmPLA1 encoding a putative phospholipase a generates haploid induction in maize</article-title>. <source>Mol. Plant</source> <volume>10</volume>, <fpage>520</fpage>&#x2013;<lpage>522</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.molp.2017.01.011</pub-id>
</citation>
</ref>
<ref id="B52">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Liu</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Li</surname> <given-names>W.</given-names>
</name>
<name>
<surname>Zhong</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Li</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Zhang</surname> <given-names>D.</given-names>
</name>
<name>
<surname>Mi</surname> <given-names>X.</given-names>
</name>
<etal/>
</person-group>. (<year>2015</year>). <article-title>Fine mapping of qhir8 affecting in <italic>vivo</italic> haploid induction in maize</article-title>. <source>Theor. Appl. Genet.</source> <volume>128</volume>, <fpage>2507</fpage>&#x2013;<lpage>2515</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1007/s00122-015-2605-y</pub-id>
</citation>
</ref>
<ref id="B53">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Marrocco</surname> <given-names>K.</given-names>
</name>
<name>
<surname>Zhou</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Bury</surname> <given-names>E.</given-names>
</name>
<name>
<surname>Dieterle</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Funk</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Genschik</surname> <given-names>P.</given-names>
</name>
<etal/>
</person-group>. (<year>2006</year>). <article-title>Functional analysis of EID1, an F-box protein involved in phytochrome A-dependent light signal transduction</article-title>. <source>Plant J.</source> <volume>45</volume>, <fpage>423</fpage>&#x2013;<lpage>438</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1111/j.1365-313X.2005.02635.x</pub-id>
</citation>
</ref>
<ref id="B54">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Milet</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Courtin</surname> <given-names>D.</given-names>
</name>
<name>
<surname>Garcia</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Perdry</surname> <given-names>H.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Mixed logistic regression in genome-wide association studies</article-title>. <source>BMC Bioinfo.</source> <volume>21</volume>, <fpage>536</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1186/s12859-020-03862-2</pub-id>
</citation>
</ref>
<ref id="B55">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Moral</surname> <given-names>R. A.</given-names>
</name>
<name>
<surname>Hinde</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Dem&#xe9;trio</surname> <given-names>C. G. B.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Half-normal plots and overdispersed models in R: The hnp Package</article-title>. <source>J. Stat.</source> <volume>81</volume>, <fpage>1</fpage>&#x2013;<lpage>23</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.18637/jss.v081.i10</pub-id>
</citation>
</ref>
<ref id="B56">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Moreno-Gonzalez</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Dudley</surname> <given-names>J. W.</given-names>
</name>
<name>
<surname>Lambert</surname> <given-names>R. J.</given-names>
</name>
</person-group> (<year>1975</year>). <article-title>Design III study of linkage disequilibrium for percent oil in maize1</article-title>. <source>Crop Sci.</source> <volume>15</volume>, <fpage>840</fpage>&#x2013;<lpage>843</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.2135/cropsci1975.0011183X001500060027x</pub-id>
</citation>
</ref>
<ref id="B57">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>M&#xf4;ro</surname> <given-names>G. V.</given-names>
</name>
<name>
<surname>Santos</surname> <given-names>M. F.</given-names>
</name>
<name>
<surname>Bento</surname> <given-names>D. A. V.</given-names>
</name>
<name>
<surname>Aguiar</surname> <given-names>A. M.</given-names>
</name>
<name>
<surname>de Souza</surname> <given-names>C. L.</given-names> <suffix>Jr.</suffix>
</name>
</person-group> (<year>2012</year>). <article-title>Genetic analysis of kernel oil content in tropical maize with design III and QTL mapping</article-title>. <source>Euphytica</source> <volume>185</volume>, <fpage>419</fpage>&#x2013;<lpage>428</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1007/s10681-011-0604-x</pub-id>
</citation>
</ref>
<ref id="B58">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Nair</surname> <given-names>S. K.</given-names>
</name>
<name>
<surname>Chaikam</surname> <given-names>V.</given-names>
</name>
<name>
<surname>Gowda</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Hindu</surname> <given-names>V.</given-names>
</name>
<name>
<surname>Melchinger</surname> <given-names>A. E.</given-names>
</name>
<name>
<surname>Boddupalli</surname> <given-names>P. M.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Genetic dissection of maternal influence on in <italic>vivo</italic> haploid induction in maize</article-title>. <source>Crop J.</source> <volume>8</volume> (<issue>2</issue>), <fpage>287</fpage>&#x2013;<lpage>298</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.cj.2019.09.00</pub-id>
</citation>
</ref>
<ref id="B59">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Nair</surname> <given-names>S. K.</given-names>
</name>
<name>
<surname>Molenaar</surname> <given-names>W.</given-names>
</name>
<name>
<surname>Melchinger</surname> <given-names>A. E.</given-names>
</name>
<name>
<surname>Boddupalli</surname> <given-names>P. M.</given-names>
</name>
<name>
<surname>Martinez</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Lopez</surname> <given-names>L. A.</given-names>
</name>
<etal/>
</person-group>. (<year>2017</year>). <article-title>Dissection of a major QTL qhir1 conferring maternal haploid induction ability in maize</article-title>. <source>Theor. Appl. Genet.</source> <volume>130</volume> (<issue>6</issue>), <fpage>1113</fpage>&#x2013;<lpage>1122</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1007/s00122-017-2873-9</pub-id>
</citation>
</ref>
<ref id="B60">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Peiffer</surname> <given-names>J. A.</given-names>
</name>
<name>
<surname>Romay</surname> <given-names>M. C.</given-names>
</name>
<name>
<surname>Gore</surname> <given-names>M. A.</given-names>
</name>
<name>
<surname>Flint-Garcia</surname> <given-names>S. A.</given-names>
</name>
<name>
<surname>Zhang</surname> <given-names>Z.</given-names>
</name>
<name>
<surname>Millard</surname> <given-names>M. J.</given-names>
</name>
<etal/>
</person-group>. (<year>2014</year>). <article-title>The genetic architecture of maize height</article-title>. <source>Genetics</source> <volume>196</volume>, <fpage>1337</fpage>&#x2013;<lpage>1356</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1534/genetics.113.159152</pub-id>
</citation>
</ref>
<ref id="B61">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Peterson</surname> <given-names>R. A.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Finding optimal normalizing transformations via best Normalize</article-title>. <source>R J.</source> <volume>13</volume>, <fpage>310</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.32614/RJ-2021-041</pub-id>
</citation>
</ref>
<ref id="B62">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Prigge</surname> <given-names>V.</given-names>
</name>
<name>
<surname>S&#xe1;nchez</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Dhillon</surname> <given-names>B. S.</given-names>
</name>
<name>
<surname>Schipprack</surname> <given-names>W.</given-names>
</name>
<name>
<surname>Araus</surname> <given-names>J. L.</given-names>
</name>
<name>
<surname>B&#xe4;nziger</surname> <given-names>M.</given-names>
</name>
<etal/>
</person-group>. (<year>2011</year>). <article-title>Doubled haploids in tropical maize: I. effects of inducers and source germplasm on in <italic>vivo</italic> haploid induction rates</article-title>. <source>Crop Sci.</source> <volume>51</volume>, <fpage>1498</fpage>&#x2013;<lpage>1506</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.2135/cropsci2010.10.0568</pub-id>
</citation>
</ref>
<ref id="B63">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Prigge</surname> <given-names>V.</given-names>
</name>
<name>
<surname>Schipprack</surname> <given-names>W.</given-names>
</name>
<name>
<surname>Mahuku</surname> <given-names>G.</given-names>
</name>
<name>
<surname>Atlin</surname> <given-names>G. N.</given-names>
</name>
<name>
<surname>Melchinger</surname> <given-names>A. E.</given-names>
</name>
</person-group> (<year>2012</year>a). <article-title>Development of in <italic>vivo</italic> haploid inducers for tropical maize breeding programs</article-title>. <source>Euphytica</source> <volume>185</volume>, <fpage>481</fpage>&#x2013;<lpage>490</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1007/s10681-012-0657-5</pub-id>
</citation>
</ref>
<ref id="B64">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Prigge</surname> <given-names>V.</given-names>
</name>
<name>
<surname>Xu</surname> <given-names>X.</given-names>
</name>
<name>
<surname>Li</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Babu</surname> <given-names>R.</given-names>
</name>
<name>
<surname>Chen</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Atlin</surname> <given-names>G. N.</given-names>
</name>
<etal/>
</person-group>. (<year>2012</year>b). <article-title>New insights into the genetics of in <italic>vivo</italic> induction of maternal haploids, the backbone of doubled haploid technology in maize</article-title>. <source>Genetics</source> <volume>190</volume>, <fpage>781</fpage>&#x2013;<lpage>793</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1534/genetics.111.133066</pub-id>
</citation>
</ref>
<ref id="B65">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Randolph</surname> <given-names>L. F.</given-names>
</name>
</person-group> (<year>1940</year>). <article-title>Note on haploid frequencies</article-title>. <source>Maize Genet. Coop. Newsl.</source> <volume>14</volume>, <fpage>23</fpage>&#x2013;<lpage>24</lpage>.</citation>
</ref>
<ref id="B66">
<citation citation-type="book">
<person-group person-group-type="author">
<collab>R Core Team</collab>
</person-group>. (<year>2022</year>). <source>R: A language and environment for statistical computing</source> (<publisher-loc>Vienna, Austria</publisher-loc>: <publisher-name>R Foundation for Statistical Computing</publisher-name>). Available at: <uri xlink:href="http://www.r-project.org">http://www.r-project.org</uri>.</citation>
</ref>
<ref id="B67">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ren</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Zhang</surname> <given-names>X.</given-names>
</name>
<name>
<surname>Li</surname> <given-names>Z.</given-names>
</name>
<name>
<surname>Wu</surname> <given-names>P.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Genetic analysis of maternal haploid inducibility for in <italic>vivo</italic> haploid induction in maize</article-title>. <source>Agriculture</source> <volume>12</volume>, <fpage>845</fpage>. doi: <pub-id pub-id-type="doi">10.3390/agriculture12060845</pub-id>
</citation>
</ref>
<ref id="B68">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>R&#xf6;ber</surname> <given-names>F.</given-names>
</name>
</person-group> (<year>1999</year>). <source>Fortpflanzungsbiologische und genetische Untersuchungen mit RFLPMarkern zur in-vivo-Haploideninduktion bei Mais</source> (<publisher-loc>Stuttgart, Germany</publisher-loc>: <publisher-name>Universit&#xe4;t Hohenheim</publisher-name>). Ph.D. Thesis.</citation>
</ref>
<ref id="B69">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>R&#xf6;ber</surname> <given-names>F. K.</given-names>
</name>
<name>
<surname>Gordillo</surname> <given-names>G. A.</given-names>
</name>
<name>
<surname>Geiger</surname> <given-names>H. H.</given-names>
</name>
</person-group> (<year>2005</year>). <article-title>
<italic>In vivo</italic> haploid induction in maize - Performance of new inducers and significance of doubled haploid lines in hybrid breeding</article-title>. <source>Maydica</source> <volume>50</volume>, <fpage>275</fpage>&#x2013;<lpage>283</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1534/genetics.111.133066</pub-id>
</citation>
</ref>
<ref id="B70">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Roest</surname> <given-names>H. P.</given-names>
</name>
<name>
<surname>Baarends</surname> <given-names>W. M.</given-names>
</name>
<name>
<surname>de Wit</surname> <given-names>J.</given-names>
</name>
<name>
<surname>van Klaveren</surname> <given-names>J. W.</given-names>
</name>
<name>
<surname>Wassenaar</surname> <given-names>E.</given-names>
</name>
<name>
<surname>Hoogerbrugge</surname> <given-names>J. W.</given-names>
</name>
<etal/>
</person-group>. (<year>2004</year>). <article-title>The ubiquitin-conjugating DNA repair enzyme HR6A is a maternal factor essential for early embryonic development in mice</article-title>. <source>Mol. Cell. Biol.</source> <volume>24</volume>, <fpage>5485</fpage>&#x2013;<lpage>5495</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1128/MCB.24.12.5485-5495.2004</pub-id>
</citation>
</ref>
<ref id="B71">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ron</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Alandete Saez</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Eshed Williams</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Fletcher</surname> <given-names>J. C.</given-names>
</name>
<name>
<surname>McCormick</surname> <given-names>S.</given-names>
</name>
</person-group> (<year>2010</year>). <article-title>Proper regulation of a sperm-specific cis-nat-siRNA is essential for double fertilization in Arabidopsis</article-title>. <source>Gen. Dev.</source> <volume>24</volume>, <fpage>1010</fpage>&#x2013;<lpage>1021</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1101/gad.1882810</pub-id>
</citation>
</ref>
<ref id="B72">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Rotarenco</surname> <given-names>V.</given-names>
</name>
<name>
<surname>Georgeta</surname> <given-names>D.</given-names>
</name>
<name>
<surname>State</surname> <given-names>D.</given-names>
</name>
<name>
<surname>Fuia</surname> <given-names>S.</given-names>
</name>
</person-group> (<year>2010</year>). <article-title>New inducers of maternal haploids in maize</article-title>. <source>Maize Genet. Coop. News.</source> <volume>84</volume>, <fpage>36</fpage>&#x2013;<lpage>50</lpage>.</citation>
</ref>
<ref id="B73">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Santos</surname> <given-names>I.G.d.</given-names>
</name>
<name>
<surname>Verzegnazzi</surname> <given-names>A. L.</given-names>
</name>
<name>
<surname>Edwards</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Fletcher</surname> <given-names>J. C.</given-names>
</name>
<name>
<surname>McCormick</surname> <given-names>S.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Usefulness of temperate-adapted maize lines developed by doubled haploid and single-seed descent methods</article-title>. <source>Theor. Appl. Genet.</source> <volume>135</volume>, <fpage>1829</fpage>&#x2013;<lpage>1841</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1007/s00122-022-04075-2</pub-id>
</citation>
</ref>
<ref id="B74">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Satoh-Nagasawa</surname> <given-names>N.</given-names>
</name>
<name>
<surname>Nagasawa</surname> <given-names>N.</given-names>
</name>
<name>
<surname>Malcomber</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Sakai</surname> <given-names>H.</given-names>
</name>
<name>
<surname>Jackson</surname> <given-names>D.</given-names>
</name>
</person-group> (<year>2006</year>). <article-title>A trehalose metabolic enzyme controls inflorescence architecture in maize</article-title>. <source>Nature</source> <volume>441</volume>, <fpage>227</fpage>&#x2013;<lpage>230</lpage>. doi: <pub-id pub-id-type="doi">10.1038/nature04725</pub-id>
</citation>
</ref>
<ref id="B75">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Schwarz</surname> <given-names>G.</given-names>
</name>
</person-group> (<year>1978</year>). <article-title>Estimating the dimension of a model</article-title>. <source>Ann. Stat.</source> <volume>6</volume>, <fpage>461</fpage>&#x2013;<lpage>464</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1214/aos/1176344136</pub-id>
</citation>
</ref>
<ref id="B76">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Segura</surname> <given-names>V.</given-names>
</name>
<name>
<surname>Vilhj&#xe1;lmsson</surname> <given-names>B. J.</given-names>
</name>
<name>
<surname>Platt</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Korte</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Seren</surname> <given-names>U.</given-names>
</name>
<name>
<surname>Long</surname> <given-names>Q.</given-names>
</name>
<etal/>
</person-group>. (<year>2012</year>). <article-title>An efficient multi-locus mixed-model approach for genome-wide association studies in structured populations</article-title>. <source>Nat. Genet.</source> <volume>44</volume>, <fpage>825</fpage>&#x2013;<lpage>830</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1038/ng.2314</pub-id>
</citation>
</ref>
<ref id="B77">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Silva</surname> <given-names>A.J.da.</given-names>
</name>
</person-group> (<year>2017</year>). <source>Modelos para an&#xe1;lise de dados superdispersos de indu&#xe7;&#xe3;o de haploidia em milho [disserta&#xe7;&#xe3;o]</source> (<publisher-loc>Piracicaba</publisher-loc>: <publisher-name>Universidade de S&#xe3;o Paulo, Escola Superior de Agricultura Luiz de Queiroz</publisher-name>). doi:&#xa0;<pub-id pub-id-type="doi">10.11606/D.11.2017.tde-20072017-133448</pub-id>
</citation>
</ref>
<ref id="B78">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Stelpflug</surname> <given-names>S. C.</given-names>
</name>
<name>
<surname>Sekhon</surname> <given-names>R. S.</given-names>
</name>
<name>
<surname>Vaillancourt</surname> <given-names>B.</given-names>
</name>
<name>
<surname>Hirsch</surname> <given-names>C. N.</given-names>
</name>
<name>
<surname>Buell</surname> <given-names>C. R.</given-names>
</name>
<name>
<surname>de Leon</surname> <given-names>N.</given-names>
</name>
<etal/>
</person-group>. (<year>2016</year>). <article-title>An expanded maize gene expression atlas based on RNA sequencing and its use to explore root development</article-title>. <source>Plant Genome</source> <volume>9</volume> (<issue>1</issue>). doi:&#xa0;<pub-id pub-id-type="doi">10.3835/plantgenome2015.04.0025</pub-id>
</citation>
</ref>
<ref id="B79">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Suzumori</surname> <given-names>N.</given-names>
</name>
<name>
<surname>Burns</surname> <given-names>K. H.</given-names>
</name>
<name>
<surname>Yan</surname> <given-names>W.</given-names>
</name>
<name>
<surname>Matzuk</surname> <given-names>M. M.</given-names>
</name>
</person-group> (<year>2003</year>). <article-title>RFPL4 interacts with oocyte proteins of the ubiquitin-proteasome degradation pathway</article-title>. <source>Proc. Natl. Acad. Sci.</source> <volume>100</volume>, <fpage>550</fpage>&#x2013;<lpage>555</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1073/pnas.0234474100</pub-id>
</citation>
</ref>
<ref id="B80">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Tembo</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Asea</surname> <given-names>G.</given-names>
</name>
<name>
<surname>Gibson</surname> <given-names>P. T.</given-names>
</name>
<name>
<surname>Okori</surname> <given-names>P.</given-names>
</name>
</person-group> (<year>2014</year>). <article-title>Quantitative Trait Loci for Resistance to <italic>Stenocarpella maydis</italic> and <italic>Fusarium graminearum</italic> Cob Rots in Tropical Maize</article-title>. <source>J. Crop Improv.</source> <volume>28</volume>, <fpage>214</fpage>&#x2013;<lpage>228</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1080/15427528.2013.872750</pub-id>
</citation>
</ref>
<ref id="B81">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Trentin</surname> <given-names>H.</given-names>
</name>
<name>
<surname>Bat&#xee;ru</surname> <given-names>G.</given-names>
</name>
<name>
<surname>Frei</surname> <given-names>U. K.</given-names>
</name>
<name>
<surname>Dutta</surname> <given-names>S.</given-names>
</name>
<name>
<surname>L&#xfc;bberstedt</surname> <given-names>T.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Investigating the effect of the interaction of maize inducer and donor backgrounds on haploid induction rates</article-title>. <source>Plants (Basel Switzerland)</source> <volume>11</volume>, <fpage>1527</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.3390/plants11121527</pub-id>
</citation>
</ref>
<ref id="B82">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Trentin</surname> <given-names>H. U.</given-names>
</name>
<name>
<surname>Yavuz</surname> <given-names>R.</given-names>
</name>
<name>
<surname>Dermail</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Frei</surname> <given-names>U. K.</given-names>
</name>
<name>
<surname>Dutta</surname> <given-names>S.</given-names>
</name>
<name>
<surname>L&#xfc;bberstedt</surname> <given-names>T.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>Comparison between inbred and hybrid maize haploid inducers</article-title>. <source>Plants</source> <volume>12</volume>, <elocation-id>1095</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.3390/plants12051095</pub-id>
</citation>
</ref>
<ref id="B83">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>VanRaden</surname> <given-names>P. M.</given-names>
</name>
</person-group> (<year>2008</year>). <article-title>Efficient methods to compute genomic predictions</article-title>. <source>J. Dairy Sci.</source> <volume>91</volume>, <fpage>4414</fpage>&#x2013;<lpage>4423</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.3168/jds.2007-0980</pub-id>
</citation>
</ref>
<ref id="B84">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Verzegnazzi</surname> <given-names>A. L.</given-names>
</name>
<name>
<surname>dos Santos</surname> <given-names>I. G.</given-names>
</name>
<name>
<surname>Krause</surname> <given-names>M. D.</given-names>
</name>
<name>
<surname>Hufford</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Frei</surname> <given-names>U. K.</given-names>
</name>
<name>
<surname>Campbell</surname> <given-names>J.</given-names>
</name>
<etal/>
</person-group>. (<year>2021</year>). <article-title>Major locus for spontaneous haploid genome doubling detected by a case&#x2013;control GWAS in exotic maize germplasm</article-title>. <source>Theor. Appl. Genet.</source> <volume>134</volume>, <fpage>1423</fpage>&#x2013;<lpage>1434</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1007/s00122-021-03780-8</pub-id>
</citation>
</ref>
<ref id="B85">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Walley</surname> <given-names>J. W.</given-names>
</name>
<name>
<surname>Sartor</surname> <given-names>R. C.</given-names>
</name>
<name>
<surname>Shen</surname> <given-names>Z.</given-names>
</name>
<name>
<surname>Schmitz</surname> <given-names>R. J.</given-names>
</name>
<name>
<surname>Wu</surname> <given-names>K. J.</given-names>
</name>
<name>
<surname>Urich</surname> <given-names>M. A.</given-names>
</name>
<etal/>
</person-group>. (<year>2016</year>). <article-title>Integration of omic networks in a developmental atlas of maize</article-title>. <source>Science</source> <volume>353</volume>, <fpage>814</fpage>&#x2013;<lpage>818</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1126/science.aag1125</pub-id>
</citation>
</ref>
<ref id="B86">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname> <given-names>H.</given-names>
</name>
<name>
<surname>Lu</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Jiang</surname> <given-names>T.</given-names>
</name>
<name>
<surname>Berg</surname> <given-names>H.</given-names>
</name>
<name>
<surname>Li</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Xia</surname> <given-names>Y.</given-names>
</name>
</person-group> (<year>2013</year>). <article-title>The Arabidopsis U&#x2013;box/ARM repeat E3 ligase AtPUB4 influences growth and degeneration of tapetal cells, and its mutation leads to conditional male sterility</article-title>. <source>Plant J.</source> <volume>74</volume>, <fpage>511</fpage>&#x2013;<lpage>523</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1111/tpj.12146</pub-id>
</citation>
</ref>
<ref id="B87">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname> <given-names>D.</given-names>
</name>
<name>
<surname>Zhong</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Feng</surname> <given-names>B.</given-names>
</name>
<name>
<surname>Qi</surname> <given-names>X.</given-names>
</name>
<name>
<surname>Yan</surname> <given-names>T.</given-names>
</name>
<name>
<surname>Liu</surname> <given-names>J.</given-names>
</name>
<etal/>
</person-group>. (<year>2023</year>). <article-title>The RUBY reporter enables efficient haploid identification in maize and tomato</article-title>.. <source>Plant Biotechnol. J.</source> <volume>21</volume>, <fpage>1707</fpage>&#x2013;<lpage>1715</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1111/pbi.14071</pub-id>
</citation>
</ref>
<ref id="B88">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Xu</surname> <given-names>X.</given-names>
</name>
<name>
<surname>Li</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Dong</surname> <given-names>X.</given-names>
</name>
<name>
<surname>Jin</surname> <given-names>W.</given-names>
</name>
<name>
<surname>Melchinger</surname> <given-names>A. E.</given-names>
</name>
<name>
<surname>Chen</surname> <given-names>S.</given-names>
</name>
</person-group> (<year>2013</year>). <article-title>Gametophytic and zygotic selection leads to segregation distortion through in <italic>vivo</italic> induction of a maternal haploid in maize</article-title>. <source>J. Exp. Bot.</source> <volume>64</volume>, <fpage>1083</fpage>&#x2013;<lpage>1096</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1093/jxb/ers393</pub-id>
</citation>
</ref>
<ref id="B89">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Cui</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Zhang</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Zhang</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Li</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Kan</surname> <given-names>X.</given-names>
</name>
<etal/>
</person-group>. (<year>2016</year>). <article-title>Confirmation and fine mapping of a major QTL for aflatoxin resistance in maize using a combination of linkage and association mapping</article-title>. <source>Toxins (Basel)</source> <volume>8</volume>, <fpage>1</fpage>&#x2013;<lpage>15</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.3390/toxins8090258</pub-id>
</citation>
</ref>
<ref id="B90">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhong</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Liu</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Qi</surname> <given-names>X.</given-names>
</name>
<name>
<surname>Jiao</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>D.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>Y.</given-names>
</name>
<etal/>
</person-group>. (<year>2019</year>). <article-title>Mutation of ZmDMP enhances haploid induction in maize</article-title>. <source>Nat. Plants</source> <volume>5</volume>, <fpage>575</fpage>&#x2013;<lpage>580</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1038/s41477-019-0443-7</pub-id>
</citation>
</ref>
</ref-list>
</back>
</article>