Single-step genome-wide association analyses for milk urea concentration in Walloon Holstein cows

The aims of this study were to conduct a single-step genome-wide association to identify genomic regions associated with milk urea (MU) and to estimate genetic correlations between MU and milk yield (MY), milk composition (calcium content [CC], fat percentage [FP], protein percentage [PP], and casein percentage [CNP]), and cheese-making properties (CMP; coagulation time [CT], curd firmness after 30 min from rennet addition [a30], and titratable acidity [TA]). The used data have been collected from 2015 to 2020 on 78,073 first-parity (485,218 test-day records) and 48,766 second-parity (284,942 test-day records) Holstein cows distributed in 671 herds in the Walloon Region of Belgium. Data of 565,533 SNP located on 29 BTA of 6,617 animals (1,712 males) were used. Random regression test-day models were used to estimate genetic parameters through the Bayesian Gibbs sampling method. The SNP solutions were estimated using a single-step genomic BLUP approach. The proportion of the total additive genetic variance explained by windows of 50 consecutive SNPs (with an average size of ∼216 kb) was calculated, and the top-3 genomic regions explaining the largest rate of the genetic variance were considered promising regions and used to identify potential candidate genes. Mean (SD) MU was 25.38 (8.02) mg/dL and 25.03 (8.06) mg/dL in the first and second lactation, respectively. Mean heritability estimates for daily MU were 0.21 and 0.23 for the first and second lactation, respectively. The genetic correlations estimated between MU and MY, milk composition, and CMP were quite low (ranged from -0.10 [CC] to 0.10 [TA] and -0.05 [CT] to 0.13 [TA] for the first and second lactations, respectively). The top-3 regions associated with MU were located from 80.61 to 80.74 Mb on BTA6, 103.26 to 103.41 Mb on BTA11, and 1.59 to 2.15 Mb on BTA14. Genes including KCNT1, MROH1, SHARPIN, TSSK5, CPSF1, HSF1, TONSL, DGAT1, SCX, and MAF1 were identified as positional candidate genes for MU. The findings of this study can be used for a better understanding of the genomic architecture underlying MU in Holstein cattle.

Genome-wide association study for body weight and feed consumption traits in Japanese quail using Bayesian approaches

The aim of this study was to perform a genome-wide association study (GWAS) based on Bayes A and Bayes B statistical methods to identify genomic loci and candidate genes associated with body weight gain, feed intake, and feed conversion ratio in Japanese quail. For this purpose, genomic data obtained from Illumina iSelect 4K quail SNP chip were utilized. After implementing various quality control steps, genotype data from a total of 875 birds for 2,015 SNP markers were used for subsequent analyses. The Bayesian analyses were performed using hibayes package in R (version 4.3.1) and Gibbs sampling algorithm. The results of the analyses showed that Bayes A accounted for 11.43, 11.65, and 11.39% of the phenotypic variance for body weight gain, feed intake, and feed conversion ratio, respectively, while the variance explained by Bayes B was 7.02, 8.61, and 6.48%, respectively. Therefore, in the current study, results obtained from Bayes A were used for further analyses. In order to perform the gene enrichment analysis and to identify the functional pathways and classes of genes that are over-represented in a large set of genes associated with each trait, all markers that accounted for more than 0.1% of the phenotypic variance for each trait were used. The results of this analysis revealed a total of 23, 38, and 14 SNP markers associated with body weight gain, feed intake, and feed conversion ratio in Japanese quail, respectively. The results of the gene enrichment analysis led to the identification of biological pathways (and candidate genes) related to lipid phosphorylation (TTC7A gene) and cell junction (FGFR4 and FLRT2 genes) associated with body weight gain, calcium signaling pathway (ADCY2 and CAMK1D genes) associated with feed intake, and glycerolipid metabolic process (LIPC gene), lipid metabolic process (ADGRF5 and ESR1 genes), and glutathione transferase (GSTK1 gene) associated with feed conversion ratio. Overall, the findings of this study can provide valuable insights into the genetic architecture of growth and feed consumption traits in Japanese quail.

Single-step genome-wide association analyses for selected infrared-predicted cheese-making traits in Walloon Holstein cows

This study aimed to perform genome-wide association study to identify genomic regions associated with milk production and cheese-making properties (CMP) in Walloon Holstein cows. The studied traits were milk yield, fat percentage, protein percentage, casein percentage (CNP), calcium content, somatic cell score (SCS), coagulation time, curd firmness after 30 min from rennet addition, and titratable acidity. The used data have been collected from 2014 to 2020 on 78,073 first-parity (485,218 test-day records), 48,766 second-parity (284,942 test-day records), and 21,948 third-parity (105,112 test-day records) Holstein cows distributed in 671 herds in the Walloon Region of Belgium. Data of 565,533 single nucleotide polymorphisms (SNP), located on 29 Bos taurus autosomes (BTA) of 6,617 animals (1,712 males), were used. Random regression test-day models were used to estimate genetic parameters through the Bayesian Gibbs sampling method. The SNP solutions were estimated using a single-step genomic BLUP approach. The proportion of the total additive genetic variance explained by windows of 50 consecutive SNPs (with an average size of ∼216 KB) was calculated, and regions accounting for at least 1.0% of the total additive genetic variance were used to search for positional candidate genes. Heritability estimates for the studied traits ranged from 0.10 (SCS) to 0.53 (CNP), 0.10 (SCS) to 0.50 (CNP), and 0.12 (SCS) to 0.49 (CNP) in the first, second, and third parity, respectively. Genome-wide association analyses identified 6 genomic regions (BTA1, BTA14 [4 regions], and BTA20) associated with the considered traits. Genes including the SLC37A1 (BTA1), SHARPIN, MROH1, DGAT1, FAM83H, TIGD5, MROH6, NAPRT, ADGRB1, GML, LYPD2, JRK (BTA14), and TRIO (BTA20) were identified as positional candidate genes for the studied CMP. The findings of this study help to unravel the genomic background of a cow's ability for cheese production and can be used for the future implementation and use of genomic evaluation to improve the cheese-making traits in Walloon Holstein cows.

JustRNA: a database of plant long noncoding RNA expression profiles and functional network

Long noncoding RNAs (lncRNAs) are regulatory RNAs involved in numerous biological processes. Many plant lncRNAs have been identified, but their regulatory mechanisms remain largely unknown. A resource that enables the investigation of lncRNA activity under various conditions is required because the co-expression between lncRNAs and protein-coding genes may reveal the effects of lncRNAs. This study developed JustRNA, an expression profiling resource for plant lncRNAs. The platform currently contains 1 088 565 lncRNA annotations for 80 plant species. In addition, it includes 3692 RNA-seq samples derived from 825 conditions in six model plants. Functional network reconstruction provides insight into the regulatory roles of lncRNAs. Genomic association analysis and microRNA target prediction can be employed to depict potential interactions with nearby genes and microRNAs, respectively. Subsequent co-expression analysis can be employed to strengthen confidence in the interactions among genes. Chromatin immunoprecipitation sequencing data of transcription factors and histone modifications were integrated into the JustRNA platform to identify the transcriptional regulation of lncRNAs in several plant species. The JustRNA platform provides researchers with valuable insight into the regulatory mechanisms of plant lncRNAs. JustRNA is a free platform that can be accessed at http://JustRNA.itps.ncku.edu.tw.

Functional colour genes and signals of selection in colour-polymorphic salamanders

Coloration has been associated with multiple biologically relevant traits that drive adaptation and diversification in many taxa. However, despite the great diversity of colour patterns present in amphibians the underlying molecular basis is largely unknown. Here, we use insight from a highly colour-variable lineage of the European fire salamander (Salamandra salamandra bernardezi) to identify functional associations with striking variation in colour morph and pattern. The three focal colour morphs-ancestral black-yellow striped, fully yellow and fully brown-differed in pattern, visible coloration and cellular composition. From population genomic analyses of up to 4,702 loci, we found no correlations of neutral population genetic structure with colour morph. However, we identified 21 loci with genotype-phenotype associations, several of which relate to known colour genes. Furthermore, we inferred response to selection at up to 142 loci between the colour morphs, again including several that relate to coloration genes. By transcriptomic analysis across all different combinations, we found 196 differentially expressed genes between yellow, brown and black skin, 63 of which are candidate genes involved in animal coloration. The concordance across different statistical approaches and 'omic data sets provide several lines of evidence for loci linked to functional differences between colour morphs, including TYR, CAMK1 and PMEL. We found little association between colour morph and the metabolomic profile of its toxic compounds from the skin secretions. Our research suggests that current ecological and evolutionary hypotheses for the origins and maintenance of these striking colour morphs may need to be revisited.

Lineage-specific expansions of TET/JBP genes and a new class of DNA transposons shape fungal genomic and epigenetic landscapes

TET/JBP dioxygenases oxidize methylpyrimidines in nucleic acids and are implicated in generation of epigenetic marks and potential intermediates for DNA demethylation. We show that TET/JBP genes are lineage-specifically expanded in all major clades of basidiomycete fungi, with the majority of copies predicted to encode catalytically active proteins. This pattern differs starkly from the situation in most other organisms that possess just a single or a few copies of the TET/JBP family. In most basidiomycetes, TET/JBP genes are frequently linked to a unique class of transposons, KDZ (Kyakuja, Dileera, and Zisupton) and appear to have dispersed across chromosomes along with them. Several of these elements typically encode additional proteins, including a divergent version of the HMG domain. Analysis of their transposases shows that they contain a previously uncharacterized version of the RNase H fold with multiple distinctive Zn-chelating motifs and a unique insert, which are predicted to play roles in structural stabilization and target sequence recognition, respectively. We reconstruct the complex evolutionary history of TET/JBPs and associated transposons as involving multiple rounds of expansion with concomitant lineage sorting and loss, along with several capture events of TET/JBP genes by different transposon clades. On a few occasions, these TET/JBP genes were also laterally transferred to certain Ascomycota, Glomeromycota, Viridiplantae, and Amoebozoa. One such is an inactive version, calnexin-independence factor 1 (Cif1), from Schizosaccharomyces pombe, which has been implicated in inducing an epigenetically transmitted prion state. We argue that this unique transposon-TET/JBP association is likely to play important roles in speciation during evolution and epigenetic regulation.