Proteomic profiles of Leptospira borgpetersenii serovar Hardjo strains JB197 and HB203 cultured at different temperatures

Leptospirosis is a global zoonotic disease affecting humans, domestic, and wild animals. Leptospira are typically shed in the urine of reservoir hosts which persist in suitable environments where incidental host transmission occurs after direct contact with infected urine or contaminated environments. Interestingly, serologically identical L. borgpetersenii serovar Hardjo strains JB197 and HB203 show divergent disease severity in the hamster model; JB197 causes severe acute infection while HB203 causes persistent chronic infection. Historically, serovar Hardjo was limited to culture at 29 °C, but utilization of HAN media allows propagation from host tissues at 37 °C. Here, the proteome of strains JB197 and HB203 were characterized after culture from experimentally challenged hamsters at 29 °C and 37 °C. Comparative analyses of JB197 and HB203 samples cultured at 29 °C yielded 425 significantly differentially expressed (DE) proteins, while strains at 37 °C yielded 613 DE proteins including prominent outer membrane proteins and known virulence factors. In agreement, membrane protein GO terms were identified by STRING network analyses along with numerous metabolic KEGG pathways consistent with condition differences. Within strain, JB197 cultured at 29 °C vs 37 °C identified 529 DE proteins, while HB203 identified 524 DE proteins. Investigating differential protein profiles provide insights into strain specific behaviors with implications for better understanding host-pathogen interactions, disease transmission, and response to environmental conditions which can contribute to vaccine development, diagnostic improvement, and ultimately leptospirosis control. SIGNIFICANCE: Leptospirosis is a devastating zoonotic disease affecting humans, wild and domestic animals around the globe. Different species and serovars of Leptospira can affect various animal host species differently; for instance, a serovar that is asymptomatic in the rat may cause severe disease in a dog or human. These differences in host response are not only found at the species and serovar level for Leptospira, but also at the strain level. A prime example comes from strains JB197 and HB203, both species L. borgpetersenii, both serovar Hardjo. Interestingly, JB197 causes a severe acute infection in the hamster while HB203 causes an asymptomatic chronic infection. Understanding these unique relationships between pathogen and host species is important, especially in the context of prevention technologies such as vaccine design, where the strain of Leptospira used as a bacterin might have different efficiencies in different hosts. In this study, proteomic profiles of strains JB197 and HB203 were analyzed, and results revealed diverse protein expression profiles of outer membrane proteins, as well as proteins functioning in motility and growth.

Effect of selection genotype on immune response to Brucella abortus RB51 in Holstein cattle

Genetic selection for milk production traits in US Holsteins has affected numerous genes associated with reproduction and immunity. This study compares the transcriptomic response of peripheral blood mononuclear cells to an in vitro Brucella abortus strain RB51 (RB51) bacterial challenge between contemporary Holsteins and Holsteins that have not been selected for milk production traits since the mid-1960s. Total RNA was extracted from peripheral blood mononuclear cells from four contemporary and four unselected lactating, primiparous cows following 24-h incubation with or without stimulation with RB51 bacteria. RNA was sequenced and reads analyzed using tools from galaxy.scinet.usda.gov. A total of 412 differentially expressed genes (false discovery rate p < 0.05, log fold change > |1|) were identified. The upregulated genes (genes with higher expression in contemporary than unselected cattle) were enriched for 19 terms/pathways, including alanine, aspartate, and glutamate metabolism, indicating a cellular stress response. Downregulated genes (genes with higher expression in unselected than contemporary cows) were enriched for 37 terms/pathways, representing diverse immune responses, including natural killer cell-mediated immunity, interferon-γ production, negative regulation of interleukin-10 production, and cytokine receptor activity indicating a broad immune response with an emphasis on immune defense. These results provide evidence that differences exist between the two genotypes in response to in vitro bacterial challenge. This suggests that contemporary cows, genetically selected for milk production, may have reduced immune function, including limitations in response to intracellular bacteria.

Fecal and vaginal microbiota of vaccinated and non-vaccinated pregnant elk challenged with Brucella abortus

Introduction

Brucella abortus is the causative agent of brucellosis in cattle and in humans, resulting in economic losses in the agricultural sector and representing a major threat to public health. Elk populations in the American Northwest are reservoirs for this bacterium and transmit the agent to domestic cattle herds. One potential strategy to mitigate the transmission of brucellosis by elk is vaccination of elk populations against B. abortus; however, elk appear to be immunologically distinct from cattle in their responses to current vaccination strategies. The differences in host response to B. abortus between cattle and elk could be attributed to differences between the cattle and elk innate and adaptive immune responses. Because species-specific interactions between the host microbiome and the immune system are also known to affect immunity, we sought to investigate interactions between the elk microbiome and B. abortus infection and vaccination.

Methods

We analyzed the fecal and vaginal microbial communities of B. abortus-vaccinated and unvaccinated elk which were challenged with B. abortus during the periparturient period.

Results

We observed that the elk fecal and vaginal microbiota are similar to those of other ruminants, and these microbial communities were affected both by time of sampling and by vaccination status. Notably, we observed that taxa representing ruminant reproductive tract pathogens tended to increase in abundance in the elk vaginal microbiome following parturition. Furthermore, many of these taxa differed significantly in abundance depending on vaccination status, indicating that vaccination against B. abortus affects the elk vaginal microbiota with potential implications for animal reproductive health.

Discussion

This study is the first to analyze the vaginal microbiota of any species of the genus Cervus and is also the first to assess the effects of B. abortus vaccination and challenge on the vaginal microbiome.

Wildlife Immune Responses to Mycobacterium bovis and to Bacille of Calmette-Guerin

Bovine tuberculosis (bTB) is a zoonotic bacterial disease presenting public health, veterinary, and economic threats around the globe. Although cattle producers rely on regular testing and management practices to minimize domestic herd exposure, wildlife species around the world continue to be the main reservoirs for disease. Wildlife reservoirs for bTB include the Eurasian badger (Meles meles) in Great Britain and Ireland, the brushtail possum (Trichosurus vulpecula) in New Zealand, wild boar (Sus scrofa) in Spain, as well as white-tailed deer (Odocoileus virginianus) in the United States and red deer (Cervus elaphus) in Spain. Although all reservoir species share the ability to infect cattle, they differ in transmission capability, disease pathogenesis, diagnostic detection, and vaccination strategies. In this review, bTB interactions with these wildlife reservoirs are discussed, illustrating the need to address bTB disease in wildlife hosts to achieve eradication in domestic livestock.

Vaccination of White-Tailed Deer with Mycobacterium bovis Bacillus Calmette-Guérin (BCG): Effect of Mycobacterium avium ssp. paratuberculosis Infection

In many parts of the world, bovine tuberculosis eradication efforts are hampered by wildlife reservoirs of Mycobacterium bovis, which serve as a constant source of M. bovis for nearby cattle. The human tuberculosis vaccine, M. bovis BCG has been investigated for use in several wildlife species, including deer. In the US, white-tailed deer in Michigan have been the source of infection for over 82 cattle herds since M. bovis was discovered in free-ranging deer in 1995. The efficacy of BCG may be influenced by many factors, including prior exposure or infection with non-tuberculous mycobacteria, that is, species other than members of the M. tuberculosis complex. M. avium subspecies paratuberculosis (Map) infection is not uncommon in ruminants such as deer. Using natural exposure to Map and experimental infection with M. bovis, we demonstrate that Map infection increased BCG vaccine efficacy as measured by lesion severity scores.

Bovine Leptospirosis Due to Persistent Renal Carriage of Leptospira borgpetersenii Serovar Tarassovi

Leptospirosis is a global zoonotic disease that causes significant morbidity and mortality in human and animal populations. Leptospira interrogans is a leading cause of human disease, and L. borgpetersenii is a leading cause of animal disease. Cattle are reservoir hosts of L. borgpetersenii serovar Hardjo, which is transmitted via urine, semen, and uterine discharges resulting in abortion and poor reproductive performance. Bovine bacterin vaccines can only protect against those serovars included in vaccine formulations and typically include serovar Hardjo among others. Genotyping and serotyping represent two different and unique methods for classifying leptospires that do not always correlate well; comprehensive characterization using either method requires recovery of isolates from infected animals. In this study, we report for the first time, isolation of L. borgpetersenii serovar Tarassovi from the urine of a dairy cow in the U.S. The classification of the isolate, designated strain MN900, was confirmed by whole-genome sequencing, serotyping with reference antisera and monoclonal antibodies, Matrix Assisted Laser Desorption/Ionization (MALDI), and immunoblotting with reference antisera. Strain MN900 was excreted in urine samples for 18 weeks even as the cow was seronegative for serovar Tarassovi. Strain MN900 has an unusual morphology since it is not as motile as other leptospires and lacks hooked ends. Serovar Tarassovi is not included in U.S. bacterin vaccines. These results demonstrate the importance of culture and concomitant genotyping and serotyping to accurately classify leptospires, and as required to design efficacious vaccine and diagnostic strategies to not only limit animal disease but reduce zoonotic risk.

Evaluation of LipL32 and LigA/LigB Knockdown Mutants in Leptospira interrogans Serovar Copenhageni: Impacts to Proteome and Virulence

Leptospirosis is a worldwide zoonosis caused by pathogenic species of the genus Leptospira. The recent application of CRISPR interference (CRISPRi) to Leptospira facilitates targeted gene silencing and provides a new tool to investigate pathogenic mechanisms of leptospirosis. CRISPRi relies on the expression of a catalytically “dead” Cas9 (dCas9) and a single-guide RNA (sgRNA). Previously, our group generated a LipL32 and a double LigA/LigB (LigAB) mutant, which, in the current study, are characterized by whole-cell proteomics in comparison with control leptospires harboring plasmid expressing dCas9 alone. Comparison of control and LigAB mutant leptospires identified 46 significantly differentially expressed (DE) proteins, including 27 proteins that were less abundant and 19 proteins that were more abundant in the LigAB mutant compared with the control. Comparison of the control and LipL32 mutant leptospires identified 243 DE proteins, of which 84 proteins were more abundant and 159 were less abundant in the LipL32 mutant strain. Significantly increased amounts of known virulence impactors and surface membrane receptors, including LipL45, LipL31, LigB, and LipL41, were identified. The virulence of LipL32 and LigAB mutants were evaluated in the hamster model of leptospirosis; the LigAB mutant was unable to cause acute disease although mutant leptospires could still be recovered from target organs, albeit at a significantly lower bacterial burden (<850 and <16-fold in liver and kidney, respectively, in comparison with control), indicating attenuation of virulence and a shift to chronic bacterial persistence. Notably, the LipL32 mutant displayed augmented virulence as evidenced by early onset of clinical symptoms and increased numbers of circulating foamy macrophages. Validation of LipL32 and LigAB mutants recovered from liver and kidney in the presence or absence of antibiotic selection revealed high plasmid stability and, by extension, gene silencing in vivo. Collectively, this work emphasizes the advantages and feasibility of using CRISPRi technology to evaluate and characterize virulence factors of leptospires and their respective host–pathogen interactions in animal models of leptospirosis. Importantly, it also provides insight into the requirements of LigA and LigB for acute disease and explores the impact of silencing expression of lipL32, which resulted in substantial changes in amounts of outer membrane proteins.

Complete Genome Sequence of Four Strains of Leptospira borgpetersenii serovar Hardjo isolated from Cattle in the Central United States

Pathogenic species of Leptospira cause leptospirosis, a global zoonotic disease affecting humans and all major livestock species. Cattle act as a reservoir host for L. borgpetersenii serovar Hardjo which colonize the kidneys and reproductive tract from which they are excreted and transmitted to other cattle via urine, semen or uterine discharges. Bovine leptospirosis results in reproductive failure, abortion, stillbirth and loss of milk production, and is an occupational risk for those working with infected animals. A recent study determined that 7.2% of cattle from an abattoir in the central United States were actively shedding pathogenic Leptospira. Here, we report and compare the complete genome sequence of four recent isolates of L. borgpetersenii serovar Hardjo designated strain TC112, TC147, TC129, and TC273.

Circulating Foamy Macrophages in the Golden Syrian Hamster (Mesocricetus auratus) Model of Leptospirosis

Leptospirosis is a world-wide zoonotic disease caused by pathogenic Leptospira and can be asymptomatic or can cause clinical signs ranging from influenza-like to multi-organ failure and death in severe cases. While species and strain specificity can play a major role in disease presentation, the hamster is susceptible to most leptospiral infections and is the model of choice for vaccine efficacy testing. During evaluation of blood smears from hamsters challenged with different species and strains of Leptospira, a circulating population of large, mononuclear, lipid-filled cells, most similar to foamy macrophages (FMs), was detected. Circulating FMs were identified by Giemsa staining and verified by scanning and transmission electron microscopy. FMs were found in the circulating blood of all Leptospira-challenged hamsters, indicating that the finding was not species or strain specific, although higher numbers of FMs tended to correlate with severity of disease. The unique finding of circulating FMs in the hamster model of leptospirosis can yield additional insights into the pathogenesis of leptospirosis and other diseases that induce circulating FMs.

Distinct transcriptional profiles of Leptospira borgpetersenii serovar Hardjo strains JB197 and HB203 cultured at different temperatures

Background

Leptospirosis is a zoonotic, bacterial disease, posing significant health risks to humans, livestock, and companion animals around the world. Symptoms range from asymptomatic to multi-organ failure in severe cases. Complex species-specific interactions exist between animal hosts and the infecting species, serovar, and strain of pathogen. Leptospira borgpetersenii serovar Hardjo strains HB203 and JB197 have a high level of genetic homology but cause different clinical presentation in the hamster model of infection; HB203 colonizes the kidney and presents with chronic shedding while JB197 causes severe organ failure and mortality. This study examines the transcriptome of L. borgpetersenii and characterizes differential gene expression profiles of strains HB203 and JB197 cultured at temperatures during routine laboratory conditions (29°C) and encountered during host infection (37°C).

Methodology/Principal findings

L. borgpetersenii serovar Hardjo strains JB197 and HB203 were isolated from the kidneys of experimentally infected hamsters and maintained at 29°C and 37°C. RNAseq revealed distinct gene expression profiles; 440 genes were differentially expressed (DE) between JB197 and HB203 at 29°C, and 179 genes were DE between strains at 37°C. Comparison of JB197 cultured at 29°C and 37°C identified 135 DE genes while 41 genes were DE in HB203 with those same culture conditions. The consistent differential expression of ligB, which encodes the outer membrane virulence factor LigB, was validated by immunoblotting and 2D-DIGE. Differential expression of lipopolysaccharide was also observed between JB197 and HB203.

Conclusions/Significance

Investigation of the L. borgpetersenii JB197 and HB203 transcriptome provides unique insight into the mechanistic differences between acute and chronic disease. Characterizing the nuances of strain to strain differences and investigating the environmental sensitivity of Leptospira to temperature is critical to the development and progress of leptospirosis prevention and treatment technologies, and is an important consideration when serovars are selected and propagated for use as bacterin vaccines as well as for the identification of novel therapeutic targets.

Leptospirosis is a global zoonotic, neglected tropical disease. Interestingly, a high level of species specificity (both bacteria and host) plays a major role in the severity of disease presentation which can vary from asymptomatic to multi-organ failure. Pathogenic Leptospira colonize the kidneys of infected individuals and are shed in urine into the environment where they can survive until they are contracted by another host. This study looks at two strains of L. borgpetersenii, HB203 and JB197 which are genetically very similar, and identical by serotyping as serovar Hardjo, yet HB203 causes a chronic infection in the hamster while JB197 causes organ failure and mortality. To better characterize bacterial factors causing different disease outcomes, we examined the gene expression profile of these strains in the context of temperatures that would reflect natural Leptospira life cycles (environmentally similar 29°C and 37°C which is more indicative of host environment). We found vast differences in gene expression both between the strains and within strains between temperatures. Characterization of the transcriptome of L. borgpetersenii serovar Hardjo strains JB197 and HB203 provides insights into factors that can determine acute versus chronic disease in the hamster model of infection. Additionally, these studies highlight strain to strain variability within the same species, and serovar, at different growth temperatures, which needs to be considered when serovars are selected and propagated for use as bacterin vaccines used to immunize domestic animal species.

Case report: characterization of a persistent, treatment-resistant, novel Staphylococcus aureus infection causing chronic mastitis in a Holstein dairy cow

Background

Mastitis is the most common health concern plaguing the modern dairy cow and costs dairy producers estimates of two billion dollars annually. Staphylococcus aureus infections are prevalent, displaying varied disease presentation and markedly low cure rates. Neutrophils are considered the first line of defense against mastitis causing bacteria and are frequently targeted in the development of treatment and prevention technologies. We describe a case of naturally occurring, chronic mastitis in a Holstein cow (1428), caused by a novel strain of S. aureus that was not able to be cleared by antibiotic treatment.

Case presentation

The infection was identified in a single quarter, 2 months into the cow’s first lactation. The infection persisted for the following 20 months, including through dry off, and a second calving and lactation. This case of mastitis was associated with a consistently high somatic cell count, however presented with no other clinical signs. This cow was unsuccessfully treated with antibiotics commonly used to treat mastitis, consisting of two rounds of treatment during lactation and an additional round at the beginning of dry off. The chronic infection was also unchanged through an experimental mid-lactation treatment with pegylated granulocyte-colony stimulating factor (PEG-gCSF) and an additional periparturient treatment with PEG-gCSF. We isolated milk neutrophils from 1428 and compared them to two cows challenged with experimental S. aureus, strain Newbould 305. Neutrophils from 1428’s milk had higher surface expression of myeloperoxidase compared to experimental Newbould challenged animals, as well as increased presence of Neutrophil Extracellular Traps. This suggests a heightened activation state of neutrophils sourced from 1428’s naturally occurring infection. Upon postmortem examination, the affected quarter revealed multifocal abscesses separated by fibrous connective tissues. Abscesses were most common in the gland cistern and collecting duct region. Microscopically, the inflammatory reaction was pyogranulomatous to granulomatous and consistent with botryomycosis. Colonies of Gram-positive cocci were found within the eosinophilic matrix of the Splendore-Hoeppli reaction within granulomas and intracellularly within the acinar epithelium.

Conclusions

Collectively, we describe a unique case of chronic mastitis, the characterization of which provides valuable insight into the mechanics of S. aureus treatment resistance and immune escape.

Investigating the Immunological and Biological Equilibrium of Reservoir Hosts and Pathogenic Leptospira: Balancing the Solution to an Acute Problem?

Leptospirosis is a devastating zoonotic disease affecting people and animals across the globe. Pathogenic leptospires are excreted in urine of reservoir hosts which directly or indirectly leads to continued disease transmission, via contact with mucous membranes or a breach of the skin barrier of another host. Human fatalities approach 60,000 deaths per annum; though most vertebrates are susceptible to leptospirosis, complex interactions between host species and serovars of Leptospira can yield disease phenotypes that vary from asymptomatic shedding in reservoir hosts, to multi-organ failure in incidental hosts. Clinical symptoms of acute leptospirosis reflect the diverse range of pathogenic species and serovars that cause infection, the level of exposure, and the relationship of the pathogen with the given host. However, in all cases, pathogenic Leptospira are excreted into the environment via urine from reservoir hosts which are uniformly recognized as asymptomatic carriers. Therefore, the reservoir host serves as the cornerstone of persistent disease transmission. Although bacterin vaccines can be used to abate renal carriage and excretion in domestic animal species, there is an urgent need to advance our understanding of immune-mediated host–pathogen interactions that facilitate persistent asymptomatic carriage. This review summarizes the current understanding of host–pathogen interactions in the reservoir host and prioritizes research to unravel mechanisms that allow for colonization but not destruction of the host. This information is required to understand, and ultimately control, the transmission of pathogenic Leptospira.

Novel Engraftment and T Cell Differentiation of Human Hematopoietic Cells in ART-/-IL2RG-/Y SCID Pigs

Pigs with severe combined immunodeficiency (SCID) are an emerging biomedical animal model. Swine are anatomically and physiologically more similar to humans than mice, making them an invaluable tool for preclinical regenerative medicine and cancer research. One essential step in further developing this model is the immunological humanization of SCID pigs. In this work we have generated T- B- NK- SCID pigs through site directed CRISPR/Cas9 mutagenesis of IL2RG within a naturally occurring DCLRE1C (ARTEMIS)-/- genetic background. We confirmed ART-/-IL2RG-/Y pigs lacked T, B, and NK cells in both peripheral blood and lymphoid tissues. Additionally, we successfully performed a bone marrow transplant on one ART-/-IL2RG-/Y male SCID pig with bone marrow from a complete swine leukocyte antigen (SLA) matched donor without conditioning to reconstitute porcine T and NK cells. Next, we performed in utero injections of cultured human CD34+ selected cord blood cells into the fetal ART-/-IL2RG-/Y SCID pigs. At birth, human CD45+ CD3ε+ cells were detected in cord and peripheral blood of in utero injected SCID piglets. Human leukocytes were also detected within the bone marrow, spleen, liver, thymus, and mesenteric lymph nodes of these animals. Taken together, we describe critical steps forwards the development of an immunologically humanized SCID pig model.

MicroRNA profiles of dry secretions through the first three weeks of the dry period from Holstein cows

In dairy cows, the period from the end of lactation through the dry period and into the transition period, requires vast physiological and immunological changes critical to mammary health. The dry period is important to the success of the next lactation and intramammary infections during the dry period will adversely alter mammary function, health and milk production for the subsequent lactation. MicroRNAs (miRNAs) are small non-coding RNAs that can post transcriptionally regulate gene expression. We sought to characterize the miRNA profile in dry secretions from the last day of lactation to 3, 10, and 21 days post dry-off. We identified 816 known and 80 novel miRNAs. We found 46 miRNAs whose expression significantly changed (q-value < 0.05) over the first three weeks of dry-off. Additionally, we examined the slopes of random regression models of log transformed normalized counts and cross analyzed the 46 significantly upregulated and downregulated miRNAs. These miRNAs were found to be associated with important components of pregnancy, lactation, as well as inflammation and disease. Detailing the miRNA profile of dry secretions through the dry-off period provides insight into the biology at work, possible means of regulation, components of resistance and/or susceptibility, and outlets for targeted therapy development.

Differential phenotype of immune cells in blood and milk following pegylated granulocyte colony-stimulating factor therapy during a chronic Staphylococcus aureus infection in lactating Holsteins

Neutrophils are principal host innate immune cell responders to mastitis infections. Thus, therapies have been developed that target neutrophil expansion. This includes the neutrophil-stimulating cytokine granulocyte colony-stimulating factor (gCSF). Pegylated gCSF (PEG-gCSF; Imrestor, Elanco Animal Health, Greenfield, IN) has been shown to reduce the natural incidence of mastitis in periparturient cows in commercial settings and reduce severity of disease against experimental mastitis challenge. Pegylated gCSF stimulates neutrophil expansion but also induces changes in monocyte and lymphocyte circulating numbers, surface protein expression changes, or both. We hypothesized that PEG-gCSF modulates surface expression of monocytes and neutrophils and facilitates their migration to the mammary gland. We challenged 8 mid-lactation Holsteins with approximately 150 cfu of Staphylococcus aureus (Newbould 305) in a single quarter via intramammary infusion. All animals developed chronic infections as assessed by bacteria counts and somatic cell counts (SCC). Ten to 16 wk postchallenge, 4 of the animals were treated with 2 subcutaneous injections of PEG-gCSF 7 d apart. Complete blood counts, SCC, bacterial counts, milk yield, feed intake, neutrophils extracellular trap analysis, and flow cytometric analyses of milk and blood samples were performed at indicated time points for 14 d after the first PEG-gCSF injection. The PEG-gCSF-treated cows had significantly increased numbers of blood neutrophils and lymphocytes compared with control cows. Flow cytometric analyses revealed increased surface expression of myeloperoxidase (MPO) on neutrophils and macrophages in milk but not in blood of treated cows. Neutrophils isolated from blood of PEG-gCSF-treated cows had decreased surface expression of CD62L (L-selectin) in blood, consistent with cell activation. Surprisingly, CD62L cell surface expression was increased on neutrophils and macrophages sourced from milk from treated animals compared with cells isolated from controls. The PEG-gCSF-treated cows did not clear the S. aureus infection, nor did they significantly differ in SCC from controls. These findings provide evidence that PEG-gCSF therapy modifies cell surface expression of neutrophils and monocytes. However, although surface MPO⁺ cells accumulate in the mammary gland, the lack of bacterial control from these milk-derived cells suggests an incomplete role for PEG-gCSF treatment against chronic S. aureus infection and possibly chronic mammary infections in general.

Identifying constraints on the higher-order structure of RNA: continued development and application of comparative sequence analysis methods

Comparative sequence analysis addresses the problem of RNA folding and RNA structural diversity, and is responsible for determining the folding of many RNA molecules, including 5S, 16S, and 23S rRNAs, tRNA, RNAse P RNA, and Group I and II introns. Initially this method was utilized to fold these sequences into their secondary structures. More recently, this method has revealed numerous tertiary correlations, elucidating novel RNA structural motifs, several of which have been experimentally tested and verified, substantiating the general application of this approach. As successful as the comparative methods have been in elucidating higher-order structure, it is clear that additional structure constraints remain to be found. Deciphering such constraints requires more sensitive and rigorous protocols, in addition to RNA sequence datasets that contain additional phylogenetic diversity and an overall increase in the number of sequences. Various RNA databases, including the tRNA and rRNA sequence datasets, continue to grow in number as well as diversity. Described herein is the development of more rigorous comparative analysis protocols. Our initial development and applications on different RNA datasets have been very encouraging. Such analyses on tRNA, 16S and 23S rRNA are substantiating previously proposed associations and are now beginning to reveal additional constraints on these molecules. A subset of these involve several positions that correlate simultaneously with one another, implying units larger than a basepair can be under a phylogenetic constraint.