Gene editing of the wheat homologs of TONNEAU1‐recruiting motif encoding gene affects grain shape and weight in wheat

Grain size and weight are important components of a suite of yield‐related traits in crops. Here, we showed that the CRISPR‐Cas9 gene editing of TaGW7, a homolog of rice OsGW7 encoding a TONNEAU1‐recruiting motif (TRM) protein, affects grain shape and weight in allohexaploid wheat. By editing the TaGW7 homoeologs in the B and D genomes, we showed that mutations in either of the two or both genomes increased the grain width and weight but reduced the grain length. The effect sizes of mutations in the TaGW7 gene homoeologs coincided with the relative levels of their expression in the B and D genomes. The effects of gene editing on grain morphology and weight traits were dosage dependent with the double‐copy mutant showing larger effect than the respective single copy mutants. The TaGW7‐centered gene co‐expression network indicated that this gene is involved in the pathways regulating cell division and organ growth, also confirmed by the cellular co‐localization of TaGW7 with α‐ and β‐tubulin proteins, the building blocks of microtubule arrays. The analyses of exome capture data in tetraploid domesticated and wild emmer, and hexaploid wheat revealed the loss of diversity around TaGW7‐associated with domestication selection, suggesting that TaGW7 is likely to play an important role in the evolution of yield component traits in wheat. Our study showed how integrating CRISPR‐Cas9 system with cross‐species comparison can help to uncover the function of a gene fixed in wheat for allelic variants targeted by domestication selection and select targets for engineering new gene variants for crop improvement.     

ADP‐glucose pyrophosphorylase genes,associated with kernel weight,underwent selection during wheat domestication and breeding

ADP‐glucose pyrophosphorylase, comprising two small subunits and two large subunits, is considered a key enzyme in the endosperm starch synthesis pathway in wheat (Triticum aestivum L.). Two genes, TaAGP‐S1‐7A and TaAGP‐L‐1B, were investigated in this study. Haplotypes of these genes were associated with thousand kernel weight (TKW) in different populations. Mean TKWs of favoured haplotypes were significantly higher than those of nonfavoured ones. Two molecular markers developed to distinguish these haplotypes could be used in molecular breeding. Frequencies of favoured haplotypes were dramatically increased in cultivars released in China after the 1940s. These favoured haplotypes were also positively selected in six major wheat production regions globally. Selection of AGP‐S1 and AGP‐L‐1B in wheat mainly occurred during and after hexaploidization. Strong additive effects of the favoured haplotypes of with other genes for starch synthesis were also detected in different populations.     

Polymorphism in the neurofibromin gene,Nf1, is associated with antagonistic selection on wing size and development time in Drosophila melanogaster

In many invertebrates, body size shows genetically based clines, with size increasing in colder climates. Large body size is typically associated with prolonged development times. We consider variation in the CNS‐specific gene neurofibromin 1 (Nf1) and its association with body size and development time. We identified two major Nf1 haplotypes in natural populations, Nf1‐insertion‐A and Nf1‐deletion‐G. These haplotypes are characterized by a 45‐base insertion/deletion (INDEL) in Nf1 intron 2 and an A/G synonymous substitution (locus L17277). Linkage disequilibrium (LD) between the INDEL and adjacent sites is high but appears to be restricted within the Nf1 gene interval. In Australia, the frequency of the Nf1‐insertion‐A haplotype increases with latitude where wing size is larger, independent of the chromosomal inversion In(3R)Payne. Unexpectedly, the Nf1‐insertion‐A haplotype is negatively associated with wing size. We found that the Nf1‐insertion‐A haplotype is enriched in females with shorter development time. This suggests that the Nf1 haplotype cline may be driven by selection for development time rather than size; females from southern (higher latitude) D. melanogaster populations maintain a rapid development time despite being relatively larger, and the higher incidence of Nf1‐insertion‐A in Southern Australia may contribute to this pattern, whereas the effects of the Nf1 haplotypes on size may be countered by other loci with antagonistic effects on size and development time. Our results point to the potential complexity involved in identifying selection on genetic variants exhibiting pleiotropic effects when studies are based on spatial patterns or association studies.     

Analysis of the functions of TaGW2 homoeologs in wheat grain weight and protein content traits

GW2 is emerging as a key genetic determinant of grain weight in cereal crops; it has three homoeologs (TaGW2‐A1, ‐B1 and ‐D1) in hexaploid common wheat (Triticum aestivum L.). Here, by analyzing the gene editing mutants that lack one (B1 or D1), two (B1 and D1) or all three (A1, B1 and D1) homoeologs of TaGW2, several insights are gained into the functions of TaGW2‐B1 and ‐D1 in common wheat grain traits. First, both TaGW2‐B1 and ‐D1 affect thousand‐grain weight (TGW) by influencing grain width and length, but the effect conferred by TaGW2‐B1 is stronger than that of TaGW2‐D1. Second, there exists functional interaction between TaGW2 homoeologs because the TGW increase shown by a double mutant (lacking B1 and D1) was substantially larger than that of their single mutants. Third, both TaGW2‐B1 and ‐D1 modulate cell number and length in the outer pericarp of developing grains, with TaGW2‐B1 being more potent. Finally, TaGW2 homoeologs also affect grain protein content as this parameter was generally increased in the mutants, especially in the lines lacking two or three homoeologs. Consistent with this finding, two wheat end‐use quality‐related parameters, flour protein content and gluten strength, were considerably elevated in the mutants. Collectively, our data shed light on functional difference between and additive interaction of TaGW2 homoeologs in the genetic control of grain weight and protein content traits in common wheat, which may accelerate further research on this important gene and its application in wheat improvement.     

The retromer protein ZmVPS29 regulates maize kernel morphology likely through an auxin‐dependent process(es)

Kernel size and morphology are two important yield‐determining traits in maize, but their molecular and genetic mechanisms are poorly characterized. Here, we identified a major QTL, qKM4.08, which explains approximately 24.20% of the kernel morphology variance in a recombinant population derived from two elite maize inbred lines, Huangzaosi (HZS, round kernel) and LV28 (slender kernel). Positional cloning and transgenic analysis revealed that qKM4.08 encodes ZmVPS29, a retromer complex component. Compared with the ZmVPS29 HZS allele, the ZmVPS29 LV28 allele showed higher expression in developing kernels. Overexpression of ZmVPS29 conferred a slender kernel morphology and increased the yield per plant in different maize genetic backgrounds. Sequence analysis revealed that ZmVPS29 has been under purifying selection during maize domestication. Association analyses identified two significant kernel morphology‐associated polymorphic sites in the ZmVPS29 promoter region that were significantly enriched in modern maize breeding lines. Further study showed that ZmVPS29 increased auxin accumulation during early kernel development by enhancing auxin biosynthesis and transport and reducing auxin degradation and thereby improved kernel development. Our results suggest that ZmVPS29 regulates kernel morphology, most likely through an auxin‐dependent process(es).     

A multiparental cross population for mapping QTL for agronomic traits in durum wheat (Triticum turgidum ssp. durum)

Multiparental cross designs for mapping quantitative trait loci (QTL) provide an efficient alternative to biparental populations because of their broader genetic basis and potentially higher mapping resolution. We describe the development and deployment of a recombinant inbred line (RIL) population in durum wheat (Triticum turgidum ssp. durum) obtained by crossing four elite cultivars. A linkage map spanning 2664 cM and including 7594 single nucleotide polymorphisms (SNPs) was produced by genotyping 338 RILs. QTL analysis was carried out by both interval mapping on founder haplotype probabilities and SNP bi‐allelic tests for heading date and maturity date, plant height and grain yield from four field experiments. Sixteen QTL were identified across environments and detection methods, including two yield QTL on chromosomes 2BL and 7AS, with the former mapped independently from the photoperiod response gene Ppd‐B1, while the latter overlapped with the vernalization locus VRN‐A3. Additionally, 21 QTL with environment‐specific effects were found. Our results indicated a prevalence of environment‐specific QTL with relatively small effect on the control of grain yield. For all traits, functionally different QTL alleles in terms of direction and size of genetic effect were distributed among parents. We showed that QTL results based on founder haplotypes closely matched functional alleles at known heading date loci. Despite the four founders, only 2.1 different functional haplotypes were estimated per QTL, on average. This durum wheat population provides a mapping resource for detailed genetic dissection of agronomic traits in an elite background typical of breeding programmes.     

Identification of Novel SNP in Promoter Sequence of TaGW2-6A Associated with Grain Weight and Other Agronomic Traits in Wheat (Triticum aestivum L.)

TaGW2 is an orthologue of rice gene OsGW2, which encodes E3 RING ubiquitin ligase and controls the grain size in rice. In wheat, three copies of TaGW2 have been identified and mapped on wheat homoeologous group 6 viz. TaGW2-6A, TaGW2-6B and TaGW2-6D. In the present study, using as many as 207 Indian wheat genotypes, we identified four SNPs including two novel SNPs (SNP-988 and SNP-494) in the promoter sequence of TaGW2-6A. All the four SNPs were G/A or A/G substitutions (transitions). Out of the four SNPs, SNP-494 was causal, since it was found associated with grain weight. The mean TGW (41.1 g) of genotypes with the allele SNP-494_A was significantly higher than mean TGW (38.6 g) of genotypes with the allele SNP-494_G. SNP-494 also regulates the expression of TaGW2-6A so that the wheat genotypes with SNP-494_G have higher expression and lower TGW and the genotypes with SNP-494_A have lower expression but higher TGW. Besides, SNP-494 was also found associated with grain length-width ratio, awn length, spike length, grain protein content, peduncle length and plant height. This suggested that gene TaGW2-6A not only controls grain size, but also controls other agronomic traits. In the promoter region, SNP-494 was present in ‘CGCG’ motif that plays an important role in Ca²⁺/calmodulin mediated regulation of genes. A user-friendly CAPS marker was also developed to identify the desirable allele of causal SNP (SNP-494) for use in marker-assisted selection for improvement of grain weight in wheat. Using four SNPs, five haplotypes were identified; of these, Hap_5 (G_A_G_A) was found to be a desirable haplotype having significantly higher grain weight (41.13g) relative to other four haplotypes (36.33-39.16 g).     

Expression of TaCYP78A3, a gene encoding cytochrome P450 CYP78A3 protein in wheat (Triticum aestivum L.), affects seed size

Several studies have described quantitative trait loci (QTL) for seed size in wheat, but the relevant genes and molecular mechanisms remain largely unknown. Here we report the functional characterization of the wheat TaCYP78A3 gene and its effect on seed size. TaCYP78A3 encoded wheat cytochrome P450 CYP78A3, and was specifically expressed in wheat reproductive organs. TaCYP78A3 activity was positively correlated with the final seed size. Its silencing caused a reduction of cell number in the seed coat, resulting in an 11% decrease in wheat seed size, whereas TaCYP78A3 over‐expression induced production of more cells in the seed coat, leading to an 11–48% increase in Arabidopsis seed size. In addition, the cell number in the final seed coat was determined by the TaCYP78A3 expression level, which affected the extent of integument cell proliferation in the developing ovule and seed. Unfortunately, TaCYP78A3 over‐expression in Arabidopsis caused a reduced seed set due to an ovule developmental defect. Moreover, TaCYP78A3 over‐expression affected embryo development by promoting embryo integument cell proliferation during seed development, which also ultimately affected the final seed size in Arabidopsis. In summary, our results indicated that TaCYP78A3 plays critical roles in influencing seed size by affecting the extent of integument cell proliferation. The present study provides direct evidence that TaCYP78A3 affects seed size in wheat, and contributes to an understanding of the cellular basis of the gene influencing seed development.     

Multiscale resistant kernel surfaces derived from inferred gene flow: An application with vernal pool breeding salamanders

The importance of assessing spatial data at multiple scales when modelling species–environment relationships has been highlighted by several empirical studies. However, no landscape genetics studies have optimized landscape resistance surfaces by evaluating relevant spatial predictors at multiple spatial scales. Here, we model multiscale/layer landscape resistance surfaces to estimate resistance to inferred gene flow for two vernal pool breeding salamander species, spotted (Ambystoma maculatum) and marbled (A. opacum) salamanders. Multiscale resistance surface models outperformed spatial layers modelled at their original spatial scale. A resistance surface with forest land cover at a 500‐m Gaussian kernel bandwidth and normalized vegetation index at a 100‐m Gaussian kernel bandwidth was the top optimized resistance surface for A. maculatum, while a resistance surface with traffic rate and topographic curvature, both at a 500‐m Gaussian kernel bandwidth, was the top optimized resistance surface for A. opacum. Species‐specific resistant kernels were fit at all vernal pools in our study area with the optimized multiscale/layer resistance surface controlling kernel spread. Vernal pools were then evaluated and scored based on surrounding upland habitat (local score) and connectivity with other vernal pools on the landscape, with resistant kernels driving vernal pool connectivity scores. As expected, vernal pools that scored highest were in areas within forested habitats and with high vernal pool densities and low species‐specific landscape resistance. Our findings highlight the success of using a novel analytical approach in a multiscale framework with applications beyond vernal pool amphibian conservation.     

Copy number variations of CBF genes at the Fr‐A2 locus are essential components of winter hardiness in wheat

Winter hardiness is important for the adaptation of wheat to the harsh winter conditions in temperate regions and is thus also an important breeding goal. Here, we employed a panel of 407 European winter wheat cultivars to dissect the genetic architecture of winter hardiness. We show that copy number variation (CNV) of CBF (C‐repeat Binding Factor) genes at the Fr‐A2 locus is the essential component for winter survival, with CBF‐A14 CNV being the most likely causal polymorphism, accounting for 24.3% of the genotypic variance. Genome‐wide association mapping identified several markers in the Fr‐A2 chromosomal region, which even after accounting for the effects of CBF‐A14 copy number explained approximately 15% of the genotypic variance. This suggests that additional, as yet undiscovered, polymorphisms are present at the Fr‐A2 locus. Furthermore, CNV of Vrn‐A1 explained an additional 3.0% of the genotypic variance. The allele frequencies of all loci associated with winter hardiness were found to show geographic patterns consistent with their role in adaptation. Collectively, our results from the candidate gene analysis, association mapping and genome‐wide prediction show that winter hardiness in wheat is a quantitative trait, but with a major contribution of the Fr‐A2 locus.     

Paternal phylogeography and genetic diversity of East Asian goats

This study was a first analysis of paternal genetic diversity for extensive Asian domestic goats using SRY gene sequences. Sequencing comparison of the SRY 3′‐untranslated region among 210 Asian goats revealed four haplotypes (Y1A, Y1B, Y2A and Y2B) derived from four variable sites including a novel substitution detected in this study. In Asian goats, the predominant haplotype was Y1A (62%) and second most common was Y2B (30%). Interestingly, the Y2B was a unique East Asian Y chromosomal variant, which differentiates eastern and western Eurasian goats. The SRY geographic distribution in Myanmar and Cambodia indicated predominant the haplotype Y1A in plains areas and a high frequency of Y2B in mountain areas. The results suggest recent genetic infiltration of modern breeds into South‐East Asian goats and an ancestral SRY Y2B haplotype in Asian native goats.     

Beta‐globin gene evolution in the ruminants: evidence for an ancient origin of sheep haplotype B

Domestic sheep (Ovis aries) can be divided into two groups with significantly different responses to hypoxic environments, determined by two allelic beta‐globin haplotypes. Haplotype A is very similar to the goat beta‐globin locus, whereas haplotype B has a deletion spanning four globin genes, including beta‐C globin, which encodes a globin with high oxygen affinity. We surveyed the beta‐globin locus using resequencing data from 70 domestic sheep from 42 worldwide breeds and three Ovis canadensis and two Ovis dalli individuals. Haplotype B has an allele frequency of 71.4% in O. aries and was homozygous (BB) in all five wild sheep. This shared ancestry indicates haplotype B is at least 2–3 million years old. Approximately 40 kb of the sequence flanking the ~37‐kb haplotype B deletion had unexpectedly low identity between haplotypes A and B. Phylogenetic analysis showed that the divergent region of sheep haplotype B is remarkably distinct from the beta‐globin loci in goat and cattle but still groups with the Ruminantia. We hypothesize that this divergent ~40‐kb region in haplotype B may be from an unknown ancestral ruminant and was maintained in the lineage to O. aries, but not other Bovidae, evolving independently of haplotype A. Alternatively, the ~40‐kb sequence in haplotype B was more recently acquired by an ancestor of sheep from an unknown non‐Bovidae ruminant, replacing part of haplotype A. Haplotype B has a lower nucleotide diversity than does haplotype A, suggesting a recent bottleneck, whereas the higher frequency of haplotype B suggests a subsequent spread through the global population of O. aries.     

Molecular diversity and landscape genomics of the crop wild relative Triticum urartu across the Fertile Crescent

Modern plant breeding can benefit from the allelic variation that exists in natural populations of crop wild relatives that evolved under natural selection in varying pedoclimatic conditions. In this study, next‐generation sequencing was used to generate 1.3 million genome‐wide single nucleotide polymorphisms (SNPs) on ex situ collections of Triticum urartu L., the wild donor of the A^u subgenome of modern wheat. A set of 75 511 high‐quality SNPs were retained to describe 298 T. urartu accessions collected throughout the Fertile Crescent. Triticum urartu showed a complex pattern of genetic diversity, with two main genetic groups distributed sequentially from west to east. The incorporation of geographical information on sampling points showed that genetic diversity was correlated to the geographical distance (R² = 0.19) separating samples from Jordan and Lebanon, from Syria and southern Turkey, and from eastern Turkey, Iran and Iraq. The wild emmer genome was used to derive the physical positions of SNPs on the seven chromosomes of the A^u subgenome, allowing us to describe a relatively slow decay of linkage disequilibrium in the collection. Outlier loci were described on the basis of the geographic distribution of the T. urartu accessions, identifying a hotspot of directional selection on chromosome 4A. Bioclimatic variation was derived from grid data and related to allelic variation using a genome‐wide association approach, identifying several marker–environment associations (MEAs). Fifty‐seven MEAs were associated with altitude and temperature measures while 358 were associated with rainfall measures. The most significant MEAs and outlier loci were used to identify genomic loci with adaptive potential (some already reported in wheat), including dormancy and frost resistance loci. We advocate the application of genomics and landscape genomics on ex situ collections of crop wild relatives to efficiently identify promising alleles and genetic materials for incorporation into modern crop breeding.     

Resistance to Asparagus virus 1 in the Wild Relative Asparagus amarus

Five asparagus cultivars, three breeding lines and the wild relative Asparagus amarus were tested for natural infection by Asparagus virus 1 (AV‐1) in experimental fields at two locations over 3 and 4 years, respectively. In the first year after re‐planting the annual crowns in the field, more than 90% of tested plants of cultivars were infected by AV‐1. In the third and fourth year, 100% of tested plants of cultivars were AV‐1 infected. In comparison, all plants of the wild relative A. amarus were completely free of AV‐1, suggesting a high level of resistance. Additionally, 1‐year‐old glasshouse‐cultivated plants of A. officinalis and A. amarus were placed in an AV‐1 provocation cabin under field conditions. Seven months later, 100% of the A. officinalis plants showed a high virus concentration in ELISA, whereas no AV‐1 was detectable in the A. amarus plants. This result was confirmed by highly sensitive AV‐1‐specific RT‐PCR. To exclude vector resistance, the feeding behaviour of green peach aphid Myzus persicae was tested over 12 h using the electrical penetration graph method. Both asparagus genotypes were accepted by the aphids as potential hosts, but the feeding time was significantly longer on A. amarus. A genetic distance analysis of the various cultivars of Asparagus officinalis and selected wild relatives of the JKI collection was carried out, resulting in a clear discrimination of cultivars and wild relatives, especially A. amarus. The potential breeding value of the putative resistance carrier is discussed.     

Identification and development of a functional marker of TaGW2 associated with grain weight in bread wheat (Triticum aestivum L.)

The OsGW2 gene is involved in rice grain development, influencing grain width and weight. Its ortholog in wheat, TaGW2, was considered as a candidate gene related to grain development. We found that TaGW2 is constitutively expressed, with three orthologs expressing simultaneously. The coding sequence (CDS) of TaGW2 is 1,275?bp encoding a protein with 424 amino acids, and has a functional domain shared with OsGW2. No divergence was detected within the CDS sequences in the same locus in ten varieties. Genome-specific primers were designed based on the sequence divergence of the promoter regions in the three orthologous genes, and TaGW2 was located in homologous group 6 chromosomes through CS nulli-tetrasomic (NT). Two SNPs were detected in the promoter region of TaGW2-6A, forming two haplotypes: Hap-6A-A (?593A and ?739G) and Hap-6A-G (?593G and ?769A). A cleaved amplified polymorphic sequence (CAPS) marker was developed based on the ?593 A-G polymorphism to distinguish the two haplotypes in TaGW2-6A. This gene was fine mapped 0.6?cM from marker cfd80.2 near the centromere in a recombinant inbred line (RIL) population. Two hundred sixty-five Chinese wheat varieties were genotyped and association analysis revealed that Hap-6A-A was significantly associated with wider grains and higher one-thousand grain weight (TGW) in two crop seasons. qRT-PCR revealed a negative relationship between TaGW2 expression level and grain width. The Hap-6A-A frequencies in Chinese varieties released at different periods showed that it had been strongly positively selected in breeding. In landraces, Hap-6A-A is mainly distributed in southern Chinese wheat regions. Association analysis also indicated that Hap-6A-A not only increased TGW by more than 3?g, but also had earlier heading and maturity. In contrast to Chinese varieties, Hap-6A-G was the predominant haplotype in European varieties; Hap-6A-A was mainly present in varieties released in the former Yugoslavia, Italy, Bulgaria, Hungary and Portugal.