Springer Science and Business Media LLC

Meiotic maps are a key tool for comparative genomics and association mapping studies. Next-generation sequencing and genotyping by sequencing are speeding the processes of SNP discovery and the development of new genetic tools, including meiotic maps for numerous species. Currently there are limited genetic resources for sockeye salmon, Oncorhynchus nerka. We develop the first dense meiotic map for sockeye salmon using a combination of novel SNPs found in restriction site associated DNA (RAD tags) and SNPs available from existing expressed sequence tag (EST) based assays. We discovered and genotyped putative SNPs in 3,430 RAD tags. We removed paralogous sequence variants leaving 1,672 SNPs; these were combined with 53 EST-based SNP genotypes for linkage mapping. The map contained 29 male and female linkage groups, consistent with the haploid chromosome number expected for sockeye salmon. The female map contains 1,057 loci spanning 4,896 cM, and the male map contains 1,118 loci spanning 4,220 cM. Regions of conservation with rainbow trout and synteny between the RAD based rainbow trout map and the sockeye salmon map were established. Using RAD sequencing and EST-based SNP assays we successfully generated the first high density linkage map for sockeye salmon.

Hessian fly (Mayetiola destructor) is one of the most destructive pests of wheat. The genes encoding 12-oxo-phytodienoic acid reductase (OPR) and lipoxygenase (LOX) play critical roles in insect resistance pathways in higher plants, but little is known about genes controlling resistance to Hessian fly in wheat. In this study, 154 F6:8 recombinant inbred lines (RILs) generated from a cross between two cultivars, ‘Jagger’ and ‘2174’ of hexaploid wheat (2n = 6 × =42; AABBDD), were used to map genes associated with resistance to Hessian fly. Two QTLs were identified. The first one was a major QTL on chromosome 1A (QHf.osu-1A), which explained 70% of the total phenotypic variation. The resistant allele at this locus in cultivar 2174 could be orthologous to one or more of the previously mapped resistance genes (H9, H10, H11, H16, and H17) in tetraploid wheat. The second QTL was a minor QTL on chromosome 2A (QHf.osu-2A), which accounted for 18% of the total phenotypic variation. The resistant allele at this locus in 2174 is collinear to an Yr17-containing-fragment translocated from chromosome 2N of Triticum ventricosum (2n = 4 × =28; DDNN) in Jagger. Genetic mapping results showed that two OPR genes, TaOPR1-A and TaOPR2-A, were tightly associated with QHf.osu-1A and QHf.osu-2A, respectively. Another OPR gene and three LOX genes were mapped but not associated with Hessian fly resistance in the segregating population. This study has located two major QTLs/genes in bread wheat that can be directly used in wheat breeding programs and has also provided insights for the genetic association and disassociation of Hessian fly resistance with OPR and LOX genes in wheat.

Resistance inducers have been used in annual crops as an alternative for disease control. Wood perennial fruit trees, such as those of the citrus species, are candidates for treatment with resistance inducers, such as salicylic acid (SA) and chitosan (CHI). However, the involved mechanisms in resistance induced by elicitors in citrus are currently few known. In the present manuscript, we report information regarding the transcriptional changes observed in sweet orange in response to exogenous applications of SA and CHI using RNA-seq technology. More genes were induced by SA treatment than by CHI treatment. In total, 1,425 differentially expressed genes (DEGs) were identified following treatment with SA, including the important genes WRKY50, PR2, and PR9, which are known to participate in the salicylic acid signaling pathway, and genes involved in ethylene/Jasmonic acid biosynthesis (ACS12, AP2 domain-containing transcription factor, and OPR3). In addition, SA treatment promoted the induction of a subset of genes involved in several metabolic processes, such as redox states and secondary metabolism, which are associated with biotic stress. For CHI treatment, there were 640 DEGs, many of them involved in secondary metabolism. For both SA and CHI treatments, the auxin pathway genes were repressed, but SA treatment promoted induction in the ethylene and jasmonate acid pathway genes, in addition to repressing the abscisic acid pathway genes. Chitosan treatment altered some hormone metabolism pathways. The DEGs were validated by quantitative Real-Time PCR (qRT-PCR), and the results were consistent with the RNA-seq data, with a high correlation between the two analyses. We expanded the available information regarding induced defense by elicitors in a species of Citrus that is susceptible to various diseases and identified the molecular mechanisms by which this defense might be mediated.

Basic helix-loop-helix (bHLH) proteins are transcription factors (TFs) that have been shown to regulate anthocyanin biosynthesis in many plant species. However, the bHLH gene family in Populus deltoids has not yet been reported. In this study, 185 PdbHLH genes were identified in the Populus deltoids genome and were classified into 15 groups based on their sequence similarity and phylogenetic relationships. Analysis of the gene structure, chromosome location and conserved motif of the PdbHLH genes were performed by bioinformatic methods. Gene duplication analyses revealed that 114 PdbHLH were expanded and retained after WGD/segmental and proximal duplication. Investigation of cis-regulatory elements of PdbHLH genes indicated that many PdbHLH genes are involved in the regulation of anthocyanin biosynthesis. The expression patterns of PdbHLHs were obtained from previous data in two colored-leaf poplar (QHP and JHP) and green leaf poplar (L2025). Further analysis revealed that 12 candidate genes, including 3 genes (PdbHLH57, PdbHLH143, and PdbHLH173) from the subgroup III(f) and 9 gene from other groups, were positively associated with anthocyanin biosynthesis. In addition, 4 genes (PdbHLH4, PdbHLH1, PdbHLH18, and PdbHLH164) may be involved in negatively regulating the anthocyanin biosynthesis. These results provide a basis for the functional characterization of bHLH genes and investigations on the molecular mechanisms of anthocyanin biosynthesis in colored-leaf poplar.

Retinal neovascularization (RNV) is a leading cause of blindness worldwide. Long non-coding RNA (lncRNA) and competing endogenous RNA (ceRNA) regulatory networks play vital roles in angiogenesis. The RNA-binding protein galectin-1 (Gal-1) participates in pathological RNV in oxygen-induced retinopathy mouse models. However, the molecular associations between Gal-1 and lncRNAs remain unclear. Herein, we aimed to explore the potential mechanism of action of Gal-1 as an RNA-binding protein. A comprehensive network of Gal-1, ceRNAs, and neovascularization-related genes was constructed based on transcriptome chip data and bioinformatics analysis of human retinal microvascular endothelial cells (HRMECs). We also conducted functional enrichment and pathway enrichment analyses. Fourteen lncRNAs, twenty-nine miRNAs, and eleven differentially expressed angiogenic genes were included in the Gal-1/ceRNA network. Additionally, the expression of six lncRNAs and eleven differentially expressed angiogenic genes were validated by qPCR in HRMECs with or without siLGALS1. Several hub genes, such as NRIR, ZFPM2-AS1, LINC0121, apelin, claudin-5, and C-X-C motif chemokine ligand 10, were found to potentially interact with Gal-1 via the ceRNA axis. Furthermore, Gal-1 may be involved in regulating biological processes related to chemotaxis, chemokine-mediated signaling, the immune response, and the inflammatory response. The Gal-1/ceRNA axis identified in this study may play a vital role in RNV. This study provides a foundation for the continued exploration of therapeutic targets and biomarkers associated with RNV.

Alternative polyadenylation (APA) has emerged as a pervasive mechanism that contributes to the transcriptome complexity and dynamics of gene regulation. The current tsunami of whole genome poly(A) site data from various conditions generated by 3′ end sequencing provides a valuable data source for the study of APA-related gene expression. Cluster analysis is a powerful technique for investigating the association structure among genes, however, conventional gene clustering methods are not suitable for APA-related data as they fail to consider the information of poly(A) sites (e.g., location, abundance, number, etc.) within each gene or measure the association among poly(A) sites between two genes. Here we proposed a computational framework, named PASCCA, for clustering genes from replicated or unreplicated poly(A) site data using canonical correlation analysis (CCA). PASCCA incorporates multiple layers of gene expression data from both the poly(A) site level and gene level and takes into account the number of replicates and the variability within each experimental group. Moreover, PASCCA characterizes poly(A) sites in various ways including the abundance and relative usage, which can exploit the advantages of 3′ end deep sequencing in quantifying APA sites. Using both real and synthetic poly(A) site data sets, the cluster analysis demonstrates that PASCCA outperforms other widely-used distance measures under five performance metrics including connectivity, the Dunn index, average distance, average distance between means, and the biological homogeneity index. We also used PASCCA to infer APA-specific gene modules from recently published poly(A) site data of rice and discovered some distinct functional gene modules. We have made PASCCA an easy-to-use R package for APA-related gene expression analyses, including the characterization of poly(A) sites, quantification of association between genes, and clustering of genes. By providing a better treatment of the noise inherent in repeated measurements and taking into account multiple layers of poly(A) site data, PASCCA could be a general tool for clustering and analyzing APA-specific gene expression data. PASCCA could be used to elucidate the dynamic interplay of genes and their APA sites among various biological conditions from emerging 3′ end sequencing data to address the complex biological phenomenon.

With the emphasis on analysing genotype-by-environment interactions within the framework of genomic selection and genome-wide association analysis, there is an increasing demand for reliable tools that can be used to simulate large-scale genomic data in order to assess related approaches. We proposed a theory to simulate large-scale genomic data on genotype-by-environment interactions and added this new function to our developed tool GPOPSIM. Additionally, a simulated threshold trait with large-scale genomic data was also added. The validation of the simulated data indicated that GPOSPIM2.0 is an efficient tool for mimicking the phenotypic data of quantitative traits, threshold traits, and genetically correlated traits with large-scale genomic data while taking genotype-by-environment interactions into account. This tool is useful for assessing genotype-by-environment interactions and threshold traits methods.

Cold hardiness is an important agronomic trait and can significantly affect grape production and quality. Until now, there are no reports focusing on cold hardiness quantitative trait loci (QTL) mapping. In this study, grapevine interspecific hybridisation was carried out with the maternal parent ‘Cabernet sauvignon’ and paternal parent ‘Zuoyouhong’. A total of 181 hybrid offspring and their parents were used as samples for restriction-site associated DNA sequencing (RAD). Grapevine cane phloem and xylem cold hardiness of the experimental material was detected using the low-temperature exotherm method in 2016, 2017 and 2018. QTL mapping was then conducted based on the integrated map. We constructed a high-density genetic linkage map with 16,076, 11,643, and 25,917 single-nucleotide polymorphism (SNP) markers anchored in the maternal, paternal, and integrated maps, respectively. The average genetic distances of adjacent markers in the maps were 0.65 cM, 0.77 cM, and 0.41 cM, respectively. Colinearity analysis was conducted by comparison with the grape reference genome and showed good performance. Six QTLs were identified based on the phenotypic data of 3 years and they were mapped on linkage group (LG) 2, LG3, and LG15. Based on QTL results, candidate genes which may be involved in grapevine cold hardiness were selected. High-density linkage maps can facilitate grapevine fine QTL mapping, genome comparison, and sequence assembly. The cold hardiness QTL mapping and candidate gene discovery performed in this study provide an important reference for molecular-assisted selection in grapevine cold hardiness breeding.