PigBiobank: a valuable resource for understanding genetic and biological mechanisms of diverse complex traits in pigs.

To fully unlock the potential of pigs as both agricultural species for animal-based protein food and biomedical models for human biology and disease, a comprehensive understanding of molecular and cellular mechanisms underlying various complex phenotypes in pigs and how the findings can be translated to other species, especially humans, are urgently needed. Here, within the Farm animal Genotype-Tissue Expression (FarmGTEx) project, we build the PigBiobank (http://pigbiobank.farmgtex.org) to systematically investigate the relationships among genomic variants, regulatory elements, genes, molecular networks, tissues and complex traits in pigs. This first version of the PigBiobank curates 71 885 pigs with both genotypes and phenotypes from over 100 pig breeds worldwide, covering 264 distinct complex traits. The PigBiobank has the following functions: (i) imputed sequence-based genotype-phenotype associations via a standardized and uniform pipeline, (ii) molecular and cellular mechanisms underlying trait-associations via integrating multi-omics data, (iii) cross-species gene mapping of complex traits via transcriptome-wide association studies, and (iv) high-quality results display and visualization. The PigBiobank will be updated timely with the development of the FarmGTEx-PigGTEx project, serving as an open-access and easy-to-use resource for genetically and biologically dissecting complex traits in pigs and translating the findings to other species.

Resting cell formation in the marine diatom Thalassiosira pseudonana.

Resting cells represent a survival strategy employed by diatoms to endure prolonged periods of unfavourable conditions. In the oceans, many diatoms sink at the end of their blooming season and therefore need to endure cold and dark conditions in the deeper layers of the water column. How they survive these conditions is largely unknown. We conducted an integrative analysis encompassing methods from histology, physiology, biochemistry, and genetics to reveal the biological mechanism of resting-cell formation in the model diatom Thalassiosira pseudonana. Resting-cell formation was triggered by a decrease in light and temperature with subsequent catabolism of storage compounds. Resting cells were characterised by an acidic and viscous cytoplasm and altered morphology of the chloroplast ultrastructure. The formation of resting cells in T. pseudonana is an energy demanding process required for a biophysical alteration of the cytosol and chloroplasts to endure the unfavourable conditions of the deeper ocean as photosynthetic organisms. However, most resting cells (> 90%) germinate upon return to favorable growth conditions.

Real-time in situ fluorescence imaging of telomerase and miR378 in living cells using a two-color DNA tetrahedron nanoprobe combined with molecular beacons.

A biosensor that can detect biomarkers accurately, quickly, and conveniently is important for the diagnosis of various diseases. However, most of the existing detection methods require sample extraction, which makes it difficult to detect and image intracellular molecules or to detect two different types of biomarkers simultaneously. In this study, we constructed a DNA tetrahedral nanoprobe (DTP) capable of detecting both miR378 and telomerase, both of which are tumor markers. In the presence of miR378, FAM on the molecular beacon of DTP fluoresced via Förster resonance energy transfer (FRET), and the limit of detection was 476 pM with excellent specificity. When present, telomerase binds to telomerase substrate (TS) primers, extending the repeat sequence (TTAGGG)n to trigger Cy3 fluorescence. A strong linear relationship existed between the fluorescence intensity of Cy3 and the number of HeLa cells. The limit of detection was 800 HeLa cells. In addition, DTP was less cytotoxic to and biocompatible with HeLa cells and fluoresced only in cancer cells, which can help to sensitively distinguish between normal and cancer cells. In conclusion, DTP can simultaneously detect the content of miR378 and activity of telomerase and realize intracellular imaging, which has broad application prospects in early cancer diagnosis and treatment.

Detection ofmiRNA-378based on a catalytic hairpin self-assembly reaction combined with gold nanoparticle colorimetry.

Recent studies have shown that abnormalmiRNA-378expression is a rule, rather than an exception, in cervical cancer and can be used as a diagnostic and prognostic biomarker to assess tumor initiation. In this study, we developed a general, sensitive strategy for detectingmiRNA-378using catalytic hairpin self-assembly (CHA) combined with gold nanoparticles (AuNP) colorimetry. The presence ofmiRNA-378triggers the repeated self-assembly of two designed hairpin DNAs (H1 and H2) into dsDNA polymers, which leads to changes in the surface plasmon resonance absorption band and the macroscopic color of the AuNP colloids due to the formation of nanoparticle-DNA conjugates. This experimental phenomenon can be observed by ultraviolet-visible spectrometry or even with the naked eye. Using this method,miRNA-378could be quantitatively detected at the picomolar level (as low as 20.7 pM). Compared with traditional methods, such as quantitative polymerase chain reaction and RNA blotting, this strategy has a simple operation, low cost, and high sensitivity and selectivity, and thus, exhibits significant potential for miRNA detection.

Evaluating three strategies of genome-wide association analysis for integrating data from multiple populations.

In livestock, genome-wide association studies (GWAS) are usually conducted in a single population (single-GWAS) with limited sample size and detection power. To enhance the detection power of GWAS, meta-analysis of GWAS (meta-GWAS) and mega-analysis of GWAS (mega-GWAS) have been proposed to integrate data from multiple populations at the level of summary statistics or individual data, respectively. However, there is a lack of comparison for these different strategies, which makes it difficult to guide the best practice of GWAS integrating data from multiple study populations. To maximize the comparison of different association analysis strategies across multiple populations, we conducted single-GWAS, meta-GWAS, and mega-GWAS for the backfat thickness of 100 kg (BFT_100) and days to 100 kg (DAYS_100) within each of the three commercial pig breeds (Duroc, Yorkshire, and Landrace). Based on controlling the genome inflation factor to one, we calculated corrected p-values (pC ). In Yorkshire, with the largest sample size, mega-GWAS, meta-GWAS and single-GWAS detected 149, 38 and 20 significant SNPs (pC < 1E-5) associated with BFT_100, as well as 26, four, and one QTL, respectively. Among them, pC of SNPs from mega-GWAS was the lowest, followed by meta-GWAS and single-GWAS. The correlation of pC among the three GWAS strategies ranged from 0.60 to 0.75 and the correlation of SNP effect values between meta-GWAS and mega-GWAS was 0.74, all showing good agreement. Collectively, even though there are differences in the integration of individual data or summary statistics, integrating data from multiple populations is an effective means of genetic argument for complex traits, especially mega-GWAS versus single-GWAS.

The allelopathic potential of red macroalga Pyropia haitanensis solvent extracts on controlling bloom-forming microalgae: Insights into the inhibitory compounds.

Various strategies have been explored to mitigate the impact of harmful algal blooms (HABs). While chemical and physical methods have traditionally been employed to regulate microalgal growth, their prolonged adverse effects on the ecosystem are a cause for concern. Recognizing the integral role of macroalgae within the ecosystem, this study reveals the anti-algal properties of solvent-based extracts derived from the red macroalga Pyropia haitanensis as a means of preventing microalgal blooms. In our investigation, we initially assessed the growth-inhibitory effects of methanol and acetone extracts from P. haitanensis on five microalgae known to contribute to bloom-formation. Significantly reduced growth was observed in all microalgal species when inoculated with both methanol and acetone extracts. Further analysis revealed the effectiveness of the methanol extract (ME), and further fractionation with petroleum ether (PE), ethyl acetate (EA), and n-butanol (NB) for testing against Skeletonema costatum and Pseudo-nitzschia pungens. The methanol fractions exhibited strong inhibition, resulting in the complete elimination of both microalgae after 96 hours of exposure to PE, EA, and NB extracts. Gas Chromatography-Mass Spectroscopy (GC-MS) analysis of the ME and its solvent fractions identified 49 confirmed compounds. These compounds are likely potential contributors to the observed inhibition of microalgal growth. In conclusion, our findings suggest that solvent extracts from P. haitanensis possess substantial potential for the control of HABs, offering a promising avenue for further research and application in ecosystem management.

Fluorescent DNA tetrahedral probe with catalytic hairpin self-assembly reaction for imaging of miR-21 and miR-155 in living cells.

A CHA-based fluorescent DNA tetrahedral probe (FDTp) has been designed to detect the microRNAs miR-21 and miR-155 sensitively and specifically in living cells. The design consisted of functional elements (H1, H2, and Protector) connected to a DNA tetrahedron modified with two pairs of fluorophores and quenching groups. In the presence of miR-21, the chain displacement effect was triggered and Cy3 fluorescence was emitted. In the presence of miR-155, the signal of the catalytic hairpin assembly (CHA) between H1 and H2 on FDTp was amplified, making the fluorescence of FAM sensitive to miR-155. Using this method, the detection limit for miR-155 was 5 pM. The FDTp successfully imaged miR-21 and miR-155 in living cells and distinguished a variety of cell lines based on their expression levels of miR-21 and miR-155. The detection and imaging of dual targets in this design ensured the accuracy of tumor diagnosis and provided a new method for early tumor diagnosis.

Integrating large-scale meta-analysis of genome-wide association studies improve the genomic prediction accuracy for combined pig populations.

The strategy of combining reference populations has been widely recognized as an effective way to enhance the accuracy of genomic prediction (GP). This study investigated the efficiency of genomic prediction using prior information and combined reference population. In total, prior information considering trait-associated single nucleotide polymorphisms (SNPs) obtained from meta-analysis of genome-wide association studies (GWAS meta-analysis) was incorporated into three models to assess the performance of GP using combined reference populations. Two different Yorkshire populations with imputed whole genome sequence (WGS) data (9,741,620 SNPs), named as P1 (1259 individuals) and P2 (1018 individuals), were used to predict genomic estimated breeding values for three live carcass traits, including backfat thickness, loin muscle area, and loin muscle depth. A 10 × 5 fold cross-validation was used to evaluate the prediction accuracy of 203 randomly selected candidate pigs from the P2 population and the reference population consisted of the remaining pigs from P2 and the stepwise added pigs from P1. By integrating SNPs with different p-value thresholds from GWAS meta-analysis downloaded from PigGTEx Project, the prediction accuracy of GBLUP, genomic feature BLUP (GFBLUP) and GBLUP given genetic architecture (BLUP|GA) were compared. Moreover, we explored effects of reference population size and heritability enrichment of genomic features on the prediction accuracy improvement of GFBLUP and BLUP|GA relative to GBLUP. The prediction accuracy of GBLUP using all WGS markers showed average improvement of 4.380% using the P1 + P2 reference population compared with the P2 reference population. Using the combined reference population, GFBLUP and BLUP|GA yielded 6.179% and 5.525% higher accuracies than GBLUP using all SNPs based on the single reference population, respectively. Positive regression coefficients were estimated in relation to the improvement in prediction accuracy (between GFBLUP/BLUP|GA and GBLUP) and the size of the reference as well as the heritability enrichment of genomic features. Compared to the classic GBLUP model, GFBLUP and BLUP|GA models integrating GWAS meta-analysis information increase the prediction accuracy and using combined populations with enlarged reference population size further enhances prediction accuracy of the two approaches. The heritability enrichment of genomic features can be used as an indicator to reflect weather prior information is accurately presented.

Integrating Multiple Database Resources to Elucidate the Gene Flow in Southeast Asian Pig Populations.

The domestic pig (Sus scrofa) and its subfamilies have experienced long-term and extensive gene flow, particularly in Southeast Asia. Here, we analyzed 236 pigs, focusing on Yunnan indigenous, European commercial, East Asian, and Southeast Asian breeds, using the Pig Genomics Reference Panel (PGRP v1) of Pig Genotype-Tissue Expression (PigGTEx) to investigate gene flow and associated complex traits by integrating multiple database resources. In this study, we discovered evidence of admixtures from European pigs into the genome of Yunnan indigenous pigs. Additionally, we hypothesized that a potential conceptual gene flow route that may have contributed to the genetic composition of the Diannan small-ear pig is a gene exchange from the Vietnamese pig. Based on the most stringent gene introgression scan using the fd statistic, we identified three specific loci on chromosome 8, ranging from 51.65 to 52.45 Mb, which exhibited strong signatures of selection and harbored the NAF1, NPY1R, and NPY5R genes. These genes are associated with complex traits, such as fat mass, immunity, and litter weight, in pigs, as supported by multiple bio-functionalization databases. We utilized multiple databases to explore the potential dynamics of genetic exchange in Southeast Asian pig populations and elucidated specific gene functionalities.

Genome-wide association analysis of heifer livability and early first calving in Holstein cattle.

BACKGROUND: The survival and fertility of heifers are critical factors for the success of dairy farms. The mortality of heifers poses a significant challenge to the management and profitability of the dairy industry. In dairy farming, achieving early first calving of heifers is also essential for optimal productivity and sustainability. Recently, Council on Dairy Cattle Breeding (CDCB) and USDA have developed new evaluations of heifer health and fertility traits. However, the genetic basis of these traits has yet to be thoroughly studied. RESULTS: Leveraging the extensive U.S dairy genomic database maintained at CDCB, we conducted large-scale GWAS analyses of two heifer traits, livability and early first calving. Despite the large sample size, we found no major QTL for heifer livability. However, we identified a major QTL in the bovine MHC region associated with early first calving. Our GO analysis based on nearby genes detected 91 significant GO terms with a large proportion related to the immune system. This QTL in the MHC region was also confirmed in the analysis of 27 K bull with imputed sequence variants. Since these traits have few major QTL, we evaluated the genome-wide distribution of GWAS signals across different functional genomics categories. For heifer livability, we observed significant enrichment in promotor and enhancer-related regions. For early calving, we found more associations in active TSS, active Elements, and Insulator. We also identified significant enrichment of CDS and conserved variants in the GWAS results of both traits. By linking GWAS results and transcriptome data from the CattleGTEx project via TWAS, we detected four and 23 significant gene-trait association pairs for heifer livability and early calving, respectively. Interestingly, we discovered six genes for early calving in the Bovine MHC region, including two genes in lymph node tissue and one gene each in blood, adipose, hypothalamus, and leukocyte. CONCLUSION: Our large-scale GWAS analyses of two heifer traits identified a major QTL in the bovine MHC region for early first calving. Additional functional enrichment and TWAS analyses confirmed the MHC QTL with relevant biological evidence. Our results revealed the complex genetic basis of heifer health and fertility traits and indicated a potential connection between the immune system and reproduction in cattle.

Comprehensive analyses of 723 transcriptomes enhance genetic and biological interpretations for complex traits in cattle.

By uniformly analyzing 723 RNA-seq data from 91 tissues and cell types, we built a comprehensive gene atlas and studied tissue specificity of genes in cattle. We demonstrated that tissue-specific genes significantly reflected the tissue-relevant biology, showing distinct promoter methylation and evolution patterns (e.g., brain-specific genes evolve slowest, whereas testis-specific genes evolve fastest). Through integrative analyses of those tissue-specific genes with large-scale genome-wide association studies, we detected relevant tissues/cell types and candidate genes for 45 economically important traits in cattle, including blood/immune system (e.g., CCDC88C) for male fertility, brain (e.g., TRIM46 and RAB6A) for milk production, and multiple growth-related tissues (e.g., FGF6 and CCND2) for body conformation. We validated these findings by using epigenomic data across major somatic tissues and sperm. Collectively, our findings provided novel insights into the genetic and biological mechanisms underlying complex traits in cattle, and our transcriptome atlas can serve as a primary source for biological interpretation, functional validation, studies of adaptive evolution, and genomic improvement in livestock.

Long-term changes in Noctiluca scintillans blooms along the Chinese coast from 1933 to 2020.

Noctiluca scintillans is one of the most common harmful algal species and widely known due to its bioluminescence. In this study, the spatial distribution, seasonal variations, and long-term trends of N. scintillans blooms in China and the related drivers were analyzed and discussed. From 1933 to 2020, a total of 265 events of N. scintillans blooms were recorded in Chinese coastal waters, with a total duration of 1052 days. The first N. scintillans bloom occurred in Zhejiang in 1933, and only three events were recorded before 1980. From 1981 to 2020, N. scintillans caused harmful algal blooms (HABs) almost every year, both the average duration and the proportion of multiphase HABs showed an increasing trend. 1986-1992, 2002-2004, and 2009-2016 were the three peak periods with a frequency of no less than five events of N. scintillans blooms per year. In terms of spatial distribution, N. scintillans blooms spread from the Southeast China Sea to the Bohai Sea after 2000, Guangdong, Fujian, and Hebei were the three provinces with the highest numbers of recorded events of N. scintillans blooms. Moreover, 86.8% of the events of N. scintillans blooms occurred in spring (March, April, and May) and summer (June, July, and August). Among environmental factors, the dissolved inorganic phosphate, dissolved silicate and chemical oxygen demand were significantly correlated with the cell density of N. scintillans during N. scintillans blooms, and most of N. scintillans blooms were recorded in the temperature range of 18.0-25.0°C. Precipitation, hydrodynamics, water temperature, and food availability might be the main factors affecting the spatial-temporal distribution of N. scintillans blooms along the Chinese coast.

Phytoplankton community structure and environmental factors during the outbreak of Crown-of-Thorns Starfish in Xisha Islands, South China Sea.

The "larval starvation hypothesis" proposed that the growing frequency of Crown-of-Thorns Starfish (CoTS) outbreaks could be attributed to increased availability of phytoplankton. However, comprehensive field investigation on the living environment of CoTS larvae and the availability of phytoplankton are still lacking. A cruise was conducted in June 2022 in Xisha Islands, South China Sea, to study the interaction between environmental conditions and phytoplankton communities during CoTS outbreak period. The average concentrations of dissolved inorganic phosphorus (0.05 ± 0.01 µmol L-1), dissolved inorganic nitrogen (0.66 ± 0.8 µmol L-1) and chlorophyll a (0.05 ± 0.05 µg L-1) suggested that phytoplankton may be limited for CoTS larvae in Xisha Islands. Microscopic observation and high-throughput sequencing were used to study the composition and structure of the phytoplankton communities. Bacillariophyta predominated in phytoplankton communities with the highest abundance and species richness. 29 dominant species, including 4 species with size-range preferred by CoTS larvae, were identified in Xisha Islands. The diversity index of all stations indicated a species-rich and structure-stable phytoplankton community in Xisha Islands during the period of CoTS outbreak, which may contribute to CoTS outbreak. These findings revealed the structure of phytoplankton community and environmental factors in the study area during CoTS outbreak, providing the groundwork for future research into the causes and processes of CoTS outbreak.

An artificial multienzyme cascade for the whole-cell synthesis of rare ketoses from glycerol.

PURPOSE: In our previous study, we constructed a one-pot multi-enzyme system for rare ketoses synthesis based on L-rhamnulose-1-phosphate aldolase (RhaD) from accessible glycerol in vitro. To eliminate tedious purification of enzymes, a facile Escherichia coli whole-cell cascade platform was established in this study. METHODS: To enhance the conversion rate, the reaction conditions, substrate concentrations and expressions of related enzymes were extensively optimized. RESULTS: The biosynthetic route for the cascade synthesis of rare ketoses in whole cells was successfully constructed and three rare ketoses including D-allulose, D-sorbose and L-fructose were produced using glycerol and D/L-glyceraldehyde (GA). Under optimized conditions, the conversion rates of rare ketoses were 85.0% and 93.0% using D-GA and L-GA as the receptor, respectively. Furthermore, alditol oxidase (AldO) was introduced to the whole-cell system to generate D-GA from glycerol, and the total production yield of D-sorbose and D-allulose was 8.2 g l-1 only from the sole carbon source glycerol. CONCLUSION: This study demonstrates a feasible and cost-efficient method for rare sugars synthesis and can also be applied to the green synthesis of other value-added chemicals from glycerol.

Functional annotation of regulatory elements in cattle genome reveals the roles of extracellular interaction and dynamic change of chromatin states in rumen development during weaning.

We profiled landscapes of bovine regulatory elements and explored dynamic changes of chromatin states in rumen development during weaning. The regulatory elements (15 chromatin states) and their coordinated activities in cattle were defined through genome-wide profiling of four histone modifications, CTCF-binding, DNA accessibility, DNA methylation, and transcriptome in rumen epithelial tissues. Each chromatin state presented specific enrichment for sequence ontology, methylation, trait-associated variants, transcription, gene expression-associated variants, selection signatures, and evolutionarily conserved elements. During weaning, weak enhancers and flanking active transcriptional start sites (TSS) were the most dynamic chromatin states and occurred in tandem with significant variations in gene expression and DNA methylation, significantly associated with stature, production, and reproduction economic traits. By comparing with in vitro cultured epithelial cells and in vivo rumen tissues, we showed the commonness and uniqueness of these results, especially the roles of cell interactions and mitochondrial activities in tissue development.

Ion fluxes Involved in the Adaptation of the Estuarine Diatom Coscinodiscus centralis Ehrenberg to Salinity Stress.

Although estuarine diatoms have a wide range of salt tolerance, they are often severely stressed by elevated salt concentrations. It remains poorly understood how estuarine diatoms maintain ionic homeostasis under high-salinity conditions. Using a scanning ion-selective electrode technique, this study determined the fluxes of H+, Na+, and K+ involved in the acclimatization of the estuarine diatom Coscinodiscus centralis Ehrenberg after an elevation in salinity from 15 psu to 35 psu. The C. centralis cells exhibited marked H+ effluxes after a transient treatment (TT, 30 min) and short-term treatment (ST, 24 h). However, a drastic shift of H+ efflux toward an influx was induced in the long-term treatment (LT, 10 days). The Na+ flux under TT, ST, and LT salinity conditions was found to accelerate the Na+ efflux. More pronounced effects were observed under the ST and LT salinity conditions compared to the TT salinity condition. The K+ influx showed a significant increase under the LT salinity condition. However, the salinity-induced Na+/H+ exchange in the estuarine diatom was inhibited by amiloride and sodium orthovanadate. These results indicate that the Na+ extrusion in salt-stressed cells is mainly the result of an active Na+/H+ antiport across the plasma membrane. The pattern of ion fluxes under the TT and ST salinity conditions were different from those under the LT salinity conditions, suggesting an incomplete regulation of the acclimation process in the estuarine diatom under short-term salinity stress.

Investigation of rumen long noncoding RNA before and after weaning in cattle.

BACKGROUND: This study aimed to identify long non-coding RNA (lncRNA) from the rumen tissue in dairy cattle, explore their features including expression and conservation levels, and reveal potential links between lncRNA and complex traits that may indicate important functional impacts of rumen lncRNA during the transition to the weaning period. RESULTS: A total of six cattle rumen samples were taken with three replicates from before and after weaning periods, respectively. Total RNAs were extracted and sequenced with lncRNA discovered based on size, coding potential, sequence homology, and known protein domains. As a result, 404 and 234 rumen lncRNAs were identified before and after weaning, respectively. However, only nine of them were shared under two conditions, with 395 lncRNAs found only in pre-weaning tissues and 225 only in post-weaning samples. Interestingly, none of the nine common lncRNAs were differentially expressed between the two weaning conditions. LncRNA averaged shorter length, lower expression, and lower conservation scores than the genome overall, which is consistent with general lncRNA characteristics. By integrating rumen lncRNA before and after weaning with large-scale GWAS results in cattle, we reported significant enrichment of both pre- and after-weaning lncRNA with traits of economic importance including production, reproduction, health, and body conformation phenotypes. CONCLUSIONS: The majority of rumen lncRNAs are uniquely expressed in one of the two weaning conditions, indicating a functional role of lncRNA in rumen development and transition of weaning. Notably, both pre- and post-weaning lncRNA showed significant enrichment with a variety of complex traits in dairy cattle, suggesting the importance of rumen lncRNA for cattle performance in the adult stage. These relationships should be further investigated to better understand the specific roles lncRNAs are playing in rumen development and cow performance.

Genome-wide association studies for immunoglobulin concentrations in colostrum and serum in Chinese Holstein.

BACKGROUND: The early death and health problems of calves caused substantial economic losses in the dairy industry. As the immune system of neonates has not been fully developed, the absorption of maternal immunoglobulin (Ig) from colostrum is essential in protecting newborn calves against common disease organisms in their early life. The overwhelming majority of Ig in bovine whey is transported from the serum. Therefore, Ig concentration in the colostrum and serum of dairy cows are critical traits when estimating the potential disease resistance of its offspring. RESULTS: Colostrum, blood, and hair follicle samples were collected from 588 Chinese Holstein cows within 24 h after calving. The concentration of total IgG, IgG1, IgG2, IgA and IgM in both colostrum and serum were detected via ELISA methods. With GCTA software, genome-wide association studies (GWASs) were performed with 91,620 SNPs genotyped by GeneSeek 150 K (140,668 SNPs) chips. As a result, 1, 5, 1 and 29 significant SNPs were detected associated with the concentrations of colostrum IgG1, IgG2, IgA IgM, and serum IgG2 at the genome-wide level (P < 3.08E-6); 11, 2, 13, 2, 12, 8, 2, 27, 1 and 4 SNPs were found significantly associated with total IgG, IgG1, IgG2, IgA and IgM in colostrum and serum at the suggestive level (P < 6.15E-5). Such SNPs located in or proximate to (±1 Mb) 423 genes, which were functionally implicated in biological processes and pathways, such as immune response, B cell activation, inflammatory response and NF-kappaB signaling pathways. By combining the biological functions and the known QTL data for immune traits in bovine, 14 promising candidate functional genes were identified for immunoglobulin concentrations in colostrum and serum in dairy cattle, they were FGFR4, FGFR2, NCF1, IKBKG, SORBS3, IGHV1S18, KIT, PTGS2, BAX, GRB2, TAOK1, ICAM1, TGFB1 and RAC3. CONCLUSIONS: In this study, we identified 14 candidate genes related to concentrations of immunoglobulins in colostrum and serum in dairy cattle by performing GWASs. Our findings provide a groundwork for unraveling the key genes and causal mutations affecting immunoglobulin concentrations in colostrum and important information for genetic improvement of such traits in dairy cattle.

Single-cell transcriptomic and chromatin accessibility analyses of dairy cattle peripheral blood mononuclear cells and their responses to lipopolysaccharide.

BACKGROUND: Gram-negative bacteria are important pathogens in cattle, causing severe infectious diseases, including mastitis. Lipopolysaccharides (LPS) are components of the outer membrane of Gram-negative bacteria and crucial mediators of chronic inflammation in cattle. LPS modulations of bovine immune responses have been studied before. However, the single-cell transcriptomic and chromatin accessibility analyses of bovine peripheral blood mononuclear cells (PBMCs) and their responses to LPS stimulation were never reported. RESULTS: We performed single-cell RNA sequencing (scRNA-seq) and single-cell sequencing assay for transposase-accessible chromatin (scATAC-seq) in bovine PBMCs before and after LPS treatment and demonstrated that seven major cell types, which included CD4 T cells, CD8 T cells, and B cells, monocytes, natural killer cells, innate lymphoid cells, and dendritic cells. Bioinformatic analyses indicated that LPS could increase PBMC cell cycle progression, cellular differentiation, and chromatin accessibility. Gene analyses further showed significant changes in differential expression, transcription factor binding site, gene ontology, and regulatory interactions during the PBMC responses to LPS. Consistent with the findings of previous studies, LPS induced activation of monocytes and dendritic cells, likely through their upregulated TLR4 receptor. NF-κB was observed to be activated by LPS and an increased transcription of an array of pro-inflammatory cytokines, in agreement that NF-κB is an LPS-responsive regulator of innate immune responses. In addition, by integrating LPS-induced differentially expressed genes (DEGs) with large-scale GWAS of 45 complex traits in Holstein, we detected trait-relevant cell types. We found that selected DEGs were significantly associated with immune-relevant health, milk production, and body conformation traits. CONCLUSION: This study provided the first scRNAseq and scATAC-seq data for cattle PBMCs and their responses to the LPS stimulation to the best of our knowledge. These results should also serve as valuable resources for the future study of the bovine immune system and open the door for discoveries about immune cell roles in complex traits like mastitis at single-cell resolution.

Towards the detection of copy number variation from single sperm sequencing in cattle.

BACKGROUND: Copy number variation (CNV) has been routinely studied using bulk-cell sequencing. However, CNV is not well studied on the single-cell level except for humans and a few model organisms. RESULTS: We sequenced 143 single sperms of two Holstein bulls, from which we predicted CNV events using 14 single sperms with deep sequencing. We then compared the CNV results derived from single sperms with the bulk-cell sequencing of one bull's family trio of diploid genomes. As a known CNV hotspot, segmental duplications were also predicted using the bovine ARS-UCD1.2 genome. Although the trio CNVs validated only some single sperm CNVs, they still showed a distal chromosomal distribution pattern and significant associations with segmental duplications and satellite repeats. CONCLUSION: Our preliminary results pointed out future research directions and highlighted the importance of uniform whole genome amplification, deep sequence coverage, and dedicated software pipelines for CNV detection using single cell sequencing data.