BackgroundLittle is known about the roles of myeloid cell subsets in kidney injury and in the limited ability of the organ to repair itself. Characterizing these cells based only on surface markers using flow cytometry might not provide a full phenotypic picture. Defining these cells at the single-cell, transcriptomic level could reveal myeloid heterogeneity in the progression and regression of kidney disease.MethodsIntegrated droplet– and plate-based single-cell RNA sequencing were used in the murine, reversible, unilateral ureteric obstruction model to dissect the transcriptomic landscape at the single-cell level during renal injury and the resolution of fibrosis. Paired blood exchange tracked the fate of monocytes recruited to the injured kidney.ResultsA single-cell atlas of the kidney generated using transcriptomics revealed marked changes in the proportion and gene expression of renal cell types during injury and repair. Conventional flow cytometry markers would not have identified the 12 myeloid cell subsets. Monocytes recruited to the kidney early after injury rapidly adopt a proinflammatory, profibrotic phenotype that expresses Arg1, before transitioning to become Ccr2+ macrophages that accumulate in late injury. Conversely, a novel Mmp12+ macrophage subset acts during repair.ConclusionsComplementary technologies identified novel myeloid subtypes, based on transcriptomics in single cells, that represent therapeutic targets to inhibit progression or promote regression of kidney disease.

Highlights d Primary senescence and secondary senescence are distinct molecular endpoints d Secondary Ras-induced senescence has a composite SASP, Notch-induced signature d Notch signaling is an essential driver of secondary senescence d Notch blunts the senescence-associated secretory phenotype in secondary senescence

Water buffalo is the second largest resource of milk supply around the world, and it is well known for its distinctive milk quality in terms of fat, protein, lactose, vitamin, and mineral contents. Understanding the genetic architecture of milk production traits is important for future improvement by the buffalo breeding industry. The advance of genome-wide association studies (GWAS) provides an opportunity to identify potential genetic variants affecting important economical traits. In the present study, GWAS was performed for 489 buffaloes with 1,424 lactation records using the 90K Affymetrix Buffalo SNP Array (Affymetrix/Thermo Fisher Scientific, Santa Clara, CA). Collectively, 4 candidate single nucleotide polymorphisms (SNP) in 2 genomic regions were found to associate with buffalo milk production traits. One region affecting milk fat and protein percentage was located on the equivalent of Bos taurus autosome (BTA)3, spanning 43.3 to 43.8 Mb, which harbored the most likely candidate genes MFSD14A, SLC35A3, and PALMD. The other region on the equivalent of BTA14 at 66.5 to 67.0 Mb contained candidate genes RGS22 and VPS13B and influenced buffalo total milk yield, fat yield, and protein yield. Interestingly, both of the regions were reported to have quantitative trait loci affecting milk performance in dairy cattle. Furthermore, we suggest that buffaloes with the C allele at AX-85148558 and AX-85073877 loci and the G allele at AX-85106096 locus can be selected to improve milk fat yield in this buffalo-breeding program. Meanwhile, the G allele at AX-85063131 locus can be used as the favorable allele for improving milk protein percentage. Genomic prediction showed that the reliability of genomic estimated breeding values (GEBV) of 6 milk production traits ranged from 0.06 to 0.22, and the correlation between estimated breeding values and GEBV ranged from 0.23 to 0.35. These findings provide useful information to understand the genetic basis of buffalo milk properties and may play a role in accelerating buffalo breeding programs using genomic approaches.