Divergent and convergent modes of interaction between wheat and Puccinia graminis f. sp tritici isolates revealed by the comparative gene co-expression network and genome analyses

Citation: Rutter, W. B., Salcedo, A., Akhunova, A., He, F., Wang, S. C., Liang, H. Q., . . . Akhunov, E. (2017). Divergent and convergent modes of interaction between wheat and Puccinia graminis f. sp tritici isolates revealed by the comparative gene co-expression network and genome analyses. Bmc Genomics, 18, 20. doi:10.1186/s12864-017-3678-6Background: Two opposing evolutionary constraints exert pressure on plant pathogens: one to diversify virulence factors in order to evade plant defenses, and the other to retain virulence factors critical for maintaining a compatible interaction with the plant host. To better understand how the diversified arsenals of fungal genes promote interaction with the same compatible wheat line, we performed a comparative genomic analysis of two North American isolates of Puccinia graminis f. sp. tritici (Pgt). Results: The patterns of inter-isolate divergence in the secreted candidate effector genes were compared with the levels of conservation and divergence of plant-pathogen gene co-expression networks (GCN) developed for each isolate. Comprative genomic analyses revealed substantial level of interisolate divergence in effector gene complement and sequence divergence. Gene Ontology (GO) analyses of the conserved and unique parts of the isolate-specific GCNs identified a number of conserved host pathways targeted by both isolates. Interestingly, the degree of inter-isolate sub-network conservation varied widely for the different host pathways and was positively associated with the proportion of conserved effector candidates associated with each sub- network. While different Pgt isolates tended to exploit similar wheat pathways for infection, the mode of plant-pathogen interaction varied for different pathways with some pathways being associated with the conserved set of effectors and others being linked with the diverged or isolate-specific effectors. Conclusions: Our data suggest that at the intra-species level pathogen populations likely maintain divergent sets of effectors capable of targeting the same plant host pathways. This functional redundancy may play an important role in the dynamic of the "arms-race" between host and pathogen serving as the basis for diverse virulence strategies and creating conditions where mutations in certain effector groups will not have a major effect on the pathogen's ability to infect the host

A haplotype map of allohexaploid wheat reveals distinct patterns of selection on homoeologous genomes

Citation: Jordan, K. W., Wang, S., Lun, Y., Gardiner, L. J., MacLachlan, R., Hucl, P., . . . Akhunov, E. (2015). A haplotype map of allohexaploid wheat reveals distinct patterns of selection on homoeologous genomes. Genome Biology, 16(1). https://doi.org/10.1186/s13059-015-0606-4Background: Bread wheat is an allopolyploid species with a large, highly repetitive genome. To investigate the impact of selection on variants distributed among homoeologous wheat genomes and to build a foundation for understanding genotype-phenotype relationships, we performed population-scale re-sequencing of a diverse panel of wheat lines. Results: A sample of 62 diverse lines was re-sequenced using the whole exome capture and genotyping-by-sequencing approaches. We describe the allele frequency, functional significance, and chromosomal distribution of 1.57 million single nucleotide polymorphisms and 161,719 small indels. Our results suggest that duplicated homoeologous genes are under purifying selection. We find contrasting patterns of variation and inter-variant associations among wheat genomes; this, in addition to demographic factors, could be explained by differences in the effect of directional selection on duplicated homoeologs. Only a small fraction of the homoeologous regions harboring selected variants overlapped among the wheat genomes in any given wheat line. These selected regions are enriched for loci associated with agronomic traits detected in genome-wide association studies. Conclusions: Evidence suggests that directional selection in allopolyploids rarely acted on multiple parallel advantageous mutations across homoeologous regions, likely indicating that a fitness benefit could be obtained by a mutation at any one of the homoeologs. Additional advantageous variants in other homoelogs probably either contributed little benefit, or were unavailable in populations subjected to directional selection. We hypothesize that allopolyploidy may have increased the likelihood of beneficial allele recovery by broadening the set of possible selection targets. © 2015 Jordan et al.; licensee BioMed Central.Additional Authors: Talbert, L.;Bansal, U. K.;Bariana, H. S.;Hayden, M. J.;Pozniak, C.;Jeddeloh, J. A.;Hall, A

QTL mapping of yield components and kernel traits in wheat cultivars TAM 112 and Duster

In the Southern Great Plains, wheat cultivars have been selected for a combination of outstanding yield and drought tolerance as a long-term breeding goal. To understand the underlying genetic mechanisms, this study aimed to dissect the quantitative trait loci (QTL) associated with yield components and kernel traits in two wheat cultivars `TAM 112' and `Duster' under both irrigated and dryland environments. A set of 182 recombined inbred lines (RIL) derived from the cross of TAM 112/Duster were planted in 13 diverse environments for evaluation of 18 yield and kernel related traits. High-density genetic linkage map was constructed using 5,081 single nucleotide polymorphisms (SNPs) from genotyping-by-sequencing (GBS). QTL mapping analysis detected 134 QTL regions on all 21 wheat chromosomes, including 30 pleiotropic QTL regions and 21 consistent QTL regions, with 10 QTL regions in common. Three major pleiotropic QTL on the short arms of chromosomes 2B (57.5 - 61.6 Mbps), 2D (37.1 - 38.7 Mbps), and 7D (66.0 - 69.2 Mbps) colocalized with genes Ppd-B1, Ppd-D1, and FT-D1, respectively. And four consistent QTL associated with kernel length (KLEN), thousand kernel weight (TKW), plot grain yield (YLD), and kernel spike-1 (KPS) (Qklen.tamu.1A.325, Qtkw.tamu.2B.137, Qyld.tamu.2D.3, and Qkps.tamu.6A.113) explained more than 5% of the phenotypic variation. QTL Qklen.tamu.1A.325 is a novel QTL with consistent effects under all tested environments. Marker haplotype analysis indicated the QTL combinations significantly increased yield and kernel traits. QTL and the linked markers identified in this study will facilitate future marker-assisted selection (MAS) for pyramiding the favorable alleles and QTL map-based cloning.Horticulture and Landscape Architectur

Molecular Dynamics Analysis Reveals Structural Insights into Mechanism of Nicotine N-Demethylation Catalyzed by Tobacco Cytochrome P450 Mono-Oxygenase

CYP82E4, a cytochrome P450 monooxygenase, has nicotine N-demethylase (NND) activity, which mediates the bioconversion of nicotine into nornicotine in senescing tobacco leaves. Nornicotine is a precursor of the carcinogen, tobacco-specific nitrosamine. CYP82E3 is an ortholog of CYP82E4 with 95% sequence identity, but it lacks NND activity. A recent site-directed mutagenesis study revealed that a single amino acid substitution, i.e., cysteine to tryptophan at the 330 position in the middle of protein, restores the NND activity of CYP82E3 entirely. However, the same amino acid change caused the loss of the NND activity of CYP82E4. To determine the mechanism of the functional turnover of the two molecules, four 3D structures, i.e., the two molecules and their corresponding cys–trp mutants were modeled. The resulting structures exhibited that the mutation site is far from the active site, which suggests that no direct interaction occurs between the two sites. Simulation studies in different biological scenarios revealed that the mutation introduces a conformation drift with the largest change at the F-G loop. The dynamics trajectories analysis using principal component analysis and covariance analysis suggests that the single amino acid change causes the opening and closing of the transfer channels of the substrates, products, and water by altering the motion of the F-G and B-C loops. The motion of helix I is also correlated with the motion of both the F-G loop and the B-C loop and; the single amino acid mutation resulted in the curvature of helix I. These results suggest that the single amino acid mutation outside the active site region may have indirectly mediated the flexibility of the F-G and B-C loops through helix I, causing a functional turnover of the P450 monooxygenase

Single-nucleotide polymorphism discovery by high-throughput sequencing in sorghum

<p>Abstract</p> <p>Background</p> <p>Eight diverse sorghum (<it>Sorghum bicolor </it>L. Moench) accessions were subjected to short-read genome sequencing to characterize the distribution of single-nucleotide polymorphisms (SNPs). Two strategies were used for DNA library preparation. Missing SNP genotype data were imputed by local haplotype comparison. The effect of library type and genomic diversity on SNP discovery and imputation are evaluated.</p> <p>Results</p> <p>Alignment of eight genome equivalents (6 Gb) to the public reference genome revealed 283,000 SNPs at ≥82% confirmation probability. Sequencing from libraries constructed to limit sequencing to start at defined restriction sites led to genotyping 10-fold more SNPs in all 8 accessions, and correctly imputing 11% more missing data, than from semirandom libraries. The SNP yield advantage of the reduced-representation method was less than expected, since up to one fifth of reads started at noncanonical restriction sites and up to one third of restriction sites predicted <it>in silico </it>to yield unique alignments were not sampled at near-saturation. For imputation accuracy, the availability of a genomically similar accession in the germplasm panel was more important than panel size or sequencing coverage.</p> <p>Conclusions</p> <p>A sequence quantity of 3 million 50-base reads per accession using a <it>Bsr</it>FI library would conservatively provide satisfactory genotyping of 96,000 sorghum SNPs. For most reliable SNP-genotype imputation in shallowly sequenced genomes, germplasm panels should consist of pairs or groups of genomically similar entries. These results may help in designing strategies for economical genotyping-by-sequencing of large numbers of plant accessions.</p

A haplotype map of allohexaploid wheat reveals distinct patterns of selection on homoeologous genomes

BACKGROUND: Bread wheat is an allopolyploid species with a large, highly repetitive genome. To investigate the impact of selection on variants distributed among homoeologous wheat genomes and to build a foundation for understanding genotype-phenotype relationships, we performed population-scale re-sequencing of a diverse panel of wheat lines. RESULTS: A sample of 62 diverse lines was re-sequenced using the whole exome capture and genotyping-by-sequencing approaches. We describe the allele frequency, functional significance, and chromosomal distribution of 1.57 million single nucleotide polymorphisms and 161,719 small indels. Our results suggest that duplicated homoeologous genes are under purifying selection. We find contrasting patterns of variation and inter-variant associations among wheat genomes; this, in addition to demographic factors, could be explained by differences in the effect of directional selection on duplicated homoeologs. Only a small fraction of the homoeologous regions harboring selected variants overlapped among the wheat genomes in any given wheat line. These selected regions are enriched for loci associated with agronomic traits detected in genome-wide association studies. CONCLUSIONS: Evidence suggests that directional selection in allopolyploids rarely acted on multiple parallel advantageous mutations across homoeologous regions, likely indicating that a fitness benefit could be obtained by a mutation at any one of the homoeologs. Additional advantageous variants in other homoelogs probably either contributed little benefit, or were unavailable in populations subjected to directional selection. We hypothesize that allopolyploidy may have increased the likelihood of beneficial allele recovery by broadening the set of possible selection targets

Mapping QTL for Yield and Its Component Traits Using Wheat (<i>Triticum aestivum</i> L.) RIL Mapping Population from TAM 113 × Gallagher

Understanding genetic architectures of yield and yield-related traits of wheat (Triticum aestivum L.) grown under dryland or irrigated conditions is pivotal for developing modern high-yielding germplasm and cultivars. The objectives of this study were to detect quantitative trait loci (QTL) linked to yield and yield components using a mapping population derived from ‘TAM 113’/‘Gallagher’, including 191 recombinant inbred lines (RILs). The population was grown in McGregor, College Station, and Bushland, Texas, for three consecutive years from 2019 to 2021. A high-density genetic map covering all 21 chromosomes was constructed using a set of 8,075 single nucleotide polymorphisms (SNPs) obtained using genotyping-by-sequencing (GBS). A total of 147 QTLs for 16 yield-related traits were identified, which included 16 QTLs consistently detected in multiple experiments and 8 QTLs that showed pleiotropic effects. Of them, five pleiotropic QTLs overlapped with the consistent QTL. They increased grain yield (YLD) up to 37.64 g m−2, thousand kernel weight (TKW) up to 1.33 g, harvest (HI) up to 0.97%, kernel length up to 0.08 mm, and kernel width up to 0.04 mm with Gallagher alleles and increased YLD up to 22.21 g m−2, kernels spike−1 up to 1.77, TKW up to 1.14 g, and HI up to 3.72% with TAM 113 alleles. One major and consistent QTL on chromosome 2D at 34.4 Mbp overlapped with the major photoperiod gene Ppd-D1 and was affected by multiple traits, including kernel diameter (DIAM), TKW, kernel hardness index (KHI), heading date (HD), and plant height (PH). Another QTL cluster region on 7D between 52 and 66 Mbp, encompassing one consistent and three pleiotropic QTLs. One of the pleiotropic QTLs at 52 Mbp increased YLD up to 24.16 g m−2, HI up to 1%, and DIAM up to 0.03 mm. This study dissected genetic loci associated with yield and yield-related traits, providing valuable information on wheat improvement using marker-assisted selection (MAS)