Combining genetical genomics and bulked segregant analysis-based differential expression: an approach to gene localization

Investigation of genetic markers for intramuscular fat in the hybrid Wagyu cattle with bulked segregant analysis

ulked Segregant Analysis (BSA) is a rapid strategy for identifying genetic markers in specific regions of the phenotypical population and it has been widely used for QTLs mapping in smaller mixed F2 and F3 populations. We applied a modified BSA method to assessed genome-wide homozygous and heterozygous linkage patterns in the Chinese Wagyu Beef Cattle F2/F3 mixed population. Two overlapped regions from F2 and F3 populations on autosomes were found with high-density heterozygote alleles between high and low intramuscular fat groups. Regions from 24.8 M ~ 29.6 M of chromosome 23 were identified as most significantly correlated to the intramuscular fat in our samples. We also identified other 4 potential loci on chromosomes 5, 9, 15, and 21 correlated with Intramuscular fat. This study provided a novel low-cost method for QTLs mapping and identify molecular markers of phenotypical changes in a small mixed population.

Mapping of the male sterile mutant gene ftms in Brassica rapa L. ssp. pekinensis via BSR-Seq combined with whole-genome resequencing

A male sterile mutant was created by ⁶⁰Co γ-rays of microspores isolated from Chinese cabbage DH line 'FT'. A candidate gene for the male sterile trait was identified as Bra010198. Male sterility is used for hybrid seed production in Chinese cabbage. In this study, we derived a male sterile mutant (ftms) from Chinese cabbage DH line 'FT' by irradiating microspores with ⁶⁰Co γ-rays and realized the rapid trait transformation from male fertility to sterility for creating valuable breeding materials. Genetic analysis indicated that the male sterile trait is controlled by a single recessive nuclear gene, ftms. Microspore development in mutant ftms was aborted at the tetrad stage and associated with severely retarded degeneration and vacuolation of tapetum. Using BSR-seq analysis, the candidate region for ftms was mapped on chromosome A05. A large F₂ population was created, and the region was narrowed to approximately 1.7-Mb between markers Indel20 and Indel14 via linkage analysis. The recombination frequency was extremely suppressed because the region was located on the chromosome A05 centromere. Whole-genome resequencing of mutant ftms and wild-type 'FT' aligned only one nonsynonymous SNP to Bra010198; this gene is a homolog of Arabidopsis KNS4/UPEX1, which encodes a putative β-(1,3)-galactosyltransferase that controls pollen exine development. Comparative sequencing verified the SNP position on the fifth exon of Bra010198 in mutant ftms. Further genotyping revealed that the male sterile phenotype was fully co-segregated with this SNP. Quantitative real-time PCR indicated that Bra0101918 specifically expressed in stamen. The data presented herein suggested that Bra010198 is a strong candidate gene for ftms. Hence, we developed a male sterile line for potential application in breeding and expanded the knowledge about the molecular mechanism underlying male sterility in Chinese cabbage.

A comparative analysis of nonhost resistance across the two Triticeae crop species wheat and barley

BACKGROUND

Nonhost resistance (NHR) protects plants against a vast number of non-adapted pathogens which implicates a potential exploitation as source for novel disease resistance strategies. Aiming at a fundamental understanding of NHR a global analysis of transcriptome reprogramming in the economically important Triticeae cereals wheat and barley, comparing host and nonhost interactions in three major fungal pathosystems responsible for powdery mildew (Blumeria graminis ff. ssp.), cereal blast (Magnaporthe sp.) and leaf rust (Puccinia sp.) diseases, was performed.

RESULTS

In each pathosystem a significant transcriptome reprogramming by adapted- or non-adapted pathogen isolates was observed, with considerable overlap between Blumeria, Magnaporthe and Puccinia. Small subsets of these general pathogen-regulated genes were identified as differentially regulated between host and corresponding nonhost interactions, indicating a fine-tuning of the general pathogen response during the course of co-evolution. Additionally, the host- or nonhost-related responses were rather specific for each pair of adapted and non-adapted isolates, indicating that the nonhost resistance-related responses were to a great extent pathosystem-specific. This pathosystem-specific reprogramming may reflect different resistance mechanisms operating against non-adapted pathogens with different lifestyles, or equally, different co-option of the hosts by the adapted isolates to create an optimal environment for infection. To compare the transcriptional reprogramming between wheat and barley, putative orthologues were identified. Within the wheat and barley general pathogen-regulated genes, temporal expression profiles of orthologues looked similar, indicating conserved general responses in Triticeae against fungal attack. However, the comparison of orthologues differentially expressed between host and nonhost interactions revealed fewer commonalities between wheat and barley, but rather suggested different host or nonhost responses in the two cereal species.

CONCLUSIONS

Taken together, our results suggest independent co-evolutionary forces acting on host pathosystems mirrored by barley- or wheat-specific nonhost responses. As a result of evolutionary processes, at least for the pathosystems investigated, NHR appears to rely on rather specific plant responses.

Altered Expression of Genes Implicated in Xylan Biosynthesis Affects Penetration Resistance against Powdery Mildew

Heteroxylan has recently been identified as an important component of papillae, which are formed during powdery mildew infection of barley leaves. Deposition of heteroxylan near the sites of attempted fungal penetration in the epidermal cell wall is believed to enhance the physical resistance to the fungal penetration peg and hence to improve pre-invasion resistance. Several glycosyltransferase (GT) families are implicated in the assembly of heteroxylan in the plant cell wall, and are likely to work together in a multi-enzyme complex. Members of key GT families reported to be involved in heteroxylan biosynthesis are up-regulated in the epidermal layer of barley leaves during powdery mildew infection. Modulation of their expression leads to altered susceptibility levels, suggesting that these genes are important for penetration resistance. The highest level of resistance was achieved when a GT43 gene was co-expressed with a GT47 candidate gene, both of which have been predicted to be involved in xylan backbone biosynthesis. Altering the expression level of several candidate heteroxylan synthesis genes can significantly alter disease susceptibility. This is predicted to occur through changes in the amount and structure of heteroxylan in barley papillae.

Assessment of Functional EST-SSR Markers (Sugarcane) in Cross-Species Transferability, Genetic Diversity among Poaceae Plants, and Bulk Segregation Analysis

Expressed sequence tags (ESTs) are important resource for gene discovery, gene expression and its regulation, molecular marker development, and comparative genomics. We procured 10000 ESTs and analyzed 267 EST-SSRs markers through computational approach. The average density was one SSR/10.45 kb or 6.4% frequency, wherein trinucleotide repeats (66.74%) were the most abundant followed by di- (26.10%), tetra- (4.67%), penta- (1.5%), and hexanucleotide (1.2%) repeats. Functional annotations were done and after-effect newly developed 63 EST-SSRs were used for cross transferability, genetic diversity, and bulk segregation analysis (BSA). Out of 63 EST-SSRs, 42 markers were identified owing to their expansion genetics across 20 different plants which amplified 519 alleles at 180 loci with an average of 2.88 alleles/locus and the polymorphic information content (PIC) ranged from 0.51 to 0.93 with an average of 0.83. The cross transferability ranged from 25% for wheat to 97.22% for Schlerostachya, with an average of 55.86%, and genetic relationships were established based on diversification among them. Moreover, 10 EST-SSRs were recognized as important markers between bulks of pooled DNA of sugarcane cultivars through BSA. This study highlights the employability of the markers in transferability, genetic diversity in grass species, and distinguished sugarcane bulks.

An LRR/Malectin Receptor-Like Kinase Mediates Resistance to Non-adapted and Adapted Powdery Mildew Fungi in Barley and Wheat

Pattern recognition receptors (PRRs) belonging to the multigene family of receptor-like kinases (RLKs) are the sensing devices of plants for microbe- or pathogen-associated molecular patterns released from microbial organisms. Here we describe Rnr8 (for Required for non-host resistance 8) encoding HvLEMK1, a LRR-malectin domain-containing transmembrane RLK that mediates non-host resistance of barley to the non-adapted wheat powdery mildew fungus Blumeria graminis f.sp. tritici. Transgenic barley lines with silenced HvLEMK1 allow entry and colony growth of the non-adapted pathogen, although sporulation was reduced and final colony size did not reach that of the adapted barley powdery mildew fungus B. graminis f.sp. hordei. Transient expression of the barley or wheat LEMK1 genes enhanced resistance in wheat to the adapted wheat powdery mildew fungus while expression of the same genes did not protect barley from attack by the barley powdery mildew fungus. The results suggest that HvLEMK1 is a factor mediating non-host resistance in barley and quantitative host resistance in wheat to the wheat powdery mildew fungus.

Discovery of genes affecting resistance of barley to adapted and non-adapted powdery mildew fungi

Background

Non-host resistance, NHR, to non-adapted pathogens and quantitative host resistance, QR, confer durable protection to plants and are important for securing yield in a longer perspective. However, a more targeted exploitation of the trait usually possessing a complex mode of inheritance by many quantitative trait loci, QTLs, will require a better understanding of the most important genes and alleles.

Results

Here we present results from a transient-induced gene silencing, TIGS, approach of candidate genes for NHR and QR in barley against the powdery mildew fungus Blumeria graminis. Genes were selected based on transcript regulation, multigene-family membership or genetic map position. Out of 1,144 tested RNAi-target genes, 96 significantly affected resistance to the non-adapted wheat- or the compatible barley powdery mildew fungus, with an overlap of four genes. TIGS results for QR were combined with transcript regulation data, allele-trait associations, QTL co-localization and copy number variation resulting in a meta-dataset of 51 strong candidate genes with convergent evidence for a role in QR.

Conclusions

This study represents an initial, functional inventory of approximately 3% of the barley transcriptome for a role in NHR or QR against the powdery mildew pathogen. The discovered candidate genes support the idea that QR in this Triticeae host is primarily based on pathogen-associated molecular pattern-triggered immunity, which is compromised by effector molecules produced by the compatible pathogen. The overlap of four genes with significant TIGS effects both in the NHR and QR screens also indicates shared components for both forms of durable pathogen resistance.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-014-0518-8) contains supplementary material, which is available to authorized users.

Regulation of flavonol content and composition in (Syrah×Pinot Noir) mature grapes: integration of transcriptional profiling and metabolic quantitative trait locus analyses

Flavonols are a ubiquitous class of flavonoids that accumulate preferentially in flowers and mature berries. Besides their photo-protective function, they play a fundamental role during winemaking, stabilizing the colour by co-pigmentation with anthocyanins and contributing to organoleptic characteristics. Although the general flavonol pathway has been genetically and biochemically elucidated, the genetic control of flavonol content and composition at harvest is still not clear. To this purpose, the grapes of 170 segregating F1 individuals from a 'Syrah'×'Pinot Noir' population were evaluated at the mature stage for the content of six flavonol aglycons in four seasons. Metabolic data in combination with genetic data enabled the identification of 16 mQTLs (metabolic quantitative trait loci). For the first time, major genetic control by the linkage group 2 (LG 2)/MYBA region on flavonol variation, in particular of tri-hydroxylated flavonols, is demonstrated. Moreover, seven regions specifically associated with the fine control of flavonol biosynthesis are identified. Gene expression profiling of two groups of individuals significantly divergent for their skin flavonol content identified a large set of differentially modulated transcripts. Among these, the transcripts coding for MYB and bZIP transcription factors, methyltranferases, and glucosyltranferases specific for flavonols, proteins, and factors belonging to the UV-B signalling pathway and co-localizing with the QTL regions are proposed as candidate genes for the fine regulation of flavonol content and composition in mature grapes.

The quest for tolerant varieties: the importance of integrating "omics" techniques to phenotyping

The primary objective of crop breeding is to improve yield and/or harvest quality while minimizing inputs. Global climate change and the increase in world population are significant challenges for agriculture and call for further improvements to crops and the development of new tools for research. Significant progress has been made in the molecular and genetic analysis of model plants. However, is science generating false expectations? Are 'omic techniques generating valuable information that can be translated into the field? The exploration of crop biodiversity and the correlation of cellular responses to stress tolerance at the plant level is currently a challenge. This viewpoint reviews concisely the problems one encounters when working on a crop and provides an outline of possible workflows when initiating cellular phenotyping via "-omic" techniques (transcriptomics, proteomics, metabolomics).

Identification and characterization of microsatellites in expressed sequence tags and their cross transferability in different plants

Expressed sequence tags (EST) are potential source for the development of genic microsatellite markers, gene discovery, comparative genomics, and other genomic studies. In the present study, 7630 ESTs were examined from NCBI for SSR identification and characterization. A total of 263 SSRs were identified with an average density of one SSR/4.2 kb (3.4% frequency). Analysis revealed that trinucleotide repeats (47.52%) were most abundant followed by tetranucleotide (19.77%), dinucleotide (19.01%), pentanucleotide (9.12%), and hexanucleotide repeats (4.56%). Functional annotation was done through homology search and gene ontology, and 35 EST-SSRs were selected. Primer pairs were designed for evaluation of cross transferability and polymorphism among 11 plants belonging to five different families. Total 402 alleles were generated at 155 loci with an average of 2.6 alleles/locus and the polymorphic information content (PIC) ranged from 0.15 to 0.92 with an average of 0.75. The cross transferability ranged from 34.84% to 98.06% in different plants, with an average of 67.86%. Thus, the validation study of annotated 35 EST-SSR markers which correspond to particular metabolic activity revealed polymorphism and evolutionary nature in different families of Angiospermic plants.

A first step toward the development of a barley NAM population and its utilization to detect QTLs conferring leaf rust seedling resistance

We suggest multi-parental nested association mapping as a valuable innovation in barley genetics, which increases the power to map quantitative trait loci and assists in extending genetic diversity of the elite barley gene pool. Plant genetic resources are a key asset to further improve crop species. The nested association mapping (NAM) approach was introduced to identify favorable genes in multi-parental populations. Here, we report toward the development of the first explorative barley NAM population and demonstrate its usefulness in a study on mapping quantitative trait loci (QTLs) for leaf rust resistance. The NAM population HEB-5 was developed from crossing and backcrossing five exotic barley donors with the elite barley cultivar 'Barke,' resulting in 295 NAM lines in generation BC1S1. HEB-5 was genetically characterized with 1,536 barley SNPs. Across HEB-5 and within the NAM families, no deviation from the expected genotype and allele frequencies was detected. Genetic similarity between 'Barke' and the NAM families ranged from 78.6 to 83.1 %, confirming the backcrossing step during population development. To explore its usefulness, a screen for leaf rust (Puccinia hordei) seedling resistance was conducted. Resistance QTLs were mapped to six barley chromosomes, applying a mixed model genome-wide association study. In total, four leaf rust QTLs were detected across HEB-5 and four QTLs within family HEB-F23. Favorable exotic QTL alleles reduced leaf rust symptoms on two chromosomes by 33.3 and 36.2 %, respectively. The located QTLs may represent new resistance loci or correspond to new alleles of known resistance genes. We conclude that the exploratory population HEB-5 can be applied to mapping and utilizing exotic QTL alleles of agronomic importance. The NAM concept will foster the evaluation of the genetic diversity, which is present in our primary barley gene pool.

A reference consensus genetic map for molecular markers and economically important traits in faba bean (Vicia faba L.)

BACKGROUND

Faba bean (Vicia faba L.) is among the earliest domesticated crops from the Near East. Today this legume is a key protein feed and food worldwide and continues to serve an important role in culinary traditions throughout Middle East, Mediterranean region, China and Ethiopia. Adapted to a wide range of soil types, the main faba bean breeding objectives are to improve yield, resistance to biotic and abiotic stresses, seed quality and other agronomic traits. Genomic approaches aimed at enhancing faba bean breeding programs require high-quality genetic linkage maps to facilitate quantitative trait locus analysis and gene tagging for use in a marker-assisted selection. The objective of this study was to construct a reference consensus map in faba bean by joining the information from the most relevant maps reported so far in this crop.

RESULTS

A combination of two approaches, increasing the number of anchor loci in diverse mapping populations and joining the corresponding genetic maps, was used to develop a reference consensus map in faba bean. The map was constructed from three main recombinant inbreed populations derived from four parental lines, incorporates 729 markers and is based on 69 common loci. It spans 4,602 cM with a range from 323 to 1041 loci in six main linkage groups or chromosomes, and an average marker density of one locus every 6 cM. Locus order is generally well maintained between the consensus map and the individual maps.

CONCLUSION

We have constructed a reliable and fairly dense consensus genetic linkage map that will serve as a basis for genomic approaches in faba bean research and breeding. The core map contains a larger number of markers than any previous individual map, covers existing gaps and achieves a wider coverage of the large faba bean genome as a whole. This tool can be used as a reference resource for studies in different genetic backgrounds, and provides a framework for transferring genetic information when using different marker technologies. Combined with syntenic approaches, the consensus map will increase marker density in selected genomic regions and will be useful for future faba bean molecular breeding applications.

Design of new genome- and gene-sourced primers and identification of QTL for seed oil content in a specially high-oil Brassica napus cultivar

Rapeseed (Brassica napus L.) is one of most important oilseed crops in the world. There are now various rapeseed cultivars in nature that differ in their seed oil content because they vary in oil-content alleles and there are high-oil alleles among the high-oil rapeseed cultivars. For these experiments, we generated doubled haploid (DH) lines derived from the cross between the specially high-oil cultivar zy036 whose seed oil content is approximately 50% and the specially low-oil cultivar 51070 whose seed oil content is approximately 36%. First, to address the deficiency in polymorphic markers, we designed 5944 pairs of newly developed genome-sourced primers and 443 pairs of newly developed primers related to oil-content genes to complement the 2244 pairs of publicly available primers. Second, we constructed a new DH genetic linkage map using 527 molecular markers, consisting of 181 publicly available markers, 298 newly developed genome-sourced markers and 48 newly developed markers related to oil-content genes. The map contained 19 linkage groups, covering a total length of 2,265.54 cM with an average distance between markers of 4.30 cM. Third, we identified quantitative trait loci (QTL) for seed oil content using field data collected at three sites over 3 years, and found a total of 12 QTL. Of the 12 QTL associated with seed oil content identified, 9 were high-oil QTL which derived from the specially high-oil cultivar zy036. Two high-oil QTL on chromosomes A2 and C9 co-localized in two out of three trials. By QTL mapping for seed oil content, we found four candidate genes for seed oil content related to four gene markers: GSNP39, GSSR161, GIFLP106 and GIFLP046. This information will be useful for cloning functional genes correlated with seed oil content in the future.

Allele mining in barley genetic resources reveals genes of race-non-specific powdery mildew resistance

Race-non-specific, or quantitative, pathogen resistance is of high importance to plant breeders due to its expected durability. However, it is usually controlled by multiple quantitative trait loci (QTL) and therefore difficult to handle in practice. Knowing the genes that underlie race-non-specific resistance (NR) would allow its exploitation in a more targeted manner. Here, we performed an association-genetic study in a customized worldwide collection of spring barley accessions for candidate genes of race-NR to the powdery mildew fungus Blumeria graminis f. sp. hordei (Bgh) and combined data with results from QTL mapping as well as functional-genomics approaches. This led to the identification of 11 associated genes with converging evidence for an important role in race-NR in the presence of the Mlo gene for basal susceptibility. Outstanding in this respect was the gene encoding the transcription factor WRKY2. The results suggest that unlocking plant genetic resources and integrating functional-genomic with genetic approaches can accelerate the discovery of genes underlying race-NR in barley and other crop plants.