Bovine Mastitis disease is a silent killer of the dairy industry in the whole world. The overuse of antibiotics to treat the disease has rendered the disease causing pathogens to be antibiotic resistant and can also induce antibiotic resistance in the human gut flora and microbiome. We realised that a non-antibiotic treatment method is the need of the hour to reduce the uncontrolled antibiotic usage in the dairy industry. To implement our non-antibiotic treatment technique it is also extremely important to devise an early subclinical bovine mastitis detection kit so that the disease gets detected in a stage where the disease is curable.
Our genetic circuit has been judiciously designed to suffice all our needs as discussed on the Description page. Our main aim of designing the GMO was such that it doesn’t harm the cows and also is a specific and effective therapeutic for this undetectable fatal disease.
Before deciding on our chassis organism we listed down the features that our chassis should fulfill such that it becomes a more effective tool to deliver our planned therapeutic with specificity and effectivity. The features that we required in our chassis are -
- It should be a gram-positive bacteria - We plan to incorporate the agrC-agrA quorum-sensing unit of Staphylococcus aureus into our chassis such that it can sense the presence of AIP-I molecules. Staphylococcus aureus itself is a gram-positive bacteria. Its quorum-sensing unit will be more easily expressed on another gram-positive bacteria than a gram-negative bacteria due to the absence of an outer membrane like gram-negative bacteria.
-
It should be a Nisin producing bacteria - The nisin biosynthetic gene cluster is essential for posttranslational modification of the nisin precursor molecule, self Immunity against Nisin, translocation, and cleaving of the leader peptide. All these processes are essential for the production of active nisin. Other than cloning a gene cluster would be difficult as that would make our designed plasmid unstable. So we need a Nisin-producing bacterial strain.
- It should be Nis T mutant - As per our mode of targeting the causative pathogens (degrade the biofilm and kill them using bacteriocins), it is needed to release the Nisin PV and the DNaseI at the same time such that we can overcome the protective effect of biofilm on the producing bacteria.
- It should be safe for the cows - We emphasized finding a probiotic strain that can be injected into the udder tissue without causing any harm to the cow.
Based on all these requirements, we chose Lactococcus Lactis LMG 7930 as our chassis organism. The reasons for selecting this strain are as follows -
- This is a probiotic strain and is already in use for treating bovine mastitis through intramammary injection.
- This is a gram-positive Nisin-producing strain but with NisT gene mutation.
Based on our directed motive, method, and approach of treating subclinical bovine mastitis our treatment section can be divided into two parts for better understanding -
- Sensing
- Killing
The end effects of broad-spectrum antibiotics usage for treating bovine mastitis has directed us to design our therapeutic to be narrow-spectrum and specific for the pathogens. We plan to incorporate a particular segment of agr quorum-sensing system. Staphylococcus aureus uses AIP-I molecules for intraspecies communication. Hence our main goal is to design a genetic circuit so that our GMO can sense the presence of AIP-I molecules in its surroundings and then regulate the downstream signaling. We designed a regulatory module that consists of agrC and agrA under a constitutive expression using pBAD promoter with a strong ribosomal binding site [B0034] and a double terminator [B0015].
We added a sensing module downstream of the regulatory module consisting of promoter P2 which is activated by agrA.
To test this sensing module, we incorporated GFP as a reporter gene downstream of the P2 promoter such that in the presence of AIP-I molecules P2 promoter can express GFP. So, we externally added AIP-I molecules to observe if the module is sensing its presence or not (details).
After analyzing the fluorescence data, we observed that our sensing module showed high fluorescence levels even in the absence of AIP-I molecules. Based on this observation, we inferred that this could be due to the leaky expression of GFP.
Now to tackle this problem, we designed a bi-directional sensing module such that it can sense the AIP-I molecules and reduce the leaky gene expression of GFP. Post these changes we observed a significant decrease in gfp expression. See results here.
Hence we decided to use this bi-directional sensing module in our proposed genetic circuit.
We have not only designed the killing of the causative pathogens but have also targeted their biofilm. The design has been such that the Nisin PV and the DNase I are released at once when the GMO is lysed by the lysis proteins.
Bacteriocin of Choice -
Bacteriocins are narrow-spectrum antimicrobials that will suffice our need to kill the targeted pathogens effectively. In the initial stages of designing the treatment part of our project, we chose NisinA bacteriocin to target Staphylococcus aureus specifically. Nisin A is a lantibiotic belonging to class I bacteriocin. It consists of 34 amino acids. It kills the pathogenic bacterial cells by pore formation. Initially, it docks to lipid II, a membrane-anchored cell wall precursor that is essential for bacterial cell wall biosynthesis. After binding, it prevents the peptidoglycan monomer from incorporation into the growing peptidoglycan network. This finally enables the C-terminal segment of nisinA to insert into the cell membrane. These nisin-lipid II complexes assemble to form a stable pore in the cell membrane of target cells. Once a pore is formed, it causes an increase in membrane permeability which leads to the dissipation of the membrane potential. Consequently, these cells cannot synthesize any energy, which ultimately leads to their death [2][3]. Therefore nisinA can very effectively kill the antibiotic-resistant S.aureus strains [4].
After receiving inputs from several experts in this field of study, we understood that to make our treatment approach more effective, we have to be more inclusive in terms of targetting more than one causative pathogen. Literature studies showed that subclinical bovine mastitis is majorly caused by Staphylococcus aureus followed by Streptococcus uberis [5]. Streptococcus uberis expresses a protein called Nisin Resistance Protein (NSR), which proteolytically cleaves Nisin A at the bond between MeLan28 and Ser29 resulting in a truncated nisin (nisin \( ^{1–28}\) ). This cleavage inactivates the Nisin A. Hence, to target both we chose Nisin PV, a bioengineered variant of Nisin A. Mutation at 29th position of Nisin A from Serine to Proline and Isoleucine to Valine at the 30th position. These mutations make this variant of bacteriocin resistant to cleavage by NSR. [6]
Though Nisin A and Ninsin PV have structurally similar stability (details) but Nisin PV is more effective with respect to our designed project. Nisin PV also has phenomenal capabilities of inhibiting biofilm formation and decreasing its viability up to 56% and 85% more than nisin A, respectively [7]. Nisin PV has been placed downstream of the P2 promoter in our genetic circuit.
Targeting the Biofilm -
Our targeted pathogens colonize the host tissues by forming an extracellular matrix, slowly helping in the development of microcolony. A biofilm is further formed by cell division and the accumulation of biomass. The biofilm protects the pathogenic cells from the bacteriocin as the biofilm shows tolerance. So we incorporated DNaseI in our genetic circuit. It was placed downstream of Nisin PV under the same P2 promoter. We observed significant degradation of the biofilm by DNase I. More details are on the Results page.
Lysis of the GMO
In the end, we need to lyse our GMO for two main reasons -
- After producing the required amount of our therapeutic molecules, the GMO will lyse itself and die. Thus serving the purpose of a kill switch.
- Both the therapeutic molecules are needed to be released at the same time to increase the effectiveness of our therapeutic approach.
These criteria can be achieved by incorporating the Lysis E7 downstream of the P3 promoter. We designed it to put under the weak P3 promoter, which gets activated by phosphorylated agrA. But the interaction between the promoter and agrA is weak compared to the interaction of P2 [8]. This weak interaction results in delayed production of Lysis proteins. Hence allowing the accumulation of Nisin PV and DNaseI. This lysis of our GMO after being used to its full potential acts as a kill switch too as it ensures the safety of the environment.
Other than that, we have incorporated another kill switch to resist damages pertaining to the accidental release of our GMO into the environment. This kill switch is based on a thymidine mutation which has made our GMO’s life span depend on the amount of thymidine available in its surrounding. Both the kill switches ensure that our proposed therapeutic is safe. More details can be found on the Safety page.
When the AIP-I concentration is above a certain threshold due to the presence of Staphylococcus aureus cells, the AIP-I molecules bind to the agrC. The agrC-AIP-I complex then acts as kinase and phosphorylates the agrA. The phosphorylated agrA then binds to P2 and P3 and activates the production of Nisin PV and DNase I downstream of P2 whereas Lysis E7 downstream of P3. But the production rate of Nisin PV and DNase I will be higher than the production rate of Lysis E7 due to the stronger affinity of phosphorylated agrA towards P2 than P3 promoter. So, our GMO will only lyse after a significant accumulation of Nisin PV and DNaseI.
If the AIP-I concentration is not above a certain threshold, then our GMO will not be switched on.
Due to Lab restrictions and time constraints, we could not clone the whole genetic circuit neither could we test the whole system. We plan to continue with this project and clone the whole genetic circuit and test it.
Detection of subclinical bovine mastitis plays an important role in laying the foundation for the treatment of the disease. Subclinical bovine mastitis is undetectable as it is symptomless, early and systematic detection of milk to identify diseased cows are a must to continue the treatment phase of the afflicted cows. Our method of detection employs the activity of SHERLOCKv2 which consists of Cas13a and Csm6. This is combined with an RNA hydrogel to give a colorimetric readout. We have tried to propose an optimisation of this method in a way that it can be used to detect any disease with viable miRNA biomarkers by coupling SHERLOCKv2 with two kinds of RNA hydrogels, PEG and polyacrylamide, depending on the initial concentrations of the targeted miRNA for the disease. Our detection pathway consists of the following steps:
- Detection of the target miRNA biomarkers by the usage of Cas13a and Csm6 constituting SHERLOCKv2[9].
- RNA-based hydrogel to give a colorimetric readout depending on the activity of the SHERLOCKv2 enzymes.[10]
- Optimisation of the hydrogel to detect multiple diseases having miRNA biomarkers.
Analysis
Our aim was to make a cheap, easy to use and real-time paper based test for the detection of subclinical bovine mastitis as it is one of the most debilitating diseases for cattle in India and hence one of the major reasons for crippling the Indian economy. Literature study showed that detection of bovine mastitis is done by measuring the somatic cell count in a test called the California Mastitis Test[11]. Some tests also monitored the levels of the enzyme N-acetyl-b-D-glucosaminidase (NAGase) in milk[12].
However these biomarkers are very vague and can be present in milk in considerable quantities during lactation periods or during any other infection[3]. These tests are therefore imparted with a relatively larger amount of subjective character with the test results being misread by various cohorts as well as the tests not being a pinpoint indicator of Subclinical Bovine Mastitis. Therefore, we decided to narrow down the most specific biomarker of bovine mastitis - microRNA that is present in the milk. These microRNAs that we intend to target are upregulated significantly in mastitis positive milk and are zero/not present in normal or non-mastitic milk. After a lot of analysis we decided on using SHERLOCK or CRISPR Cas13 to target these specific microRNA molecules. However the cleavage activity of Cas13 was enhanced by preamplification that involved recombinase polymerase amplification (RPA). We wanted to create a test kit without the process of amplification of the sample ( to reduce number and complication of steps of detection), hence we switched to SHERLOCKv2[9] involving the use of orthogonal CRISPR enzymes Cas13a and Csm6. Also we created an RNA based hydrogel structure that would allow us to have a visual colorimetric detector of the disease. This colorimetric indicator of the hydrogel is dependent on the activity of the enzymes of SHERLOCKv2 which allows us to create a rapid paper based detector of the subclinical bovine mastitis. Later we optimised the structure of the hydrogel, to detect any disease with an upregulation of a miRNA biomarker.
SHERLOCKv2
SHERLOCKv2 is an advancement that was integrated into the old SHERLOCK[13](specific high-sensitivity enzymatic reporter unlocking) by the Zhang lab. Orthogonal CRISPR enzymes Cas13a and Csm6 are multiplexed together to detect RNA molecules. Cas13a belongs to the type VI CRISPR-Cas system[14] whereas Csm6 belongs to the type III CRISPR-Cas system[7]. Cas13a is specifically activated when it encounters an RNA sequence complementary to its guide RNA[13].Csm6 is mainly activated by hex-adenylate and sometimes by tetra-adenylate sequences[15][16]. Both Cas13a and Csm6 carry out collateral cleavage of RNA[9], cleaving whatever RNA they can find in their surroundings, even when it is not complementary to their guide RNA sequence.
SHERLOCKv2 acts by the method of initial activation of cas13a which then goes on to generate cleavage residues/overhangs of 2’-3’ cyclic phosphate hex adenylate by collateral cleavage of RNA in the sample. These overhangs activate Csm6 which also carries out collateral cleavage. According to literature, it was found that this can be used to detect as low as 2 attomolar input RNA, however this process involved preamplification of the sample using Recombinase Polymerase Amplification(RPA)[9].
miRNA biomarkers
Our method of detection for subclinical bovine mastitis relies on the stark upregulation of certain miRNA biomarkers when the disease is caused. The miRNAs targeted are bta mir-184[17] and Bta-miR-222[18]. The miRNA sequences are as follows:
miRNA | Sequence |
---|---|
bta-miR-184 | UGGACGGAGAACUGAUAAGGGU |
bta-miR-222 | AGCUACAUCUGGCUACUGGGU |
Bta-miR-184 is seen to be zero in the non-mastitic milk sample. However this miRNA is expressed only in mastitis that is caused by Staphylococcus aureus[19]. Bta-miR-222 however is seen in very negligible amounts in a non-mastitic sample of milk. This miRNA is upregulated greatly in cases of mastitis that is caused by both Staphylococcus aureus and Streptococcus uberis with a p value <0.001. It was also seen that the latter miRNA had a sensitivity of 94% and a specificity of 94%[20].
The guide RNA sequences for CRISPR Cas13a are made complementary to these miRNA sequences. In a cow with subclinical bovine mastitis these miRNAs are present in the milk[19][20]. This activates Cas13a which in turn activates the Csm6. Both Cas13a and Csm6 together collaterally cleave whatever RNA they can find in the sample. For a non-mastitic cow these target miRNAs are not present(zero concentration)[19][20] in the milk and hence Cas13a and Csm6 do not get activated. Thus no cleavage of the RNA present in the milk sample occurs.
RNA HYDROGEL
In order to have a colorimetric readout without the RNA preamplification step, we have decided to create a RNA based hydrogel. The hydrogel has linker structures that are made of single stranded RNA. If the CRISPR enzymes are activated due to a mastitis positive milk sample, the linker RNAs holding the gel together are cut due to the collateral RNA cleavage activity of Cas13a and Csm6. If the linker RNAs holding the hydrogel together are cut, the gel loses its consistency and becomes watery.
This allows the flow of a colorimetric buffer solution to denote a positive test result. On the contrary, a mastitic negative sample not containing the target miRNAs will not trigger the activation of the CRISPR enzymes and the linker RNAs are intact. They hold the gel together and due to the consistency of the gel, the buffer cannot pass through[21].
Paper Based Test strip - μpad
The μPAD or microfluidic test pad is made by Prof. James Collins’ lab at MIT[13]. The pad consists of a filter paper coated with hydrophobic wax in such a manner that four spots and one channel is left without wax. These areas are hydrophilic and allow water based coloured buffers to pass through. The strip is then folded in a manner to line the hydrophilic spots one above the other by folding. This creates a stack and a continuous flow channel that ends in the final flow channel[21]
The buffer is loaded onto the first spot of the μPAD. The sample is added to the second spot. The gel goes into the third spot and the fourth spot contains coloured dyes. If the sample contains the target miRNAs, the CRISPR enzymes are activated. This cleaves the linker RNAs holding the gel together which makes the hydrogel lose its consistency. The buffer then flows through the gel and takes the coloured dyes with it. This gives a positive readout. On the other hand if the sample is negative and does not contain the target miRNA, the Cas enzymes are never activated. This leaves the linker RNA intact and holds the gel together. Hence the buffer cannot flow through and this gives a negative readout.
Optimisation of the test:
In order to make the test more inclusive, we decided to optimise it for all diseases that have miRNA biomarkers. In this case we decided to create two kinds of RNA based hydrogel - PEG and polyacrylamide based. They are used as follows:
- When the target miRNA biomarker for a disease is present in very preliminary amounts in normal samples but is upregulated multiple folds during the disease. In such cases, even normal samples will give a positive result as the cas enzymes will get activated. Hence in these cases we will use the PEG hydrogel synthesis method.
- When the target miRNA is zero in normal samples but is greatly upregulated in diseased samples. In that case, polyacrylamide RNA based hydrogel is used.
As the initial concentration of the miRNAs that we have targeted for our detection process is zero, therefore for the testing of subclinical bovine mastitis, we have employed the polyacrylamide RNA based hydrogel. However we have also explained the working reactions and method of usage of the PEG hydrogel as well. Weplan to experimentally test our proposal in due time.
Polyacrylamide based hydrogel [10]
Polyacrylamide RNA based hydrogel is used when the initial concentration of the targeted miRNA in normal condition is zero. However these miRNAs are significantly upregulated during diseased conditions.
It is made by polymerising acrylamide using APS/TEMED. Two such polyacrylamide backbones are made and each of these backbones are functionalized with a different RNA motif. These RNA motifs contain an initial methacryl group such that it can be incorporated into the polymerisation reaction during the formation of the polyacrylamide gel backbones. The reaction is as follows:
The whole Polyacrylamide RNA based hydrogel synthesis process can be accessed here.
The two single stranded RNA motifs are functionalized to the two polyacrylamide backbones denoted in yellow and green. Since RNA is unstable, the methacryl functionalized RNA motifs are very small,
RNA sequence | Nucleotide sequence |
Green | CCGCC |
Yellow | UUAUU |
A linker single stranded RNA is synthesised containing fifteen nucleotides. Five nucleotides from each end are complementary to the sequences that the RNA motifs polymerised to the backbone contain.
Usage
The linker ssRNAs are incubated with the cas enzymes in an eppendorf tube and the separated backbones are plated on the μpad. This mixture is poured into the second spot of the μpad. If the target miRNA is not present in the sample, the Cas molecules are not activated and the ssRNA linkers are intact. The non degraded linker RNA holds the two backbones together due to complementarity. This gives consistency to the hydrogel and the colored buffer does not pass through giving a negative result.
However if the target miRNA is present, the Cas proteins are activated which leads to collateral cleavage of the linker RNA residues. On adding this to the polyacrylamide spearted backbones, the linker RNA is no longer able to hold the two backbones together leading to the flow out of the coloured buffer and a positive result.
Polyethylene glycol(PEG) Based hydrogel [10]
Polyethylene glycol(PEG) RNA based hydrogel is used when the target miRNA biomarker for a disease is present in very preliminary amounts in normal samples but is upregulated multiple folds during the disease. In such cases, even normal samples will give a positive result as the cas enzymes will get activated. Hence in these cases we will use the PEG hydrogel synthesis method.
The method of synthesis involves initial functionalization of single stranded RNA molecules to disulphide bonds which is then cleaved using (tris(2-carboxyethyl)phosphine) to form single stranded RNA attached to two thiol moieties. This is then functionalized to a vinyl sulphone backbone to make a mesh like structure. This whole PEG hydrogel is then put on the third spot of the μpad.
The synthesis reaction in short is as follows:
The full synthesis reaction can be found here.
Usage
The PEG gel is already synthesized with the ssRNA hybridised to the vinyl sulphone backbone. When this is plated on the μpad, the flow of the buffer will not only depend on the activity of Cas13a and Csm6 but also on the thickness of the gel and the time elapsed. Hence in presence of minor quantities of miRNA, only a few Cas13a molecules will become active, in turn activating fewer Csm6 molecules. This activity of the cas enzymes will not be enough to cleave the whole plated hydrogel and hence within a stipulated threshold time, the buffer will not flow through.
This threshold time as well as the thickness of the gel plated needs to be in order to devise an effective interval where minor base level miRNA presence will not give a false positive result. However large upregulations of the same miRNA will activate many more Cas molecules leading to the cleavage of the whole plated hydrogel within the same threshold time leading to a positive result.
To conclude, we expect that our chosen path of detection and treatment of a subclinical bovine mastitis spurs the much needed change and brings about a decrease in the quantity of antibiotics used. We also intend to bring our detection kit to experimental reality along with the chalked out optimisations so that detection of diseases can be done easily in a paper-based manner.
References
- Lubelski J, Rink R, Khusainov R, Moll GN, Kuipers OP. Biosynthesis, immunity, regulation, mode of action and engineering of the model lantibiotic nisin. Cell Mol Life Sci. 2008 Feb;65(3):455-76. doi: 10.1007/s00018-007-7171-2. PMID: 17965835.
- Sun, Z., Zhong, J., Liang, X., Liu, J., Chen, X., & Huan, L. (2009). Novel Mechanism for Nisin Resistance via Proteolytic Degradation of Nisin by the Nisin Resistance Protein NSR. Antimicrobial Agents and Chemotherapy, 53(5), 1964–1973. https://doi.org/10.1128/aac.01382-08
- Simons, A., Alhanout, K., & Duval, R. E. (2020). Bacteriocins, Antimicrobial Peptides from Bacterial Origin: Overview of Their Biology and Their Impact against Multidrug-Resistant Bacteria. Microorganisms, 8(5), 639. https://doi.org/10.3390/microorganisms8050639
- Okuda K, Zendo T, Sugimoto S, Iwase T, Tajima A, Yamada S, Sonomoto K, Mizunoe Y. Effects of bacteriocins on methicillin-resistant Staphylococcus aureus biofilm. Antimicrob Agents Chemother. 2013 Nov;57(11):5572-9. doi: 10.1128/AAC.00888-13. Epub 2013 Aug 26. PMID: 23979748; PMCID: PMC3811281.
- Birhanu M, Leta S, Mamo G, Tesfaye S. Prevalence of bovine subclinical mastitis and isolation of its major causes in Bishoftu Town, Ethiopia. BMC Res Notes. 2017 Dec 21;10(1):767. doi: 10.1186/s13104-017-3100-0. PMID: 29268785; PMCID: PMC5740909.
- Field D, Blake T, Mathur H, O' Connor PM, Cotter PD, Paul Ross R, Hill C. Bioengineering nisin to overcome the nisin resistance protein. Mol Microbiol. 2019 Mar;111(3):717-731. doi: 10.1111/mmi.14183. Epub 2019 Feb 12. PMID: 30537404.
- Pérez-Ibarreche M, Field D, Ross RP, Hill C. A Bioengineered Nisin Derivative To Control Streptococcus uberis Biofilms. Appl Environ Microbiol. 2021 Jul 27;87(16):e0039121. doi: 10.1128/AEM.00391-21. Epub 2021 Jul 27. PMID: 34105992; PMCID: PMC8315181.
- Conformational features of the Staphylococcus aureus AgrA-promoter interactions rationalize quorum-sensing triggered gene expression Pages 124-134, ISSN 2405-5808,https://doi.org/10.1016/j.bbrep.2016.03.012.(https://www.sciencedirect.com/science/article/pii/S2405580816300334)
- Gootenberg, Jonathan S., et al. “Multiplexed and Portable Nucleic Acid Detection Platform with Cas13, Cas12a, and Csm6.” Science, vol. 360, no. 6387, 2018, pp. 439–444., doi:10.1126/science.aaq0179.
- Gayet, Raphael V., et al. “Creating CRISPR-Responsive Smart Materials for Diagnostics and Programmable Cargo Release.” Nature Protocols, vol. 15, no. 9, 2020, pp. 3030–3063., doi:10.1038/s41596-020-0367-8.
- Moroni, Paolo, et al. “Diseases of the Teats and Udder.” Rebhun's Diseases of Dairy Cattle, 2018, pp. 389–465., doi:10.1016/b978-0-323-39055-2.00008-5.
- Hovinen, Mari et al. “N-acetyl -β-D-glucosaminidase activity in cow milk as an indicator of mastitis.” The Journal of dairy research vol. 83,2 (2016): 219-27. doi:10.1017/S0022029916000224
- Abudayyeh, Omar O., et al. “RNA Targeting with CRISPR–Cas13.” Nature, vol. 550, no. 7675, 2017, pp. 280–284., doi:10.1038/nature24049.
- O'connell, Mitchell R. “Molecular Mechanisms of RNA Targeting by Cas13-Containing Type VI CRISPR–Cas Systems.” Journal of Molecular Biology, vol. 431, no. 1, 2019, pp. 66–87., doi:10.1016/j.jmb.2018.06.029.
- Smalakyte, Dalia, et al. “Type III-A CRISPR-Associated Protein Csm6 Degrades Cyclic Hexa-Adenylate Activator Using Both CARF and HEPN Domains.” Nucleic Acids Research, vol. 48, no. 16, 2020, pp. 9204–9217., doi:10.1093/nar/gkaa634.
- Garcia-Doval, Carmela, et al. “Activation and Self-Inactivation Mechanisms of the Cyclic Oligoadenylate-Dependent CRISPR Ribonuclease Csm6.” Nature Communications, vol. 11, no. 1, 2020, doi:10.1038/s41467-020-15334-5.
- MiRBase, www.mirbase.org/cgi-bin/mirna_entry.pl?acc=MI0009757
- MiRBase, www.mirbase.org/cgi-bin/mirna_entry.pl?acc=MI0004744.
- Lai, Yu-Chang, et al. “Inflammation-Related MicroRNA Expression Level in the Bovine Milk Is Affected by Mastitis.” Plos One, vol. 12, no. 5, 2017, doi:10.1371/journal.pone.0177182.
- Ma, Shaoyang, et al. “Identification and Characterization of Differentially Expressed Exosomal MicroRNAs in Bovine Milk Infected with Staphylococcus Aureus.” BMC Genomics, vol. 20, no. 1, 2019, doi:10.1186/s12864-019-6338-1
- English, Max A., et al. “Programmable CRISPR-Responsive Smart Materials.” Science, vol. 365, no. 6455, 2019, pp. 780–785., doi:10.1126/science.aaw5122.