Conserved globulin gene across eight grass genomes identify fundamental units of the loci encoding seed storage proteins

The wheat high molecular weight (HMW) glutenins are important seed storage proteins that determine bread-making quality in hexaploid wheat (Triticum aestivum). In this study, detailed comparative sequence analyses of large orthologous HMW glutenin genomic regions from eight grass species, representing a wide evolutionary history of grass genomes, reveal a number of lineage-specific sequence changes. These lineage-specific changes, which resulted in duplications, insertions, and deletions of genes, are the major forces disrupting gene colinearity among grass genomes. Our results indicate that the presence of the HMW glutenin gene in Triticeae genomes was caused by lineage-specific duplication of a globulin gene. This tandem duplication event is shared by Brachypodium and Triticeae genomes, but is absent in rice, maize, and sorghum, suggesting the duplication occurred after Brachypodium and Triticeae genomes diverged from the other grasses ~35 Ma ago. Aside from their physical location in tandem, the sequence similarity, expression pattern, and conserved cis-acting elements responsible for endosperm-specific expression further support the paralogous relationship between the HMW glutenin and globulin genes. While the duplicated copy in Brachypodium has apparently become nonfunctional, the duplicated copy in wheat has evolved to become the HMW glutenin gene by gaining a central prolamin repetitive domain.     

The limited role of differential fractionation in genome content variation and function in maize (Zea mays L.) inbred lines

Maize is a diverse paleotetraploid species with considerable presence/absence variation and copy number variation. One mechanism through which presence/absence variation can arise is differential fractionation. Fractionation refers to the loss of duplicate gene pairs from one of the maize subgenomes during diploidization. Differential fractionation refers to non‐shared gene loss events between individuals following a whole‐genome duplication event. We investigated the prevalence of presence/absence variation resulting from differential fractionation in the syntenic portion of the genome using two whole‐genome de novo assemblies of the inbred lines B73 and PH207. Between these two genomes, syntenic genes were highly conserved with less than 1% of syntenic genes being subject to differential fractionation. The few variably fractionated syntenic genes that were identified are unlikely to contribute to functional phenotypic variation, as there is a significant depletion of these genes in annotated gene sets. In further comparisons of 60 diverse inbred lines, non‐syntenic genes were six times more likely to be variable than syntenic genes, suggesting that comparisons among additional genome assemblies are not likely to result in the discovery of large‐scale presence/absence variation among syntenic genes.     

Rye Pm8 and wheat Pm3 are orthologous genes and show evolutionary conservation of resistance function against powdery mildew

The improvement of wheat through breeding has relied strongly on the use of genetic material from related wild and domesticated grass species. The 1RS chromosome arm from rye was introgressed into wheat and crossed into many wheat lines, as it improves yield and fungal disease resistance. Pm8 is a powdery mildew resistance gene on 1RS which, after widespread agricultural cultivation, is now widely overcome by adapted mildew races. Here we show by homology‐based cloning and subsequent physical and genetic mapping that Pm8 is the rye orthologue of the Pm3 allelic series of mildew resistance genes in wheat. The cloned gene was functionally validated as Pm8 by transient, single‐cell expression analysis and stable transformation. Sequence analysis revealed a complex mosaic of ancient haplotypes among Pm3‐ and Pm8‐like genes from different members of the Triticeae. These results show that the two genes have evolved independently after the divergence of the species 7.5 million years ago and kept their function in mildew resistance. During this long time span the co‐evolving pathogens have not overcome these genes, which is in strong contrast to the breakdown of Pm8 resistance since its introduction into commercial wheat 70 years ago. Sequence comparison revealed that evolutionary pressure acted on the same subdomains and sequence features of the two orthologous genes. This suggests that they recognize directly or indirectly the same pathogen effectors that have been conserved in the powdery mildews of wheat and rye.     

Nucleotide‐binding resistance gene signatures in sugar beet,insights from a new reference genome

Nucleotide‐binding (NB‐ARC), leucine‐rich‐repeat genes (NLRs) account for 60.8% of resistance (R) genes molecularly characterized from plants. NLRs exist as large gene families prone to tandem duplication and transposition, with high sequence diversity among crops and their wild relatives. This diversity can be a source of new disease resistance, but difficulty in distinguishing specific sequences from homologous gene family members hinders characterization of resistance for improving crop varieties. Current genome sequencing and assembly technologies, especially those using long‐read sequencing, are improving resolution of repeat‐rich genomic regions and clarifying locations of duplicated genes, such as NLRs. Using the conserved NB‐ARC domain as a model, 231 tentative NB‐ARC loci were identified in a highly contiguous genome assembly of sugar beet, revealing diverged and truncated NB‐ARC signatures as well as full‐length sequences. The NB‐ARC‐associated proteins contained NLR resistance gene domains, including TIR, CC and LRR, as well as other integrated domains. Phylogenetic relationships of partial and complete domains were determined, and patterns of physical clustering in the genome were evaluated. Comparison of sugar beet NB‐ARC domains to validated R‐genes from monocots and eudicots suggested extensive Beta vulgaris‐specific subfamily expansions. The NLR landscape in the rhizomania resistance conferring Rz region of Chromosome 3 was characterized, identifying 26 NLR‐like sequences spanning 20 MB. This work presents the first detailed view of NLR family composition in a member of the Caryophyllales, builds a foundation for additional disease resistance work in B. vulgaris, and demonstrates an additional nucleic‐acid‐based method for NLR prediction in non‐model plant species.     

Micro-colinearity between rice, <Emphasis Type="Italic">Brachypodium</Emphasis>, and <Emphasis Type="Italic">Triticum monococcum</Emphasis> at the wheat domestication locus <Emphasis Type="Italic">Q</Emphasis>

Brachypodium, a wild temperate grass with a small genome, was recently proposed as a new model organism for the large-genome grasses. In this study, we evaluated gene content and microcolinearity between diploid wheat (Triticum monococcum), Brachypodium sylvaticum, and rice at a local genomic region harboring the major wheat domestication gene Q. Gene density was much lower in T. monococcum (one per 41 kb) because of gene duplication and an abundance of transposable elements within intergenic regions as compared to B. sylvaticum (one per 14 kb) and rice (one per 10 kb). For the Q gene region, microcolinearity was more conserved between wheat and rice than between wheat and Brachypodium because B. sylvaticum contained two genes apparently not present within the orthologous regions of T. monococcum and rice. However, phylogenetic analysis of Q and leukotriene A-4 hydrolase-like gene orthologs, which were colinear among the three species, showed that Brachypodium is more closely related to wheat than rice, which agrees with previous studies. We conclude that Brachypodium will be a useful tool for gene discovery, comparative genomics, and the study of evolutionary relationships among the grasses but will not preclude the need to conduct large-scale genomics experiments in the Triticeae.     

A New Class of Wheat Gliadin Genes and Proteins

The utility of mining DNA sequence data to understand the structure and expression of cereal prolamin genes is demonstrated by the identification of a new class of wheat prolamins. This previously unrecognized wheat prolamin class, given the name δ-gliadins, is the most direct ortholog of barley γ3-hordeins. Phylogenetic analysis shows that the orthologous δ-gliadins and γ3-hordeins form a distinct prolamin branch that existed separate from the γ-gliadins and γ-hordeins in an ancestral Triticeae prior to the branching of wheat and barley. The expressed δ-gliadins are encoded by a single gene in each of the hexaploid wheat genomes. This single δ-gliadin/γ3-hordein ortholog may be a general feature of the Triticeae tribe since examination of ESTs from three barley cultivars also confirms a single γ3-hordein gene. Analysis of ESTs and cDNAs shows that the genes are expressed in at least five hexaploid wheat cultivars in addition to diploids Triticum monococcum and Aegilops tauschii. The latter two sequences also allow assignment of the δ-gliadin genes to the A and D genomes, respectively, with the third sequence type assumed to be from the B genome. Two wheat cultivars for which there are sufficient ESTs show different patterns of expression, i.e., with cv Chinese Spring expressing the genes from the A and B genomes, while cv Recital has ESTs from the A and D genomes. Genomic sequences of Chinese Spring show that the D genome gene is inactivated by tandem premature stop codons. A fourth δ-gliadin sequence occurs in the D genome of both Chinese Spring and Ae. tauschii, but no ESTs match this sequence and limited genomic sequences indicates a pseudogene containing frame shifts and premature stop codons. Sequencing of BACs covering a 3 Mb region from Ae. tauschii locates the δ-gliadin gene to the complex Gli-1 plus Glu-3 region on chromosome 1.     

The highly divergent Jekyll genes,required for sexual reproduction,are lineage specific for the related grass tribes Triticeae and Bromeae

Phylogenetically related groups of species contain lineage‐specific genes that exhibit no sequence similarity to any genes outside the lineage. We describe here that the Jekyll gene, required for sexual reproduction, exists in two much diverged allelic variants, Jek1 and Jek3. Despite low similarity, the Jek1 and Jek3 proteins share identical signal peptides, conserved cysteine positions and direct repeats. The Jek1/Jek3 sequences are located at the same chromosomal locus and inherited in a monogenic Mendelian fashion. Jek3 has a similar expression as Jek1 and complements the Jek1 function in Jek1‐deficient plants. Jek1 and Jek3 allelic variants were almost equally distributed in a collection of 485 wild and domesticated barley accessions. All domesticated barleys harboring the Jek1 allele belong to single haplotype J1‐H1 indicating a genetic bottleneck during domestication. Domesticated barleys harboring the Jek3 allele consisted of three haplotypes. Jekyll‐like sequences were found only in species of the closely related tribes Bromeae and Triticeae but not in other Poaceae. Non‐invasive magnetic resonance imaging revealed intrinsic grain structure in Triticeae and Bromeae, associated with the Jekyll function. The emergence of Jekyll suggests its role in the separation of the Bromeae and Triticeae lineages within the Poaceae and identifies the Jekyll genes as lineage‐specific.     

Rickettsial evolution in the light of comparative genomics

Rickettsia are best known as strictly intracellular vector‐borne bacteria that cause mild to severe diseases in humans and other animals. Recent advances in molecular tools and biological experiments have unveiled a wide diversity of Rickettsia spp. that include species with a broad host range and some species that act as endosymbiotic associates. Molecular phylogenies of Rickettsia spp. contain some ambiguities, such as the position of R. canadensis and relationships within the spotted fever group. In the modern era of genomics, with an ever‐increasing number of sequenced genomes, there is enhanced interest in the use of whole‐genome sequences to understand pathogenesis and assess evolutionary relationships among rickettsial species. Rickettsia have small genomes (1.1–1.5 Mb) as a result of reductive evolution. These genomes contain split genes, gene remnants and pseudogenes that, owing to the colinearity of some rickettsial genomes, may represent different steps of the genome degradation process. Genomics reveal extreme genome reduction and massive gene loss in highly vertebrate‐pathogenic Rickettsia compared to less virulent or endosymbiotic species. Information gleaned from rickettsial genomics challenges traditional concepts of pathogenesis that focused primarily on the acquisition of virulence factors. Another intriguing phenomenon about the reduced rickettsial genomes concerns the large fraction of non‐coding DNA and possible functionality of these “non‐coding” sequences, because of the high conservation of these regions. Despite genome streamlining, Rickettsia spp. contain gene families, selfish DNA, repeat palindromic elements and genes encoding eukaryotic‐like motifs. These features participate in sequence and functional diversity and may play a crucial role in adaptation to the host cell and pathogenesis. Genome analyses have identified a large fraction of mobile genetic elements, including plasmids, suggesting the possibility of lateral gene transfer in these intracellular bacteria. Phylogenetic analyses have identified several candidates for horizontal gene acquisition among Rickettsia spp. including tra, pat2, and genes encoding for the type IV secretion system and ATP/ADP translocase that may have been acquired from bacteria living in amoebae. Gene loss, gene duplication, DNA repeats and lateral gene transfer all have shaped rickettsial genome evolution. A comprehensive analysis of the entire genome, including genes and non‐coding DNA, will help to unlock the mysteries of rickettsial evolution and pathogenesis.     

Uncovering the evolutionary origin of blue anthocyanins in cereal grains

Functional divergence after gene duplication plays a central role in plant evolution. Among cereals, only Hordeum vulgare (barley), Triticum aestivum (wheat) and Secale cereale (rye) accumulate delphinidin‐derived (blue) anthocyanins in the aleurone layer of grains, whereas Oryza sativa (rice), Zea mays (maize) and Sorghum bicolor (sorghum) do not. The underlying genetic basis for this natural occurrence remains elusive. Here, we mapped the barley Blx1 locus involved in blue aleurone to an approximately 1.13 Mb genetic interval on chromosome 4HL, thus identifying a trigenic cluster named MbHF35 (containing HvMYB4H, HvMYC4H and HvF35H). Sequence and expression data supported the role of these genes in conferring blue‐coloured (blue aleurone) grains. Synteny analyses across monocot species showed that MbHF35 has only evolved within distinct Triticeae lineages, as a result of dispersed gene duplication. Phylogeny analyses revealed a shared evolution pattern for MbHF35 in Triticeae, suggesting that these genes have co‐evolved together. We also identified a Pooideae‐specific flavonoid 3′,5′‐hydroxylase (F3′5′H) lineage, termed here Mo_F35H2, which has a higher amino acid similarity with eudicot F3′5′Hs, demonstrating a scenario of convergent evolution. Indeed, selection tests identified 13 amino acid residues in Mo_F35H2 that underwent positive selection, possibly driven by protein thermostablility selection. Furthermore, through the interrogation of barley germplasm there is evidence that HvMYB4H and HvMYC4H have undergone human selection. Collectively, our study favours blue aleurone as a recently evolved trait resulting from environmental adaptation. Our findings provide an evolutionary explanation for the absence of blue anthocyanins in other cereals and highlight the importance of gene functional divergence for plant diversity and environmental adaptation.     

NBS-LRR sequence family is associated with leaf and stripe rust resistance on the end of homoeologous chromosome group 1S of wheat

A detailed RFLP map was constructed of the distal end of the short arm of chromosome 1D of Aegilops tauschii and wheat. At least two unrelated resistance-gene analogs (RGAs) mapped close to known leaf rust resistance genes (Lr21 and Lr40) located distal to seed storage protein genes on chromosome 1DS. One of the two RGA clones, which was previously shown to be part of a candidate gene for stripe rust resistance (Yr10) located within the homoeologous region on 1BS, identified at least three gene family members on chromosome 1DS of Ae. tauschii. One of the gene members co-segregated with the leaf rust resistance genes, Lr21 and Lr40, in Ae. tauschii and wheat segregating families. Hence, a RGA clone derived from a candidate gene for stripe rust resistance located on chromosome 1BS detected candidate genes for leaf rust resistance located in the corresponding region on 1DS of wheat. Received: 10 January 2000 / Accepted: 25 March 2000     

Sequencing chromosome 5D of Aegilops tauschii and comparison with its allopolyploid descendant bread wheat (Triticum aestivum)

Flow cytometric sorting of individual chromosomes and chromosome‐based sequencing reduces the complexity of large, repetitive Triticeae genomes. We flow‐sorted chromosome 5D of Aegilops tauschii, the D genome donor of bread wheat and sequenced it by Roche 454 GS FLX platform to approximately 2.2x coverage. Repetitive sequences represent 81.09% of the survey sequences of this chromosome, and Class I retroelements are the prominent type, with a particular abundance of LTR/Gypsy superfamily. Nonrepetitive sequences were assembled to cover 17.76% of the total chromosome regions. Up to 6188 nonrepetitive gene loci were predicted to be encoded by the 5D chromosome. The numbers and chromosomal distribution patterns of tRNA genes suggest abundance in tRNA^Lys and tRNA^Met species, while the nonrepetitive assembly reveals tRNA^Ala species as the most abundant type. A comparative analysis of the genomic sequences of bread wheat and Aegilops chromosome 5D indicates conservation of gene content. Orthologous unique genes, matching Aegilops 5D sequences, numbered 3730 in barley, 5063 in Brachypodium, 4872 in sorghum and 4209 in rice. In this study, we provide a chromosome‐specific view into the structure and organization of the 5D chromosome of Ae. tauschii, the D genome ancestor of bread wheat. This study contributes to our understanding of the chromosome‐level evolution of the wheat genome and presents a valuable resource in wheat genomics due to the recent hybridization of Ae. tauschii genome with its tetraploid ancestor.     

Genome‐wide association study identifies an NLR gene that confers partial resistance to Magnaporthe oryzae in rice

Because of the frequent breakdown of major resistance (R) genes, identification of new partial R genes against rice blast disease is an important goal of rice breeding. In this study, we used a core collection of the Rice Diversity Panel II (C‐RDP‐II), which contains 584 rice accessions and are genotyped with 700 000 single‐nucleotide polymorphism (SNP) markers. The C‐RDP‐II accessions were inoculated with three blast strains collected from different rice‐growing regions in China. Genome‐wide association study identified 27 loci associated with rice blast resistance (LABRs). Among them, 22 LABRs were not associated with any known blast R genes or QTLs. Interestingly, a nucleotide‐binding site leucine‐rich repeat (NLR) gene cluster exists in the LABR12 region on chromosome 4. One of the NLR genes is highly conserved in multiple partially resistant rice cultivars, and its expression is significantly up‐regulated at the early stages of rice blast infection. Knockout of this gene via CRISPR‐Cas9 in transgenic plants partially reduced blast resistance to four blast strains. The identification of this new non‐strain specific partial R gene, tentatively named rice blast Partial Resistance gene 1 (PiPR1), provides genetic material that will be useful for understanding the partial resistance mechanism and for breeding durably resistant cultivars against blast disease of rice.     

Structural variation and rates of genome evolution in the grass family seen through comparison of sequences of genomes greatly differing in size

Homology was searched with genes annotated in the Aegilops tauschii pseudomolecules against genes annotated in the pseudomolecules of tetraploid wild emmer wheat, Brachypodium distachyon, sorghum and rice. Similar searches were performed with genes annotated in the rice pseudomolecules. Matrices of collinear genes and rearrangements in their order were constructed. Optical BioNano genome maps were constructed and used to validate rearrangements unique to the wild emmer and Ae. tauschii genomes. Most common rearrangements were short paracentric inversions and short intrachromosomal translocations. Intrachromosomal translocations outnumbered segmental intrachromosomal duplications. The densities of paracentric inversion lengths were approximated by exponential distributions in all six genomes. Densities of collinear genes along the Ae. tauschii chromosomes were highly correlated with meiotic recombination rates but those of rearrangements were not, suggesting different causes of the erosion of gene collinearity and evolution of major chromosome rearrangements. Frequent rearrangements sharing breakpoints suggested that chromosomes have been rearranged recurrently at some sites. The distal 4 Mb of the short arms of rice chromosomes Os11 and Os12 and corresponding regions in the sorghum, B. distachyon and Triticeae genomes contain clusters of interstitial translocations including from 1 to 7 collinear genes. The rates of acquisition of major rearrangements were greater in the large wild emmer wheat and Ae. tauschii genomes than in the lineage preceding their divergence or in the B. distachyon, rice and sorghum lineages. It is suggested that synergy between large quantities of dynamic transposable elements and annual growth habit have been the primary causes of the fast evolution of the Triticeae genomes.     

The bottle gourd genome provides insights into Cucurbitaceae evolution and facilitates mapping of a Papaya ring‐spot virus resistance locus

Bottle gourd (Lagenaria siceraria) is an important vegetable crop as well as a rootstock for other cucurbit crops. In this study, we report a high‐quality 313.4‐Mb genome sequence of a bottle gourd inbred line, USVL1VR‐Ls, with a scaffold N50 of 8.7 Mb and the longest of 19.0 Mb. About 98.3% of the assembled scaffolds are anchored to the 11 pseudomolecules. Our comparative genomic analysis identifies chromosome‐level syntenic relationships between bottle gourd and other cucurbits, as well as lineage‐specific gene family expansions in bottle gourd. We reconstructed the genome of the most recent common ancestor of Cucurbitaceae, which revealed that the ancestral Cucurbitaceae karyotypes consisted of 12 protochromosomes with 18 534 protogenes. The 12 protochromosomes are largely retained in the modern melon genome, while have undergone different degrees of shuffling events in other investigated cucurbit genomes. The 11 bottle gourd chromosomes derive from the ancestral Cucurbitaceae karyotypes followed by 19 chromosomal fissions and 20 fusions. The bottle gourd genome sequence has facilitated the mapping of a dominant monogenic locus, Prs, conferring Papaya ring‐spot virus (PRSV) resistance in bottle gourd, to a 317.8‐kb region on chromosome 1. We have developed a cleaved amplified polymorphic sequence (CAPS) marker tightly linked to the Prs locus and demonstrated its potential application in marker‐assisted selection of PRSV resistance in bottle gourd. This study provides insights into the paleohistory of Cucurbitaceae genome evolution, and the high‐quality genome sequence of bottle gourd provides a useful resource for plant comparative genomics studies and cucurbit improvement.     

GABI‐DUPLO: a collection of double mutants to overcome genetic redundancy in Arabidopsis thaliana

Owing to duplication events in its progenitor, more than 90% of the genes in the Arabidopsis thaliana genome are members of multigene families. A set of 2108 gene families, each consisting of precisely two unlinked paralogous genes, was identified in the nuclear genome of A. thaliana on the basis of sequence similarity. A systematic method for the creation of double knock‐out lines for such gene pairs, designated as DUPLO lines, was established and 200 lines are now publicly available. Their initial phenotypic characterisation led to the identification of seven lines with defects that emerge only in the adult stage. A further six lines display seedling lethality and 23 lines were lethal before germination. Another 14 lines are known to show phenotypes under non‐standard conditions or at the molecular level. Knock‐out of gene pairs with very similar coding sequences or expression profiles is more likely to produce a mutant phenotype than inactivation of gene pairs with dissimilar profiles or sequences. High coding sequence similarity and highly similar expression profiles are only weakly correlated, implying that promoter and coding regions of these gene pairs display different degrees of diversification.