Structural and functional partitioning of bread wheat chromosome 3B

In Planta Delivery of Chemical Compounds into Barley Meiocytes: EdU as Compound Example

Here, we describe a protocol for in planta delivery of chemical compounds into meiocytes of different barley genotypes not impacting plant fertility allowing to harvest seeds from treated plants. Compound uptake into meiocytes is assessed by determining 5-ethynyl-2'-deoxyuridine (EdU) incorporation. Similar to EdU, other compounds being soluble in an aqueous solution can be delivered in planta before/during meiosis to decipher their impact on meiosis and meiotic recombination.We give practical advice on how to deliver EdU as compound example (delivery via injection or needle and thread, addition of detergents or surfactants to increase compound uptake), how in planta compound delivery can be established for your plant material under specific growing conditions, how to generate and characterize barley hybrid plants, and how to conduct a meiotic cytological study of (treated) barley plants.

Genome-wide characterization and expression and co-expression analysis suggested diverse functions of WOX genes in bread wheat

WUSCHEL-related homeobox (WOX) genes belong to the homeobox superfamily, are plant-specific and play vital functions in the growth and development. Herein, we identified a total of 43 TaWOX genes in the allohexaploid (AABBDD) genome of Triticum aestivum L. These genes were distributed on the various chromosomes of each subgenome (A, B and D). The phylogenetic analysis showed the clustering of TaWOXs into three clades: ancient, intermediate and modern or WUS. The gene and protein structures including exon/intron organization, intron phases, and domain and motif distribution were found to be conserved in each phylogenetic clade. The subcellular localization was predicted as nuclear. The Ka/Ks analyses suggested the purifying selection of paralogous genes. The differential expression profiling of various TaWOXs in numerous tissue developmental stages and different layers of grains suggested their role in growth and development. Moreover, a few genes exhibited modulated expression during abiotic and biotic stress conditions, which revealed their roles in stress response. The occurrence of various cis-acting regulatory elements further confirmed their role in plant development and stress tolerance. The co-expression analyses suggested the interactions of these genes with other genes, involved in various processes including plant development, signalling and stress responses. The present study reported several characteristic features of TaWOXs genes that can be useful for further characterization in future studies.

A haplotype-led approach to increase the precision of wheat breeding

Crop productivity must increase at unprecedented rates to meet the needs of the growing worldwide population. Exploiting natural variation for the genetic improvement of crops plays a central role in increasing productivity. Although current genomic technologies can be used for high-throughput identification of genetic variation, methods for efficiently exploiting this genetic potential in a targeted, systematic manner are lacking. Here, we developed a haplotype-based approach to identify genetic diversity for crop improvement using genome assemblies from 15 bread wheat (Triticum aestivum) cultivars. We used stringent criteria to identify identical-by-state haplotypes and distinguish these from near-identical sequences (~99.95% identity). We showed that each cultivar shares ~59 % of its genome with other sequenced cultivars and we detected the presence of extended haplotype blocks containing hundreds to thousands of genes across all wheat chromosomes. We found that genic sequence alone was insufficient to fully differentiate between haplotypes, as were commonly used array-based genotyping chips due to their gene centric design. We successfully used this approach for focused discovery of novel haplotypes from a landrace collection and documented their potential for trait improvement in modern bread wheat. This study provides a framework for defining and exploiting haplotypes to increase the efficiency and precision of wheat breeding towards optimising the agronomic performance of this crucial crop.

Brinton, Uauy and colleagues utilize genomic data from the 10+ Wheat Genome Project to develop a useful tool for studying and generating new wheat cultivars. This framework uses advanced exploitation of wheat haplotypes to bring newfound precision and efficiency to wheat breeding.

Genome-Wide Characterization of OFP Family Genes in Wheat (Triticum aestivum L.) Reveals That TaOPF29a-A Promotes Drought Tolerance

OVATE family proteins (OFPs) are plant-specific transcription factors that play important roles in plant development. Although common wheat (Triticum aestivum L.) is a major staple food worldwide, OFPs have not been systematically analyzed in this important crop. Here, we performed a genome-wide survey of OFP genes in wheat and identified 100 genes belonging to 34 homoeologous groups. Arabidopsis thaliana, rice (Oryza sativa), and wheat OFP genes were divided into four subgroups based on their phylogenetic relationships. Structural analysis indicated that only four TaOFPs contain introns. We mapped the TaOFP genes onto the wheat chromosomes and determined that TaOFP17 was duplicated in this crop. A survey of cis-acting elements along the promoter regions of TaOFP genes suggested that subfunctionalization of homoeologous genes might have occurred during evolution. The TaOFPs were highly expressed in wheat, with tissue- or organ-specific expression patterns. In addition, these genes were induced by various hormone and stress treatments. For instance, TaOPF29a-A was highly expressed in roots in response to drought stress. Wheat plants overexpressing TaOPF29a-A had longer roots and higher dry weights than nontransgenic plants under drought conditions, suggesting that this gene improves drought tolerance. Our findings provide a starting point for further functional analysis of this important transcription factor family and highlight the potential of using TaOPF29a-A to genetically engineer drought-tolerant crops.

A CNL protein in wild emmer wheat confers powdery mildew resistance

Powdery mildew, a fungal disease caused by Blumeria graminis f. sp. tritici (Bgt), has a serious impact on wheat production. Loss of resistance in cultivars prompts a continuing search for new sources of resistance. Wild emmer wheat (Triticum turgidum ssp. dicoccoides, WEW), the progenitor of both modern tetraploid and hexaploid wheats, harbors many powdery mildew resistance genes. We report here the positional cloning and functional characterization of Pm41, a powdery mildew resistance gene derived from WEW, which encodes a coiled-coil, nucleotide-binding site and leucine-rich repeat protein (CNL). Mutagenesis and stable genetic transformation confirmed the function of Pm41 against Bgt infection in wheat. We demonstrated that Pm41 was present at a very low frequency (1.81%) only in southern WEW populations. It was absent in other WEW populations, domesticated emmer, durum, and common wheat, suggesting that the ancestral Pm41 was restricted to its place of origin and was not incorporated into domesticated wheat. Our findings emphasize the importance of conservation and exploitation of the primary WEW gene pool, as a valuable resource for discovery of resistance genes for improvement of modern wheat cultivars.

Genome-wide characterization and expression analysis suggested diverse functions of the mechanosensitive channel of small conductance-like (MSL) genes in cereal crops

Mechanosensitive ion channels are pore-forming transmembrane proteins that allow ions to move down their electrochemical gradient in response to mechanical stimuli. They participate in many plant developmental processes including the maintenance of plastid shape, pollen tube growth, etc. Herein, a total of 11, 10, 6, 30, 9, and 8 MSL genes were identified in Aegilops tauschii, Hordeum vulgare, Sorghum bicolor, Triticum aestivum, Triticum urartu, and Zea mays, respectively. These genes were located on various chromosomes of their respective cereal, while MSLs of T. urartu were found on scaffolds. The phylogenetic analysis, subcellular localization, and sequence homology suggested clustering of MSLs into two classes. These genes consisted of cis-regulatory elements related to growth and development, responsive to light, hormone, and stress. Differential expression of various MSL genes in tissue developmental stages and stress conditions revealed their precise role in development and stress responses. Altered expression during CaCl2 stress suggested their role in Ca2+ homeostasis and signaling. The co-expression analysis suggested their interactions with other genes involved in growth, defense responses etc. A comparative expression profiling of paralogous genes revealed either retention of function or pseudo-functionalization. The present study unfolded various characteristics of MSLs in cereals, which will facilitate their in-depth functional characterization in future studies.

MSH2 shapes the meiotic crossover landscape in relation to interhomolog polymorphism in Arabidopsis

During meiosis, DNA double-strand breaks undergo interhomolog repair to yield crossovers between homologous chromosomes. To investigate how interhomolog sequence polymorphism affects crossovers, we sequenced multiple recombinant populations of the model plant Arabidopsis thaliana. Crossovers were elevated in the diverse pericentromeric regions, showing a local preference for polymorphic regions. We provide evidence that crossover association with elevated diversity is mediated via the Class I crossover formation pathway, although very high levels of diversity suppress crossovers. Interhomolog polymorphism causes mismatches in recombining molecules, which can be detected by MutS homolog (MSH) mismatch repair protein heterodimers. Therefore, we mapped crossovers in a msh2 mutant, defective in mismatch recognition, using multiple hybrid backgrounds. Although total crossover numbers were unchanged in msh2 mutants, recombination was remodelled from the diverse pericentromeres towards the less-polymorphic sub-telomeric regions. Juxtaposition of megabase heterozygous and homozygous regions causes crossover remodelling towards the heterozygous regions in wild type Arabidopsis, but not in msh2 mutants. Immunostaining showed that MSH2 protein accumulates on meiotic chromosomes during prophase I, consistent with MSH2 regulating meiotic recombination. Our results reveal a pro-crossover role for MSH2 in regions of higher sequence diversity in A. thaliana.

TIP_finder: An HPC Software to Detect Transposable Element Insertion Polymorphisms in Large Genomic Datasets

Transposable elements (TEs) are non-static genomic units capable of moving indistinctly from one chromosomal location to another. Their insertion polymorphisms may cause beneficial mutations, such as the creation of new gene function, or deleterious in eukaryotes, e.g., different types of cancer in humans. A particular type of TE called LTR-retrotransposons comprises almost 8% of the human genome. Among LTR retrotransposons, human endogenous retroviruses (HERVs) bear structural and functional similarities to retroviruses. Several tools allow the detection of transposon insertion polymorphisms (TIPs) but fail to efficiently analyze large genomes or large datasets. Here, we developed a computational tool, named TIP_finder, able to detect mobile element insertions in very large genomes, through high-performance computing (HPC) and parallel programming, using the inference of discordant read pair analysis. TIP_finder inputs are (i) short pair reads such as those obtained by Illumina, (ii) a chromosome-level reference genome sequence, and (iii) a database of consensus TE sequences. The HPC strategy we propose adds scalability and provides a useful tool to analyze huge genomic datasets in a decent running time. TIP_finder accelerates the detection of transposon insertion polymorphisms (TIPs) by up to 55 times in breast cancer datasets and 46 times in cancer-free datasets compared to the fastest available algorithms. TIP_finder applies a validated strategy to find TIPs, accelerates the process through HPC, and addresses the issues of runtime for large-scale analyses in the post-genomic era. TIP_finder version 1.0 is available at https://github.com/simonorozcoarias/TIP_finder.

SPO11.2 is essential for programmed double-strand break formation during meiosis in bread wheat (Triticum aestivum L.)

Meiotic recombination is initiated by formation of DNA double-strand breaks (DSBs). This involves a protein complex that includes in plants the two similar proteins, SPO11-1 and SPO11-2. We analysed the sequences of SPO11-2 in hexaploid bread wheat (Triticum aestivum), as well as in its diploid and tetraploid progenitors. We investigated its role during meiosis using single, double and triple mutants. The three homoeologous SPO11-2 copies of hexaploid wheat exhibit high nucleotide and amino acid similarities with those of the diploids, tetraploids and Arabidopsis. Interestingly, however, two nucleotides deleted in exon-2 of the A copy lead to a premature stop codon and suggest that it encodes a non-functional protein. Remarkably, the mutation was absent from the diploid A-relative Triticum urartu, but present in the tetraploid Triticum dicoccoides and in different wheat cultivars indicating that the mutation occurred after the first polyploidy event and has since been conserved. We further show that triple mutants with all three copies (A, B, D) inactivated are sterile. Cytological analyses of these mutants show synapsis defects, accompanied by severe reductions in bivalent formation and numbers of DMC1 foci, thus confirming the essential role of TaSPO11-2 in meiotic recombination in wheat. In accordance with its 2-nucleotide deletion in exon-2, double mutants for which only the A copy remained are also sterile. Notwithstanding, some DMC1 foci remain visible in this mutant, suggesting a residual activity of the A copy, albeit not sufficient to restore fertility.

Structural Variations Affecting Genes and Transposable Elements of Chromosome 3B in Wheats

Structural variations (SVs) such as copy number and presence–absence variations are polymorphisms that are known to impact genome composition at the species level and are associated with phenotypic variations. In the absence of a reference genome sequence, their study has long been hampered in wheat. The recent production of new wheat genomic resources has led to a paradigm shift, making possible to investigate the extent of SVs among cultivated and wild accessions. We assessed SVs affecting genes and transposable elements (TEs) in a Triticeae diversity panel of 45 accessions from seven tetraploid and hexaploid species using high-coverage shotgun sequencing of sorted chromosome 3B DNA and dedicated bioinformatics approaches. We showed that 23% of the genes are variable within this panel, and we also identified 330 genes absent from the reference accession Chinese Spring. In addition, 60% of the TE-derived reference markers were absent in at least one accession, revealing a high level of intraspecific and interspecific variability affecting the TE space. Chromosome extremities are the regions where we observed most of the variability, confirming previous hypotheses made when comparing wheat with the other grasses. This study provides deeper insights into the genomic variability affecting the complex Triticeae genomes at the intraspecific and interspecific levels and suggests a phylogeny with independent hybridization events leading to different hexaploid species.

Comprehensive analysis of polygalacturonase gene family highlights candidate genes related to pollen development and male fertility in wheat (Triticum aestivum L.)

Four polygalacturonase gene family members were highlighted that contribute to elucidate the roles of polygalacturonase during the fertility conversion process in male-sterile wheat. Polygalacturonase (PG) belongs to a large family of hydrolases with important functions in cell separation during plant growth and development via the degradation of pectin. Specific expressed PGs in anthers may be significant for male sterility research and hybrid wheat breeding, but they have not been characterized in wheat (Triticum aestivum L.). In this study, we systematically studied the PG gene family using the latest published wheat reference genomic information. In total, 113 wheat PG genes were identified, which could be classified into six categories A-F according to their structure characteristics and phylogenetic comparisons with Arabidopsis and rice. Polyploidy and segmental duplications in wheat were proved to be mainly responsible for the expansion of the wheat PG gene family. RNA-seq showed that TaPGs have specific temporal and spatial expression characteristics, in which 12 TaPGs with spike-specific expression patterns were detected by qRT-PCR in different fertility anthers of KTM3315A, a thermo-sensitive cytoplasmic male-sterile wheat. Four of them specific upregulated (TaPG09, TaPG95, and TaPG93) or downregulated (TaPG87) at trinucleate stage of fertile anthers, and further aligning with the homologous in Arabidopsis revealed that they may undertake functions such as anther dehiscence, separation of pollen, pollen development, and pollen tube elongation, thereby inducing male fertility conversion in KTM3315A. These findings facilitate function investigations of the wheat PG gene family and provide new insights into the fertility conversion mechanism in male-sterile wheat.

Phylogeny and gene expression of the complete NITRATE TRANSPORTER 1/PEPTIDE TRANSPORTER FAMILY in Triticum aestivum

NPF genes encode membrane transporters involved in the transport of a large variety of substrates including nitrate and peptides. The NPF gene family has been described for many plants, but the whole NPF gene family for wheat has not been completely identified. The release of the wheat reference genome has enabled the identification of the entire wheat NPF gene family. A systematic analysis of the whole wheat NPF gene family was performed, including responses of specific gene expression to development and nitrogen supply. A total of 331 NPF genes (113 homoeologous groups) have been identified in wheat. The chromosomal location of the NPF genes is unevenly distributed, with predominant occurrence in the long arms of the chromosomes. The phylogenetic analysis indicated that wheat NPF genes are closely clustered with Arabidopsis, Brachypodium, and rice orthologues, and subdivided into eight subfamilies. The expression profiles of wheat NPF genes were examined using RNA-seq data, and a subset of 44 NPF genes (homoeologous groups) with contrasting expression responses to nitrogen and/or development in different tissues were identified. The systematic identification of gene composition, chromosomal locations, evolutionary relationships, and expression profiles contributes to a better understanding of the roles of the wheat NPF genes and lays the foundation for further functional analysis in wheat.

Illegitimate Recombination Between Homeologous Genes in Wheat Genome

Polyploidies produce a large number of duplicated regions and genes in genomes, which have a long-term impact and stimulate genetic innovation. The high similarity between homeologous chromosomes, forming different subgenomes, or homologous regions after genome repatterning, may permit illegitimate DNA recombination. Here, based on gene colinearity, we aligned the (sub)genomes of common wheat (Triticum aestivum, AABBDD genotype) and its relatives, including Triticum urartu (AA), Aegilops tauschii (DD), and T. turgidum ssp. dicoccoides (AABB) to detect the homeologous (paralogous or orthologous) colinear genes within and between (sub)genomes. Besides, we inferred more ancient paralogous regions produced by a much ancient grass-common tetraploidization. By comparing the sequence similarity between paralogous and orthologous genes, we assumed abnormality in the topology of constructed gene trees, which could be explained by gene conversion as a result of illegitimate recombination. We found large numbers of inferred converted genes (>2,000 gene pairs) suggested long-lasting genome instability of the hexaploid plant, and preferential donor roles by DD genes. Though illegitimate recombination was much restricted, duplicated genes produced by an ancient whole-genome duplication, which occurred millions of years ago, also showed evidence of likely gene conversion. As to biological function, we found that ~40% catalytic genes in colinearity, including those involved in starch biosynthesis, were likely affected by gene conversion. The present study will contribute to understanding the functional and structural innovation of the common wheat genome.

Differential chromatin accessibility landscape reveals structural and functional features of the allopolyploid wheat chromosomes

Background

Our understanding of how the complexity of the wheat genome influences the distribution of chromatin states along the homoeologous chromosomes is limited. Using a differential nuclease sensitivity assay, we investigate the chromatin states of the coding and repetitive regions of the allopolyploid wheat genome.

Results

Although open chromatin is found to be significantly enriched around genes, the majority of MNase-sensitive regions are located within transposable elements (TEs). Chromatin of the smaller D genome is more accessible than that of the larger A and B genomes. Chromatin states of different TEs vary among families and are influenced by the TEs’ chromosomal position and proximity to genes. While the chromatin accessibility of genes is influenced by proximity to TEs, and not by their position on the chromosomes, we observe a negative chromatin accessibility gradient along the telomere-centromere axis in the intergenic regions, positively correlated with the distance between genes. Both gene expression levels and homoeologous gene expression bias are correlated with chromatin accessibility in promoter regions. The differential nuclease sensitivity assay accurately predicts previously detected centromere locations. SNPs located within more accessible chromatin explain a higher proportion of genetic variance for a number of agronomic traits than SNPs located within more closed chromatin.

Conclusions

Chromatin states in the wheat genome are shaped by the interplay of repetitive and gene-encoding regions that are predictive of the functional and structural organization of chromosomes, providing a powerful framework for detecting genomic features involved in gene regulation and prioritizing genomic variation to explain phenotypes.

Molecular characterization of ascorbate peroxidase (APX) and APX-related (APX-R) genes in Triticum aestivum L

Ascorbate peroxidases (APXs) are heme-dependent H₂O₂ scavenging enzymes involved in myriad biological processes. Herein, a total of 21 TaAPX and six TaAPX-R genes were identified from the A, B and D sub-genomes of Triticum aestivum L. The occurrence of three paralogous gene pairs with unequal evolutionary rate suggested functional divergence. The phylogenetic analysis formed four distinct clades having conserved gene and protein architecture, and sub-cellular localization. The tertiary structure analysis revealed the presence of helices and coils and residues involved in ligand binding. Transcriptional profiling of each TaAPX and TaAPX-R gene suggested their specific role during development and stress response. Modulated transcript expression and APX enzyme activity during various stress conditions indicated their role in stress response. Interaction analyses suggested their association with other genes, miRNAs and various legends. The present study reported numerous features of these genes, and may provide a platform for their detailed functional characterization in future studies.

AXR1 affects DNA methylation independently of its role in regulating meiotic crossover localization

Meiotic crossovers (COs) are important for reshuffling genetic information between homologous chromosomes and they are essential for their correct segregation. COs are unevenly distributed along chromosomes and the underlying mechanisms controlling CO localization are not well understood. We previously showed that meiotic COs are mis-localized in the absence of AXR1, an enzyme involved in the neddylation/rubylation protein modification pathway in Arabidopsis thaliana. Here, we report that in axr1^-/-, male meiocytes show a strong defect in chromosome pairing whereas the formation of the telomere bouquet is not affected. COs are also redistributed towards subtelomeric chromosomal ends where they frequently form clusters, in contrast to large central regions depleted in recombination. The CO suppressed regions correlate with DNA hypermethylation of transposable elements (TEs) in the CHH context in axr1^-/- meiocytes. Through examining somatic methylomes, we found axr1^-/- affects DNA methylation in a plant, causing hypermethylation in all sequence contexts (CG, CHG and CHH) in TEs. Impairment of the main pathways involved in DNA methylation is epistatic over axr1^-/- for DNA methylation in somatic cells but does not restore regular chromosome segregation during meiosis. Collectively, our findings reveal that the neddylation pathway not only regulates hormonal perception and CO distribution but is also, directly or indirectly, a major limiting pathway of TE DNA methylation in somatic cells.

In sexually reproducing organisms, each parent transmits one and only one copy of each chromosome to their progeny via their packaging in haploid gametes. To ensure the proper transmission of the chromosomes, pairs of homologous chromosomes must associate and exchange genetic information (also called reciprocal recombination) during a special division called meiosis that lead to the formation of the gametes. The recombination process is highly controlled in terms of number and localization of the events along the chromosomes. Disruption of this control may cause an inappropriate transmission of the chromosomes in the gametes leading to abnormal chromosome numbers in the offspring which is usually deleterious. In the plant Arabidopis thaliana, we show that when the pathway modifying proteins through ubiquitination/neddylation is impaired, the number of reciprocal recombination events is maintained but they are delocalized toward the ends of the chromosomes and some chromosomes do not exchange material. We also detected changes of patterns for DNA methylation, an epigenetic modification localised on DNA cytosines. Furthermore, we demonstrate that the methylation of cytosines is not causal to the localization change of meiotic recombination events.

The Battle to Sequence the Bread Wheat Genome: A Tale of the Three Kingdoms

In the year 2018, the world witnessed the finale of the race to sequence the genome of the world's most widely grown crop, the common wheat. Wheat has been known to bear a notoriously large and complicated genome of a polyploidy nature. A decade competition to sequence the wheat genome initiated with a single consortium of multiple countries, taking a conventional strategy similar to that for sequencing Arabidopsis and rice, became ferocious over time as both sequencing technologies and genome assembling methodologies advanced. At different stages, multiple versions of genome sequences of the same variety (e.g., Chinese Spring) were produced by several groups with their special strategies. Finally, 16 years after the rice genome was finished and 9 years after that of maize, the wheat research community now possesses its own reference genome. Armed with these genomics tools, wheat will reestablish itself as a model for polyploid plants in studying the mechanisms of polyploidy evolution, domestication, genetic and epigenetic regulation of homoeolog expression, as well as defining its genetic diversity and breeding on the genome level. The enhanced resolution of the wheat genome should also help accelerate development of wheat cultivars that are more tolerant to biotic and/or abiotic stresses with better quality and higher yield.

ASY1 acts as a dosage-dependent antagonist of telomere-led recombination and mediates crossover interference in Arabidopsis

During meiosis, interhomolog recombination produces crossovers and noncrossovers to create genetic diversity. Meiotic recombination frequency varies at multiple scales, with high subtelomeric recombination and suppressed centromeric recombination typical in many eukaryotes. During recombination, sister chromatids are tethered as loops to a polymerized chromosome axis, which, in plants, includes the ASY1 HORMA domain protein and REC8-cohesin complexes. Using chromatin immunoprecipitation, we show an ascending telomere-to-centromere gradient of ASY1 enrichment, which correlates strongly with REC8-cohesin ChIP-seq data. We mapped crossovers genome-wide in the absence of ASY1 and observe that telomere-led recombination becomes dominant. Surprisingly, asy1/+ heterozygotes also remodel crossovers toward subtelomeric regions at the expense of the pericentromeres. Telomeric recombination increases in asy1/+ occur in distal regions where ASY1 and REC8 ChIP enrichment are lowest in wild type. In wild type, the majority of crossovers show interference, meaning that they are more widely spaced along the chromosomes than expected by chance. To measure interference, we analyzed double crossover distances, MLH1 foci, and fluorescent pollen tetrads. Interestingly, while crossover interference is normal in asy1/+, it is undetectable in asy1 mutants, indicating that ASY1 is required to mediate crossover interference. Together, this is consistent with ASY1 antagonizing telomere-led recombination and promoting spaced crossover formation along the chromosomes via interference. These findings provide insight into the role of the meiotic axis in patterning recombination frequency within plant genomes.

Measuring Performance Metrics of Machine Learning Algorithms for Detecting and Classifying Transposable Elements

Because of the promising results obtained by machine learning (ML) approaches in several fields, every day is more common, the utilization of ML to solve problems in bioinformatics. In genomics, a current issue is to detect and classify transposable elements (TEs) because of the tedious tasks involved in bioinformatics methods. Thus, ML was recently evaluated for TE datasets, demonstrating better results than bioinformatics applications. A crucial step for ML approaches is the selection of metrics that measure the realistic performance of algorithms. Each metric has specific characteristics and measures properties that may be different from the predicted results. Although the most commonly used way to compare measures is by using empirical analysis, a non-result-based methodology has been proposed, called measure invariance properties. These properties are calculated on the basis of whether a given measure changes its value under certain modifications in the confusion matrix, giving comparative parameters independent of the datasets. Measure invariance properties make metrics more or less informative, particularly on unbalanced, monomodal, or multimodal negative class datasets and for real or simulated datasets. Although several studies applied ML to detect and classify TEs, there are no works evaluating performance metrics in TE tasks. Here, we analyzed 26 different metrics utilized in binary, multiclass, and hierarchical classifications, through bibliographic sources, and their invariance properties. Then, we corroborated our findings utilizing freely available TE datasets and commonly used ML algorithms. Based on our analysis, the most suitable metrics for TE tasks must be stable, even using highly unbalanced datasets, multimodal negative class, and training datasets with errors or outliers. Based on these parameters, we conclude that the F1-score and the area under the precision-recall curve are the most informative metrics since they are calculated based on other metrics, providing insight into the development of an ML application.

WFhb1-1 plays an important role in resistance against Fusarium head blight in wheat

Fusarium head blight (FHB) is a severe disease of wheat (Triticum aestivum L.). Qfhb1 is the most important quantitative trait locus (QTL) for FHB resistance. We previously identified wheat gene WFhb1-1 (aka WFhb1-c1) as a candidate for FHB resistance gene. Here we report that WFhb1-1 has been cloned. The gene (GenBank # KU304333.1) consists of a single exon, encoding a putative membrane protein of 127 amino acids. WFhb1-1 protein produced in Pichia pastoris inhibits growth of both F. graminearum and P. pastoris in culture. Western Blotting with anti- WFhb1-1 antibody revealed a significant decrease (p < 0.01) in WFhb1-1 accumulation, 12 hours post Fusarium inoculation in non-Qfhb1-carrier wheat but not in Qfhb1-carrier wheat. Overexpressing WFhb1-1 in non-Qfhb1-carrier wheat led to a significant decrease (p < 0.01) in Fusarium-damaged rachis rate, Fusarium-diseased kernel rate and DON content in harvested kernels, while silencing WFhb1-1 in Qfhb1-carrier wheat resulted in a significant increase (p < 0.01) in FHB severity. Therefore, WFhb1-1 is an important FHB resistance gene with a potential antifungal function and probably a key functional component of Qfhb1 in wheat. A model regarding how WFhb1-1 functions in FHB resistance/susceptibility is hypothesized and discussed.

Genome-Scale Analysis of Homologous Genes among Subgenomes of Bread Wheat (Triticum aestivum L.)

Determining the distribution and correspondence of genome-scale homologous genes in wheat are effective ways to uncover chromosome rearrangement that has occurred during crop evolution and domestication, which can contribute to improvements in crop breeding. High-resolution and comprehensive analysis of the wheat genome by the International Wheat Genome Sequencing Consortium (IWGSC) revealed a total of 88,733 high-confidence homologous genes of four major types (1:1:1, 1:1:0, 0:1:1 and 1:0:1) among the A, B and D subgenomes of wheat. This data was used to compare homologous gene densities among chromosomes, clarify their distribution and correspondence relationship, and compare their functional enrichment. The average density of 1:1:1 homologous genes was about 10 times more than the density of the other three types of homologous genes, although the homologous gene densities of the various chromosomes were similar within each homologous type. Three regions of exceptional density were detected in 1:1:1 homologous genes, the isolate peak on the tail of chromosome 4A, and the desert regions at the start of chromosome 7A and 7D. The correspondence between homologous genes of the wheat subgenomes demonstrated translocation between the tail segments of chromosome 4A and 5A, and the inversion of the segment of original 5A and 7B into the tail of 4A. The homologous genes on the inserting segments of 5A and 7B to 4A were highly enriched in nitrogen, primary metabolite and small molecular metabolism processes, compared with genes on other regions of the original 4A chromosome. This study provides a refined genome-scale reference of homologous genes for wheat molecular research and breeding, which will help to broaden the application of the wheat genome and can be used as a template for research on other polyploid plants.

Partitioning and physical mapping of wheat chromosome 3B and its homoeologue 3E in Thinopyrum elongatum by inducing homoeologous recombination

We performed homoeologous recombination-based partitioning and physical mapping of wheat chromosome 3B and Th. elongatum chromosome 3E, providing a unique physical framework of this homoeologous pair for genome studies. The wheat (Triticum aestivum, 2n = 6x = 42, AABBDD) and Thinopyrum elongatum (2n = 2x = 14, EE) genomes can be differentiated from each other by fluorescent genomic in situ hybridization (FGISH) as well as molecular markers. This has facilitated homoeologous recombination-based partitioning and engineering of their genomes for physical mapping and alien introgression. Here, we constructed a special wheat genotype, which was double monosomic for wheat chromosome 3B and Th. elongatum chromosome 3E and homozygous for the ph1b mutant, to induce 3B-3E homoeologous recombination. Totally, 81 3B-3E recombinants were recovered and detected in the primary, secondary, and tertiary homoeologous recombination cycles by FGISH. Comparing to the primary recombination, the secondary and tertiary recombination shifted toward the proximal regions due to the increase in homology between the pairing partners. The 3B-3E recombinants were genotyped by high-throughput wheat 90-K single nucleotide polymorphism (SNP) arrays and their recombination breakpoints physically mapped based on the FGISH patterns and SNP results. The 3B-3E recombination physically partitioned chromosome 3B into 38 bins, and 429 SNPs were assigned to the distinct bins. Integrative analysis of FGISH and SNP results led to the construction of a composite bin map for chromosome 3B. Additionally, we developed 22 SNP-derived semi-thermal asymmetric reverse PCR markers specific for chromosome 3E and constructed a comparative map of homoeologous chromosomes 3E, 3B, 3A, and 3D. In summary, this work provides a unique physical framework for further studies of the 3B-3E homoeologous pair and diversifies the wheat genome for wheat improvement.

The improved assembly of 7DL chromosome provides insight into the structure and evolution of bread wheat

Wheat is one of the most important staple crops worldwide and also an excellent model species for crop evolution and polyploidization studies. The breakthrough of sequencing the bread wheat genome and progenitor genomes lays the foundation to decipher the complexity of wheat origin and evolutionary process as well as the genetic consequences of polyploidization. In this study, we sequenced 3286 BACs from chromosome 7DL of bread wheat cv. Chinese Spring and integrated the unmapped contigs from IWGSC v1 and available PacBio sequences to close gaps present in the 7DL assembly. In total, 8043 out of 12 825 gaps, representing 3 491 264 bp, were closed. We then used the improved assembly of 7DL to perform comparative genomic analysis of bread wheat (Ta7DL) and its D donor, Aegilops tauschii (At7DL), to identify domestication signatures. Results showed a strong syntenic relationship between Ta7DL and At7DL, although some small rearrangements were detected at the distal regions. A total of 53 genes appear to be lost genes during wheat polyploidization, with 23% (12 genes) as RGA (disease resistance gene analogue). Furthermore, 86 positively selected genes (PSGs) were identified, considered to be domestication-related candidates. Finally, overlapping of QTLs obtained from GWAS analysis and PSGs indicated that TraesCS7D02G321000 may be one of the domestication genes involved in grain morphology. This study provides comparative information on the sequence, structure and organization between bread wheat and Ae. tauschii from the perspective of the 7DL chromosome, which contribute to better understanding of the evolution of wheat, and supports wheat crop improvement.

SNP markers for early identification of high molecular weight glutenin subunits (HMW-GSs) in bread wheat

A set of eight SNP markers was developed to facilitate the early selection of HMW-GS alleles in breeding programmes. In bread wheat (Triticum aestivum), the high molecular weight glutenin subunits (HMW-GSs) are the most important determinants of technological quality. Known to be very diverse, HMW-GSs are encoded by the tightly linked genes Glu-1-1 and Glu-1-2. Alleles that improve the quality of dough have been identified. Up to now, sodium dodecyl sulphate-polyacrylamide gel electrophoresis (SDS-PAGE) of grain proteins is the most widely used for their identification. To facilitate the early selection of HMW-GS alleles in breeding programmes, we developed DNA-based molecular markers. For each accession of a core collection (n = 364 lines) representative of worldwide bread wheat diversity, HMW-GSs were characterized by both genotyping and SDS-PAGE. Based on electrophoresis, we observed at least 8, 22 and 9 different alleles at the Glu-A1, Glu-B1 and Glu-D1 loci, respectively, including new variants. We designed a set of 17 single-nucleotide polymorphism (SNP) markers that were representative of the most frequent SDS-PAGE alleles at each locus. At Glu-A1 and Glu-D1, two and three marker-based haplotypes, respectively, captured the diversity of the SDS-PAGE alleles rather well. Discrepancies were found mainly for the Glu-B1 locus. However, statistical tests revealed that two markers at each Glu-B1 gene and their corresponding haplotypes were more significantly associated with the rheological properties of the dough than were the relevant SDS-PAGE alleles. To conclude, this study demonstrates that the SNP markers developed provide additional information on HMW-GS diversity. Two markers at Glu-A1, four at Glu-B1 and two at Glu-D1 constitute a useful toolbox for breeding wheat to improve end-use value.

CNVmap: A Method and Software To Detect and Map Copy Number Variants from Segregation Data

Single nucleotide polymorphisms (SNPs) are used widely for detecting quantitative trait loci, or for searching for causal variants of diseases. Nevertheless, structural variations such as copy-number variants (CNVs) represent a large part of natural genetic diversity, and contribute significantly to trait variation. Numerous methods and softwares based on different technologies (amplicons, CGH, tiling, or SNP arrays, or sequencing) have already been developed to detect CNVs, but they bypass a wealth of information such as genotyping data from segregating populations, produced, e.g., for QTL mapping. Here, we propose an original method to both detect and genetically map CNVs using mapping panels. Specifically, we exploit the apparent heterozygous state of duplicated loci: peaks in appropriately defined genome-wide allelic profiles provide highly specific signatures that identify the nature and position of the CNVs. Our original method and software can detect and map automatically up to 33 different predefined types of CNVs based on segregation data only. We validate this approach on simulated and experimental biparental mapping panels in two maize populations and one wheat population. Most of the events found correspond to having just one extra copy in one of the parental lines, but the corresponding allelic value can be that of either parent. We also find cases with two or more additional copies, especially in wheat, where these copies locate to homeologues. More generally, our computational tool can be used to give additional value, at no cost, to many datasets produced over the past decade from genetic mapping panels.

TaCYP81D5, one member in a wheat cytochrome P450 gene cluster, confers salinity tolerance via reactive oxygen species scavenging

As one of the largest gene families in plants, the cytochrome P450 monooxygenase genes (CYPs) are involved in diverse biological processes including biotic and abiotic stress response. Moreover, P450 genes are prone to expanding due to gene tandem duplication during evolution, resulting in generations of novel alleles with the neo-function or enhanced function. Here, the bread wheat (Triticum aestivum) gene TaCYP81D5 was found to lie within a cluster of five tandemly arranged CYP81D genes, although only a single such gene (BdCYP81D1) was present in the equivalent genomic region of the wheat relative Brachypodium distachyon. The imposition of salinity stress could up-regulate TaCYP81D5, but the effect was abolished in plants treated with an inhibitor of reactive oxygen species synthesis. In SR3, a wheat cultivar with an elevated ROS content, the higher expression and the rapider response to salinity of TaCYP81D5 were related to the chromatin modification. Constitutively expressing TaCYP81D5 enhanced the salinity tolerance both at seedling and reproductive stages of wheat via accelerating ROS scavenging. Moreover, an important component of ROS signal transduction, Zat12, was proven crucial in this process. Though knockout of solely TaCYP81D5 showed no effect on salinity tolerance, knockdown of BdCYP81D1 or all TaCYP81D members in the cluster caused the sensitivity to salt stress. Our results provide the direct evidence that TaCYP81D5 confers salinity tolerance in bread wheat and this gene is prospective for crop improvement.

Genome-wide and SNP network analyses reveal genetic control of spikelet sterility and yield-related traits in wheat

Revealing the genetic factors underlying yield and agronomic traits in wheat are an imperative need for covering the global food demand. Yield boosting requires a deep understanding of the genetic basis of grain yield-related traits (e.g., spikelet fertility and sterility). Here, we have detected much natural variation among ancient hexaploid wheat accessions in twenty-two agronomic traits collected over eight years of field experiments. A genome-wide association study (GWAS) using 15 K single nucleotide polymorphisms (SNPs) was applied to detect the genetic basis of studied traits. Subsequently, the GWAS output was reinforced via other statistical and bioinformatics analyses to detect putative candidate genes. Applying the genome-wide SNP-phenotype network defined the most decisive SNPs underlying the traits. Six pivotal SNPs, co-located physically within the genes encoding enzymes, hormone response, metal ion transport, and response to oxidative stress have been identified. Of these, metal ion transport and Gibberellin 2-oxidases (GA2oxs) genes showed strong involvement in controlling the spikelet sterility, which had not been reported previously in wheat. SNP-gene haplotype analysis confirmed that these SNPs influence spikelet sterility, especially the SNP co-located on the exon of the GA2ox gene. Interestingly, these genes were highly expressed in the grain and spike, demonstrating their pivotal role in controlling the trait. The integrative analysis strategy applied in this study, including GWAS, SNP-phenotype network, SNP-gene haplotype, expression analysis, and genome-wide prediction (GP), empower the identification of functional SNPs and causal genes. GP outputs obtained in this study are encouraging for the implementation of the traits to accelerate yield improvement by making an early prediction of complex yield-related traits in wheat. Our findings demonstrate the usefulness of the ancient wheat material as a valuable resource for yield-boosting. This is the first comprehensive genome-wide analysis for spikelet sterility in wheat, and the results provide insights into yield improvement.

Characterization of wheat homeodomain-leucine zipper family genes and functional analysis of TaHDZ5-6A in drought tolerance in transgenic Arabidopsis

BACKGROUND

Many studies in Arabidopsis and rice have demonstrated that HD-Zip transcription factors play important roles in plant development and responses to abiotic stresses. Although common wheat (Triticum aestivum L.) is one of the most widely cultivated and consumed food crops in the world, the function of the HD-Zip proteins in wheat is still largely unknown.

RESULTS

To explore the potential biological functions of HD-Zip genes in wheat, we performed a bioinformatics and gene expression analysis of the HD-Zip family. We identified 113 HD-Zip members from wheat and classified them into four subfamilies (I-IV) based on phylogenic analysis against proteins from Arabidopsis, rice, and maize. Most HD-Zip genes are represented by two to three homeoalleles in wheat, which are named as TaHDZX_ZA, TaHDZX_ZB, or TaHDZX_ZD, where X denotes the gene number and Z the wheat chromosome on which it is located. TaHDZs in the same subfamily have similar protein motifs and intron/exon structures. The expression profiles of TaHDZ genes were analysed in different tissues, at different stages of vegetative growth, during seed development, and under drought stress. We found that most TaHDZ genes, especially those in subfamilies I and II, were induced by drought stress, suggesting the potential importance of subfamily I and II TaHDZ members in the responses to abiotic stress. Compared with wild-type (WT) plants, transgenic Arabidopsis plants overexpressing TaHDZ5-6A displayed enhanced drought tolerance, lower water loss rates, higher survival rates, and higher proline content under drought conditions. Additionally, the transcriptome analysis identified a number of differentially expressed genes between 35S::TaHDZ5-6A transgenic and wild-type plants, many of which are involved in stress response.

CONCLUSIONS

Our results will facilitate further functional analysis of wheat HD-Zip genes, and also indicate that TaHDZ5-6A may participate in regulating the plant response to drought stress. Our experiments show that TaHDZ5-6A holds great potential for genetic improvement of abiotic stress tolerance in crops.

Examining the Effects of Temperature on Recombination in Wheat

Meiotic recombination plays a crucial role in the generation of new varieties. The effectiveness of recombination is limited by the distribution of crossover events, which in wheat and many other crops is skewed toward the distal regions of the chromosomes. Whole-genome sequencing of wheat has revealed that there are numerous important genes in the pericentromeric regions, which are inaccessible to manipulation due to the lack of crossover events. Studies in barley have shown that the distribution of recombination events can be shifted toward the centromeres by increasing temperature during meiosis. Here we present an analysis of the effects of temperature on the distribution and frequency of recombination events in wheat. Our data show that although increased temperature during meiosis does cause an inward shift in recombination distribution for some chromosomes, its overall utility is limited, with many genes remaining highly linked.

Thaumatin-like protein kinases: Molecular characterization and transcriptional profiling in five cereal crops

Thaumatin-like protein kinases (TLPKs) are defense related proteins having antimicrobial property. Herein, we identified two TLPKs in the genome of Brachypodium distachyon and Oryza sativa, four in Hordeum vulgare and Sorghum bicolor, and 16 in Triticum aestivum. All the TLPKs were located at only one chromosome in each plant except T. aestivum, where they were located on chromosome 2 and chromosome 3. Paralogous analysis suggested the occurrence of one duplication event (DE) in each B. distachyon and O. sativa, two in H. vulgare while four DEs in T. aestivum genome during the evolution of TLPKs. The majority of TLPKs were intron less, while a few contains one or two introns. The introns were found in each 0, 1 and 2 phase. Protein structure analysis suggested the occurrence of a thaumatin and a kinase domain with a transmembrane (TM) helix in each TLPK. Further, a thaumatin family signature motif "GX[GF]XCXT[GA]DCX(1,2)GX(2,3)C", a "REDDD" motif and 16 cysteine residues were found conserved in the majority of TLPKs. Expression analysis indicated variable expression of TLPKs in various tissues of different cereal crops. They were high expressing in reproductive tissues in B. distachyon, while in leaves in T. aestivum. Modulated expression of TaTLPKs in the presence of fungal pathogen, and heat, drought and salt stress in T. aestivum suggested their roles in stress response. Co-expression analysis showed interaction of TLPKs with various development and stress related genes. The results indicated diverse roles of TLPKs, which can be utilized for the development of eco-friendly pest resistant crops in future.

A unique chromosome translocation disrupting ClWIP1 leads to gynoecy in watermelon

To understand sex determination in watermelon (Citrullus lanatus), a spontaneous gynoecious watermelon mutant, XHBGM, was selected from the monoecious wild type XHB. Using map-based cloning, resequencing and fluorescence in situ hybridization analysis, a unique chromosome translocation between chromosome 2 and chromosome 3 was found in XHBGM. Based on the breakpoint location in chromosome 2, a putative C2H2 zinc finger transcription factor gene, ClWIP1 (gene ID Cla008537), an orthologue of the melon gynoecy gene CmWIP1, was disrupted. Using clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated system 9 to edit ClWIP1, we obtained gynoecious watermelon lines. Functional studies showed that ClWIP1 is expressed specifically in carpel primordia and is related to the abortion of carpel primordia in early floral development. To identify the cellular and metabolic processes associated with ClWIP1, we compared the shoot apex transcriptomes of two gynoecious mutants and their corresponding wild types. Transcriptome analysis showed that differentially expressed genes related to the ethylene and cytokinin pathways were upregulated in the gynoecious mutants. This study explores the molecular mechanism of sex determination in watermelon and provides a theoretical and technical basis for breeding elite gynoecious watermelon lines.

Thaumatin-like protein kinases: Molecular characterization and transcriptional profiling in five cereal crops

Thaumatin-like protein kinases (TLPKs) are defense related proteins having antimicrobial property. Herein, we identified two TLPKs in the genome of Brachypodium distachyon and Oryza sativa, four in Hordeum vulgare and Sorghum bicolor, and 16 in Triticum aestivum. All the TLPKs were located at only one chromosome in each plant except T. aestivum, where they were located on chromosome 2 and chromosome 3. Paralogous analysis suggested the occurrence of one duplication event (DE) in each B. distachyon and O. sativa, two in H. vulgare while four DEs in T. aestivum genome during the evolution of TLPKs. The majority of TLPKs were intron less, while a few contains one or two introns. The introns were found in each 0, 1 and 2 phase. Protein structure analysis suggested the occurrence of a thaumatin and a kinase domain with a transmembrane (TM) helix in each TLPK. Further, a thaumatin family signature motif "GX[GF]XCXT[GA]DCX(1,2)GX(2,3)C", a "REDDD" motif and 16 cysteine residues were found conserved in the majority of TLPKs. Expression analysis indicated variable expression of TLPKs in various tissues of different cereal crops. They were high expressing in reproductive tissues in B. distachyon, while in leaves in T. aestivum. Modulated expression of TaTLPKs in the presence of fungal pathogen, and heat, drought and salt stress in T. aestivum suggested their roles in stress response. Co-expression analysis showed interaction of TLPKs with various development and stress related genes. The results indicated diverse roles of TLPKs, which can be utilized for the development of eco-friendly pest resistant crops in future.

Redistribution of Meiotic Crossovers Along Wheat Chromosomes by Virus-Induced Gene Silencing

Meiotic recombination is the main driver of genetic diversity in wheat breeding. The rate and location of crossover (CO) events are regulated by genetic and epigenetic factors. In wheat, most COs occur in subtelomeric regions but are rare in centromeric and pericentric areas. The aim of this work was to increase COs in both "hot" and "cold" chromosomal locations. We used Virus-Induced gene Silencing (VIGS) to downregulate the expression of recombination-suppressing genes XRCC2 and FANCM and of epigenetic maintenance genes MET1 and DDM1 during meiosis. VIGS suppresses genes in a dominant, transient and non-transgenic manner, which is convenient in wheat, a hard-to-transform polyploid. F1 hybrids of a cross between two tetraploid lines whose genome was fully sequenced (wild emmer and durum wheat), were infected with a VIGS vector ∼ 2 weeks before meiosis. Recombination was measured in F2 seedlings derived from F1-infected plants and non-infected controls. We found significant up and down-regulation of CO rates along subtelomeric regions as a result of silencing either MET1, DDM1 or XRCC2 during meiosis. In addition, we found up to 93% increase in COs in XRCC2-VIGS treatment in the pericentric regions of some chromosomes. Silencing FANCM showed no effect on CO. Overall, we show that CO distribution was affected by VIGS treatments rather than the total number of COs which did not change. We conclude that transient silencing of specific genes during meiosis can be used as a simple, fast and non-transgenic strategy to improve breeding abilities in specific chromosomal regions.

Genome-wide analysis of MIKC-type MADS-box genes in wheat: pervasive duplications, functional conservation and putative neofunctionalization

Wheat (Triticum aestivum) is one of the most important crops worldwide. Given a growing global population coupled with increasingly challenging cultivation conditions, facilitating wheat breeding by fine-tuning important traits is of great importance. MADS-box genes are prime candidates for this, as they are involved in virtually all aspects of plant development. Here, we present a detailed overview of phylogeny and expression of 201 wheat MIKC-type MADS-box genes. Homoeolog retention is significantly above the average genome-wide retention rate for wheat genes, indicating that many MIKC-type homoeologs are functionally important and not redundant. Gene expression is generally in agreement with the expected subfamily-specific expression pattern, indicating broad conservation of function of MIKC-type genes during wheat evolution. We also found extensive expansion of some MIKC-type subfamilies, especially those potentially involved in adaptation to different environmental conditions like flowering time genes. Duplications are especially prominent in distal telomeric regions. A number of MIKC-type genes show novel expression patterns and respond, for example, to biotic stress, pointing towards neofunctionalization. We speculate that conserved, duplicated and neofunctionalized MIKC-type genes may have played an important role in the adaptation of wheat to a diversity of conditions, hence contributing to the importance of wheat as a global staple food.

A systems biology approach uncovers a gene co-expression network associated with cell wall degradability in maize

Understanding the mechanisms triggering variation of cell wall degradability is a prerequisite to improving the energy value of lignocellulosic biomass for animal feed or biorefinery. Here, we implemented a multiscale systems approach to shed light on the genetic basis of cell wall degradability in maize. We demonstrated that allele replacement in two pairs of near-isogenic lines at a region encompassing a major quantitative trait locus (QTL) for cell wall degradability led to phenotypic variation of a similar magnitude and sign to that expected from a QTL analysis of cell wall degradability in the F271 × F288 recombinant inbred line progeny. Using DNA sequences within the QTL interval of both F271 and F288 inbred lines and Illumina RNA sequencing datasets from internodes of the selected near-isogenic lines, we annotated the genes present in the QTL interval and provided evidence that allelic variation at the introgressed QTL region gives rise to coordinated changes in gene expression. The identification of a gene co-expression network associated with cell wall-related trait variation revealed that the favorable F288 alleles exploit biological processes related to oxidation-reduction, regulation of hydrogen peroxide metabolism, protein folding and hormone responses. Nested in modules of co-expressed genes, potential new cell-wall regulators were identified, including two transcription factors of the group VII ethylene response factor family, that could be exploited to fine-tune cell wall degradability. Overall, these findings provide new insights into the regulatory mechanisms by which a major locus influences cell wall degradability, paving the way for its map-based cloning in maize.

Functional mutation allele mining of plant architecture and yield-related agronomic traits and characterization of their effects in wheat

BACKGROUND

Wheat mutant resources with phenotypic variation have been developed in recent years. These mutants might carry favorable mutation alleles, which have the potential to be utilized in the breeding process. Plant architecture and yield-related features are important agronomic traits for wheat breeders and mining favorable alleles of these traits will improve wheat characteristics.

RESULTS

Here we used 190 wheat phenotypic mutants as material and by analyzing their SNP variation and phenotypic data, mutation alleles for plant architecture and yield-related traits were identified, and the genetic effects of these alleles were evaluated. In total, 32 mutation alleles, including three pleiotropic alleles, significantly associated with agronomic traits were identified from the 190 wheat mutant lines. The SNPs were distributed on 12 chromosomes and were associated with plant height (PH), tiller number, flag leaf angle (FLA), thousand grain weight (TGW), and other yield-related traits. Further phenotypic analysis of multiple lines carrying the same mutant allele was performed to determine the effect of the allele on the traits of interest. PH-associated SNPs on chromosomes 2BL, 3BS, 3DL, and 5DL might show additive effects, reducing PH by 10.0 cm to 31.3 cm compared with wild type, which means that these alleles may be favorable for wheat improvement. Only unfavorable mutation alleles that reduced TGW and tiller number were identified. A region on chromosome 5DL with mutation alleles for PH and TGW contained several long ncRNAs, and their sequences shared more than 90% identity with cytokinin oxidase/dehydrogenase genes. Some of the mutation alleles we mined were colocalized with previously reported QTLs or genes while others were novel; these novel alleles could also result in phenotypic variation.

CONCLUSION

Our results demonstrate that favorable mutation alleles are present in mutant resources, and the region between 409.5 to 419.8 Mb on chromosome 5DL affects wheat plant height and thousand grain weight.

A systematic review of the application of machine learning in the detection and classification of transposable elements

Background

Transposable elements (TEs) constitute the most common repeated sequences in eukaryotic genomes. Recent studies demonstrated their deep impact on species diversity, adaptation to the environment and diseases. Although there are many conventional bioinformatics algorithms for detecting and classifying TEs, none have achieved reliable results on different types of TEs. Machine learning (ML) techniques can automatically extract hidden patterns and novel information from labeled or non-labeled data and have been applied to solving several scientific problems.

Methodology

We followed the Systematic Literature Review (SLR) process, applying the six stages of the review protocol from it, but added a previous stage, which aims to detect the need for a review. Then search equations were formulated and executed in several literature databases. Relevant publications were scanned and used to extract evidence to answer research questions.

Results

Several ML approaches have already been tested on other bioinformatics problems with promising results, yet there are few algorithms and architectures available in literature focused specifically on TEs, despite representing the majority of the nuclear DNA of many organisms. Only 35 articles were found and categorized as relevant in TE or related fields.

Conclusions

ML is a powerful tool that can be used to address many problems. Although ML techniques have been used widely in other biological tasks, their utilization in TE analyses is still limited. Following the SLR, it was possible to notice that the use of ML for TE analyses (detection and classification) is an open problem, and this new field of research is growing in interest.

TRITEX: chromosome-scale sequence assembly of Triticeae genomes with open-source tools

Chromosome-scale genome sequence assemblies underpin pan-genomic studies. Recent genome assembly efforts in the large-genome Triticeae crops wheat and barley have relied on the commercial closed-source assembly algorithm DeNovoMagic. We present TRITEX, an open-source computational workflow that combines paired-end, mate-pair, 10X Genomics linked-read with chromosome conformation capture sequencing data to construct sequence scaffolds with megabase-scale contiguity ordered into chromosomal pseudomolecules. We evaluate the performance of TRITEX on publicly available sequence data of tetraploid wild emmer and hexaploid bread wheat, and construct an improved annotated reference genome sequence assembly of the barley cultivar Morex as a community resource.

TRITEX: chromosome-scale sequence assembly of Triticeae genomes with open-source tools

Chromosome-scale genome sequence assemblies underpin pan-genomic studies. Recent genome assembly efforts in the large-genome Triticeae crops wheat and barley have relied on the commercial closed-source assembly algorithm DeNovoMagic. We present TRITEX, an open-source computational workflow that combines paired-end, mate-pair, 10X Genomics linked-read with chromosome conformation capture sequencing data to construct sequence scaffolds with megabase-scale contiguity ordered into chromosomal pseudomolecules. We evaluate the performance of TRITEX on publicly available sequence data of tetraploid wild emmer and hexaploid bread wheat, and construct an improved annotated reference genome sequence assembly of the barley cultivar Morex as a community resource.

Assessing by Modeling the Consequences of Increased Recombination in Recurrent Selection of Oryza sativa and Brassica rapa

Meiotic recombination generates genetic diversity but in most species the number of crossovers per meiosis is limited. Previous modeling studies showed that increasing recombination can enhance response to selection. However, such studies did not assume a specific method of modifying recombination. Our objective was to test whether two methods used to increase recombination in plants could increase genetic gain in a population undergoing recurrent selection such as in genomic selection programs. The first method, in Oryza sativa, used a mutant of anti-crossover genes, increasing global recombination without affecting the recombination landscape shape. The second one used the ploidy level of a cross between Brassica rapa and Brassica napus, increasing recombination especially in pericentromeric regions. Our modeling framework used these recombination landscapes and sampled quantitative trait loci positions from the actual gene distributions. We simulated selection programs with initially a cross between two inbred lines, for two species. Increased recombination enhanced the response to selection. The amount of enhancement in the cumulative gain largely depended on the species and the number of quantitative trait loci (2, 10, 20, 50, 200 or 1000 per chromosome). Genetic gains were increased up to 30% after 20 generations. Furthermore, increasing recombination in cold regions was the most effective: the gain was larger by 25% with the first method and 34% with the second one in B. rapa, and 12% compared to 16% in O. sativa In summary, increased recombination enhances the genetic gain in long-term selection programs, with visible effects after four to five generations.

Identification of the quantitative trait loci controlling spike-related traits in hexaploid wheat (Triticum aestivum L.)

Totally, 48 loci responsible for six spike-related traits were identified in wheat, and a major locus QGl-4A for grain length was mapped and validated for marker-assisted selection in breeding. Wheat yield is determined by the number of spikes, number of grains per spike (GN), and one-thousand kernel weight (TKW), among which GN and TKW are greatly related to the spike development and thus the spike-related traits, including spike length (SL), number of spikelet per spike (SN), grain length (GL) and grain width (GW). To identify the key loci governing the spike-related traits (SL, SN, GN, TKW, GL and GW), we conducted the quantitative trait loci (QTL) analysis combined with wheat 660K SNP chip and Kompetitive allele-specific PCR (KASP) assay, using the F₂ and F_2:3 populations derived from Luohan6 (LH6) with big spike and grain and Zhengmai366 with small spike and grain, and identified a total of 48 QTLs on 18 chromosomes. Moreover, a major stable QTL for GL on chromosome 4A, designated as QGl-4A, was mapped into a 0.37 cM interval between KASP markers Xib4A-10 and Xib4A-12, corresponding to 20 Mb physical region in the Chinese Spring genome. This QTL explained 17.30% and 5.12% of the phenotypic variation for GL in the F₂ and F_2:3 populations. Further association analysis of flanking markers Xib4A-10 and Xib4A-12 in 192 wheat varieties showed that these two markers could be used for marker-assisted selection in breeding. These results provide valuable information for map-based cloning of the target genes involved in the regulation of spike-related traits in common wheat.

Engineering meiotic recombination pathways in rice

In the last 15 years, outstanding progress has been made in understanding the function of meiotic genes in the model dicot and monocot plants Arabidopsis and rice (Oryza sativa L.), respectively. This knowledge allowed to modulate meiotic recombination in Arabidopsis and, more recently, in rice. For instance, the overall frequency of crossovers (COs) has been stimulated 2.3- and 3.2-fold through the inactivation of the rice FANCM and RECQ4 DNA helicases, respectively, two genes involved in the repair of DNA double-strand breaks (DSBs) as noncrossovers (NCOs) of the Class II crossover pathway. Differently, the programmed induction of DSBs and COs at desired sites is currently explored by guiding the SPO11-1 topoisomerase-like transesterase, initiating meiotic recombination in all eukaryotes, to specific target regions of the rice genome. Furthermore, the inactivation of 3 meiosis-specific genes, namely PAIR1, OsREC8 and OsOSD1, in the Mitosis instead of Meiosis (MiMe) mutant turned rice meiosis into mitosis, thereby abolishing recombination and achieving the first component of apomixis, apomeiosis. The successful translation of Arabidopsis results into a crop further allowed the implementation of two breakthrough strategies that triggered parthenogenesis from the MiMe unreduced clonal egg cell and completed the second component of diplosporous apomixis. Here, we review the most recent advances in and future prospects of the manipulation of meiotic recombination in rice and potentially other major crops, all essential for global food security.

The wheat SnRK1α family and its contribution to Fusarium toxin tolerance

Deoxynivalenol (DON) is a mycotoxin produced by phytopathogenic Fusarium fungi in cereal grain and plays a role as a disease virulence factor. TaFROG (Triticum aestivum Fusarium Resistance Orphan Gene) enhances wheat resistance to DON and it interacts with a sucrose non-fermenting-1 (SNF1)-related protein kinase 1 catalytic subunit α (SnRK1α). This protein kinase family is central integrator of stress and energy signalling, regulating plant metabolism and growth. Little is known regarding the role of SnRK1α in the biotic stress response, especially in wheat. In this study, 15 wheat (Triticum aestivum) SnRK1α genes (TaSnRK1αs) belonging to four homoeologous groups were identified in the wheat genome. TaSnRK1αs are expressed ubiquitously in all organs and developmental stages apart from two members predominantly detected in grain. While DON treatment had either no effect or downregulated the transcription of TaSnRK1αs, it increased both the kinase activity associated with SnRK1α and the level of active (phosphorylated) SnRK1α. Down-regulation of two TaSnRK1αs homoeolog groups using virus induced gene silencing (VIGS) increased the DON-induced damage of wheat spikelets. Thus, we demonstrate that TaSnRK1αs contribute positively to wheat tolerance of DON and conclude that this gene family may provide useful tools for the improvement of crop biotic stress resistance.

TaCER1-1A is involved in cuticular wax alkane biosynthesis in hexaploid wheat and responds to plant abiotic stresses

To protect above-ground plant organs from excessive water loss, their surfaces are coated by waxes. The genes involved in wax formation have been investigated in detail in Arabidopsis but scarcely in crop species. Here, we aimed to isolate and characterize a CER1 enzyme responsible for formation of the very long-chain alkanes present in high concentrations especially during late stages of wheat development. On the basis of comparative wax and transcriptome analyses of various wheat organs, we selected TaCER1-1A as a primary candidate and demonstrated that it was located to the endoplasmic reticulum, the subcellular compartment for wax biosynthesis. A wheat nullisomic-tetrasomic substitution line lacking TaCER1-1A had significantly reduced amounts of C₃₃ alkane, whereas rice plants overexpressing TaCER1-1A showed substantial increases of C₂₅ -C₃₃ alkanes relative to wild type control. Similarly, heterologous expression of TaCER1-1A in Arabidopsis wild type and the cer1 mutant resulted in increased levels of unbranched alkanes, iso-branched alkanes and alkenes. Finally, the expression of TaCER1-1A was found activated by abiotic stresses and abscisic acid treatment, resulting in increased production of alkanes in wheat. Taken together, our results demonstrate that TaCER1-1A plays an important role in wheat wax alkane biosynthesis and involved in responding to drought and other environmental stresses.

A Co-Expression Network in Hexaploid Wheat Reveals Mostly Balanced Expression and Lack of Significant Gene Loss of Homeologous Meiotic Genes Upon Polyploidization

Polyploidization has played an important role in plant evolution. However, upon polyploidization, the process of meiosis must adapt to ensure the proper segregation of increased numbers of chromosomes to produce balanced gametes. It has been suggested that meiotic gene (MG) duplicates return to a single copy following whole genome duplication to stabilize the polyploid genome. Therefore, upon the polyploidization of wheat, a hexaploid species with three related (homeologous) genomes, the stabilization process may have involved rapid changes in content and expression of MGs on homeologous chromosomes (homeologs). To examine this hypothesis, sets of candidate MGs were identified in wheat using co-expression network analysis and orthology informed approaches. In total, 130 RNA-Seq samples from a range of tissues including wheat meiotic anthers were used to define co-expressed modules of genes. Three modules were significantly correlated with meiotic tissue samples but not with other tissue types. These modules were enriched for GO terms related to cell cycle, DNA replication, and chromatin modification and contained orthologs of known MGs. Overall, 74.4% of genes within these meiosis-related modules had three homeologous copies which was similar to other tissue-related modules. Amongst wheat MGs identified by orthology, rather than co-expression, the majority (93.7%) were either retained in hexaploid wheat at the same number of copies (78.4%) or increased in copy number (15.3%) compared to ancestral wheat species. Furthermore, genes within meiosis-related modules showed more balanced expression levels between homeologs than genes in non-meiosis-related modules. Taken together, our results do not support extensive gene loss nor changes in homeolog expression of MGs upon wheat polyploidization. The construction of the MG co-expression network allowed identification of hub genes and provided key targets for future studies.

Reduced stomatal density in bread wheat leads to increased water-use efficiency

Wheat is a staple crop, frequently cultivated in water-restricted environments. Improving crop water-use efficiency would be desirable if grain yield can be maintained. We investigated whether a decrease in wheat stomatal density via the manipulation of epidermal patterning factor (EPF) gene expression could improve water-use efficiency. Our results show that severe reductions in stomatal density in EPF-overexpressing wheat plants have a detrimental outcome on yields. However, wheat plants with a more moderate reduction in stomatal density (i.e. <50% reduction in stomatal density on leaves prior to tillering) had yields indistinguishable from controls, coupled with an increase in intrinsic water-use efficiency. Yields of these moderately reduced stomatal density plants were also comparable with those of control plants under conditions of drought and elevated CO2. Our data demonstrate that EPF-mediated control of wheat stomatal development follows that observed in other grasses, and we identify the potential of stomatal density as a tool for breeding wheat plants that are better able to withstand water-restricted environments without yield loss.

Genome-wide variation patterns between landraces and cultivars uncover divergent selection during modern wheat breeding

Genetic diversity, population structure, LD decay, and selective sweeps in 687 wheat accessions were analyzed, providing relevant guidelines to facilitate the use of the germplasm in wheat breeding. Common wheat (Triticum aestivum L.) is one of the most widely grown crops in the world. Landraces were subjected to strong human-mediated selection in developing high-yielding, good quality, and widely adapted cultivars. To investigate the genome-wide patterns of allelic variation, population structure and patterns of selective sweeps during modern wheat breeding, we tested 687 wheat accessions, including landraces (148) and cultivars (539) mainly from China and Pakistan in a wheat 90 K single nucleotide polymorphism array. Population structure analysis revealed that cultivars and landraces from China and Pakistan comprised three relatively independent genetic clusters. Cultivars displayed lower nucleotide diversity and a wider average LD decay across whole genome, indicating allelic erosion and a diversity bottleneck due to the modern breeding. Analysis of genetic differentiation between landraces and cultivars from China and Pakistan identified allelic variants subjected to selection during modern breeding. In total, 477 unique genome regions showed signatures of selection, where 109 were identified in both China and Pakistan germplasm. The majority of genomic regions were located in the B genome (225), followed by the A genome (175), and only 77 regions were located in the D genome. EigenGWAS was further used to identify key selection loci in modern wheat cultivars from China and Pakistan by comparing with global winter wheat and spring wheat diversity panels, respectively. A few known functional genes or loci found within these genome regions corresponded to known phenotypes for disease resistance, vernalization, quality, adaptability and yield-related traits. This study uncovered molecular footprints of modern wheat breeding and explained the genetic basis of polygenic adaptation in wheat. The results will be useful for understanding targets of modern wheat breeding, and in devising future breeding strategies to target beneficial alleles currently not pursued.

Highly Efficient and Heritable Targeted Mutagenesis in Wheat via the Agrobacterium tumefaciens-Mediated CRISPR/Cas9 System

The CRISPR/Cas9 system has been successfully used in hexaploid wheat. Although it has been reported that the induced mutations can be passed to the next generation, gene editing and transmission patterns in later generations still need to be studied. In this study, we demonstrated that the CRISPR/Cas9 system could achieve efficient mutagenesis in five wheat genes via Agrobacterium-mediated transformation of an sgRNA targeting the D genome, an sgRNA targeting both the A and B homologues and three tri-genome guides targeting the editing of all three homologues. High mutation rates and putative homozygous or biallelic mutations were observed in the T0 plants. The targeted mutations could be stably inherited by the next generation, and the editing efficiency of each mutant line increased significantly across generations. The editing types and inheritance of targeted mutagenesis were similar, which were not related to the targeted subgenome number. The presence of Cas9/sgRNA could cause new mutations in subsequent generations, while mutated lines without Cas9/sgRNA could retain the mutation type. Additionally, off-target mutations were not found in sequences that were highly homologous to the selected sgRNA sequences. Overall, the results suggested that CRISPR/Cas9-induced gene editing via Agrobacterium-mediated transformation plays important roles in wheat genome engineering.

Draft Genome Sequence of Tubulinosema ratisbonensis, a Microsporidian Species Infecting the Model Organism Drosophila melanogaster

We present the draft genome sequence of Tubulinosema ratisbonensis, a microsporidium species infecting Drosophila melanogaster A total of 3,013 protein-encoding genes and an array of transposable elements were identified. This work represents a necessary step to develop a novel model of host-parasite relationships using the highly tractable genetic model D. melanogaster.