The maturation and aging trajectory of Marchantia polymorpha at single-cell resolution.

Bryophytes represent a sister to the rest of land plants. Despite their evolutionary importance and relatively simple body plan, a comprehensive understanding of the cell types and transcriptional states that underpin the temporal development of bryophytes has not been achieved. Using time-resolved single-cell RNA sequencing, we define the cellular taxonomy of Marchantia polymorpha across asexual reproduction phases. We identify two maturation and aging trajectories of the main plant body of M. polymorpha at single-cell resolution: the gradual maturation of tissues and organs along the tip-to-base axis of the midvein and the progressive decline of meristem activities in the tip along the chronological axis. Specifically, we observe that the latter aging axis is temporally correlated with the formation of clonal propagules, suggesting an ancient strategy to optimize allocation of resources to producing offspring. Our work thus provides insights into the cellular heterogeneity that underpins the temporal development and aging of bryophytes.

Multiple introductions and recombination events underlie the emergence of a hyper-transmissible Cryptosporidium hominis subtype in the USA.

The parasite Cryptosporidium hominis is a leading cause of the diarrheal disease cryptosporidiosis, whose incidence in the United States has increased since 2005. Here, we show that the newly emerged and hyper-transmissible subtype IfA12G1R5 is now dominant in the United States. In a comparative analysis of 127 newly sequenced and 95 published C. hominis genomes, IfA12G1R5 isolates from the United States place into three of the 14 clusters (Pop6, Pop13, and Pop14), indicating that this subtype has multiple ancestral origins. Pop6 (IfA12G1R5a) has an East Africa origin and has recombined with autochthonous subtypes after its arrival. Pop13 (IfA12G1R5b) is imported from Europe, where it has recombined with the prevalent local subtype, whereas Pop14 (IfA12G1R5c) is a progeny of secondary recombination between Pop6 and Pop13. Selective sweeps in invasion-associated genes have accompanied the emergence of the dominant Pop14. These observations offer insights into the emergence and evolution of hyper-transmissible pathogens.

FindIT2: an R/Bioconductor package to identify influential transcription factor and targets based on multi-omics data.

BACKGROUND: Transcription factors (TFs) play central roles in regulating gene expression. With the rapid growth in the use of high-throughput sequencing methods, there is a need to develop a comprehensive data processing and analyzing framework for inferring influential TFs based on ChIP-seq/ATAC-seq datasets. RESULTS: Here, we introduce FindIT2 (Find Influential TFs and Targets), an R/Bioconductor package for annotating and processing high-throughput multi-omics data. FindIT2 supports a complete framework for annotating ChIP-seq/ATAC-seq peaks, identifying TF targets by the combination of ChIP-seq and RNA-seq datasets, and inferring influential TFs based on different types of data input. Moreover, benefited from the annotation framework based on Bioconductor, FindIT2 can be applied to any species with genomic annotations, which is particularly useful for the non-model species that are less well-studied. CONCLUSION: FindIT2 provides a user-friendly and flexible framework to generate results at different levels according to the richness of the annotation information of user's species. FindIT2 is compatible with all the operating systems and is released under Artistic-2.0 License. The source code and documents are freely available through Bioconductor ( https://bioconductor.org/packages/devel/bioc/html/FindIT2.html ).

Age and episode-associated occurrence of Cryptosporidium species and subtypes in a birth-cohort of dairy calves.

The role of species-specific immunity in infection patterns of Cryptosporidium spp. in humans and farm animals is not well understood. In the present study, the dynamics of Cryptosporidium infections in a natural cryptosporidiosis model was examined using genotyping, subtyping and whole genome sequencing tools. In a cross-sectional survey of Cryptosporidium spp. in 934 dairy cattle on one farm, marked age-associated differences in the distribution of Cryptosporidium species and C. bovis subtypes were observed. In a closely followed longitudinal birth cohort study of 81 calves over a 9-month period, shedding of C. parvum oocysts by the IIdA19G1 subtype started at 4 days, peaked at 2 weeks and ended mostly by 4 weeks. In contrast, the shedding of C. bovis oocysts started at 2 weeks, peaked initially at 6 weeks and had a second wave during 15th to 23rd weeks. For C. ryanae, calves had mostly only one episode of infection by one subtype, with accumulative infection increasing much slower than C. parvum and C. bovis. Overall, the accumulative infection rates and mean duration of oocyst shedding for calves in the cohort were 97.4% (76/78) and 2.3 weeks, 100.0% (80/80) and 3.9 weeks, and 78.7% (63/80) and 3.2 weeks for C. parvum, C. bovis and C. ryanae, respectively. The oocyst shedding intensity was much lower in C. bovis and C. ryanae infections compared with C. parvum infection, and in the second episode of C. bovis infection compared with the first episode. The two episodes of C. bovis infections were caused by different genome types that differed mostly in nine genes. Cryptosporidium parvum infection was associated with the occurrence of watery diarrhoea. Data from the natural history study of cryptosporidiosis indicate that despite the existence of acquired immunity against homologous pathogens, neonatal animals experience waves of Cryptosporidium infections by different species and genome types.

Dynamic chromatin state profiling reveals regulatory roles of auxin and cytokinin in shoot regeneration.

Shoot regeneration is mediated by the sequential action of two phytohormones, auxin and cytokinin. However, the chromatin regulatory landscapes underlying this dynamic response have not yet been studied. In this study, we jointly profiled chromatin accessibility, histone modifications, and transcriptomes to demonstrate that a high auxin/cytokinin ratio environment primes Arabidopsis shoot regeneration by increasing the accessibility of the gene loci associated with pluripotency and shoot fate determination. Cytokinin signaling not only triggers the commitment of the shoot progenitor at later stages but also allows chromatin to maintain shoot identity genes at the priming stage. Our analysis of transcriptional regulatory dynamics further identifies a catalog of regeneration cis-elements dedicated to cell fate transitions and uncovers important roles of BES1, MYC, IDD, and PIF transcription factors in shoot regeneration. Our results, thus, provide a comprehensive resource for studying cell reprogramming in plants and provide potential targets for improving future shoot regeneration efficiency.

Development and Application of a gp60-Based Subtyping Tool for Cryptosporidium bovis.

Cryptosporidium bovis is a common enteric pathogen in bovine animals. The research on transmission characteristics of the pathogen is hampered by the lack of subtyping tools. In this study, we retrieve the nucleotide sequence of the 60 kDa glycoprotein (GP60) from the whole genome sequences of C. bovis we obtained previously and analyze its sequence characteristics. Despite a typical structure of the GP60 protein, the GP60 of C. bovis had only 19.3-45.3% sequence identity to those of other Cryptosporidium species. On the basis of the gene sequence, a subtype typing tool was developed for C. bovis and used in the analysis of 486 C. bovis samples from dairy cattle, yaks, beef cattle, and water buffalos from China. Sixty-eight sequence types were identified from 260 subtyped samples, forming six subtype families, namely XXVIa to XXVIf. The mosaic sequence patterns among subtype families and the 121 potential recombination events identified among the sequences both suggest the occurrence of genetic recombination at the locus. No obvious host adaptation and geographic differences in the distribution of subtype families were observed. Most farms with more extensive sampling had more than one subtype family, and the dominant subtype families on a farm appeared to differ between pre- and post-weaned calves, indicating the likely occurrence of multiple episodes of C. bovis infections. There was an association between XXVId infection and occurrence of moderate diarrhea in dairy cattle. The subtyping tool developed and the data generated in the study might improve our knowledge of the genetic diversity and transmission of C. bovis.

Genus-level evolutionary relationships of FAR proteins reflect the diversity of lifestyles of free-living and parasitic nematodes.

BACKGROUND: Nematodes are a widespread and diverse group comprising free-living and parasitic species, some of which have major detrimental effects on crops, animals, and human health. Genomic comparisons of nematodes may help reveal the genetic bases for the evolution of parasitic lifestyles. Fatty acid and retinol-binding proteins (FARs) are thought to be unique to nematodes and play essential roles in their development, reproduction, infection, and possibly parasitism through promoting the uptake, transport, and distribution of lipid and retinol. However, the evolution of FAR family proteins across the phylum Nematoda remains elusive. RESULTS: We report here the evolutionary relationship of the FAR gene family across nematodes. No FAR was found in Trichocephalida species and Romanomermis culicivorax from Clade I, and FAR could be found in species from Clades III, IV, and V. FAR proteins are conserved in Clade III species and separated into three clusters. Tandem duplications and high divergence events lead to variable richness and low homology of FARs in Steinernema of Clade IVa, Strongyloides of Clade IVb, and intestinal parasitic nematodes from Clades Vc and Ve. Moreover, different richness and sequence variations of FARs in pine wood, root-knot, stem, and cyst nematodes might be determined by reproduction mode or parasitism. However, murine lungworm Angiostrongylus and bovine lungworm Dictyocaulus viviparus from Clade Vd have only 3-4 orthologs of FAR. RNA-seq data showed that far genes, especially far-1 and far-2, were highly expressed in most nematodes. Angiostrongylus cantonensis FAR-1 and FAR-3 have low sequence homology and distinct ligand-binding properties, leading to differences in the cavity volume of proteins. These data indicate that FAR proteins diverged early and experienced low selective pressure to form genus-level diversity. The far genes are present in endophyte or root-colonized bacteria of Streptomyces, Kitasatospora sp., Bacillus subtilis, and Lysobacter, suggesting that bacterial far genes might be derived from plant-parasitic nematodes by horizontal gene transfer. CONCLUSIONS: Data from these comparative analyses have provided insights into genus-level diversity of FAR proteins in the phylum Nematoda. FAR diversification provides a glimpse into the complicated evolution history across free-living and parasitic nematodes.