Kimma: flexible linear mixed effects modeling with kinship covariance for RNA-seq data.

MOTIVATION: The identification of differentially expressed genes (DEGs) from transcriptomic datasets is a major avenue of research across diverse disciplines. However, current bioinformatic tools do not support covariance matrices in DEG modeling. Here, we introduce kimma (Kinship In Mixed Model Analysis), an open-source R package for flexible linear mixed effects modeling including covariates, weights, random effects, covariance matrices, and fit metrics. RESULTS: In simulated datasets, kimma detects DEGs with similar specificity, sensitivity, and computational time as limma unpaired and dream paired models. Unlike other software, kimma supports covariance matrices as well as fit metrics like Akaike information criterion (AIC). Utilizing genetic kinship covariance, kimma revealed that kinship impacts model fit and DEG detection in a related cohort. Thus, kimma equals or outcompetes current DEG pipelines in sensitivity, computational time, and model complexity. AVAILABILITY AND IMPLEMENTATION: Kimma is freely available on GitHub https://github.com/BIGslu/kimma with an instructional vignette at https://bigslu.github.io/kimma_vignette/kimma_vignette.html.

Changes in the host transcriptome and microbial metatranscriptome of the ileum of dairy calves subjected to artificial dosing of exogenous rumen contents.

Development of a properly functioning gastrointestinal tract (GIT) at an early age is critical for the wellbeing and lifetime productivity of dairy cattle. The role of early microbial colonization on GIT development in neonatal cattle and the associated molecular changes remain largely unknown, particularly for the small intestine. In this study, we performed artificial dosing of exogenous rumen fluid during the early life of the calf, starting at birth through the weaning transition at 8 wk. Six calves were included in this study. At 8 wk of age, tissue from the ileum was collected and subjected to host transcriptome and microbial metatranscriptome analysis using RNA sequencing. A total of 333 genes showed significant differential expression (DE) (fold-change ≥2; adjusted P < 0.1, mean read-count ≥10) between the treated and control calves. Gene ontology analysis indicated that these DE genes are predominantly associated with processes related to the host immune response (P < 0.0001). Association analysis between the host gene expression and the microbial genus abundance identified 57 genes as having significant correlation with the ileum microbial genera (P < 0.0001). Of these, three genes showed significant association with six microbial genera: lysozyme 2 (LYZ2), fatty acid binding protein 5 (FABP5), and fucosyltransferase (FUT1). Specifically, the profound increase in expression of LYZ2 in treated calves suggests the initiation of antibacterial activity and innate response from the host. Despite the limitation of a relatively small sample size, this study sheds light on the potential impact of early introduction of microbes on the small intestine of calves.

Validating the Use of Bovine Buccal Sampling as a Proxy for the Rumen Microbiota by Using a Time Course and Random Forest Classification Approach.

Analysis of the cow microbiome, as well as host genetic influences on the establishment and colonization of the rumen microbiota, is critical for development of strategies to manipulate ruminal function toward more efficient and environmentally friendly milk production. To this end, the development and validation of noninvasive methods to sample the rumen microbiota at a large scale are required. In this study, we further optimized the analysis of buccal swab samples as a proxy for direct bacterial samples of the rumen of dairy cows. To identify an optimal time for sampling, we collected buccal swab and rumen samples at six different time points relative to animal feeding. We then evaluated several biases in these samples using a machine learning classifier (random forest) to select taxa that discriminate between buccal swab and rumen samples. Differences in the inverse Simpson's diversity, Shannon's evenness, and Bray-Curtis dissimilarities between methods were significantly less apparent when sampling was performed prior to morning feeding (P < 0.05), suggesting that this time point was optimal for representative sampling. In addition, the random forest classifier was able to accurately identify nonrumen taxa, including 10 oral and putative feed-associated taxa. Two highly prevalent (>60%) taxa in buccal and rumen samples had significant variance in relative abundances between sampling methods but could be qualitatively assessed via regular buccal swab sampling. This work not only provides new insights into the oral community of ruminants but also further validates and refines buccal swabbing as a method to assess the rumen bacterial in large herds.IMPORTANCE The gastrointestinal tracts of ruminants harbor a diverse microbial community that coevolved symbiotically with the host, influencing its nutrition, health, and performance. While the influence of environmental factors on rumen microbes is well documented, the process by which host genetics influences the establishment and colonization of the rumen microbiota still needs to be elucidated. This knowledge gap is due largely to our inability to easily sample the rumen microbiota. There are three common methods for rumen sampling but all of them present at least one disadvantage, including animal welfare, sample quality, labor, and scalability. The development and validation of noninvasive methods, such as buccal swabbing, for large-scale rumen sampling is needed to support studies that require large sample sizes to generate reliable results. The validation of buccal swabbing will also support the development of molecular tools for the early diagnosis of metabolic disorders associated with microbial changes in large herds.

Feeding modes shape the acquisition and structure of the initial gut microbiota in newborn lambs.

Early gut microbial colonization is important for postnatal metabolic and immune development. However, little is known about the effects of different feeding modes (suckling versus bottle-feeding) or microbial sources on this process in farm animals. We found that suckled and bottle-fed newborn lambs had their own distinct gut microbiota. Results from 16S rRNA gene sequencing and qPCR showed that, compared with suckling, bottle feeding significantly increased the abundances of Escherichia/Shigella, Butyricicoccus, and Clostridium XlVa, while significantly decreased the abundance of Clostridium XI. The higher levels of Escherichia/Shigella in bottle-fed lambs suggest that artificial feeding may increase the number of potential pathogens and delay the establishment of the anaerobic environment and anaerobic microbes. Feeding modes also affected the direct transmission of bacteria from the mother and the environment to newborns. The SourceTracker analysis estimated that the early gut microbes of suckled lambs were mainly derived from the mother's teats (43%) and ambient air (28%); whereas those of bottle-fed lambs were dominated by bacteria from the mother's vagina (46%), ambient air (31%), and the sheep pen floor (12%). These findings advance our understanding of gut microbiota in early life and may help design techniques to improve gut microbiota and health.

Transient changes in milk production efficiency and bacterial community composition resulting from near-total exchange of ruminal contents between high- and low-efficiency Holstein cows.

The objectives of this study were to determine if milk production efficiency (MPE) is altered by near-total exchange of ruminal contents between high- (HE) and low-MPE (LE) cows and to characterize ruminal bacterial community composition (BCC) before exchange and over time postexchange. Three pairs of ruminally cannulated, third-lactation cows were selected whose MPE (energy-corrected milk per unit of dry matter intake) differed over their first 2 lactations. Approximately 95% of ruminal contents were exchanged between cows within each pair. Ruminal pH and volatile fatty acid (VFA) profiles, along with BCC (characterized by sequencing of the variable 4 region of 16S rRNA genes), were assessed just before feeding on d -8, -7, -5, -4, -1, 1, 2, 3, 7, 10, 14, 21, 28, 35, 42, and 56, relative to the exchange date. High-MPE cows had higher total ruminal VFA concentrations, higher molar percentages of propionate and valerate, and lower molar percentages of acetate and butyrate than did LE cows, and re-established these differences 1 d after contents exchange. Across all LE cows, MPE increased during 7 d postexchange but declined thereafter. Two of the 3 HE cows displayed decreases in MPE following introduction of the ruminal contents from the corresponding LE cow, but MPE increased in the third HE cow, which was determined to be an outlier. For all 6 cows, both liquid- and solids-associated BCC differed between individuals within a pair before contents exchange. Upon exchange, BCC of both phases in all 3 pairs was more similar to that of the donor inoculum than to preexchange host BCC. For 5 of 6 cows, the solids-associated community returned within 10 d to more resemble the preexchange community of that host than that of the donor community. Individual variability before the exchange was greater in liquids than in solids, as was the variability in response of bacterial communities to the exchange. Individual cows varied in their response, but generally moved toward re-establishment of their preexchange communities by 10 d after contents exchange. By contrast, ruminal pH and VFA profiles returned to preexchange levels within 1 d. Despite the small number of cows studied, the data suggest an apparent role for the ruminal bacterial community as a determinant of MPE.

Multi-omics in nasal epithelium reveals three axes of dysregulation for asthma risk in the African Diaspora populations.

Asthma has striking disparities across ancestral groups, but the molecular underpinning of these differences is poorly understood and minimally studied. A goal of the Consortium on Asthma among African-ancestry Populations in the Americas (CAAPA) is to understand multi-omic signatures of asthma focusing on populations of African ancestry. RNASeq and DNA methylation data are generated from nasal epithelium including cases (current asthma, N = 253) and controls (never-asthma, N = 283) from 7 different geographic sites to identify differentially expressed genes (DEGs) and gene networks. We identify 389 DEGs; the top DEG, FN1, was downregulated in cases (q = 3.26 × 10-9) and encodes fibronectin which plays a role in wound healing. The top three gene expression modules implicate networks related to immune response (CEACAM5; p = 9.62 × 10-16 and CPA3; p = 2.39 × 10-14) and wound healing (FN1; p = 7.63 × 10-9). Multi-omic analysis identifies FKBP5, a co-chaperone of glucocorticoid receptor signaling known to be involved in drug response in asthma, where the association between nasal epithelium gene expression is likely regulated by methylation and is associated with increased use of inhaled corticosteroids. This work reveals molecular dysregulation on three axes - increased Th2 inflammation, decreased capacity for wound healing, and impaired drug response - that may play a critical role in asthma within the African Diaspora.

A quality improvement study on the relationship between intranasal povidone-iodine and anesthesia and the nasal microbiota of surgery patients.

INTRODUCTION: The composition of the nasal microbiota in surgical patients in the context of general anesthesia and nasal povidone-iodine decolonization is unknown. The purpose of this quality improvement study was to determine: (i) if general anesthesia is associated with changes in the nasal microbiota of surgery patients and (ii) if preoperative intranasal povidone-iodine decolonization is associated with changes in the nasal microbiota of surgery patients. MATERIALS AND METHODS: One hundred and fifty-one ambulatory patients presenting for surgery were enrolled in a quality improvement study by convenience sampling. Pre- and post-surgery nasal samples were collected from patients in the no intranasal decolonization group (control group, n = 54). Pre-decolonization nasal samples were collected from the preoperative intranasal povidone-iodine decolonization group (povidone-iodine group, n = 97). Intranasal povidone-iodine was administered immediately prior to surgery and continued for 20 minutes before patients proceeded for surgery. Post-nasal samples were then collected. General anesthesia was administered to both groups. DNA from the samples was extracted for 16S rRNA sequencing on an Illumina MiSeq. RESULTS: In the control group, there was no evidence of change in bacterial diversity between pre- and post-surgery samples. In the povidone-iodine group, nasal bacterial diversity was greater in post-surgery, relative to pre-surgery (Shannon's Diversity Index (P = 0.038), Chao's richness estimate (P = 0.02) and Inverse Simpson index (P = 0.027). Among all the genera, only the relative abundance of the genus Staphylococcus trended towards a decrease in patients after application (FDR adjusted P = 0.06). Abundant genera common to both povidone-iodine and control groups included Staphylococcus, Bradyrhizobium, Corynebacterium, Dolosigranulum, Lactobacillus, and Moraxella. CONCLUSIONS: We found general anesthesia was not associated with changes in the nasal microbiota. Povidone-iodine treatment was associated with nasal microbial diversity and decreased abundance of Staphylococcus. Future studies should examine the nasal microbiota structure and function longitudinally in surgical patients receiving intranasal povidone-iodine.

Assessing the Response of Ruminal Bacterial and Fungal Microbiota to Whole-Rumen Contents Exchange in Dairy Cows.

A major goal for the dairy industry is to improve overall milk production efficiency (MPE). With the advent of next-generation sequencing and advanced methods for characterizing microbial communities, efforts are underway to improve MPE by manipulating the rumen microbiome. Our previous work demonstrated that a near-total exchange of whole rumen contents between pairs of lactating Holstein dairy cows of disparate MPE resulted in a reversal of MPE status for ∼10 days: historically high-efficiency cows decreased in MPE, and historically low-efficiency cows increased in MPE. Importantly, this switch in MPE status was concomitant with a reversal in the ruminal bacterial microbiota, with the newly exchanged bacterial communities reverting to their pre-exchange state. However, this work did not include an in-depth analysis of the microbial community response or an interrogation of specific taxa correlating to production metrics. Here, we sought to better understand the response of rumen communities to this exchange protocol, including consideration of the rumen fungi. Rumen samples were collected from 8 days prior to, and 56 days following the exchange and were subjected to 16S rRNA and ITS amplicon sequencing to assess bacterial and fungal community composition, respectively. Our results show that the ruminal fungal community did not differ significantly between hosts of disparate efficiency prior to the exchange, and no change in community structure was observed over the time course. Correlation of microbial taxa to production metrics identified one fungal operational taxonomic unit (OTU) in the genus Neocallimastix that correlated positively to MPE, and several bacterial OTUs classified to the genus Prevotella. Within the Prevotella, Prevotella_1 was found to be more abundant in high-efficiency cows whereas Prevotella_7 was more abundant in low-efficiency cows. Overall, our results suggest that the rumen bacterial community is a primary microbial driver of host efficiency, that the ruminal fungi may not have as significant a role in MPE as previously thought, and that more work is needed to better understand the functional roles of specific ruminal microbial community members in modulating MPE.

The Effects of Artificially Dosed Adult Rumen Contents on Abomasum Transcriptome and Associated Microbial Community Structure in Calves.

This study aimed to investigate the changes in abomasum transcriptome and the associated microbial community structure in young calves with artificially dosed, adult rumen contents. Eight young bull calves were randomly dosed with freshly extracted rumen contents from an adult cow (high efficiency (HE), n = 4), or sterilized rumen content (Con, n = 4). The dosing was administered within 3 days of birth, then at 2, 4, and 6 weeks following the initial dosing. Abomasum tissues were collected immediately after sacrifice at 8 weeks of age. Five genera (Tannerella, Desulfovibrio, Deinococcus, Leptotrichia, and Eubacterium; p < 0.05) showed significant difference in abundance between the treatments. A total of 975 differentially expressed genes were identified (p < 0.05, fold-change > 1.5, mean read-counts > 5). Pathway analysis indicated that up-regulated genes were involved in immune system process and defense response to virus, while the down-regulated genes involved in ion transport, ATP biosynthetic process, and mitochondrial electron transport. Positive correlation (r > 0.7, p < 0.05) was observed between TRPM4 gene and Desulfovibrio, which was significantly higher in the HE group. TRPM4 had a reported role in the immune system process. In conclusion, the dosing of adult rumen contents to calves can alter not only the composition of active microorganisms in the abomasum but also the molecular mechanisms in the abomasum tissue, including reduced protease secretion and decreased hydrochloric acid secretion.

Transcriptomics analysis of host liver and meta-transcriptome analysis of rumen epimural microbial community in young calves treated with artificial dosing of rumen content from adult donor cow.

In mammals, microbial colonization of the digestive tract (GIT) occurs right after birth by several bacterial phyla. Numerous human and mouse studies have reported the importance of early gut microbial inhabitants on host health. However, few attempts have been undertaken to directly interrogate the role of early gut/rumen microbial colonization on GIT development or host health in neonatal ruminants through artificial manipulation of the rumen microbiome. Thus, the molecular changes associated with bacterial colonization are largely unknown in cattle. In this study, we dosed young calves with exogenous rumen fluid obtained from an adult donor cow, starting at birth, and repeated every other week until six weeks of age. Eight Holstein bull calves were included in this study and were separated into two groups of four: the first group was treated with rumen content freshly extracted from an adult cow, and the second group was treated with sterilized rumen content. Using whole-transcriptome RNA-sequencing, we investigated the transcriptional changes in the host liver, which is a major metabolic organ and vital to the calf's growth performance. Additionally, the comparison of rumen epimural microbial communities between the treatment groups was performed using the rRNA reads generated by sequencing. Liver transcriptome changes were enriched with genes involved in cell signaling and protein phosphorylation. Specifically, up-regulation of SGPL1 suggests a potential increase in the metabolism of sphingolipids, an essential molecular signal for bacterial survival in digestive tracts. Notably, eight genera, belonging to four phyla, had significant increases in abundance in treated calves. Our study provides insight into host liver transcriptome changes associated with early colonization of the microbial communities in neonatal calves. Such knowledge provides a foundation for future probiotics-based research in microbial organism mediated rumen development and nutrition in ruminants.

Author Correction: Transcriptomics analysis of host liver and meta-transcriptome analysis of rumen epimural microbial community in young calves treated with artificial dosing of rumen content from adult donor cow.

A correction to this article has been published and is linked from the HTML and PDF versions of this paper. The error has been fixed in the paper.