Use of the EM algorithm to detect QTL affecting multiple-traits in an across half-sib family analysis

QTL detection experiments in livestock species commonly use the half-sib design. Each male is mated to a number of females, each female producing a limited number of progeny. Analysis consists of attempting to detect associations between phenotype and genotype measured on the progeny. When family sizes are limiting experimenters may wish to incorporate as much information as possible into a single analysis. However, combining information across sires is problematic because of incomplete linkage disequilibrium between the markers and the QTL in the population. This study describes formulæ for obtaining MLEs via the expectation maximization (EM) algorithm for use in a multiple-trait, multiple-family analysis. A model specifying a QTL with only two alleles, and a common within sire error variance is assumed. Compared to single-family analyses, power can be improved up to fourfold with multi-family analyses. The accuracy and precision of QTL location estimates are also substantially improved. With small family sizes, the multi-family, multi-trait analyses reduce substantially, but not totally remove, biases in QTL effect estimates. In situations where multiple QTL alleles are segregating the multi-family analysis will average out the effects of the different QTL alleles.     

Population dynamic of the extinct European aurochs: genetic evidence of a north-south differentiation pattern and no evidence of post-glacial expansion

Background

Various evolutionary models have been proposed to interpret the fate of paralogous duplicates, which provides substrates on which evolution selection could act. In particular, domestication, as a special selection, has played important role in crop cultivation with divergence of many genes controlling important agronomic traits. Recent studies have indicated that a pair of duplicate genes was often sub-functionalized from their ancestral functions held by the parental genes. We previously demonstrated that the rice cell-wall invertase (CWI) gene GIF1 that plays an important role in the grain-filling process was most likely subjected to domestication selection in the promoter region. Here, we report that GIF1 and another CWI gene OsCIN1 constitute a pair of duplicate genes with differentiated expression and function through independent selection.

Results

Through synteny analysis, we show that GIF1 and another cell-wall invertase gene OsCIN1 were paralogues derived from a segmental duplication originated during genome duplication of grasses. Results based on analyses of population genetics and gene phylogenetic tree of 25 cultivars and 25 wild rice sequences demonstrated that OsCIN1 was also artificially selected during rice domestication with a fixed mutation in the coding region, in contrast to GIF1 that was selected in the promoter region. GIF1 and OsCIN1 have evolved into different expression patterns and probable different kinetics parameters of enzymatic activity with the latter displaying less enzymatic activity. Overexpression of GIF1 and OsCIN1 also resulted in different phenotypes, suggesting that OsCIN1 might regulate other unrecognized biological process.

Conclusion

How gene duplication and divergence contribute to genetic novelty and morphological adaptation has been an interesting issue to geneticists and biologists. Our discovery that the duplicated pair of GIF1 and OsCIN1 has experienced sub-functionalization implies that selection could act independently on each duplicate towards different functional specificity, which provides a vivid example for evolution of genetic novelties in a model crop. Our results also further support the established hypothesis that gene duplication with sub-functionalization could be one solution for genetic adaptive conflict.     

Sexual dimorphism in the postcranial skeleton of New World primates

This study examines sexual dimorphism in 24 dimensions of the postcranial skeleton of four platyrrhine species: Callithrix jacchus, Saguinus nigricollis, Saimiri sciureus, and Cebus albifrons. The two callitrichid species show a relatively small amount of variation in the degree of sexual dimorphism among the different dimensions. Variation is considerably higher in the two cebid species as reflected by a mosaic pattern of sexual dimorphisms with males being significantly larger than females in some dimensions, and females significantly larger than males in others. In dimensions of the pectoral girdle and limb bones, males and females in each of the two cebid species are essentially scaled versions of each other, with males being peramorphic compared to females. This pattern is primarily the result of time hypermorphosis, i.e. an extension of the growth period in time in males. Rate hypermorphosis, i.e. an increase in the rate of growth in time in males, appears to play an additional role, however, in S. sciureus. By contrast, in dimensions of the true pelvis, sex differences in shape are dissociated from those in size. They are interpreted as the result of acceleration, i.e. increase in rate of shape change in females, as an adaptation to obstetrical functions. Interspecific analyses indicate positive allometry of mean degree of postcranial dimorphism with respect to body size. This coincides with previous findings by Leutenegger and Cheverud [1982, 1985] on the scaling of sexual dimorphism in body weight and canine size, and thus supports their model which posits selection on body size as the prime mover for the evolution of sexual dimorphism.     

Structural Basis for Type VI Secretion Effector Recognition by a Cognate Immunity Protein

The type VI secretion system (T6SS) has emerged as an important mediator of interbacterial interactions. A T6SS from Pseudomonas aeruginosa targets at least three effector proteins, type VI secretion exported 1–3 (Tse1–3), to recipient Gram-negative cells. The Tse2 protein is a cytoplasmic effector that acts as a potent inhibitor of target cell proliferation, thus providing a pronounced fitness advantage for P. aeruginosa donor cells. P. aeruginosa utilizes a dedicated immunity protein, type VI secretion immunity 2 (Tsi2), to protect against endogenous and intercellularly-transferred Tse2. Here we show that Tse2 delivered by the T6SS efficiently induces quiescence, not death, within recipient cells. We demonstrate that despite direct interaction of Tsi2 and Tse2 in the cytoplasm, Tsi2 is dispensable for targeting the toxin to the secretory apparatus. To gain insights into the molecular basis of Tse2 immunity, we solved the 1.00 Å X-ray crystal structure of Tsi2. The structure shows that Tsi2 assembles as a dimer that does not resemble previously characterized immunity or antitoxin proteins. A genetic screen for Tsi2 mutants deficient in Tse2 interaction revealed an acidic patch distal to the Tsi2 homodimer interface that mediates toxin interaction and immunity. Consistent with this finding, we observed that destabilization of the Tsi2 dimer does not impact Tse2 interaction. The molecular insights into Tsi2 structure and function garnered from this study shed light on the mechanisms of T6 effector secretion, and indicate that the Tse2–Tsi2 effector–immunity pair has features distinguishing it from previously characterized toxin–immunity and toxin–antitoxin systems.     

Succession in Stream Biofilms is an Environmentally Driven Gradient of Stress Tolerance

The century-long research on succession has bestowed us with a number of theories, but little agreement on what causes species replacements through time. The majority of studies has explored the temporal trends of individual species in plant and much less so in microbial communities, arguing that interspecific interactions, especially competition, play a key role in community organization throughout succession. In this experimental investigation of periphytic succession in re-circulating laboratory streams, we examined the density and the relative abundance of diatoms and soft algae for 35 days across gradients of low to high nutrient supply (nitrogen + phosphorus) and low to intermediate current velocity (10 vs. 30 cm·s⁻¹). All algal species were classified into trophic groups and morphological guilds, both of which responded more strongly to nutrient than current velocity manipulations, as shown by regression analyses. We concluded that within the manipulated environmental ranges: (1) Succession was a gradient of stress tolerance, driven primarily by nutrient supply and secondarily, by current velocity. Nutrient supply had a qualitative effect in determining whether the contribution of species tolerant vs. sensitive to nutrient limitation would increase through time, while current velocity had a quantitative influence and affected only the rate of this increase. (2) The mechanism of algal succession at a functional level was a neutral coexistence, whereby the tolerant low profile guild maintained high density when overgrown by sensitive species, while sensitive species, constituting mostly the motile and high profile guilds, were neither facilitated nor inhibited by tolerant species but controlled by the environment. It is suggested that the mechanism of succession may depend on the level of biological organization with interspecific interactions giving way to neutral coexistence along the hierarchy from species to functional groups. Considering that the functional makeup is strictly environmentally defined, while species composition reflects local and regional species pools that may exhibit substantial geographic variability, succession is deterministic at a functional level but stochastic at a species level.     

Fate of the Benzene Ring of Linear Alkylbenzene Sulfonate in Natural Waters

The biodegradability of linear alkylbenzene sulfonate (LAS) was studied in water samples collected from a receiving stream at locations above and below the discharge of a municipal wastewater treatment plant. Rates of primary biodegradation were determined for a commercial LAS mixture by a modified methylene blue-active substance method. Rates of LAS ultimate degradation were determined by radiochemical methods, using a C₁₂ LAS homolog uniformly labeled with ¹⁴C in the benzene ring. The C₁₂ LAS was tested at low concentrations (50 and 500 μg/liter) comparable to those existing in the receiving stream. Loss of methylene blue-active substance response over time occurred rapidly in water samples containing sediment collected from below the treatment plant, with an estimated half-life for LAS of 0.23 days. Evolution of ¹⁴CO₂ during mineralization of the benzene ring occurred rapidly in the same samples, with a half-life for the benzene ring of 0.73 day. Mineralization of the benezene ring was also observed in river water containing no sediments and in river water and sediment samples collected from above the treatment plant. However, the rate of degradation was reduced in these cases, with half-lives for ring carbon ranging from 1.4 to 14 days. Although LAS degradation was enhanced in the presence of sediments, adsorption of LAS to the clay-silt fraction of river sediments was low, and most of the radioactivity was bound to biomass.     

Molecular phylogenetics of western North American frogs of the Rana boylii species group

Phylogenetic relationships among frogs of the genus Rana from western North America are investigated using 2013 aligned bases of mitochondrial DNA sequence from the genes encoding ND1 (subunit one of NADH dehydrogenase), tRNA(Ile), tRNA(Gln), tRNA(Met), ND2, tRNA(Trp), tRNA(Ala), tRNA(Asn), tRNA(Cys), tRNA(Tyr), and COI (subunit I of cytochrome c oxidase), plus the origin for light-strand replication (O(L)) between the tRNA(Asn) and tRNA(Cys) genes. The aligned sequences contain 401 phylogenetically informative characters. A well-resolved phylogenetic hypothesis in which the Rana boylii species group (R. aurora, R. boylii, R. cascadae, R. muscosa, and R. pretiosa) is monophyletic is obtained. Molecular sequence divergence suggests that the R. boylii species group is approximately 8 million years old. The traditional hypothesis showing monophyly of the yellow-legged frogs (R. boylii and R. muscosa) is statistically rejected in favor of a hypothesis in which R. aurora, R. cascadae, and R. muscosa form a clade. Reanalyses of published nuclear ribosomal DNA restriction-site data and allozymic data support a monophyletic R. boylii group, but do not effectively resolve relationships among species within this group. Eight populations of R. muscosa form two major clades separated by a biogeographic break in the Sierra Nevada of California. This biogeographic break is broadly concordant with breaks found in four other amphibian and reptilian taxa. The two major clades within R. muscosa are estimated to have diverged approximately 2.2 million years before present. Each of these major clades contains two subgroups showing approximately 1.5 million years divergence, implicating climatic effects of Pleistocene glaciation in vicariance. The four distinct subgroups of R. muscosa separated by at least 1.4 million years of evolutionary divergence are suggested as potential units for conservation.