Correcting mortality for loss to follow-up: a nomogram applied to antiretroviral treatment programmes in sub-Saharan Africa

Background The World Health Organization estimates that in sub-Saharan Africa about 4 million HIV-infected patients had started antiretroviral therapy (ART) by the end of 2008. Loss of patients to follow-up and care is an important problem for treatment programmes in this region. As mortality is high in these patients compared to patients remaining in care, ART programmes with high rates of loss to follow-up may substantially underestimate mortality of all patients starting ART. Methods and Findings We developed a nomogram to correct mortality estimates for loss to follow-up, based on the fact that mortality of all patients starting ART in a treatment programme is a weighted average of mortality among patients lost to follow-up and patients remaining in care. The nomogram gives a correction factor based on the percentage of patients lost to follow-up at a given point in time, and the estimated ratio of mortality between patients lost and not lost to follow-up. The mortality observed among patients retained in care is then multiplied by the correction factor to obtain an estimate of programme-level mortality that takes all deaths into account. A web calculator directly calculates the corrected, programme-level mortality with 95% confidence intervals (CIs). We applied the method to 11 ART programmes in sub-Saharan Africa. Patients retained in care had a mortality at 1 year of 1.4% to 12.0%; loss to follow-up ranged from 2.8% to 28.7%; and the correction factor from 1.2 to 8.0. The absolute difference between uncorrected and corrected mortality at 1 year ranged from 1.6% to 9.8%, and was above 5% in four programmes. The largest difference in mortality was in a programme with 28.7% of patients lost to follow-up at 1 year. Conclusions The amount of bias in mortality estimates can be large in ART programmes with substantial loss to follow-up. Programmes should routinely report mortality among patients retained in care and the proportion of patients lost. A simple nomogram can then be used to estimate mortality among all patients who started ART, for a range of plausible mortality rates among patients lost to follow-up.

Hibernacula Microclimate and White-nose Syndrome Susceptibility in the Little Brown Myotis (Myotis lucifugus)

The objective of this project was to determine the relationship between hibernacula microclimate and White-nose Syndrome (WNS), an emerging infectious disease in bats. Microclimate was examined on a species scale and at the level of the individual bat to determine if there was a difference in microclimate preference between healthy and WNS-affected little brown myotis (Myotis lucifugus) and to determine the role of microclimate in disease progression. There is anecdotal evidence that colder, drier hibernacula are less affected by WNS. This was tested by placing rugged temperature and humidity dataloggers in field sites throughout the eastern USA, experimentally determining the response to microclimate differences in captive bats, and testing microclimate roosting preference. This study found that microclimate significantly differed from the entrance of a hibernaculum versus where bats traditionally roost. It also found hibernaculum temperature and sex had significant impacts on survival in WNS-affected bats. Male bats with WNS had increased survivability over WNS-affected female bats and WNS bats housed below the ideal growth range of the fungus that causes WNS, Geomyces destructans, had increased survival over those housed at warmer temperatures. The results from this study are immediately applicable to (1) predict which hibernacula are more likely to be infected next winter, (2) further our understanding of WNS, and (3) determine if direct mitigation strategies, such as altering the microclimate of mines, will be effective ways to combat the spread of the fungus.

Exploration Of Survival Traits In The White-Nose Syndrome-Affected Big Brown Bat (Eptesicus Fuscus)

The widespread mortality of hibernating bats is associated with the emerging infectious disease white-nose syndrome (WNS), and has provoked a strong interest in understanding which bats will survive, and why? The ability of infected bats to resist WNS may depend upon variation in the expression of different characteristics. In a captive colony of big brown bats, I sought to characterize the phenotypic variability, repeatability, and survivability for several key ¿survival¿ traits, including: torpor patterns, microclimate preferences, and wound healing capacity. Torpor patterns were profiled using temperature sensitive dataloggers throughout the hibernation season, while microclimate preferences were quantified by using temperature-graded boxes and thermal imaging. In order to assess wound healing capacity, small wing biopsies were obtained from each bat and healing progress was tracked for one month. Individuals exhibited a wide range of phenotypes that were significantly influenced by sex and body condition. Repeatability estimates suggest that there is not a strong genetic basis for the observed variation in torpor patterns or microclimate preferences. Certain phenotypes (e.g., BMI) were associated with an increased probability of overwinter survivorship, which suggests a basis for intra-species differences in WNS susceptibility. The results from this project provide novel insight into what we know about ¿who will survive,¿ and will influence the direction and implementation of future conservation and mitigation strategies.

Paleoenvironment and Paleoecology of a Late Paleocene High-Latitude Terrestrial Succession, Arkose Ridge Formation at Box Canyon, Southern Talkeetna Mountains, Alaska

Paleogene sedimentary rocks of the Arkose Ridge Formation (Talkeetna Mountains, Alaska) preserve a record of a fluvial-lacustrine depositional environment and its forested ecosystem in an active basin among the convergent margin tectonic processes that shaped southern Alaska. An -800 m measured succession at Box Canyon indicates braid-plain deposition with predominantly gravelly deposits low in the exposure to sandy and muddy facies associations below an overlying lava flow sequence. U-Pb geochronology on zircons from a tuff and a sandstone within the measured section, as well as an Ar/Ar date from the overlying lava constrain the age of the sedimentary succession to between similar to 59 Ma and 48 Ma Fossil plant remains occur throughout the Arkose Ridge Formation as poorly-preserved coalified woody debris and fragmentary leaf impressions. At Box Canyon, however, a thin la-custrine depositional lens of rhythmically laminated mudrocks yielded fish fossils and a well-preserved floral assemblage including foliage and reproductive organs representing conifers, sphenopsids, monocots, and dicots. Leaf physiognomic methods to estimate paleoclimate were applied to the dicot leaf collection and indicate warm temperate paleotemperatures (-11-15 +/- -4 degrees C MAT) and elevated paleoprecipitation (-120 cm/yr MAP) estimates as compared to modem conditions; results that are parallel with previously published estimates from the partly coeval Chickaloon Formation deposited in more distal depositional environments in the same basin. The low abundance of leaf herbivory in the Box Canyon dicot assemblage (-9% of leaves damaged) is also similar to the results from assemblages in the meander-plain depositional systems of the Chickaloon. This new suite of data informs models of the tectonostratigraphic evolution of southern Alaska and the developing understanding of terrestrial paleoecology and paleoclimate at high latitudes during the Late Paleocene-Early Eocene greenhouse climate phase. (c) 2014 Elsevier B.V. All rights reserved.

Patch selection by wintering white-tailed deer : consequences for eastern hemlock regeneration at local and landscape scales

Tsuga canadensis (eastern hemlock) is a highly shade-tolerant, late-successional, and long-lived conifer species found throughout eastern North America. It is most often found in pure or nearly pure stands, because highly acidic and nutrient poor forest floor conditions are thought to favor T. canadensis regeneration while simultaneously limiting the establishment of some hardwood species with greater nutrient requirements. Once a common species, T. canadensis is currently experiencing widescale declines across its range. The hemlock woolly adelgid (Adelges tsugae) is decimating the population across its eastern distribution. Across the Upper Great Lakes region, where the adelgid is currently being held at bay by cold winter temperatures, T. canadensis has been experiencing failures in regeneration attributed, in part, to herbivory by white-tailed deer (Odocoileus virginianus). Deer utilize T. canadensis stands as winter habitat in areas of high snow depth. Tsuga canadensis, once a major component of these forests, currently exists at just a fraction of its pre-settlement abundance due to historic logging and contemporary forest management practices, and what remains is found in small remnant patches surrounded by second- and third-growth deciduous forests. The deer population across the region, however, is likely double that of pre-European settlement times. In this dissertation I explore the relationship between white-tailed deer use of T. canadensis as winter habitat and the effect this use is having on regeneration and forest succession. For this research I quantified stand composition and structure and abiotic variables of elevation and snow depth in 39 randomly selected T. canadensis stands from across the western Upper Peninsula of Michigan. I also quantified composition and the configuration of the landscapes surrounding these stands. I measured relative deer use of T. canadensis stands as pellet group piles deposited in each stand during each of three consecutive winters, 2005-06, 2006-07, and 2007-08. The results of this research suggest that deer use of T. canadensis stands as winter habitat is influenced primarily by snow depth, elevation, and the composition and configuration of the greater landscapes surrounding these stands. Specifically, stands with more heterogeneous landscapes surrounding them (i.e., a patchy mosaic of conifer, deciduous, and open cover) had higher relative deer use than stands surrounded by homogenous deciduous forest cover. Additionally, the intensity of use and the number of stands used was greater in years with higher average snow depth. Tsuga canadensis regeneration in these stands was negatively associated with deer use and Acer saccharum (sugar maple) basal area. Of the 39 stands, 17 and 22 stands had no T. canadensis regeneration in small and large sapling categories, respectively. Acer saccharum was the most common understory tree species, and the importance of A. saccharum in the understory (stems < 10 cm dbh) of the stands was positively associated with overstory A. saccharum dominance. Tsuga canadensis establishment was associated with high-decay coarse woody debris and moss, and deciduous leaf litter inputs in these stands may be limiting access to these important microsites. Furthermore, A. saccharum is more tolerant to the effects of deer herbivory than T. canadensis, giving A. saccharum a competitive advantage in stands being utilized as winter habitat by deer. My research suggests that limited microsite availability, in conjunction with deer herbivory, may be leading to an erosion in T. canadensis patch stability and an altered successional trajectory toward one of A. saccharum dominance, an alternately stable climax species.

Multi-node TOA-DOA cooperative LOS-NLOS localization : enabling high accuracy and reliability

This dissertation investigates high performance cooperative localization in wireless environments based on multi-node time-of-arrival (TOA) and direction-of-arrival (DOA) estimations in line-of-sight (LOS) and non-LOS (NLOS) scenarios. Here, two categories of nodes are assumed: base nodes (BNs) and target nodes (TNs). BNs are equipped with antenna arrays and capable of estimating TOA (range) and DOA (angle). TNs are equipped with Omni-directional antennas and communicate with BNs to allow BNs to localize TNs; thus, the proposed localization is maintained by BNs and TNs cooperation. First, a LOS localization method is proposed, which is based on semi-distributed multi-node TOA-DOA fusion. The proposed technique is applicable to mobile ad-hoc networks (MANETs). We assume LOS is available between BNs and TNs. One BN is selected as the reference BN, and other nodes are localized in the coordinates of the reference BN. Each BN can localize TNs located in its coverage area independently. In addition, a TN might be localized by multiple BNs. High performance localization is attainable via multi-node TOA-DOA fusion. The complexity of the semi-distributed multi-node TOA-DOA fusion is low because the total computational load is distributed across all BNs. To evaluate the localization accuracy of the proposed method, we compare the proposed method with global positioning system (GPS) aided TOA (DOA) fusion, which are applicable to MANETs. The comparison criterion is the localization circular error probability (CEP). The results confirm that the proposed method is suitable for moderate scale MANETs, while GPS-aided TOA fusion is suitable for large scale MANETs. Usually, TOA and DOA of TNs are periodically estimated by BNs. Thus, Kalman filter (KF) is integrated with multi-node TOA-DOA fusion to further improve its performance. The integration of KF and multi-node TOA-DOA fusion is compared with extended-KF (EKF) when it is applied to multiple TOA-DOA estimations made by multiple BNs. The comparison depicts that it is stable (no divergence takes place) and its accuracy is slightly lower than that of the EKF, if the EKF converges. However, the EKF may diverge while the integration of KF and multi-node TOA-DOA fusion does not; thus, the reliability of the proposed method is higher. In addition, the computational complexity of the integration of KF and multi-node TOA-DOA fusion is much lower than that of EKF. In wireless environments, LOS might be obstructed. This degrades the localization reliability. Antenna arrays installed at each BN is incorporated to allow each BN to identify NLOS scenarios independently. Here, a single BN measures the phase difference across two antenna elements using a synchronized bi-receiver system, and maps it into wireless channel’s K-factor. The larger K is, the more likely the channel would be a LOS one. Next, the K-factor is incorporated to identify NLOS scenarios. The performance of this system is characterized in terms of probability of LOS and NLOS identification. The latency of the method is small. Finally, a multi-node NLOS identification and localization method is proposed to improve localization reliability. In this case, multiple BNs engage in the process of NLOS identification, shared reflectors determination and localization, and NLOS TN localization. In NLOS scenarios, when there are three or more shared reflectors, those reflectors are localized via DOA fusion, and then a TN is localized via TOA fusion based on the localization of shared reflectors.

Comprehensive review and application of particle image velocimetry

For a fluid dynamics experimental flow measurement technique, particle image velocimetry (PIV) provides significant advantages over other measurement techniques in its field. In contrast to temperature and pressure based probe measurements or other laser diagnostic techniques including laser Doppler velocimetry (LDV) and phase Doppler particle analysis (PDPA), PIV is unique due to its whole field measurement capability, non-intrusive nature, and ability to collect a vast amount of experimental data in a short time frame providing both quantitative and qualitative insight. These properties make PIV a desirable measurement technique for studies encompassing a broad range of fluid dynamics applications. However, as an optical measurement technique, PIV also requires a substantial technical understanding and application experience to acquire consistent, reliable results. Both a technical understanding of particle image velocimetry and practical application experience are gained by applying a planar PIV system at Michigan Technological University’s Combustion Science Exploration Laboratory (CSEL) and Alternative Fuels Combustion Laboratory (AFCL). Here a PIV system was applied to non-reacting and reacting gaseous environments to make two component planar PIV as well as three component stereographic PIV flow field velocity measurements in conjunction with chemiluminescence imaging in the case of reacting flows. This thesis outlines near surface flow field characteristics in a tumble strip lined channel, three component velocity profiles of non-reacting and reacting swirled flow in a swirl stabilized lean condition premixed/prevaporized-fuel model gas turbine combustor operating on methane at 5-7 kW, and two component planar PIV measurements characterizing the AFCL’s 1.1 liter closed combustion chamber under dual fan driven turbulent mixing flow.

Computational study of internal and external condensing flows and experimental synthesis to investigate their attainability and stability in ground-based and space-based environments

This doctoral thesis presents the computational work and synthesis with experiments for internal (tube and channel geometries) as well as external (flow of a pure vapor over a horizontal plate) condensing flows. The computational work obtains accurate numerical simulations of the full two dimensional governing equations for steady and unsteady condensing flows in gravity/0g environments. This doctoral work investigates flow features, flow regimes, attainability issues, stability issues, and responses to boundary fluctuations for condensing flows in different flow situations. This research finds new features of unsteady solutions of condensing flows; reveals interesting differences in gravity and shear driven situations; and discovers novel boundary condition sensitivities of shear driven internal condensing flows. Synthesis of computational and experimental results presented here for gravity driven in-tube flows lays framework for the future two-phase component analysis in any thermal system. It is shown for both gravity and shear driven internal condensing flows that steady governing equations have unique solutions for given inlet pressure, given inlet vapor mass flow rate, and fixed cooling method for condensing surface. But unsteady equations of shear driven internal condensing flows can yield different “quasi-steady” solutions based on different specifications of exit pressure (equivalently exit mass flow rate) concurrent to the inlet pressure specification. This thesis presents a novel categorization of internal condensing flows based on their sensitivity to concurrently applied boundary (inlet and exit) conditions. The computational investigations of an external shear driven flow of vapor condensing over a horizontal plate show limits of applicability of the analytical solution. Simulations for this external condensing flow discuss its stability issues and throw light on flow regime transitions because of ever-present bottom wall vibrations. It is identified that laminar to turbulent transition for these flows can get affected by ever present bottom wall vibrations. Detailed investigations of dynamic stability analysis of this shear driven external condensing flow result in the introduction of a new variable, which characterizes the ratio of strength of the underlying stabilizing attractor to that of destabilizing vibrations. Besides development of CFD tools and computational algorithms, direct application of research done for this thesis is in effective prediction and design of two-phase components in thermal systems used in different applications. Some of the important internal condensing flow results about sensitivities to boundary fluctuations are also expected to be applicable to flow boiling phenomenon. Novel flow sensitivities discovered through this research, if employed effectively after system level analysis, will result in the development of better control strategies in ground and space based two-phase thermal systems.

USING PROBABILISTIC GRAPHICAL MODELS TO DRAW INFERENCES IN SENSOR NETWORKS WITH TRACKING APPLICATIONS

Sensor networks have been an active research area in the past decade due to the variety of their applications. Many research studies have been conducted to solve the problems underlying the middleware services of sensor networks, such as self-deployment, self-localization, and synchronization. With the provided middleware services, sensor networks have grown into a mature technology to be used as a detection and surveillance paradigm for many real-world applications. The individual sensors are small in size. Thus, they can be deployed in areas with limited space to make unobstructed measurements in locations where the traditional centralized systems would have trouble to reach. However, there are a few physical limitations to sensor networks, which can prevent sensors from performing at their maximum potential. Individual sensors have limited power supply, the wireless band can get very cluttered when multiple sensors try to transmit at the same time. Furthermore, the individual sensors have limited communication range, so the network may not have a 1-hop communication topology and routing can be a problem in many cases. Carefully designed algorithms can alleviate the physical limitations of sensor networks, and allow them to be utilized to their full potential. Graphical models are an intuitive choice for designing sensor network algorithms. This thesis focuses on a classic application in sensor networks, detecting and tracking of targets. It develops feasible inference techniques for sensor networks using statistical graphical model inference, binary sensor detection, events isolation and dynamic clustering. The main strategy is to use only binary data for rough global inferences, and then dynamically form small scale clusters around the target for detailed computations. This framework is then extended to network topology manipulation, so that the framework developed can be applied to tracking in different network topology settings. Finally the system was tested in both simulation and real-world environments. The simulations were performed on various network topologies, from regularly distributed networks to randomly distributed networks. The results show that the algorithm performs well in randomly distributed networks, and hence requires minimum deployment effort. The experiments were carried out in both corridor and open space settings. A in-home falling detection system was simulated with real-world settings, it was setup with 30 bumblebee radars and 30 ultrasonic sensors driven by TI EZ430-RF2500 boards scanning a typical 800 sqft apartment. Bumblebee radars are calibrated to detect the falling of human body, and the two-tier tracking algorithm is used on the ultrasonic sensors to track the location of the elderly people.

Long-Range Pollution Transport: Trans-Atlantic Mechanisms and Lagrangian Modeling Methods

Over the past several decades, it has become apparent that anthropogenic activities have resulted in the large-scale enhancement of the levels of many trace gases throughout the troposphere. More recently, attention has been given to the transport pathway taken by these emissions as they are dispersed throughout the atmosphere. The transport pathway determines the physical characteristics of emissions plumes and therefore plays an important role in the chemical transformations that can occur downwind of source regions. For example, the production of ozone (O3) is strongly dependent upon the transport its precursors undergo. O3 can initially be formed within air masses while still over polluted source regions. These polluted air masses can experience continued O3 production or O3 destruction downwind, depending on the air mass's chemical and transport characteristics. At present, however, there are a number of uncertainties in the relationships between transport and O3 production in the North Atlantic lower free troposphere. The first phase of the study presented here used measurements made at the Pico Mountain observatory and model simulations to determine transport pathways for US emissions to the observatory. The Pico Mountain observatory was established in the summer of 2001 in order to address the need to understand the relationships between transport and O3 production. Measurements from the observatory were analyzed in conjunction with model simulations from the Lagrangian particle dispersion model (LPDM), FLEX-PART, in order to determine the transport pathway for events observed at the Pico Mountain observatory during July 2003. A total of 16 events were observed, 4 of which were analyzed in detail. The transport time for these 16 events varied from 4.5 to 7 days, while the transport altitudes over the ocean ranged from 2-8 km, but were typically less than 3 km. In three of the case studies, eastward advection and transport in a weak warm conveyor belt (WCB) airflow was responsible for the export of North American emissions into the FT, while transport in the FT was governed by easterly winds driven by the Azores/Bermuda High (ABH) and transient northerly lows. In the fourth case study, North American emissions were lofted to 6-8 km in a WCB before being entrained in the same cyclone's dry airstream and transported down to the observatory. The results of this study show that the lower marine FT may provide an important transport environment where O3 production may continue, in contrast to transport in the marine boundary layer, where O3 destruction is believed to dominate. The second phase of the study presented here focused on improving the analysis methods that are available with LPDMs. While LPDMs are popular and useful for the analysis of atmospheric trace gas measurements, identifying the transport pathway of emissions from their source to a receptor (the Pico Mountain observatory in our case) using the standard gridded model output, particularly during complex meteorological scenarios can be difficult can be difficult or impossible. The transport study in phase 1 was limited to only 1 month out of more than 3 years of available data and included only 4 case studies out of the 16 events specifically due to this confounding factor. The second phase of this study addressed this difficulty by presenting a method to clearly and easily identify the pathway taken by only those emissions that arrive at a receptor at a particular time, by combining the standard gridded output from forward (i.e., concentrations) and backward (i.e., residence time) LPDM simulations, greatly simplifying similar analyses. The ability of the method to successfully determine the source-to-receptor pathway, restoring this Lagrangian information that is lost when the data are gridded, is proven by comparing the pathway determined from this method with the particle trajectories from both the forward and backward models. A sample analysis is also presented, demonstrating that this method is more accurate and easier to use than existing methods using standard LPDM products. Finally, we discuss potential future work that would be possible by combining the backward LPDM simulation with gridded data from other sources (e.g., chemical transport models) to obtain a Lagrangian sampling of the air that will eventually arrive at a receptor.

Cellular senescence of white blood cells in very long-term survivors after allogeneic hematopoietic stem cell transplantation: the role of chronic graft-versus-host disease and female donor sex

In this single-center, cross-sectional study, we evaluated 44 very long-term survivors with a median follow-up of 17.5 years (range, 11-26 years) after hematopoietic stem cell transplantation. We assessed the telomere length difference in human leukocyte antigen-identical donor and recipient sibling pairs and searched for its relationship with clinical factors. The telomere length (in kb, mean +/- SD) was significantly shorter in all recipient blood cells compared with their donors' blood cells (P < .01): granulocytes (6.5 +/- 0.9 vs 7.1 +/- 0.9), naive/memory T cells (5.7 +/- 1.2 vs 6.6 +/- 1.2; 5.2 +/- 1.0 vs 5.7 +/- 0.9), B cells (7.1 +/- 1.1 vs 7.8 +/- 1.1), and natural killer/natural killer T cells (4.8 +/- 1.0 vs 5.6 +/- 1.3). Chronic graft-versus-host disease (P < .04) and a female donor (P < .04) were associated with a greater difference in telomere length between donor and recipient. Critically short telomeres have been described in degenerative diseases and secondary malignancies. If this hypothesis can be confirmed, identification of recipients at risk for cellular senescence could become part of monitoring long-term survivors after hematopoietic stem cell transplantation.

Evaluation of Binocular Eye Trackers and Algorithms for 3D Gaze Interaction in Virtual Reality Environments

Tracking user’s visual attention is a fundamental aspect in novel human-computer interaction paradigms found in Virtual Reality. For example, multimodal interfaces or dialogue-based communications with virtual and real agents greatly benefit from the analysis of the user’s visual attention as a vital source for deictic references or turn-taking signals. Current approaches to determine visual attention rely primarily on monocular eye trackers. Hence they are restricted to the interpretation of two-dimensional fixations relative to a defined area of projection. The study presented in this article compares precision, accuracy and application performance of two binocular eye tracking devices. Two algorithms are compared which derive depth information as required for visual attention-based 3D interfaces. This information is further applied to an improved VR selection task in which a binocular eye tracker and an adaptive neural network algorithm is used during the disambiguation of partly occluded objects.

E-Learning Systems, Environments and Approaches

The volume consists of twenty-five chapters selected from among peer-reviewed papers presented at the CELDA (Cognition and Exploratory Learning in the Digital Age) 2013 Conference held in Fort Worth, Texas, USA, in October 2013 and also from world class scholars in e-learning systems, environments and approaches. The following sub-topics are included: Exploratory Learning Technologies (Part I), e-Learning social web design (Part II), Learner communities through e-Learning implementations (Part III), Collaborative and student-centered e-Learning design (Part IV). E-Learning has been, since its initial stages, a synonym for flexibility. While this dynamic nature has mainly been associated with time and space it is safe to argue that currently it embraces other aspects such as the learners’ profile, the scope of subjects that can be taught electronically and the technology it employs. New technologies also widen the range of activities and skills developed in e-Learning. Electronic learning environments have evolved past the exclusive delivery of knowledge. Technology has endowed e-Learning with the possibility of remotely fomenting problem solving skills, critical thinking and team work, by investing in information exchange, collaboration, personalisation and community building.

Protein changes and proteolytic degradation in red and white clover plants subjected to waterlogging

Red (Trifolium pratense L., cv. “Start”) and white clover varieties (Trifolium repens L., cv. “Debut” and cv. “Haifa”) were waterlogged for 14 days and subsequently recovered for the period of 21 days. Physiological and biochemical responses of the clover varieties were distinctive, which suggested different sensitivity toward flooding. The comparative study of morphological and biochemical parameters such as stem length, leaflet area, dry weight, protein content, protein pattern and proteolytic degradation revealed prominent changes under waterlogging conditions. Protease activity in the stressed plants increased significantly, especially in red clover cv. “Start”, which exhibited eightfold higher azocaseinolytic activity compared to the control. Changes in the protein profiles were detected by SDS-PAGE electrophoresis. The specific response of some proteins (Rubisco, Rubisco-binding protein, Rubisco activase, ClpA and ClpP protease subunits) toward the applied stress was assessed by immunoblotting. The results characterized the red clover cultivar “Start” as the most sensitive toward waterlogging, expressing reduced levels of Rubisco large and small subunits, high content of ClpP protease subunits and increased activity of protease isoforms.