Four loci differentially expressed in muscle tissue depending on water-holding capacity are associated with meat quality in commercial pig herds

Four genes, VTN, KERA, LYZ, and a non-annotated EST (Affymetrix probe set ID: Ssc.25503.1.S1_at), whose candidacy for traits related to water-holding capacity of meat arises from their trait-dependent differential expression, were selected for candidate gene analysis. Based on in silico analysis SNPs were detected, confirmed by sequencing and used to genotype animals of 4 pig populations including 3 commercial herds of Pietrain (PI), Pietrain × (German Large White × German Landrace) (PIF1), German Landrace (DL) and 1 experimental F₂ population Duroc × Pietrain (DUPI). Comparative and genetic mapping established the location of VTN on SSC12, of LYZ and KERA on SSC5 and of UN on SSC7, coinciding with QTL regions for meat quality traits. VTN showed association with pH1, pH24 and drip loss. LYZ revealed association with conductivity 24, pH1 and drip loss. KERA was associated with pH. UN showed association with pH24 and drip loss, respectively. However, none of the candidate genes showed significant associations for a particular trait across all populations. This may be due to breed specific effects that are related to the differences in meat quality of theses pig breeds. The studies revealed statistic evidence for a link of genetic variation at these loci or close to them and promoted those four candidate genes as functional and/or positional candidate genes for meat quality traits.     

High‐density linkage map reveals QTL underlying growth traits in AP13×VS16 biparental population of switchgrass

Switchgrass (Panicum virgatum L.), a native warm‐season perennial grass, is being considered as a feedstock for biofuel production in the United States. To expedite its genetic improvement and enhance genetic gain per selection cycle, application of marker‐assisted selection is indispensable. A high‐density linkage map was constructed in a pseudo‐F₁ testcross mapping population of AP13×VS16, consisting of 349 progenies. A total of 8,757 single nucleotide polymorphism (SNP) markers generated through genotype‐by‐sequencing (GBS) were used to construct the linkage map. The total map length spans up to 2,540.2 cM with the marker density of one marker in every 0.25–0.34 cM. Spring green‐up (SG), days to flowering (FL), and the vegetative growth period (VP) data were analyzed and used for quantitative trait loci (QTL) mapping. The population showed significant variations and exhibited transgressive segregation for SG, FL, and VP. QTL analyses were performed using trait mean of each year and location along with BLUP (best linear unbiased prediction) values of the traits. A total of 35, 37, and 34 QTL for SG, FL, and VP, respectively, were identified. Phenotypic variability explained by each QTL ranged from 11.29% to 27.85%. The additive genetic effects of individual QTL ranged from ?1.81 to 2.40, ?6.12 to 7.58, and ?16.01 to 6.38 for SG, FL, and VP, respectively. Comparing major QTL regions in the switchgrass genome, 20 candidate genes were identified which were reported to be involved in growth‐, development‐, and flowering‐related traits in switchgrass.     

Differential expression of genes in milk of dairy cattle during lactation

The milk fat globule (MFG) is one of the most representative of mammary gland tissues and can be utilized to study gene expression of lactating cows during lactation. In this study, RNA‐seq technology was employed to detect differential expression of genes in MFGs at day 10 and day 70 after calving between two groups of cows with extremely high (H group) and low (L group) 305‐day milk yield, milk fat yield and milk protein yield. In total, 1232, 81, 429 and 178 significantly differentially expressed genes (false discovery rate q < 0.05) were detected between H10 and L10, H70 and L70, H10 and H70, and L10 and L70 respectively. Gene Ontology enrichment and pathway analysis revealed that these differentially expressed genes were enriched in biological processes involved in mammary gland development, protein and lipid metabolism process, signal transduction, cellular process, differentiation and immune function. Among these differentially expressed genes, 178 (H10 vs. L10), 4 (H70 vs. L70), 68 (H10 vs. H70) and 22 (L10 vs. L70) were found to be located within previously reported QTL regions for milk production traits. Based on these results, some promising candidate genes for milk production traits in dairy cattle were suggested.     

Weighted gene co-expression network analysis reveals potential candidate genes affecting drip loss in pork

Drip loss is an essential evaluation indicator for pork quality. It is closely related to other meat quality indicators, including water-holding capacity, water loss rate and pH value at 45 min (pH₁) and 24 h post-mortem (pH₂), and is influenced by environmental and genetic factors and their interactions. We previously conducted differentially expressed gene analysis to identify candidate genes affecting drip loss using eight individuals with extremely high- and low-drip loss selected from 28 purebred Duroc pigs. Using 28 identical samples, in the present study, we performed weighted gene co-expression network analysis with drip loss and drip loss-related traits, including water-holding capacity, water loss rate, pH₁ and pH₂. A total of 25 modules were identified, and five of them correlated with at least two drip loss or drip loss-related traits. After functional enrichment analysis of genes in the five modules, three modules were found to be critical, as their genes were significantly involved in amino acid metabolism, immune response and apoptosis, which have potential relationships with drip loss. Furthermore, we identified five candidate genes affecting drip loss in one critical module, AASS, BCKDHB, ALDH6A1, MUT and MCCC1, as they overlapped with differentially expressed genes detected in our previous study, exhibited protein–protein interactions and had potential biological functions in affecting drip loss according to the literature. The outcomes of the present study enhance our understanding of the molecular mechanisms underlying drip loss and will aid in improving the pork quality.     

Discovery of genomic regions and candidate genes controlling shelling percentage using QTL‐seq approach in cultivated peanut (Arachis hypogaea L.)

Cultivated peanut (Arachis hypogaea L.) is an important grain legume providing high‐quality cooking oil, rich proteins and other nutrients. Shelling percentage (SP) is the 2nd most important agronomic trait after pod yield and this trait significantly affects the economic value of peanut in the market. Deployment of diagnostic markers through genomics‐assisted breeding (GAB) can accelerate the process of developing improved varieties with enhanced SP. In this context, we deployed the QTL‐seq approach to identify genomic regions and candidate genes controlling SP in a recombinant inbred line population (Yuanza 9102 × Xuzhou 68‐4). Four libraries (two parents and two extreme bulks) were constructed and sequenced, generating 456.89–790.32 million reads and achieving 91.85%–93.18% genome coverage and 14.04–21.37 mean read depth. Comprehensive analysis of two sets of data (Yuanza 9102/two bulks and Xuzhou 68‐4/two bulks) using the QTL‐seq pipeline resulted in discovery of two overlapped genomic regions (2.75 Mb on A09 and 1.1 Mb on B02). Nine candidate genes affected by 10 SNPs with non‐synonymous effects or in UTRs were identified in these regions for SP. Cost‐effective KASP (Kompetitive Allele‐Specific PCR) markers were developed for one SNP from A09 and three SNPs from B02 chromosome. Genotyping of the mapping population with these newly developed KASP markers confirmed the major control and stable expressions of these genomic regions across five environments. The identified candidate genomic regions and genes for SP further provide opportunity for gene cloning and deployment of diagnostic markers in molecular breeding for achieving high SP in improved varieties.     

Contribution of nitrogen fixation to first year Miscanthus × giganteus

Many characteristics make Miscanthus × giganteus an appealing bioenergy feedstock in temperate North America, but the degree to which this plant species interacts with nitrogen‐fixing bacteria remains understudied. Demonstration of associative nitrogen fixation in Miscanthus would support management with minimal fertilizer inputs that is demanded of long‐term biofuel sustainability. As a first step, we investigate the role of biological nitrogen fixation in nutrition of immature Miscanthus and temporal dynamics of plant‐associated nitrogen fixers. The contribution of biological nitrogen fixation to plant nitrogen acquisition in first year Miscanthus × giganteus was estimated using a yield‐dependent ¹⁵N isotope dilution model. Temporal changes in plant‐associated diazotroph relative abundance and community composition were analyzed with quantitative PCR and terminal restriction fragment length polymorphism of the nifH gene in rhizome and rhizosphere DNA extracts. We estimate 16% of new plant nitrogen was derived by nitrogen fixation during the growing season, despite non‐limiting soil nitrogen. Diazotroph communities from rhizome and rhizosphere changed with plant development and endophytic nitrogen fixers had significantly higher relative abundance and altered community composition at sampling dates in July and August. This study provides evidence for a small, but measurable, benefit of associative nitrogen fixation to first year Miscanthus × giganteus that underscores the potential and need for selection of breeding lines that maximize this trait.     

Genes with expression levels correlating to drip loss prove association of their polymorphism with water holding capacity of pork

Six genes that were known to exhibit expression levels that are correlated to drip loss BVES, SLC3A2, ZDHHC5, CS, COQ9, and EGFR have been for candidate gene analysis. Based on in silico analysis SNPs were detected, confirmed by sequencing, and used for genotyping. The SNPs were genotyped in about 1,800 animals from six pig populations including commercial herds of Pietrain (PI) and German Landrace (DL), different commercial herds of Pietrain × (German Large White × German Landrace) (PIF1(a/b/c)), and one experimental F₂-population Duroc × Pietrain (DUPI). Comparative and genetic mapping established the location of BVES on SSC1, of SLC3A2 and ZDHHC5 on SSC2, of CS on SSC5, of COQ9 on SSC6 and of EGFR on SSC9, respectively, coinciding with QTL regions for carcass and meat quality traits. BVES, SLC3A2, and CS revealed association at least with drip loss and with several other measures of water holding capacity (WHC). Moreover, COQ9 and EGFR were associated with several meat quality traits such as meat color and/or thawing loss. This study reveals statistic evidence in addition to the functional relationship of these genes to WHC previously evidenced by expression analysis. This study reveals positional and genetic statistical evidence for a link of genetic variation at these loci or close to them and promotes those six candidate genes as functional and/or positional candidate genes for meat quality traits.     

Genes influencing spinal bone mineral density in inbred F344, LEW, COP, and DA rats

Previously, we identified the regions of chromosomes 10q12–q31 and 15p16–q21 harbor quantitative trait loci (QTLs) for lumbar volumetric bone mineral density (vBMD) in female F2 rats derived from Fischer 344 (F344) × Lewis (LEW) and Copenhagen 2331 (COP) × Dark Agouti (DA) crosses. The purpose of this study is to identify the candidate genes within these QTL regions contributing to the variation in lumbar vBMD. RNA was extracted from bone tissue of F344, LEW, COP, and DA rats. Microarray analysis was performed using Affymetrix Rat Genome 230 2.0 Arrays. Genes differentially expressed among the rat strains were then ranked based on the strength of the correlation with lumbar vBMD in F2 animals derived from these rats. Quantitative PCR (qPCR) analysis was performed to confirm the prioritized candidate genes. A total of 285 genes were differentially expressed among all strains of rats with a false discovery rate less than 10%. Among these genes, 18 candidate genes were prioritized based on their strong correlation (r ² > 0.90) with lumbar vBMD. Of these, 14 genes (Akap1, Asgr2, Esd, Fam101b, Irf1, Lcp1, Ltc4s, Mdp-1, Pdhb, Plxdc1, Rabep1, Rhot1, Slc2a4, Xpo4) were confirmed by qPCR. We identified several novel candidate genes influencing spinal vBMD in rats.     

A whole genome sequence association study for puberty in a large Duroc × Erhualian F2 population

A large proportion of gilts and sows are culled from reproduction populations because of anestrus and pubertal reproductive failure. Selecting early onset of puberty gilts has a favorable effect on sows’ reproductivity. However, age at puberty is hard to be routinely measured in commercial herds. With molecular genetic predictors, identifying individuals that have a propensity for early onset of puberty can be simplified. We previously performed genome scanning and a genome‐wide association study for puberty in an F2 resource population using 183 microsatellites and 62 125 SNPs respectively. The detection power and resolution of identified quantitative trait loci were very low. Herein, we re‐sequenced 19 founders of the F2 resource population in high coverage, and whole genome sequences of F2 individuals were imputed to perform an association study for reproductive traits. A total of 2339 SNPs associated with pubertal reproductive failure were identified in the region of 30.94–40.74 Mb on SSC7, with the top one, positioned at 33.36 Mb, explaining 16% of the phenotypic variances. We improved the magnitude of the P‐value by 10E+5 to 10E+7 using the whole genome sequence rather than using low/middle density markers as in previous studies, and we narrowed down the QTL confidence interval to 5.25 Mb. Combining the annotation of gene function, RAB23 and BAK1 were perceived as the most compelling candidate genes. The identified loci may be useful in culling sows failing to show estrus by marker‐assisted selection to increase reproductive efficiency of swine herds.     

Systems genetic analysis of binge-like eating in a C57BL/6J x DBA/2J-F2 cross

Binge eating is a heritable trait associated with eating disorders and refers to the rapid consumption of a large quantity of energy-dense food that is, associated with loss of control and negative affect. Binge eating disorder is the most common eating disorder in the United States; however, the genetic basis is unknown. We previously identified robust mouse inbred strain differences between C57BL/6J and DBA/2J in binge-like eating of sweetened palatable food in an intermittent access, conditioned place preference paradigm. To map the genetic basis of changes in body weight and binge-like eating (BLE) and to identify candidate genes, we conducted quantitative trait locus (QTL) analysis in 128 C57BL/6J x DBA/2J-F2 mice combined with PheQTL and trait covariance analysis in GeneNetwork2 using legacy BXD-RI trait datasets. We identified a QTL on Chromosome 18 influencing changes in body weight across days in females (log of the odds [LOD] = 6.3; 1.5-LOD: 3–12 cM) that contains the candidate gene Zeb1. We also identified a sex-combined QTL influencing initial palatable food intake on Chromosome 5 (LOD = 5.8; 1.5-LOD: 21–28 cM) that contains the candidate gene Lcorl and a second QTL influencing escalated palatable food intake on Chromosome 6 in males (LOD = 5.4; 1.5-LOD: 50–59 cM) that contains the candidate genes Adipor2 and Plxnd1. Finally, we identified a suggestive QTL in females for slope of BLE on distal Chromosome 18 (LOD = 4.1; p = 0.055; 1.5-LOD: 23–35 cM). Future studies will use BXD-RI strains to fine map loci and support candidate gene nomination for gene editing.     

Evolution of soil carbon stocks under Miscanthus × giganteus and Miscanthus sinensis across contrasting environmental conditions

Miscanthus is a C4 bioenergy perennial crop characterized by its high potential yield. Our study aimed to compare the carbon storage capacities of Miscanthus sinensis (M. sinensis) with that of Miscanthus × giganteus (M. × giganteus) in field conditions in different types of soils in France. We set up a multi‐environment experimental network. On each trial, we tested two treatments: M. × giganteus established from rhizomes (Gr) and M. sinensis transplanted seedlings (Sp). We quantified the soil organic carbon (SOC) stock at equivalent soil mass for both genotypes in 2014 and 2019 and for two sampling depths: L1 (ca. 0–5 cm) and L1‐2 (ca. 0–30 cm). We also calculated the total and annual variation of the SOC stock and investigated factors that could explain the variation and the initial state of the SOC stock. ANOVAs were performed to compare the SOC stock, as well as the SOC stock variation rates across treatments and soil layers. Results showed that the soil bulk density did not vary significantly between 2014 and 2019 for both treatments (Gr and Sp). The SOC concentration (i.e. SOC expressed in g/kg) increased significantly between 2014 and 2019 in L1, whereas no significant evolution was found in L2 (ca. 5–30 cm). The SOC stock (i.e. SOC expressed in t/ha) increased significantly in the superficial layer L1 for M. × giganteus and M. sinensis, by 0.48 ± 0.41 and 0.54 ± 0.25 t ha^?1 year^?1 on average, respectively, although no significant change was detected in the layer L1‐2 for both genotypes. Moreover, SOC stocks in 2019 did not differ significantly between M. × giganteus and M. sinensis in the soil layers L1 and L1‐2. Lastly, our results showed that the initial SOC stock was significantly higher when miscanthus was grown after set‐aside than after annual crops.     

A new PCR-based marker on chromosome 4AL for resistance to pre-harvest sprouting in wheat (<Emphasis Type="Italic">Triticum aestivum</Emphasis> L.)

Pre-harvest sprouting (PHS) is a complex trait controlled by multiple genes with strong interaction between environment and genotype that makes it difficult to select breeding materials by phenotypic assessment. One of the most important genes for pre-harvest sprouting resistance is consistently identified on the long arm of chromosome 4A. The 4AL PHS tolerance gene has therefore been targeted by Australian white-grained wheat breeders. A new robust PCR marker for the PHS QTL on wheat chromosome 4AL based on candidate genes search was developed in this study. The new marker was mapped on 4AL deletion bin 13-0.59-0.66 using 4AL deletion lines derived from Chinese Spring. This marker is located on 4AL between molecular markers Xbarc170 and Xwg622 in the doubled-haploid wheat population Cranbrook × Halberd. It was mapped between molecular markers Xbarc170 and Xgwm269 that have been previously shown to be closely linked to grain dormancy in the doubled haploid wheat population SW95-50213 × Cunningham and was co-located with Xgwm269 in population Janz × AUS1408. This marker offers an additional efficient tool for marker-assisted selection of dormancy for white-grained wheat breeding. Comparative analysis indicated that the wheat chromosome 4AL QTL for seed dormancy and PHS resistance is homologous with the barley QTL on chromosome 5HL controlling seed dormancy and PHS resistance. This marker will facilitate identification of the gene associated with the 4A QTL that controls a major component of grain dormancy and PHS resistance.     

Chemical Composition and Antibacterial Activity of Iranian Lavandula × hybrida

Lavandin (Lavandula × hybrida) is an evergreen shrub and cultivated worldwide for its essential oil which possesses various biological activities. In this study, the essential oils were isolated from the leaves of ten lavandin populations in western Iran. The hydrodistilled essential oils were analyzed by GC‐FID/MS. Results indicated significant differences (P ≤ 0.05) among the various populations for the main essential oil constituents. The major components from different populations were 1,8‐cineole (31.64 – 47.94%), borneol (17.11 – 26.14%), and camphor (8.41 – 12.68%). In vitro antibacterial activity was evaluated against S. agalactiae, S. aureus, E. coli, and K. pneumoniae. The inhibition zones were in the range of 09.36 mm for S. aureus to 23.30 mm for E. coli. Results indicated that there was a significant correlation between essential oil composition and level of antibacterial efficacy expressed as inhibition zones.     

Transcriptome Profiling of Embryonic Development Rate in Rainbow Trout Advanced Backcross Introgression Lines

In rainbow trout (Oncorhynchus mykiss) and other fishes, embryonic development rate is an ecologically and evolutionarily important trait that is closely associated with survival and physiological performance later in life. To identify genes differentially regulated in fast and slow-developing embryos of rainbow trout, we examined gene expression across developmental time points in rainbow trout embryos possessing alleles linked to a major quantitative trait loci (QTL) for fast versus slow embryonic development rate. Whole genome expression microarray analyses were conducted using embryos from a fourth generation backcross family, whereby each backcross generation involved the introgression of the fast-developing alleles for a major development rate QTL into a slow-developing clonal line of rainbow trout. Embryos were collected at 15, 19, and 28 days post-fertilization; sex and QTL genotype were determined using molecular markers, and cDNA from 48 embryos were used for microarray analysis. A total of 183 features were identified with significant differences between embryonic development rate genotypes. Genes associated with cell cycle growth, muscle contraction and protein synthesis were expressed significantly higher in embryos with the fast-developing allele (Clearwater) than those with the slow-developing allele (Oregon State University), which may associate with fast growth and early body mass construction in embryo development. Across time points, individuals with the fast-developing QTL allele appeared to have earlier onset of these developmental processes when compared to individuals with the slow development alleles, even as early as 15 days post-fertilization. Differentially expressed candidate genes chosen for linkage mapping were localized primarily to regions outside of the major embryonic development rate QTL, with the exception of a single gene (very low-density lipoprotein receptor precursor).     

Identification of two major QTL for yellow seed color in two crosses of resynthesized <Emphasis Type="Italic">Brassica napus</Emphasis> line No. 2127-17

Yellow seed color, which results from a thinner seed coat, is associated with improved feed quality of rapeseed (Brassica napus L.) meal and increased oil and protein content. As this trait follows various genetic models under different genetic backgrounds, a study was performed in two genetic backgrounds to gain a better understanding of the genetic mechanisms underlying yellow seed color. The quantitative trait locus (QTL) analysis was undertaken using two crosses, Quantum ×No. 2127-17 (HZ-1) and No. 2127-17 × 94,570 (HZ-2). In the HZ-1 population, three putative QTL were detected in linkage groups N18, N5, and N3, respectively. For all of them, yellow seed color arose from the No. 2127-17 alleles. Of these QTL, the one in linkage group N18 (Bnsc-18a) explained more than half of the phenotypic variation. In the HZ-2 population, three QTL were found in linkage groups N9, N18, and N8, respectively. Of these QTL, that in linkage group N9 (Bnsc-9a) explained more than half of the phenotypic variation, whereas the QTL Bnsc-18a had a low seed color value and explained only 9.03–11.72% of the phenotypic variation. Bulked segregant analysis (BSA) of the extremes of a BC1 population derived from the cross of No. 2127-17 × 94,570 (HZ-3) identified one major gene that was identical with the QTL Bnsc-9a detected in the HZ-2 population. The QTL Bnsc-18a was common in the HZ-1 and HZ-2 populations, and the others were population-specific. These results suggested that different black-seeded forms had different seed color genes.