The Evolution, Expression Patterns, and Domestication Selection Analysis of the Annexin Gene Family in the Barley Pan-Genome

Plant annexins constitute a conserved protein family that plays crucial roles in regulating plant growth and development, as well as in responses to both biotic and abiotic stresses. In this study, a total of 144 annexin genes were identified in the barley pan-genome, comprising 12 reference genomes, including cultivated barley, landraces, and wild barley. Their chromosomal locations, physical-chemical characteristics, gene structures, conserved domains, and subcellular localizations were systematically analyzed to reveal the certain differences between wild and cultivated populations. Through a cis-acting element analysis, co-expression network, and large-scale transcriptome analysis, their involvement in growth, development, and responses to various stressors was highlighted. It is worth noting that HvMOREXann5 is only expressed in pistils and anthers, indicating its crucial role in reproductive development. Based on the resequencing data from 282 barley accessions worldwide, genetic variations in thefamily were investigated, and the results showed that 5 out of the 12 identified HvMOREXanns were affected by selection pressure. Genetic diversity and haplotype frequency showed notable reductions between wild and domesticated barley, suggesting that a genetic bottleneck occurred on the annexin family during the barley domestication process. Finally, qRT-PCR analysis confirmed the up-regulation of HvMOREXann7 under drought stress, along with significant differences between wild accessions and varieties. This study provides some insights into the genome organization and genetic characteristics of the annexin gene family in barley at the pan-genome level, which will contribute to better understanding its evolution and function in barley and other crops.

Genome-wide association studies of peduncle length in wheat under rain-fed and irrigating field conditions

Drought is one of the most destructive environmental factors limiting wheat production and food security globally. Peduncle length (PLE) is an important morphological trait to determine plant architecture, photosynthate transport, and yield formation, which is also considered a useful index for drought tolerance in wheat. However, the genetic basis of wheat PLE is not well studied at present. Here, a large-scale genome-wide association study (GWAS) of PLE was performed using a panel of 282 wheat accessions with the Wheat 660K SNP array genotyping under rain-fed and irrigating field conditions. Totally, 1,301 significant marker-trait associations (MTAs) were identified using the threshold of p-value < 4.16 × 10^-4, five of which were high-confidence. Furthermore, combining GWAS intervals, previously reported QTLs, expression levels, homologous genes, and selected sweep analysis, a total of 5 candidate genes were detected to associate with drought stress. Moreover, the expression levels of TraesCS2A02G082100 were significantly up-regulated under drought conditions and co-localized in the selected sweep region, suggesting it is a drought-responsive gene. Our results shed light on the genetic basis underlying wheat drought tolerance, which accelerates the marker-assistant selection and genetic improvement through genomic breeding in wheat.

A large-scale genomic association analysis identifies the candidate causal genes conferring stripe rust resistance under multiple field environments

The incorporation of resistance genes into wheat commercial varieties is the ideal strategy to combat stripe or yellow rust (YR). In a search for novel resistance genes, we performed a large-scale genomic association analysis with high-density 660K single nucleotide polymorphism (SNP) arrays to determine the genetic components of YR resistance in 411 spring wheat lines. Following quality control, 371 972 SNPs were screened, covering over 50% of the high-confidence annotated gene space. Nineteen stable genomic regions harbouring 292 significant SNPs were associated with adult-plant YR resistance across nine environments. Of these, 14 SNPs were localized in the proximity of known loci widely used in breeding. Obvious candidate SNP variants were identified in certain confidence intervals, such as the cloned gene Yr18 and the major locus on chromosome 2BL, despite a large extent of linkage disequilibrium. The number of causal SNP variants was refined using an independent validation panel and consideration of the estimated functional importance of each nucleotide polymorphism. Interestingly, four natural polymorphisms causing amino acid changes in the gene TraesCS2B01G513100 that encodes a serine/threonine protein kinase (STPK) were significantly involved in YR responses. Gene expression and mutation analysis confirmed that STPK played an important role in YR resistance. PCR markers were developed to identify the favourable TraesCS2B01G513100 haplotype for marker-assisted breeding. These results demonstrate that high-resolution SNP-based GWAS enables the rapid identification of putative resistance genes and can be used to improve the efficiency of marker-assisted selection in wheat disease resistance breeding.

Genomic and Phenotypic Divergence in Wild Barley Driven by Microgeographic Adaptation

Microgeographic adaptation is a fundamental driving force of evolution, but the underlying causes remain undetermined. Here, the phenotypic, genomic and transcriptomic variations of two wild barley populations collected from sharply divergent and adjacent micro-geographic sites to identify candidate genes associated with edaphic local adaptation are investigated. Common garden and reciprocal transplant studies show that large phenotypic differentiation and local adaptation to soils occur between these populations. Genetic, phylogenetic and admixture analyses based on population resequencing show that significant genetic divergences occur between basalt and chalk populations. These divergences are consistent with the phenotypic variations observed in the field. Genome sweep analyses reveal 162.7 Mb of selected regions driven by edaphic local adaptation, in which 445 genes identified, including genes associated with root architecture, metal transport/detoxification, and ABA signaling. When the phenotypic, genomic and transcriptomic data are combined, HvMOR, encoding an LBD transcription factor, is determined to be the vital candidate for regulating the root architecture to adapt to edaphic conditions at the microgeographic scale. This study provides new insights into the genetic basis of edaphic adaptation and demonstrates that edaphic factors may contribute to the evolution and speciation of barley.

GENOME‐WIDE IDENTIFICATION AND ANALYSIS OF THE CslF GENE FAMILY IN BARLEY (Hordeum vulgare L.)

The Cellulose synthase-like (Csl) F family has been considered as one of the most crucial genes regulating β-glucan synthesis. It is a cereal cell wall component, holding advantages for human nutrition, but disadvantages in animal nutrition, malting and brewing industries. Based on a genome-wide search method, present study identified barley CslF gene family members by considering the importance of (1,3;1,4)-β-D-glucan and newly completed barley genome. A sum of Eighteen CslF genes were recognized in the barley genome.  Then, phylogenetic analyses classified them into 3 groups, that shared conserved motif compositions. A new motif, D,D,WQxxD was also found, which was responsible for the cellulose synthase. Furthermore, using RNA-seq data,  the HvCslFs expression profiles were systematically examined in different tissues and tissue-specific candidates were found. Lastly, interaction network analysis identified 11 CslF genes involved in the interaction network. All together, prsent task provides valuable evidence about the genomic organization and evolutionary relationship of the CslF gene family in barley, and facilitate the functional surveys of CslF genes in barley and beyond.

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.

N6-methyladenosine regulatory machinery in plants: composition, function and evolution

N6-methyladenosine (m6A) RNA methylation, one of the most pivotal internal modifications of RNA, is a conserved post-transcriptional mechanism to enrich and regulate genetic information in eukaryotes. The scope and function of this modification in plants has been an intense focus of study, especially in model plant systems. The characterization of plant m6A writers, erasers and readers, as well as the elucidation of their functions, is currently one of the most fascinating hotspots in plant biology research. The functional analysis of m6A in plants will be booming in the foreseeable future, which could contribute to crop genetic improvement through epitranscriptome manipulation. In this review, we systematically analysed and summarized recent advances in the understanding of the structure and composition of plant m6A regulatory machinery, and the biological functions of m6A in plant growth, development and stress response. Finally, our analysis showed that the evolutionary relationships between m6A modification components were highly conserved across the plant kingdom.

Shifting the limits in wheat research and breeding using a fully annotated reference genome

An annotated reference sequence representing the hexaploid bread wheat genome in 21 pseudomolecules has been analyzed to identify the distribution and genomic context of coding and noncoding elements across the A, B, and D subgenomes. With an estimated coverage of 94% of the genome and containing 107,891 high-confidence gene models, this assembly enabled the discovery of tissue- and developmental stage-related coexpression networks by providing a transcriptome atlas representing major stages of wheat development. Dynamics of complex gene families involved in environmental adaptation and end-use quality were revealed at subgenome resolution and contextualized to known agronomic single-gene or quantitative trait loci. This community resource establishes the foundation for accelerating wheat research and application through improved understanding of wheat biology and genomics-assisted breeding.

Conservation analysis of long non-coding RNAs in plants

Long non-coding RNAs (lncRNAs) are gene regulators that have vital roles in development and adaptation to the environment in eukaryotes. However, the structural and evolutionary analyses of plant lncRNAs are limited. In this study, we performed an analysis of lncRNAs in five monocot and five dicot species. Our results showed that plant lncRNA genes were generally shorter and had fewer exons than protein-coding genes. The numbers of lncRNAs were positively correlated with the numbers of protein-coding genes in different plant species, despite a high range of variation. Sequence conservation analysis showed that the majority of lncRNAs had high sequence conservation at the intra-species and sub-species levels, reminiscent of protein-coding genes. At the inter-species level, a subset of lncRNAs were highly diverged at the nucleotide level, but conserved by position. Interestingly, we found that plant lncRNAs have identical splicing signals, and those which can form precursors or targets of miRNAs have a conservative identity in different species. We also revealed that most of the lowly expressed lncRNAs were tissue-specific, while those highly conserved were constitutively transcribed. Meanwhile, we characterized a subset of rice lncRNAs that were co-expressed with their adjacent protein-coding genes, suggesting they may play cis-regulatory roles. These results will contribute to understanding the biological significance and evolution of lncRNAs in plants.

Characterization of microRNAs and their targets in wild barley (Hordeum vulgare subsp. spontaneum) using deep sequencing

MicroRNAs (miRNA) are a class of small, endogenous RNAs that play a negative regulatory role in various developmental and metabolic processes of plants. Wild barley (Hordeum vulgare subsp. spontaneum), as the progenitor of cultivated barley (Hordeum vulgare subsp. vulgare), has served as a valuable germplasm resource for barley genetic improvement. To survey miRNAs in wild barley, we sequenced the small RNA library prepared from wild barley using the Illumina deep sequencing technology. A total of 70 known miRNAs and 18 putative novel miRNAs were identified. Sequence analysis revealed that all of the miRNAs identified in wild barley contained the highly conserved hairpin sequences found in barley cultivars. MiRNA target predictions showed that 12 out of 52 miRNA families were predicted to target transcription factors, including 8 highly conserved miRNA families in plants and 4 wheat-barley conserved miRNA families. In addition to transcription factors, other predicted target genes were involved in diverse physiological and metabolic processes and stress defense. Our study for the first time reported the large-scale investigation of small RNAs in wild barley, which will provide essential information for understanding the regulatory role of miRNAs in wild barley and also shed light on future practical utilization of miRNAs for barley improvement.

Global Identification of MicroRNAs and Their Targets in Barley under Salinity Stress

Salinity is a major limiting factor for agricultural production worldwide. A better understanding of the mechanisms of salinity stress response will aid efforts to improve plant salt tolerance. In this study, a combination of small RNA and mRNA degradome sequencing was used to identify salinity responsive-miRNAs and their targets in barley. A total of 152 miRNAs belonging to 126 families were identified, of which 44 were found to be salinity responsive with 30 up-regulated and 25 down-regulated respectively. The majority of the salinity-responsive miRNAs were up-regulated at the 8h time point, while down-regulated at the 3h and 27h time points. The targets of these miRNAs were further detected by degradome sequencing coupled with bioinformatics prediction. Finally, qRT-PCR was used to validate the identified miRNA and their targets. Our study systematically investigated the expression profile of miRNA and their targets in barley during salinity stress phase, which can contribute to understanding how miRNAs respond to salinity stress in barley and other cereal crops.

Preparation of high molecular weight gDNA and bacterial artificial chromosome (BAC) libraries in plants

Bacterial artificial chromosome (BAC) libraries are extremely valuable large-insert DNA libraries for physical mapping, positional cloning, comparative genomic analysis, complete genome sequencing, and evolutionary studies. Due to their stability and relative simplicity BAC libraries are most preferred over other approaches for cloning large genomic DNA fragments for large-insert libraries. Isolation of intact high molecular weight (HMW) DNA is a critical step underlying the success of large-insert genomic DNA library construction. It requires the isolation of purified nuclei, embedding them into LMP agarose plugs, restriction digestion of the plugs, and quite often size selection using PFGE and electro-elution of insert DNA. The construction of BAC libraries is complex and challenging for most molecular laboratories. To facilitate the construction of BAC libraries, we present a step-by-step protocol for isolation of HMW DNA and construction of plant BAC libraries.

Dispersion and domestication shaped the genome of bread wheat

Despite the international significance of wheat, its large and complex genome hinders genome sequencing efforts. To assess the impact of selection on this genome, we have assembled genomic regions representing genes for chromosomes 7A, 7B and 7D. We demonstrate that the dispersion of wheat to new environments has shaped the modern wheat genome. Most genes are conserved between the three homoeologous chromosomes. We found differential gene loss that supports current theories on the evolution of wheat, with greater loss observed in the A and B genomes compared with the D. Analysis of intervarietal polymorphisms identified fewer polymorphisms in the D genome, supporting the hypothesis of early gene flow between the tetraploid and hexaploid. The enrichment for genes on the D genome that confer environmental adaptation may be associated with dispersion following wheat domestication. Our results demonstrate the value of applying next-generation sequencing technologies to assemble gene-rich regions of complex genomes and investigate polyploid genome evolution. We anticipate the genome-wide application of this reduced-complexity syntenic assembly approach will accelerate crop improvement efforts not only in wheat, but also in other polyploid crops of significance.

Development of chromosome-arm-specific microsatellite markers in Triticum aestivum (Poaceae) using NGS technology

PREMISE OF THE STUDY

The aim of this study was to assess the feasibility of developing chromosome-arm-specific microsatellite markers in wheat on a large scale based on chromosome survey sequences obtained with next-generation sequencing (NGS) technology.

METHODS AND RESULTS

The Illumina Hi Seq2000 sequencing platform was used to sequence DNA of isolated wheat chromosome-arm 7DL. The data were assembled and microsatellite loci were identified computationally. In total, 16315 microsatellites were identified from 161061 assembled contigs. Thirty-three markers were randomly selected for validation across 20 diverse wheat cultivars. Two nulli-tetrasomic stocks were also screened to validate the specificity of the newly developed markers.

CONCLUSIONS

This is the first study on identification of chromosome-arm-specific microsatellite markers using NGS technology. These new chromosome-arm-specific markers will facilitate saturation of the 7DL genetic map, and their availability will support genetic mapping and positional cloning in wheat.

Identification and characterization of microRNAs from barley (Hordeum vulgare L.) by high-throughput sequencing

MicroRNAs (miRNAs) are a class of endogenous RNAs that regulates the gene expression involved in various biological and metabolic processes. Barley is one of the most important cereal crops worldwide and is a model organism for genetic and genomic studies in Triticeae species. However, the miRNA research in barley has lagged behind other model species in grass family. To obtain more information of miRNA genes in barley, we sequenced a small RNA library created from a pool of equal amounts of RNA from four different tissues using Solexa sequencing. In addition to 126 conserved miRNAs (58 families), 133 novel miRNAs belonging to 50 families were identified from this sequence data set. The miRNA* sequences of 15 novel miRNAs were also discovered, suggesting the additional evidence for existence of these miRNAs. qRT-PCR was used to examine the expression pattern of six randomly selected miRNAs. Some miRNAs involved in drought and salt stress response were also identified. Furthermore, the potential targets of these putative miRNAs were predicted using the psRNATarget tools. Our results significantly increased the number of novel miRNAs in barley, which should be useful for further investigation into the biological functions and evolution of miRNAs in barley and other species.

Genetic diversity of mitochondrial cytochrome b gene in Chinese native buffalo

The origins of the domestic water buffalo remain contentious. To better understand the origins of Chinese water buffalo, we sequenced the complete mitochondrial cytochrome b (MT-CYB) gene from 270 individuals representing 13 Chinese domestic swamp buffalo populations. We found genetic evidence of introgression of river buffalo into Chinese swamp buffalo herds. Swamp buffalo haplotypes can be divided into two highly divergent lineages (A and B), suggesting that Chinese native swamp buffalo have two maternal origins. We found that the A→G transition in the buffalo MT-CYB gene stop codon resulted in buffalo haplotypes being terminated by one of two stop codons: AGA or AGG. AGA is common to river buffalo and lineage A of swamp buffalo, while AGG is specific to lineage B of swamp buffalo. Lineage A appears to have been domesticated in China. Further genetic evidence is required to clarify the origins of lineage B.

A MAPK gene from Dead Sea fungus confers stress tolerance to lithium salt and freezing-thawing: Prospects for saline agriculture

The Dead Sea is one of the most saline lakes on earth ( approximately 340 g/liter salinity) and is approximately 10 times saltier than the oceans. Eurotium herbariorum, a common fungal species, was isolated from its water. EhHOG gene, encoding a mitogen-activated protein kinase (MAPK) that plays an essential role in the osmoregulatory pathway in yeast and many other eukaryotes, was isolated from E. herbariorum. The deduced amino acid sequences of EhHOG indicated high similarity with homologous genes from Aspergillus nidulans, Saccharomyces cerevisiae, and Schizosaccharomyces pombe and contained a TGY motif for phosphorylation by MAPK kinase. When EhHOG was expressed in S. cerevisiae hog1Delta mutant, the growth and aberrant morphology of hog1Delta mutant was restored under high osmotic stress condition. Moreover, intracellular glycerol content in the transformant increased to a much higher level than that in the mutant during salt-stress situations. hog1Delta mutant complemented by EhHOG outperformed the wild type or had higher genetic fitness under high Li(+) and freezing-thawing conditions. The present study revealed the putative presence of a high-osmolarity glycerol response (HOG) pathway in E. herbariorum and the significance of EhHOG in osmotic regulation, heat stress, freeze stress, and oxidative stress. The Dead Sea is becoming increasingly more saline while the fungi living in it evolutionarily adapt to its high-saline environment, particularly with the extraordinarily high Li(+) concentration. The Dead Sea is potentially an excellent model for studies of evolution under extreme environments and is an important gene pool for future agricultural genetic engineering prospects.

Wild barley eibi1 mutation identifies a gene essential for leaf water conservation

Drought is a major abiotic stress that limits plant growth and crop productivity. A spontaneous wilty mutant (eibi1) hypersensitive to drought was identified from wild barley (Hordeum spontaneum Koch). eibi1 showed the highest relative water loss rate among the known wilty mutants, which indicates that eibi1 is one of the most drought-sensitive mutants. eibi1 had the same abscisic acid (ABA) level, the same ability to accumulate stress-induced ABA, and the same stomatal movement in response to light, dark, drought, and exogenous ABA as the wild type, revealing that eibi1 was neither an ABA-deficient nor an ABA-insensitive mutant. The eibi1 leaves had a larger chlorophyll efflux rate in 80% ethanol than the wild-type leaves; and the transpiration rate of eibi1 was more closely related to chlorophyll efflux rate than to stomatal density, demonstrating that the cuticle of eibi1 was defective. eibi1 will be a promising candidate to study the actual barrier layer in the cuticle that limits water loss of the plant. Exogenous ABA reduced leaf length growth in eibi1 more than in the wild type, implying an interaction on plant growth of ABA signal transduction and the eibi1 product. One may infer that the eibi1 product may reverse the growth inhibition induced by ABA.

Non-gridded library: a new approach for BAC (bacterial artificial chromosome) exploitation in hexaploid wheat (Triticum aestivum)

The feasibility of exploiting non-gridded bacterial artificial chromosome (BAC) libraries and some major factors affecting the efficiency of handling such libraries were studied in hexaploid wheat. Even for a bacterial culture containing only 55% recombinants, some 2000 BAC clones with inserts ranging from 45 to 245 kb could be pooled. The pooled BAC clones could be amplified by culturing for up to 6 h without losing any target clones. These results imply that even for hexaploid wheat, which has an extremely large genome, some 250 pools are sufficient for a BAC library that should satisfy many research objectives. This non-gridded strategy would dramatically reduce the cost and make robotic equipment non-essential in exploiting BAC technology. To construct a representative library and to minimise clone competition, thawing and re-freezing ligation mixtures and bacterial cultures should be avoided in BAC library construction and application.

Polymorphisms in the α-amy1 gene of wild and cultivated barley revealed by the polymerase chain reaction

α-Amylases are the key enzymes involved in the hydrolysis of starch in plants. The polymerase chain reaction (PCR) was used to detect polymorphisms in the length of amplified sequences between the annealing sites of two primers derived from published α-amy1 gene sequences in barley. These two primers (Bsw1 and Bsw7), flanking the promoter region and the first exon, amplified two PCR fragments in barley. One of the amplified products, with the expected length of 820 bp, appeared together with another shorter PCR band of around 750 bp. This 750-bp fragment seems to be derived from an α-amylase gene not reported previously. Both of the PCR products could be amplified from the two-rowed barley varieties tested, including cv Himalaya from which the sequence information was obtained. Five of the six-rowed barley varieties also have the two PCR fragments whereas another two have only the long fragment. These two fragments seem to be unique to barley, neither of them could be amplified from other cereals; for example, wheat, rye or sorghum. These two α-amylase fragments were mapped to the long arm of 6H, the location of the α-amy1 genes, using wheat-barley addition lines. Amplification of genomic DNA from wild barley accessions with primers Bsw1 and Bsw7 indicated that both of the fragments could be present, or the long and short fragments could be present alone. The results also demonstrated that the genes specifying these two fragments could be independent from each other in barley. The conserved banding pattern of these two fragments in the two-rowed barley varieties implies that artificial selection from these genes may have played an important role in the evolution of cultivated barley from wild barley.

Identification and mapping of polymorphisms in cereals based on the polymerase chain reaction

The polymerase chain reaction (PCR) can be used to detect polymorphisms in the length of amplified sequences between the annealing sites of two synthetic DNA primers. When the distance varies between two individuals then the banding pattern generated by the PCR reaction is essentially a genetic polymorphism and can be mapped in the same way as other genetic markers. This procedure has been used in a number of eukaryotes. Here we report the use of PCR to detect genetic polymorphisms in cereals. Known gene sequences can be used to design primers and detect polymorphic PCR products. This is demonstrated with primers to the α-amylase gene family. A second approach is to use semi-random primers to target diverse regions of the genome. For this purpose the consensus sequences at the intron-exon splice junctions were used. The targeting of the intronexon splice junctions in conjunction with primers of random and defined sequences, such as α-amylase, provides a source of extensive variation in PCR products. These polymorphisms can be mapped as standard genetic markers.