genetic regulation of flowering is a key element of adapting plants to cultivation in different latitudes , daylength responses and/or purposes ( such as seed / grain or biomass production ) . in the semi - arid tropical habitats to which many of the world s most important seed / grain crops
are native , flowering during  short  days of < 12.5  h of light in a 24-h period coordinates gametogenesis and reproduction with generally favorable rainfall , temperature , and solar radiation ( harper 1977 ) .
photoperiod - insensitive (  day - neutral  ) mutants of many seed and grain crops that flower based on accumulation of heat units have become widely used in temperate latitudes , producing higher seed / grain yields than short - day genotypes which generally flower late in the growing season . in crops grown for biomass rather than seed / grain ,
flowering is generally undesirable , in sugarcane , e.g. , reducing yields ( julien and soopramanien 1976 ; long 1976 ; julien et al .
1978 ; heinz 1987 ) and increasing disease susceptibility ( ricaud et al . 1980 ) .
flowering is widely acknowledged as a significant obstacle in maximizing the cellulosic yield potential of temperate biomass crops such as miscanthus , a close relative of sugarcane ( jakob et al . 2009 ) .
its small diploid genome with minimal gene duplication and substantial homozygosity make sorghum ( sorghum bicolor l. moench . ) an attractive genomic model ( paterson et al .
2009 ) for the andropogoneae , tropical grasses in which biochemical and morphological specializations (  c4  photosynthesis ) improve carbon assimilation at high temperatures and which include complex paleo- and/or neopolyploid genomes of important crops such as sugarcane , miscanthus , maize and pearl millet ( the latter being paniceae ) .
parsimony suggests the sorghum karyotype to resemble that of the andropogoneae common ancestor ( wang et al .
my of divergence from rice , sorghum has experienced only 3% differential gene loss and little rearrangement ( paterson et al .
2009 ) . with an estimated 96  my of  abstinence  from genome duplication ( wang et al .
2015 ) , functions of sorghum genes are likely to still resemble those of the common grass ancestor . in contrast , in only 15  my , genome duplication and associated  fractionations  ( gene and cis - regulator losses ) in maize ( schnable et al . 2010 ; woodhouse et al .
sorghum evolution and improvement has yielded a remarkable diversity of morphologies , with genotypes of 6  m or more in height preferred in africa as an important dual - purpose ( grain  +  straw ) crop .
in contrast , genotypes of 1  m or less are widely utilized in mechanized agriculture .
introduced into usa 150  years ago , sorghum is grown on 810 million acres and has a farm - gate value of typically $1  billion / year ( usda - nass 2013 ) .
its drought resistance provides an important  failsafe  in the agriculture of the us southern plains that often receive too little rain for other grains .
likewise , in arid countries of northeast and west africa , sorghum contributes 2639% of calories in the human diet ( fao statistical database on crops 2011 ) .
in addition to its importance for food , feed , and straw , sorghum is the # 2  us source of fuel ethanol and the availability of  sweet sorghums 
, s. bicolor genotypes with stalk sugar concentrations approaching those of sugarcane , make it suitable for commercial - scale biofuel production from each of three feedstocks , starch , sugar and cellulose ( rooney et al .
the requirement for qualitatively different flowering regimes hinders both use of the rich diversity of tropical crop germplasm in temperate latitudes , and of high - yielding improved temperate germplasm in the tropics .
for example , wild forms of sorghum and also many tropical cultivars flower only during  short  day photoperiods of 12.5  h daylight or less ( quinby and karper 1945 ) , coinciding with generally favorable rainfall , temperature , and solar radiation ( harper 1977 ) . in temperate latitudes ,
such obligate short - day genotypes typically flower near the end of the growing season , by which time the majority of solar radiation has already been missed , and often with too little time to mature seeds before freezing temperatures arrive .
much higher grain yield is generally realized in temperate regions from photoperiod - insensitive (  day - neutral  ) genotypes that flower much earlier and make better use of season - long solar radiation .
have sought to render rich diversity of native tropical crop germplasm more accessible to temperate breeding programs by crossing to a temperate donor genotype and  converting  to ( selecting for ) day - neutral flowering while backcrossing to otherwise restore the genotype of the tropical recurrent parent . to date , in sorghum alone > 700 genotypes have been converted over 50  years of effort , in most cases involving four generations of backcrossing each followed by two generations of selfing to identify recessive day - neutral flowering segregants .
previously , we studied f2 progeny of a cross between s. bicolor btx623 , a day - neutral breeding line , and s. propinquum , a tropical and short - day flowering wild relative . most of the f2 population , with at least one copy of the s. propinquum short day flowering allele , flowered in mid - september or later in the latitude where the experiment was done ( 30.6n ; college station , tx ) , consistent with the short - day flowering of the tropical s. propinquum parent .
we showed a single quantitative trait locus ( qtl ) , flravgd1 , to account for 85.7% of this phenotypic variation in flowering time ( lin et al .
1995 ) , and to correspond approximately to the locus of the triticeae photoperiodic genetic loci , ppd1 - 3 and ppd - h1 ( paterson , lin , et al .
we show by integration of positional and association genetic evidence that the short - day flowering qtl flravgd1 locates on chromosome 6 ( formerly linkage group d ) in a 10-kb interval within the physically largest centimorgan in the sorghum genome .
this interval contains a single annotated gene that we show to partly complement the day - neutral flowering phenotype .
the gene accounting for flravgd1 is sb06g012260 , a member of the  flowering locus t  ( ft ) family of transcription factors that is highly divergent from most known floral regulators in this family .
sb06g012260 is unique to panicoids , evolving as a single gene duplication shortly after the oryzoid ( rice )  panicoid ( sorghum , maize ) divergence .
sb06g012260 suppresses flowering , although it is quite distant evolutionarily from other ft family members that are floral suppressors ( pin et al .
we discuss its relationship to named sorghum flowering genes and implications of these findings for candidacy of other genes , as well as its evolution in an invasive sorghum species that has adapted to temperate latitudes .
to conduct interval mapping of flowering time in sorghum , we previously analyzed an f2 population of s. bicolor cultivar btx623s .
propinquum using 78 rflp loci spanning 935  cm with an average distance of 14  cm between markers ( lin et al .
a qtl explaining 85.7% of phenotypic variation in flowering time , flravgd1 , was placed in the 21-cm interval between dna markers psb095 and psb428a . to more finely map flravgd1
, 34 f2 plants were selected that were putatively recombinant in the 21-cm interval between dna markers psb095 and psb428a on chromosome 6 .
an additional 27 dna markers were applied to pooled dna from 50 to 150 selfed f3 progenies that were also grown in the field near college station , texas .
about 4 of the 34 f3 families , # 10 , 187 , 191 , and 211 , were excluded because the dna marker genotypes of f2 and pooled f3 tissue were not consistent ( # 211 ) , or because the flravgd1 locus genotype of their f2 parents predicted from phenotypic segregation in f3 progenies was contradicted by both flanking markers and by virtually all other markers on the chromosome ( # 10 , 187 , 191 ) .
each such inconsistency would have required a double recombination event to explain , and 3 such events among 34 progeny in a 21-cm interval are highly improbable .
a few such incongruous plants were also observed in the f2 generation , and are an important example of the need for progeny testing to reliably fine - map flowering , which can be influenced by other genetic effects , temperature , and other factors such as some diseases ( quinby 1974 ) .
based on progeny testing , flravgd1 was placed between dna markers psb1113 and cdsr084 , estimated to range from 0.3 to 1.1  cm apart in two different mapping populations studied ( the larger distance is indicated in fig .
bac clones were identified containing each of these dna markers , but lengthy efforts to  chromosome walk  in this region failed due to repeatedly encountering repetitive dna near bac ends .
the sorghum genome sequence revealed this 1.1  cm interval to include 34  mbp , thus having 60-fold less recombination than the genome - wide average ( paterson et al .
( a ) linkage mapping of flravgd1 to sorghum chromosome 6 ( lg d ; lin et al .
1995 ) ; ( b ) progeny testing in 30 f3 families that were recombinant in the interval containing flravgd1 delineated the locus to a 1.1-cm region including 400 genes ; ( c ) analysis of 90 diverse exotic sorghums and their corresponding converted derivatives delimited flravgd1 to a 4.1-mb region including 63 genes .
( d ) association genetics implicated sb06g012260 in short - day flowering conferred by flravgd1 .
( e ) the day - neutral haplotype of sb06g012260 contains three deletions in the 5 region , one removing a caat box essential to many eukaryotic promoters . 
 
( a ) linkage mapping of flravgd1 to sorghum chromosome 6 ( lg d ; lin et al .
1995 ) ; ( b ) progeny testing in 30 f3 families that were recombinant in the interval containing flravgd1 delineated the locus to a 1.1-cm region including 400 genes ; ( c ) analysis of 90 diverse exotic sorghums and their corresponding converted derivatives delimited flravgd1 to a 4.1-mb region including 63 genes .
( d ) association genetics implicated sb06g012260 in short - day flowering conferred by flravgd1 .
( e ) the day - neutral haplotype of sb06g012260 contains three deletions in the 5 region , one removing a caat box essential to many eukaryotic promoters . to better circumscribe the location of flravgd1 in the 34-mbp target region , and to confirm that the s. propinquum - derived locus also accounts for intraspecific variation within s. bicolor ( the cultigen ) , we utilized recombination accumulated during  conversion  of diverse short - day tropical sorghums to day - neutral flowering ( stephens et al . 1967 ) .
sorghum conversion involves four backcrosses each followed by two generations of selfing ( stephens et al .
1967 ) , a breeding scheme that in regions held heterozygous is expected to impart 3.5 more recombinational information to converted genotypes than f2 or recombinant inbred genotypes ( allard 1956 ) .
ninety diverse exotic sorghums and their corresponding converted derivatives ( supplementary table s1 , supplementary material online ) were genotyped at nine simple sequence repeat ( ssr ) loci distributed approximately evenly through the 34  mbp target interval inferred for the qtl based on our s. bicolors .
1995 ) , and closely coincides with flravgd1 , strongly supporting the hypothesis that the s. propinquum - derived qtl locus also accounts for intraspecific variation in short - day vs. day - neutral photoperiod response in s. bicolor .
in 384 diverse sorghums ( hash et al . 2008 ) phenotyped for  photoperiod response index 
( pri : see materials and methods ) , we identified and genotyped dna polymorphisms in annotated genes near the conversion peak .
tassel was used to perform tests of association , employing population structure covariates and a kinship matrix for the sorghum germplasm panel based on published ssrs ( hash et al .
the most significant association ( p  <  10  ) was found with sb06g012260 ( fig .
1d , with raw data in additional data file 2 ) , a gene containing a pebp domain ( interpro ipr008914 ) shared by floral regulators in arabidopsis ( ft ) ( kardailsky et al . 1999 ) and oryza ( hd3a ) ( kojima et al .
2002 ) . by comparing s. propinquum bac sequences containing sb06g012260 to the s. bicolor draft genome and resequencing polymorphic sites in the 384 diverse sorghums , we found two major sb06g012260 haplotypes .
one haplotype is widespread across tropical africa , and is shared by short - day s. propinquum .
the other haplotype was most abundant in accessions from south africa , the most temperate part of the natural range of sorghum and where day - neutral flowering is common ( fig .
the predominantly south african haplotype was characterized by four deletions , respectively , of 423  nt starting 5901  nt upstream of the sb06g012260 transcription - start site ; 4186  nt starting 4553  nt upstream ; 3  nt starting 219  nt upstream ; and 27  nt in the second intron of sb06g012260 ( fig . 1e ) .
additional indels of 2 and 7  nt ( 5451 and 5025  nt upstream ) , and three synonymous mutations in exon 1 and two in exon 2 were not analyzed in depth . 

flowering of 384 diverse sorghums were compared in the 20072008 post - rainy ( short day ) , and 2008 rainy ( long day ) seasons in peninsular india to determine  photoperiod response index  ( as described in methods ) .
hierarchical clustering of sorghum accessions was based on distances between ssr genotypes ( number of bands not shared ) . while short - day sorghums experience delayed flowering under long days ( > 12  h ) , long days accelerated flowering by 30  days for many genotypes from south africa , the most temperate part of the natural range .
east asian ( asiae ) sorghums , the largest temperate - adapted group from the northern hemisphere , also had accelerated long day flowering . 
 
flowering of 384 diverse sorghums were compared in the 20072008 post - rainy ( short day ) , and 2008 rainy ( long day ) seasons in peninsular india to determine  photoperiod response index  ( as described in methods ) .
hierarchical clustering of sorghum accessions was based on distances between ssr genotypes ( number of bands not shared ) . while short - day sorghums experience delayed flowering under long days ( > 12  h ) , long days accelerated flowering by 30  days for many genotypes from south africa , the most temperate part of the natural range .
east asian ( asiae ) sorghums , the largest temperate - adapted group from the northern hemisphere , also had accelerated long day flowering .
s1 , supplementary material online ) containing short - day s. propinquum sb06g012260 alleles subcloned from a bac and transformed into the day - neutral accession tx430 using published methods ( howe et al .
2006 ) each delayed flowering of transgenic f2 progeny in long days ( table 1 ) . widely used for sorghum transformation because of its high efficiency ( howe et al .
2006 ) , tx430 has day - neutral flowering but has an unusual haplotype , with deletion of seven amino acids in the 4th exon of sb06g012260 .
independent t0 transformants were selfed to produce t1 segregating progenies , then 1524 plants from each t1 family were evaluated in the greenhouse under ambient long day conditions ( at 33.95 n latitude ) , planting on 17 may ( 13  h , 58  min photoperiod ) and recording the number of days to flower emergence .
plants were genotyped by pcr to determine allele state for the transgene . among 13 transformation events carrying a 5 
kb construct limited to sb06g012260 and its immediate upstream elements , two conferred statistically significant delays averaging 13.1 ( p  =  0.03 ) and 24.8  days ( p  =  0.09 ) , and one line showed unexpected accelerated flowering ( 14.1  days , p  =  0.05 ) . among 10 independent events harboring a 10 
kb construct spanning the entire haplotype ( from sb06g012260 through the 4,186  nt element ) , transgenic f2 progeny of only three showed significantly altered flowering , with delays of 4.1 ( p  =  0.002 ) , 4.2 ( p  =  0.07 ) and 5.2 ( p  =  0.008 ) days . 
 table 1regression analysis of relationship of flowering time to genotype for transgenic constructs ( transgene events with p  <  0.1 are shown).constructeventwild - type*transgenic*transgene effect ( days)ppap5.21451.7564.9113.150.0310pap5.21661.0085.8224.810.0857pap5.22183.0068.8514.150.0463pap10.2348.2952.444.160.0710pap10.2647.0051.084.080.0019pap10.2749.7555.005.250.0079*least - squares means .
regression analysis of relationship of flowering time to genotype for transgenic constructs ( transgene events with p  <  0.1 are shown ) .
the predominant day - neutral sb06g012260 haplotype includes one mutation likely to alter gene function .
the 3-bp deletion located 219  nt upstream of sb06g012260 removed a predicted caat box , an invariant dna sequence in many eukaryotic promoters required for sufficient transcription ( berg et al .
the 423-bp deletion removes a non - autonomous cacta transposon ; and the 4,186-nt deletion removes an open reading frame that matches a gene on chr .
7 of day - neutral sorghum ( sb07g008600 , blastp 93.01% i d , e - value 0 ) , but with a  stop  codon in its first exon .
sb02g029725 , a family member that is more closely related to hd3a and ft than to sb06g012260 , maps near the likelihood peak of another sorghum flowering qtl , flravgb1 ( lin et al . 1995 ) .
resequencing of this gene in the 384-member diversity panel used above ( hash et al .
2008 ) revealed two abundant haplotypes ( resembling s. propinquum and btx623 , respectively ) , which showed highly significant association with pri ( p  =  1.5310  ) .
thus , at least two members of the ft gene family are implicated in the modulation of flowering in sorghum .
several other sorghum , saccharum and maize ft - containing genes occur in a clade that includes both rice hd3a ( os06g06320.1 ) ( kojima et al .
indeed , sb10g003940.1 is a colinear ortholog of hd3a , which also formed a tandem duplicate in the rice lineage ( os06g06300 : fig .
3 ) . however , the sequence of rice hd3a is more similar to arabidopsis ft and another arabidopsis homolog ( at4g20370.1 ) than to sb06g012260 ( fig .
six homologs were found in arabidopsis ( prefixed at , including ft ; kardailsky et al .
1999 ) , 19 in rice ( os , including hd3a ; kojima et al .
2002 ) , 19 in sorghum ( sb ) , 26 in maize ( grm or ac ) and 8 in sugarcane ( put-157a - saccharum_officinarum ) , identified by blast from http://www.plantgdb.org/prj/estcluster/progress.php>plantgdb , last accessed june 22 , 2016 .
a phylogenetic tree of inferred homologous protein sequences was made at http://www.phylogeny.fr/ , last accessed june 22 , 2016 , using muscle for alignment and maximum - likelihood ( phyml ) to determine the tree .
four distinct sub - families including ftl - like , ft / ftl - like , tfl1-like and mft - like proteins , are highlighted on the tree .
the impact on flowering (  promoting  or  suppressing  ) are shown for several well studied pebp genes ( reviewed in klintens et al .
the panel on the right along the tree shows multiple alignments based on the arabidopsis ft protein at amino acid position 85 at exon 2 , positions 128141 ( p - loop domain ) , positions 150152 and position 164 that were shown to have critical regulatory roles . for full alignments across the entire length of these proteins ,
six homologs were found in arabidopsis ( prefixed at , including ft ; kardailsky et al .
1999 ) , 19 in rice ( os , including hd3a ; kojima et al .
2002 ) , 19 in sorghum ( sb ) , 26 in maize ( grm or ac ) and 8 in sugarcane ( put-157a - saccharum_officinarum ) , identified by blast from http://www.plantgdb.org/prj/estcluster/progress.php>plantgdb , last accessed june 22 , 2016 .
a phylogenetic tree of inferred homologous protein sequences was made at http://www.phylogeny.fr/ , last accessed june 22 , 2016 , using muscle for alignment and maximum - likelihood ( phyml ) to determine the tree .
four distinct sub - families including ftl - like , ft / ftl - like , tfl1-like and mft - like proteins , are highlighted on the tree .
the impact on flowering (  promoting  or  suppressing  ) are shown for several well studied pebp genes ( reviewed in klintens et al .
the panel on the right along the tree shows multiple alignments based on the arabidopsis ft protein at amino acid position 85 at exon 2 , positions 128141 ( p - loop domain ) , positions 150152 and position 164 that were shown to have critical regulatory roles . for full alignments across the entire length of these proteins ,
an sb06g012260 allele found at high frequency in weedy and invasive s. halepense may have contributed to its spread .
sorghum halepense (  johnson grass  ) , a tetraploid derived naturally from wild s. bicolor and s. propinquum progenitors ( paterson , schertz , et al .
1995 ) , has spread across much of asia , africa , europe , north and south america , and australia . its establishment in usa is post - columbian , by introduction as a prospective forage and/or contaminant of sorghum seedlots ( mcwhorter 1971 ) . with inbreeding progenitors of tropical origin and therefore being putatively homozygous for short - day flowering , the original tetraploid s. halepense is expected to have had four copies of the short - day haplotype . among the few old world s. halepense accessions in the us national plant germplasm collection , we confirmed by pcr that pi209217 from south africa and pi271616 from india are homozygous for the short - day haplotype , and confirmed by growouts that these two genotypes do not flower in long days . with maximum seed production under 10.512  h daylengths consistent with short - day flowering , us ( temperate ) s. halepense nonetheless flowers and produces ample seed in photoperiods of up to 16  h ( warwick and black 1983 ) .
we genotyped the four loci diagnostic of the short - day haplotype ( 423  nt , 4,186  nt , 3  nt , and intron indels ) in 480 s. halepense plants sampled equally from populations in ga , tx ( 2 ) , ne , and nj described previously ( morrell et al .
at the two terminal loci of the haplotype , 81.6% ( 423  nt indel ) and 88.2% ( intron indel ) were homozygous for the short - day alleles , consistent with the expected genotype of their tropical progenitors .
however , at the two internal loci of the haplotype ( 4,186 and 3  nt indels ) , only 1.1% and 8.0% of plants were homozygous for the short - day - associated alleles ( fig . 4 and supplementary table s2 , supplementary material online ) .
only 39-bp upstream from the locus of the caat box deletion in day - neutral s. bicolor , 85.9% of s. halepense plants have at least one copy of a 4-nt insertion that disrupts a tc - rich repeat , a cis - acting motif enriched in promoters of photoperiod - responsive genes ( mongkolsiriwatana et al .
we have not found this mutation in any other members of the genus ( unpublished data ) , suggesting that it is rare or novel .
further , at the 4,186-nt indel , 98.9% of s. halepense also have at least one allele from the day - neutral haplotype ( fig .
thus , it appears that an sb06g012260 mutation in naturalized us s. halepense independent of those in s. bicolor has resulted in 
4genotype distributions of five us johnson grass populations near the sb06g012260 gene . among 480 plants sampled equally from ga , tx ( 2 ) , ne , and nj populations ( morrell et al .
2005 ) , 81.6% and 88.2% were homozygous for the short - day haplotype ( blue bars ) at two terminal loci ( 423  nt , 27  bp intron indels ) , but only 1.1% and 8.0% at two internal loci ( 4,186 and 3  nt indels ) .
homozygosity for day - neutral alleles ( green ) is nearly absent from the ga sample ( from the region where johnson grass is thought to have been introduced to usa ) , but exceeds 10% in the two northerly populations ( ne , nj ) where day - neutral flowering would be most advantageous , and is intermediate in the two tx populations . 
 
genotype distributions of five us johnson grass populations near the sb06g012260 gene . among 480 plants sampled equally from ga , tx ( 2 ) , ne , and nj populations ( morrell et al .
2005 ) , 81.6% and 88.2% were homozygous for the short - day haplotype ( blue bars ) at two terminal loci ( 423  nt , 27  bp intron indels ) , but only 1.1% and 8.0% at two internal loci ( 4,186 and 3  nt indels ) .
homozygosity for day - neutral alleles ( green ) is nearly absent from the ga sample ( from the region where johnson grass is thought to have been introduced to usa ) , but exceeds 10% in the two northerly populations ( ne , nj ) where day - neutral flowering would be most advantageous , and is intermediate in the two tx populations . based on several lines of evidence , sb06g012260 appears to have evolved as a single - gene duplication ( fig .
5 ) shortly after the oryzoid ( rice )  panicoid ( sorghum , maize ) divergence .
first , its divergence from its nearest sorghum homolog , sb04g008320 is estimated at 40  ma ( ks  =  0.43 , based on a widely - used evolutionary rate for cereal genes ; gaut et al .
1996 ) , suggesting that it evolved after oryzoid ( rice )  panicoid ( sorghum , maize ) divergence ( consistent with gene tree topography , fig .
3 ) . sb04g008320 has a colinear rice ortholog ( os02g13830.1 ) but sb06g012260 lacks colinear orthologs in both rice and brachypodium distachyon , although possessing them in setaria and maize ( fig . 5 ) .
this pattern of orthology is indicative of an origin of sb06g012260 in panicoids , shortly after divergence from a common ancestor shared with oryzoids ( rice ) and pooids ( brachypodium ) . 

rectangles represent predicted gene models with colors showing relative orientations ( blue : same strand , green : opposite strand ) .
matching gene pairs are displayed as connecting shades . 
  microsynteny pattern around sb06g012260 across five grasses .
rectangles represent predicted gene models with colors showing relative orientations ( blue : same strand , green : opposite strand ) . matching gene pairs
several independent lines of evidence including fine mapping , association genetics , mutant complementation , and evolutionary analysis all implicate a single gene , sb06g012260 , as the cause of the flravgd1 qtl that accounted for 85.7% of variation in flowering time of a temperatetropical sorghum population ( lin et al .
identification of this gene was complicated by its location in the physically largest single cm in the sorghum genome , containing 5% of sorghum genomic dna and 1.3% ( 400 ) of genes ( paterson et al .
2009 ) , with enrichment of retroelements and other non - genic dna . while both its panicoid - specific origin ( fig . 5 ) and
1 ) falsify the inference that ma1 corresponds to the locus of the triticeae photoperiodic genetic loci , ppd1 - 3 and ppd - h1 ( paterson , lin , et al .
10 paralog should exist at 105  mb , near the flowering qtl ( koester et al . 1993 ) . however , the nearest ft homolog to this location ( ac217051.3_fg006 : chr .
10 , 114  mb ) has a sequence so divergent from sb06g012260 that it appears non - orthologous .
the wild - type sb06g012260 is the first floral suppressor discovered in the  ftl  clade of the ft gene family although another clade member is a floral promoter , zcn8 ( meng et al .
3 ) . despite being a floral suppressor , sb06g012260 resembles all four functionally characterized members of the  tfl1  clade of the family ( hanzawa et al .
2016 ) . among the 23 transgene events that we evaluated , 5 showed statistically significant ( 0.1 or less , with 3 at 0.05 or less ) delays in flowering , whereas only one showed the acceleration that might be expected from general overexpression of a florigen .
the candidate mutation that removed a caat box upstream of the day - neutral sb06g012260 allele would be expected to reduce gene expression , also consistent with the hypothesis that wild - type sb06g012260 is a floral suppressor .
there is growing appreciation that floral regulation and other aspects of growth and development are regulated by antagonistic actions of multiple members of ft and other gene families ( lifschitz et al .
exchange of a single amino acid between arabidopsis florigen ft and floral repressor tfl1 is sufficient to exchange the functions of these opposing genes ( hanzawa et al . 2005 ) .
the finding that mutant complementation studies of sb06g012260 only partially recapitulated the short - day phenotype , underline the need for more work to understand the function and regulation of sb06g012260 .
it is noteworthy that the overwhelming majority of transgenic lines with significant deviations from tx430 had delayed flowering , not the general promotion of flowering that might be explained simply by overexpression of many ft - like genes ( mcgarry and ayre 2012 ) .
however , it is somewhat perplexing that only one of the 23 events reached the average 24.6 ( 3.5 ) day delay between the reference genetic stock 100 m ( murphy et al .
2011 ) that contains ma-1 ( a dominant allele conferring short - day flowering in tropical sorghums ; quinby and karper 1945 ) and tx430 under our conditions .
shorter flowering delays than the ma1 reference genotype 100 m relative to putatively near - isogenic sm100 ( murphy et al . 2011 ) may indicate that some distant regulatory elements were missing from the transgene constructs and/or that its native heterochromatic chromatin environment is important to its natural function .
a key element of the additional information needed is the exact timing and location of sb06g012260 gene expression . in more
than 1.7 trillion reads from ( unpublished ) sorghum transcriptomic data pooled across leaves , stems , and flowers of 40 diverse genotypes grown by several different investigators in different environments , we found only four reads from sb06g012260  two each in two different genotypes .
a detailed study of 19 sorghum pebp genes in roots , leaves ( four stages ) , stems , shoot apices and floral heads by semiquantitative rt - pcr also detected no expression of sb06g012260 in any tissue ( wolabu et al .
further , two of the three other sorghum members of the subclade including sb06g012260 , sb04g008320 ( sbft7 ) and sb03g002500 ( sbft13 ) , also show no expression in any tissue , although the third member , sb10g021790 ( sbft9 ) shows moderate expression in leaves ( wolabu et al .
finally , the nearest maize homolog of sb06g012260 , zcn21 , shows no discernible expression in a broad sampling of tissues ( danilevskaya et al .
in contrast , sb02g029725 ( sbft6 ) , which we implicate above in the flravgb1 sorghum flowering qtl , showed strong expression in leaves across all developmental stages and was also detected in the floral head and florets , but does not induce flowering in arabidopsis ( wolabu et al .
other members of its subclade ( sbft3 , 4 , 5 , and 11 ) all show expression but in different tissues and developmental states , and like sb06g012260 do not show expression patterns or genic interactions consistent with being florigens ( wolabu et al .
where does sb06g012260 fit among the named sorghum flowering loci ? based on the high frequency at which this region has long been known ( lin et al .
1995 ) to be selected for during sorghum conversion , we inferred that the qtl in this region was ma1 , the sorghum maturity ( flowering ) locus long known to have the largest effect on phenotype ( quinby and karper 1945 ) .
moreover , based on relatively coarse resolution qtl data , we previously suggested correspondence of flravgd1 to a homoeologous series of triticeae photoperiodic ( ppd1 - 3 ) loci ( paterson , lin , et al .
others have suggested prr37 to be ma1 based on positional proximity to a flowering trait mapped in genetic populations different from ours , together with its expression pattern in response to photoperiodic cues ( murphy et al .
, the transition from the heterochromatic surroundings of sb06g012260 to euchromatin , and an associated increase in recombination , coincides with a precipitous drop in levels of  conversion  at 37  mb ( fig .
prr37 is in a euchromatic genomic location ( 40.3  mb ) experiencing only 15% conversion ( fig .
curiously , study of 390 exotic - converted pairs genotyped at 46,062 markers ( albeit with 66% missing data before imputation ) suggested an introgression peak at 42  mb , different from both our peak and prr37 ( thurber et al .
the fact that this region has 60-fold less recombination than the genome - wide average ( fig .
1 ) , suggests that limited recombination rather than multiple loci under selection ( thurber et al .
2013 ) accounts for the high frequency of introgression in this region , and indeed across much of the chromosome .
evidence in support of prr37 as ma1 is further confounded by factors that were not known to its proponents ( murphy et al . 2011 ) or their reviewers .
prr37 was not tested by mutant complementation , but was a positional candidate that showed striking expression differences based on comparison of near - isogenic lines called 100 m and sm100 that differ in prr37 alleles ( murphy et al .
s2 , supplementary material online ) that the short - day flowering 100 m line is introgressed with not only a putatively short - day prr37 allele but also with the short - day sb06g012260 haplotype , based on genotyping at both the 423 and 4,186  nt indels that are on the distal side of the gene relative to prr37 .
this indicates that the chromosome segment by which 100 m and sm100 differ contains both prr37 and sb06g012260 , as well as at least 82 intervening protein - coding genes and 6.75  mb dna ( prr37 [ sb06g014570 ] has an address of chr6 : 40,280,41440,290,602 .
phenotypic differences between these lines could therefore be explained by either of these two genes , interactions between the two , or other factors .
further , in the genotype that was the source of the short - day flowering allele that we mapped ( s. propinquum ) , prr37 is non - functional with a 2-nt insertion near the 5 end causing 19 nonsense mutations ( supplementary fig .
s3 , supplementary material online ) , effectively ruling out that it could confer a dominant phenotype in crosses with s. bicolor . finally , the fine - scale mapping that we describe herein ( fig .
in addition to ma1 , a more recently discovered flowering locus , ma6 , also maps to chromosome 6 very near the location of prr37 ( brady 2006 )  indeed , one report stated that prr37 corresponded to ma6 ( mullet et al .
however , members of the same group who implicated prr37 as ma1 have more recently implicated a new positional candidate as ma6 based on apparent functional polymorphism and expression differences between long - day and short - day conditions ( murphy et al . 2014 ) 
whereas this new candidate corresponds to orthologs in several other cereals such as the rice ghd7 gene , once again evidence from mutant complementation is lacking .
the time and cost associated with sorghum transformation is clearly a hindrance to validation of candidate genes in this promising botanical model . in partial summary
, we consider the most probable scenario to be that sb06g012260 is the cause of the ma1 short - day flowering trait in sorghum , as supported by independent lines of evidence including fine mapping , association genetics , mutant complementation , and evolutionary analysis .
evidence in support of another nearby gene ( prr37 ) as ma1 is now known to be confounded with previously unknown factors , specifically the presence of sb06g012260 and at least 82 additional annotated genes on the crucial introgressed chromosome segment .
the non - functional prr37 allele in our population rules out the possibility that the two genes are functioning in concert to confer the trait . while this new evidence does not contradict the hypotheses that each gene independently contributes to the flowering phenotype , or that prr37 may be ma6 ( mullet et al . 2012 ) , recent data implicating the sorghum ortholog of rice grain number , plant height and heading date-7 ( ghd7 ) as ma6 is more compelling ( murphy et al . 2014 ) , albeit needing support from mutant complementation
the evolution of sb06g012260 from a single gene duplicate , rather than a dosage - balanced gene from whole - genome duplication , follows a trend exemplified by key members of the c4 photosynthesis pathway ( wang et al .
wgd - duplicates are a rich potential source of genetic novelty , with longer half - lives than single gene duplicates ( lynch et al .
however , single - gene duplicates , often lacking their ancestral regulatory elements and in different chromatin environments from their ancestral gene(s ) , may have greater per - gene potential for the evolution of novelty as reflected by greater divergence of expression patterns than wgd - duplicates ( wang et al .
an interesting hypothesis for further investigation is whether the survival of single - gene duplicates such as sb06g012260 may be favored by location in recombinationally - recalcitrant heterochromatin , where neutral or slightly - deleterious mutations tend to survive longer than in euchromatin ( bowers et al .
2009 ) and offer more opportunity for a mutation that confers new function to coincide with an environment that drives it to high frequency .
recalcitrance of the ma1 region of chromosome 6 to recombination may have contributed to the evolution of a  coadapted gene complex  ( mayr 1954 ) . beyond flowering ,
qtls have been associated with this region affecting kernel weight ( paterson , lin , et al . 1995 ) , tillering , rhizomatousness and regrowth ( paterson , schertz , et al .
in particular , the ma1 region also holds dw2 , the gene of largest effect on sorghum stature ( height ) ( lin et al . 1995 ) but which has long been asserted to be separable from ma1 by recombination .
quinby suggested that ma1 and dw2 were different closely - linked genes , with ca .
8% crossing over ( quinby 1974 ) , but only evaluating 47 families based on phenotype . based on our observation that late flowering can occasionally be a result of factors other than allelic status at the ma1 locus and that progeny testing is necessary to validate it , such a small study must be considered tenuous . among the
30 validated f3 families in our study , three showed different segregation patterns for flowering time and plant height . since these 30 individuals comprised all
confirmed recombinants in the region from a population of 370 individuals , this suggests a 0.5-cm linkage distance between ma1 and dw2 ( lin 1998 ) . while the strongest association of allelic variation with plant height in the diversity panel was at sb06g012260 itself ( p  =  0.007 )
, we also found a marginally significant association at sb06g007330 ( p  =  0.023 ) , a putative cation efflux family protein .
further study is needed to determine if sb06g007330 is dw2  however an intriguing hypothesis for further testing is that increased height confers a competitive advantage in light interception , and alleles conferring this trait might have become abundant more quickly if co - transmitted with alleles for optimal flowering time in the native tropical environment of sorghum . moreover ,
further dissection of the ma1 region may reveal whether qtls for other traits that have been mapped to the region ( lin et al . 1995 ; paterson , lin , et al .
1995 ; paterson , schertz , et al . , 1995 ) are pleiotropic consequences of ma1 or dw2 , or represent additional members of a coadapted gene complex .
engineering of genotypes that silence the short - day flowering trait may render obsolete the need to laboriously  convert 
tropical grasses to day - neutral flowering by twelve ( ! ) generations of breeding ( stephens et al .
1967 ) , potentially dramatically accelerating cross - utilization of temperate and tropical germplasm for sorghum , sugarcane , and many other crops .
targeted selection or engineering of strong floral repressor alleles in biomass crops may confer consistent high yields ( julien and soopramanien 1976 ; long 1976 ; julien et al .
rflp and ssr mapping used published methods , with markers drawn from published maps ( bowers et al .
phenotyping of f3 families was based on 50150 selfed progenies per family , grown in the field near college station , texas , at ambient daylengths described in detail and with flowering dates recorded as described ( lin et al .
association genetics used a 384-member worldwide sorghum diversity panel from icrisat , previously characterized with 41 ssr markers ( hash et al .
2008 ) , evaluated in 2007 under short - day conditions ( 11.812.15  h light ) and high humidity , under which short - day sorghums are expected to initiate flowering promptly .
a 2008 planting was characterized by a transition from long to short - day ( 13.111.0  h ) photoperiod and dry conditions , and short - day sorghums would be expected to delay flowering .
flowering time was the number of days required for 50% of the plants in a single row to flower ( dfl50% ) .
photoperiod response index ( pri ) was defined as the mean difference in dfl50% between the two planting seasons ( i.e. , pri  =  dfl50%2008dfl50%2007 ) .
resequencing used bigdye terminator chemistry , and sequences were manually checked and aligned for single nucleotide polymorphism ( snp ) identification with sequencher 4.1 .
the quantity and frequency of haplotypes , and linkage disequilibrium were determined by haplotyper 1.0 , and tassel 2.1 , respectively .
tassel was used to perform tests of association , employing population structure covariates and a kinship matrix for the gcp / icrisat germplasm panel based on published ssrs ( hash et al .
the bac was sequenced to confirm the integrity of sb06g012260 and identify restriction sites , then completely digested with pvuii to recover a 10.5-kb fragment that extends 5,220-bp upstream and 234-bp downstream from sb06g012260 . likewise ,
the bac clone was completely digested with restriction enzymes stui and bstz17i to recover a 5.2-kb fragment that extends 434-bp upstream and 234-bp downstream from sb06g012260 .
both fragments were cloned into transformation vector pzp211 , and the clones ( supplementary fig .
s1 , supplementary material online ) re - sequenced to confirm the integrity of the constructs .
independent t0 transformants were selfed to produce t1 segregating progenies , then 1524 plants from each t1 family were evaluated in the greenhouse under ambient long day conditions ( at 33.95 n latitude ) , recording the number of days from planting on 17 may to flower emergence .
genotyping of the four loci diagnostic of the ma1 haplotype ( 423  nt , 4186  nt , 3  nt , and intron indels ) in s. halepense used 480 plants sampled equally from ga , tx ( 2 ) , ne , and nj populations described previously ( morrell et al .
supplementary figs . s1s4 and tables s1 and s2 are available at molecular biology and evolution online ( http://www.mbe.oxfordjournals.org/ ) .
 of central importance in adapting plants of tropical origin to temperate cultivation has been selection of daylength - neutral genotypes that flower early in the temperate summer and take full advantage of its long days . a cross between tropical and temperate sorghums [ sorghum propinquum ( kunth ) hitchc.s . 
 bicolor ( l. ) moench ] , revealed a quantitative trait locus ( qtl ) , flravgd1 , accounting for 85.7% of variation in flowering time under long days . 
 fine - scale genetic mapping placed flravgd1 on chromosome 6 within the physically largest centimorgan in the genome . 
 forward genetic data from  converted  sorghums validated the qtl . 
 association genetic evidence from a diversity panel delineated the qtl to a 10-kb interval containing only one annotated gene , sb06g012260 , that was shown by reverse genetics to complement a recessive allele . 
 sb06g012260 ( sbft12 ) contains a phosphatidylethanolamine - binding ( pebp ) protein domain characteristic of members of the  ft  family of flowering genes acting as a floral suppressor . 
 sb06g012260 appears to have evolved 40  
 ma in a panicoid ancestor after divergence from oryzoid and pooid lineages . 
 a species - specific sb06g012260 mutation may have contributed to spread to temperate regions by s. halepense (  johnsongrass  ) , one of the world s most widespread invasives . 
 alternative alleles for another family member , sb02g029725 ( sbft6 ) , mapping near another flowering qtl , also showed highly significant association with photoperiod response index ( p  =  1.5310    6 ) . 
 the evolution of sb06g012260 adds to evidence that single gene duplicates play large roles in important environmental adaptations . 
 increased knowledge of sb06g012260 opens new doors to improvement of sorghum and other grain and cellulosic biomass crops . 
Introduction
Results
Discussion
Materials and Methods
Supplementary Material