Impact of visitors and hospital staff on nosocomial transmission and spread to community.

We develop a deterministic meta-population model to qualitatively capture some key features of disease transmission between a community and its healthcare facility. We consider the disease transmission dynamics within a healthcare facility and between the healthcare facility and its community. The focus of this study is to quantify the roles of the healthcare workers at and visitors to this healthcare facility in shaping the transmission dynamics during a disease outbreak. We stratify the total population into the general population in the community and the healthcare workers and visitors in the healthcare facility, to account for nosocomial transmission in the case when an individual in the community may be exposed to an infection due to a visit to the healthcare facility. Equilibrium stability analysis is carried out to inform long-term outcomes of disease dynamics in the coupled community-health care facility system. The basic reproduction number is calculated and its dependence on the waiting time and various disease transmission rates is analyzed. Numerical simulations are performed with pertussis as the disease in question. The results show that waiting time only affects the peak number of infections in the waiting reception area. The results also indicate that transmission rate of infective residents in the community and the transmission rate of the infective visitors at the healthcare facility have decisive impact on disease eradication/persistence of the coupled system; while other modes of transmissions are less important, affecting the peak number of infections at best.

Subfunctionalization reduces the fitness cost of gene duplication in humans by buffering dosage imbalances.

BACKGROUND: Driven essentially by random genetic drift, subfunctionalization has been identified as a possible non-adaptive mechanism for the retention of duplicate genes in small-population species, where widespread deleterious mutations are likely to cause complementary loss of subfunctions across gene copies. Through subfunctionalization, duplicates become indispensable to maintain the functional requirements of the ancestral locus. Yet, gene duplication produces a dosage imbalance in the encoded proteins and thus, as investigated in this paper, subfunctionalization must be subject to the selective forces arising from the fitness bottleneck introduced by the duplication event. RESULTS: We show that, while arising from random drift, subfunctionalization must be inescapably subject to selective forces, since the diversification of expression patterns across paralogs mitigates duplication-related dosage imbalances in the concentrations of encoded proteins. Dosage imbalance effects become paramount when proteins rely on obligatory associations to maintain their structural integrity, and are expected to be weaker when protein complexation is ephemeral or adventitious. To establish the buffering effect of subfunctionalization on selection pressure, we determine the packing quality of encoded proteins, an established indicator of dosage sensitivity, and correlate this parameter with the extent of paralog segregation in humans, using species with larger population -and more efficient selection- as controls. CONCLUSIONS: Recognizing the role of subfunctionalization as a dosage-imbalance buffer in gene duplication events enabled us to reconcile its mechanistic nonadaptive origin with its adaptive role as an enabler of the evolution of genetic redundancy. This constructive role was established in this paper by proving the following assertion: If subfunctionalization is indeed adaptive, its effect on paralog segregation should scale with the dosage sensitivity of the duplicated genes. Thus, subfunctionalization becomes adaptive in response to the selection forces arising from the fitness bottleneck imposed by gene duplication.

Genome-wide detection of hybrid genes with multiple components in human.

BACKGROUND: Previous studies showed that gene hybrid is one of the principal processes for generating new genes. Although some gene hybrid events have been reported to be inter- or intra-species, there lacks a well-organized method for large scale detection of the events with multiple components. Hence in this study, we focus on building up an efficient method for exploring all candidates of gene hybrid events in human genome and provide useful results for further study. FINDINGS: We have developed a method designated Triad Comparison Algorithm (TCA) to detect all potential N-hybrid events (i.e., an N-hybrid gene and its N non-overlapping component regions derived from N different genes) in human genome. The results reveal that there are many convoluted N-hybrid events with multiple components (N > 2) and that the most complicated N-hybrid genes detected in human by TCA are composed of six component regions. Interestingly, our results show that most of the hybrid events belong to the 3-hybrid category. Furthermore, we observe that a single gene might participate in different events. Twelve genes were found to have dual identities contained in different N-hybrid events (i.e., they were identified as hybrid genes as well as component genes). This points out that to a certain extent the gene hybrid mechanism has generated new genes during the course of human genome evolutionary history. CONCLUSION: An efficient method, TCA, is developed for exploring all candidates of hybrid genes in the human genome and provides useful results for the evolutionary analysis. The advantage of TCA is its power of detecting any kinds of hybrid events in any species with a large genome size.

Improved variance estimators for one- and two-parameter models of nucleotide substitution.

The current variance estimators for Jukes and Cantor's one-parameter model and Kimura's two-parameter model tend to underestimate the true variances when the true proportion of differences between the two sequences under study is not small. In this paper, we developed improved variance estimators, using a higher-order Taylor expansion and empirical methods. The new estimators outperform the conventional estimators and provide accurate estimates of the true variances.

Opposite evolutionary effects between different alternative splicing patterns.

Alternative splicing (AS) has been recognized as a mechanism of relaxing selection pressure on protein subsequences. Here, we show that AS may also yield contrary evolutionary effects. We compare the evolutionary rates of 2 types of alternatively spliced exons (ASEs)-simple and complex. The former does not change the boundaries of its flanking exons, whereas the latter does. By analyzing over 26,000 human-mouse orthologous exons, we demonstrate that complex ASEs have lower Ka and Ka/Ks ratio and higher Ks than constitutively spliced exons (CSEs), whereas simple ASEs have evolutionary rates to the opposite of CSEs. Our results indicate that complex ASEs are subject to stronger selection pressure than CSEs at the protein level, but the trend is reversed at the RNA level. Therefore, the previous view that ASEs accelerate evolution of protein subsequences needs to be modified.

Comparison of three methods for estimating rates of synonymous and nonsynonymous nucleotide substitutions.

Three frequently used methods for estimating the synonymous and nonsynonymous substitution rates (Ks and Ka) were evaluated and compared for their accuracies; these methods are denoted by LWL85, LPB93, and GY94, respectively. For this purpose, we used a codon-evolution model to obtain the expected Ka and Ks values for the above three methods and compared the values with those obtained by the three methods. We also proposed some modifications of LWL85 and LPB93 to increase their accuracies. Our computer simulations under the codon-evolution model showed that for sequences < or =300 codons, the performance of GY94 may not be reliable. For longer sequences, GY94 is more accurate for estimating the Ka/Ks ratio than the modified LPB93 and LWL85 in the majority of the cases studied. This is particularly so when k > or = 3, which is the transition/transversion (mutation) rate ratio. However, when k is approximately 2 and when the sequence divergence is relatively large, the modified LWL85 performed better than GY94 and the modified LPB93. The inferiority of LPB93 to LWL85 is surprising because LPB93 was intended to improve LWL85. Also, it has been thought that the codon-based method of GY94 is better than the heuristic method of LWL85, but our simulation results showed that in many cases, the opposite was true, even though our simulation was based on the codon-evolution model.

Comparative analysis of the receptor-like kinase family in Arabidopsis and rice.

Receptor-like kinases (RLKs) belong to the large RLK/Pelle gene family, and it is known that the Arabidopsis thaliana genome contains >600 such members, which play important roles in plant growth, development, and defense responses. Surprisingly, we found that rice (Oryza sativa) has nearly twice as many RLK/Pelle members as Arabidopsis does, and it is not simply a consequence of a larger predicted gene number in rice. From the inferred phylogeny of all Arabidopsis and rice RLK/Pelle members, we estimated that the common ancestor of Arabidopsis and rice had >440 RLK/Pelles and that large-scale expansions of certain RLK/Pelle members and fusions of novel domains have occurred in both the Arabidopsis and rice lineages since their divergence. In addition, the extracellular domains have higher nonsynonymous substitution rates than the intracellular domains, consistent with the role of extracellular domains in sensing diverse signals. The lineage-specific expansions in Arabidopsis can be attributed to both tandem and large-scale duplications, whereas tandem duplication seems to be the major mechanism for recent expansions in rice. Interestingly, although the RLKs that are involved in development seem to have rarely been duplicated after the Arabidopsis-rice split, those that are involved in defense/disease resistance apparently have undergone many duplication events. These findings led us to hypothesize that most of the recent expansions of the RLK/Pelle family have involved defense/resistance-related genes.

Plasmid-bearing, plasmid-free organisms competing for two complementary nutrients in a chemostat.

A model of competition for two complementary nutrients between plasmid-bearing and plasmid-free organisms in a chemostat is proposed. A rigorous mathematical analysis of the global asymptotic behavior of the model is presented. The work extends the model of competition for a single-limited nutrient studied by Stephanopoulos and Lapidus [Chem. Engng. Sci. 443 (1988) 49] and Hsu, Waltman and Wolkowicz [J. Math. Biol. 32 (1994) 731].