A new approach for efficient genotype imputation using information from relatives

Background

Genotype imputation can help reduce genotyping costs particularly for implementation of genomic selection. In applications entailing large populations, recovering the genotypes of untyped loci using information from reference individuals that were genotyped with a higher density panel is computationally challenging. Popular imputation methods are based upon the Hidden Markov model and have computational constraints due to an intensive sampling process. A fast, deterministic approach, which makes use of both family and population information, is presented here. All individuals are related and, therefore, share haplotypes which may differ in length and frequency based on their relationships. The method starts with family imputation if pedigree information is available, and then exploits close relationships by searching for long haplotype matches in the reference group using overlapping sliding windows. The search continues as the window size is shrunk in each chromosome sweep in order to capture more distant relationships.

Results

The proposed method gave higher or similar imputation accuracy than Beagle and Impute2 in cattle data sets when all available information was used. When close relatives of target individuals were present in the reference group, the method resulted in higher accuracy compared to the other two methods even when the pedigree was not used. Rare variants were also imputed with higher accuracy. Finally, computing requirements were considerably lower than those of Beagle and Impute2. The presented method took 28 minutes to impute from 6 k to 50 k genotypes for 2,000 individuals with a reference size of 64,429 individuals.

Conclusions

The proposed method efficiently makes use of information from close and distant relatives for accurate genotype imputation. In addition to its high imputation accuracy, the method is fast, owing to its deterministic nature and, therefore, it can easily be used in large data sets where the use of other methods is impractical.     

The Effects of Demography and Long-Term Selection on the Accuracy of Genomic Prediction with Sequence Data

The use of dense SNPs to predict the genetic value of an individual for a complex trait is often referred to as “genomic selection” in livestock and crops, but is also relevant to human genetics to predict, for example, complex genetic disease risk. The accuracy of prediction depends on the strength of linkage disequilibrium (LD) between SNPs and causal mutations. If sequence data were used instead of dense SNPs, accuracy should increase because causal mutations are present, but demographic history and long-term negative selection also influence accuracy. We therefore evaluated genomic prediction, using simulated sequence in two contrasting populations: one reducing from an ancestrally large effective population size (N_e) to a small one, with high LD common in domestic livestock, while the second had a large constant-sized N_e with low LD similar to that in some human or outbred plant populations. There were two scenarios in each population; causal variants were either neutral or under long-term negative selection. For large N_e, sequence data led to a 22% increase in accuracy relative to ∼600K SNP chip data with a Bayesian analysis and a more modest advantage with a BLUP analysis. This advantage increased when causal variants were influenced by negative selection, and accuracy persisted when 10 generations separated reference and validation populations. However, in the reducing N_e population, there was little advantage for sequence even with negative selection. This study demonstrates the joint influence of demography and selection on accuracy of prediction and improves our understanding of how best to exploit sequence for genomic prediction.     

Population Growth in Planaria Dugesia tigrina (Gerard) : Regulation by the absolute number in the population

Planaria reproduce by transverse fission. Isolated worms increase in number exponentially, while social animals at the same density are inhibited in terms of numerical increase, but over a 25 day period undergo a larger increase in mass. Isolated posterior fission products reproduce faster than isolated anterior fission products. Regulation of population growth is independent of density over a 16-fold range and regulatory factors cannot be demonstrated in the medium. Regulation of population growth depends on direct contact between animals. Fission period varies from individual to individual and from period to period for a given individual. Doubling time is related to the absolute number of individuals comprising the population as follows: P_N = (P_M · N)/(K + N), where P_N is the doubling period of a population of N individuals, P_M is the doubling time of an infinitely large population, N is the number of individuals in the population, and K is the number of individuals in a population the period of which is one-half P_M. At 22°–24°C P_M is estimated to be 43.3 days and K is 1.87 individuals. A model system assumes that inhibitor flows through the population from animal to animal from the slowest to the fastest animal in the population thus acting to synchronize population increase as well as to determine the rate of population growth. A possible source of the inhibitor is discussed.     

Two Structurally Different Dienelactone Hydrolases (TfdEI and TfdEII) from Cupriavidus necator JMP134 Plasmid pJP4 Catalyse Cis- and Trans-Dienelactones with Similar Efficiency

In this study, dienelactone hydrolases (TfdEI and TfdEII) located on plasmid pJP4 of Cupriavidus necator JMP134 were cloned, purified, characterized and three dimensional structures were predicted. tfdEI and tfdEII genes were cloned into pET21b vector and expressed in E. coli BL21(DE3). The enzymes were purified by applying ultra-membrane filtration, anion-exchange QFF and gel-filtration columns. The enzyme activity was determined by using cis-dienelactone. The three-dimensional structure of enzymes was predicted using SWISS-MODEL workspace and the biophysical properties were determined on ExPASy server. Both TfdEI and TfdEII (M_r 25 kDa) exhibited optimum activity at 37°C and pH 7.0. The enzymes retained approximately 50% of their activity after 1 h of incubation at 50°C and showed high stability against denaturing agents. The TfdEI and TfdEII hydrolysed cis-dienelactone at a rate of 0.258 and 0.182 µMs⁻¹, with a K_m value of 87 µM and 305 µM, respectively. Also, TfdEI and TfdEII hydrolysed trans-dienelactone at a rate of 0.053 µMs⁻¹ and 0.0766 µMs⁻¹, with a K_m value of 84 µM and 178 µM, respectively. The TfdEI and TfdEII k_cat/K_m ratios were 0.12 µM⁻¹s⁻¹and 0.13 µM⁻¹s⁻¹ and 0.216 µM⁻¹s⁻¹ and 0.094 µM⁻¹s⁻¹ for for cis- and trans-dienelactone, respectively. The k_cat/K_m ratios for cis-dienelactone show that both enzymes catalyse the reaction with same efficiency even though K_m value differs significantly. This is the first report to characterize and compare reaction kinetics of purified TfdEI and TfdEII from Cupriavidus necator JMP134 and may be helpful for further exploration of their catalytic mechanisms.     

Structural and thermodynamic characterization of the recognition of the S100-binding peptides TRTK12 and p53 by calmodulin

Calmodulin (CaM) is a multifunctional messenger protein that activates a wide variety of signaling pathways in eukaryotic cells in a calcium-dependent manner. CaM has been proposed to be functionally distinct from the S100 proteins, a related family of eukaryotic calcium-binding proteins. Previously, it was demonstrated that peptides derived from the actin-capping protein, TRTK12, and the tumor-suppressor protein, p53, interact with multiple members of the S100 proteins. To test the specificity of these peptides, they were screened using isothermal titration calorimetry against 16 members of the human S100 protein family, as well as CaM, which served as a negative control. Interestingly, both the TRTK12 and p53 peptides were found to interact with CaM. These interactions were further confirmed by both fluorescence and nuclear magnetic resonance spectroscopies. These peptides have distinct sequences from the known CaM target sequences. The TRTK12 peptide was found to independently interact with both CaM domains and bind with a stoichiometry of 2:1 and dissociations constants K_d,C-term = 2 ± 1 µM and K_d,N-term = 14 ± 1 µM. In contrast, the p53 peptide was found to interact only with the C-terminal domain of CaM, K_d,C-term =2 ± 1 µM, 25°C. Using NMR spectroscopy, the locations of the peptide binding sites were mapped onto the structure of CaM. The binding sites for both peptides were found to overlap with the binding interface for previously identified targets on both domains of CaM. This study demonstrates the plasticity of CaM in target binding and may suggest a possible overlap in target specificity between CaM and the S100 proteins.     

Transporting Antitumor Drug Tamoxifen and Its Metabolites, 4-Hydroxytamoxifen and Endoxifen by Chitosan Nanoparticles

Synthetic and natural polymers are often used as drug delivery systems in vitro and in vivo. Biodegradable chitosan of different sizes were used to encapsulate antitumor drug tamoxifen (Tam) and its metabolites 4-hydroxytamoxifen (4-Hydroxytam) and endoxifen (Endox). The interactions of tamoxifen and its metabolites with chitosan 15, 100 and 200 KD were investigated in aqueous solution, using FTIR, fluorescence spectroscopic methods and molecular modeling. The structural analysis showed that tamoxifen and its metabolites bind chitosan via both hydrophilic and hydrophobic contacts with overall binding constants of K _tam-ch-15  = 8.7 (±0.5)×10³ M⁻¹, K _tam-ch-100  = 5.9 (±0.4)×10⁵ M⁻¹, K _tam-ch-200  = 2.4 (±0.4)×10⁵ M⁻¹ and K _{hydroxytam-ch-15}  = 2.6(±0.3)×10⁴ M⁻¹, K _{hydroxytam – ch-100}  = 5.2 (±0.7)×10⁶ M⁻¹ and K _{hydroxytam-ch-200}  = 5.1 (±0.5)×10⁵ M⁻¹, K _endox-ch-15  = 4.1 (±0.4)×10³ M⁻¹, K _endox-ch-100  = 1.2 (±0.3)×10⁶ M⁻¹ and K _endox-ch-200  = 4.7 (±0.5)×10⁵ M⁻¹ with the number of drug molecules bound per chitosan (n) 2.8 to 0.5. The order of binding is ch-100>200>15 KD with stronger complexes formed with 4-hydroxytamoxifen than tamoxifen and endoxifen. The molecular modeling showed the participation of polymer charged NH₂ residues with drug OH and NH₂ groups in the drug-polymer adducts. The free binding energies of −3.46 kcal/mol for tamoxifen, −3.54 kcal/mol for 4-hydroxytamoxifen and −3.47 kcal/mol for endoxifen were estimated for these drug-polymer complexes. The results show chitosan 100 KD is stronger carrier for drug delivery than chitosan-15 and chitosan-200 KD.     

Isometric scaling of above- and below-ground biomass at the individual and community levels in the understorey of a sub-tropical forest

Background and Aims Empirical studies and allometric partitioning (AP) theory indicate that plant above-ground biomass (M_A) scales, on average, one-to-one (isometrically) with below-ground biomass (M_R) at the level of individual trees and at the level of entire forest communities. However, the ability of the AP theory to predict the biomass allocation patterns of understorey plants has not been established because most previous empirical tests have focused on canopy tree species or very large shrubs.Methods In order to test the AP theory further, 1586 understorey sub-tropical forest plants from 30 sites in south-east China were harvested and examined. The numerical values of the scaling exponents and normalization constants (i.e. slopes and y-intercepts, respectively) of log–log linear M_A vs. M_R relationships were determined for all individual plants, for each site, across the entire data set, and for data sorted into a total of 19 sub-sets of forest types and successional stages. Similar comparisons of M_A/M_R were also made.Key Results The data revealed that the mean M_A/M_R of understorey plants was 2·44 and 1·57 across all 1586 plants and for all communities, respectively, and M_A scaled nearly isometrically with respect to M_R, with scaling exponents of 1·01 for all individual plants and 0·99 for all communities. The scaling exponents did not differ significantly among different forest types or successional stages, but the normalization constants did, and were positively correlated with M_A/M_R and negatively correlated with scaling exponents across all 1586 plants.Conclusions The results support the AP theory’s prediction that M_A scales nearly one-to-one with M_R (i.e. M_A ∝ M_R ^≈1·0) and that plant biomass partitioning for individual plants and at the community level share a strikingly similar pattern, at least for the understorey plants examined in this study. Furthermore, variation in environmental conditions appears to affect the numerical values of normalization constants, but not the scaling exponents of the M_A vs. M_R relationship. This feature of the results suggests that plant size is the primary driver of the M_A vs. M_R biomass allocation pattern for understorey plants in sub-tropical forests.     

Haplotypes at LBX1 Have Distinct Inheritance Patterns with Opposite Effects in Adolescent Idiopathic Scoliosis

Adolescent idiopathic scoliosis (AIS) is a clinically significant disorder with high heritability that affects 2–4% of the population. Genome-wide association studies have identified LBX1 as a strong susceptibility locus for AIS in Asian and Caucasian populations. Here we further dissect the genetic association with AIS in a Caucasian population. To identify genetic markers associated with AIS we employed a genome-wide association study (GWAS) design comparing 620 female Caucasian patients who developed idiopathic scoliosis during adolescence with 1,287 ethnically matched females who had normal spinal curves by skeletal maturity. The genomic region around LBX1 was imputed and haplotypes investigated for genetic signals under different inheritance models. The strongest signal was identified upstream of LBX1 (rs11190878, P_trend = 4.18×10^-9, OR = 0.63[0.54–0.74]). None of the remaining SNPs pass the genome-wide significance threshold. We found rs11190870, downstream of LBX1 and previously associated with AIS in Asian populations, to be in modest linkage disequilibrium (LD) with rs11190878 (r² = 0.40, D'' = 0.81). Haplotype analysis shows that rs11190870 and rs11190878 track a single risk factor that resides on the ancestral haplotype and is shared across ethnic groups. We identify six haplotypes at the LBX1 locus including two strongly associated haplotypes; a recessive risk haplotype (TTA, Control_freq = 0.52, P = 1.25×10^-9, OR = 1.56), and a co-dominant protective haplotype (CCG, Control_freq = 0.28, P = 2.75×10^-7, OR = 0.65). Together the association signals from LBX1 explain 1.4% of phenotypic variance. Our results identify two clinically relevant haplotypes in the LBX1-region with opposite effects on AIS risk. The study demonstrates the utility of haplotypes over un-phased SNPs for individualized risk assessment by more strongly delineating individuals at risk for AIS without compromising the effect size.     

Repeatability of adaptation in experimental populations of different sizes

The degree to which evolutionary trajectories and outcomes are repeatable across independent populations depends on the relative contribution of selection, chance and history. Population size has been shown theoretically and empirically to affect the amount of variation that arises among independent populations adapting to the same environment. Here, we measure the contribution of selection, chance and history in different-sized experimental populations of the unicellular alga Chlamydomonas reinhardtii adapting to a high salt environment to determine which component of evolution is affected by population size. We find that adaptation to salt is repeatable at the fitness level in medium (N_e = 5 × 10⁴) and large (N_e = 4 × 10⁵) populations because of the large contribution of selection. Adaptation is not repeatable in small (N_e = 5 × 10³) populations because of large constraints from history. The threshold between stochastic and deterministic evolution in this case is therefore between effective population sizes of 10³ and 10⁴. Our results indicate that diversity across populations is more likely to be maintained if they are small. Experimental outcomes in large populations are likely to be robust and can inform our predictions about outcomes in similar situations.     

Genetic signature of the natural gene pool of Tilia cordata Mill. in Lithuania: Compound evolutionary and anthropogenic effects

Tilia cordata Mill. is a valuable tree species enriching the ecological values of the coniferous‐dominated boreal forests in Europe. Following the historical decline, spreading of Tilia sp. is challenged by the elevated inbreeding and habitat fragmentation. We studied the geographical distribution of genetic diversity of Tilia cordata populations in Lithuania. We used 14 genomic microsatellite markers to genotype 543 individuals from 23 wild‐growing populations. We found that Tilia cordata retained high levels of genetic diversity (population F _IS = 0–0.15, H _o = 0.53–0.69, H _e = 0.56–0.75). AMOVA, Bayesian clustering, and Monmonier''s barrier detection indicate weak but significant differentiation among the populations (F _ST = 0.037***) into geographically interpretable clusters of (a) western Lithuania with high genetic heterogeneity but low genetic diversity, bottleneck effects, (b) relatively higher genetic diversity of Tilia cordata on rich and most soils of midland lowland, and (c) the most differentiated populations on poor soils of the coolest northeastern highland possessing the highest rare allele frequency but elevated inbreeding and bottleneck effects. Weak genetic differentiation among the Tilia cordata populations in Lithuania implies common ancestry, absence of strong adaptive gradients, and effective genetic exchange possible mediated via the riparian networks. A hypothesis on riparian networks as gene flow mediators in Tilia cordata was raised based on results of this study.     

The Effects of Temperature and Body Mass on Jump Performance of the Locust Locusta migratoria

Locusts jump by rapidly releasing energy from cuticular springs built into the hind femur that deform when the femur muscle contracts. This study is the first to examine the effect of temperature on jump energy at each life stage of any orthopteran. Ballistics and high-speed cinematography were used to quantify the energy, distance, and take-off angle of the jump at 15, 25, and 35°C in the locust Locusta migratoria. Allometric analysis across the five juvenile stages at 35°C reveals that jump distance (D; m) scales with body mass (M; g) according to the power equation D = 0.35M ^{0.17±0.08 (95% CI)}, jump take-off angle (A; degrees) scales as A = 52.5M ^0.00±0.06, and jump energy (E; mJ per jump) scales as E = 1.91M ^1.14±0.09. Temperature has no significant effect on the exponent of these relationships, and only a modest effect on the elevation, with an overall Q₁₀ of 1.08 for jump distance and 1.09 for jump energy. On average, adults jump 87% farther and with 74% more energy than predicted based on juvenile scaling data. The positive allometric scaling of jump distance and jump energy across the juvenile life stages is likely facilitated by the concomitant relative increase in the total length (L _f+t; mm) of the femur and tibia of the hind leg, L _f+t = 34.9M ^0.37±0.02. The weak temperature-dependence of jump performance can be traced to the maximum tension of the hind femur muscle and the energy storage capacity of the femur''s cuticular springs. The disproportionately greater jump energy and jump distance of adults is associated with relatively longer (12%) legs and a relatively larger (11%) femur muscle cross-sectional area, which could allow more strain loading into the femur''s cuticular springs. Augmented jump performance in volant adult locusts achieves the take-off velocity required to initiate flight.     

Length-dependent structure formation in Friedreich ataxia (GAA)n*(TTC)n repeats at neutral pH

More than 15 human genetic diseases have been associated with the expansion of trinucleotide DNA repeats, which may involve the formation of non-duplex DNA structures. The slipped-strand nucleation of duplex DNA within GC-rich trinucleotide repeats may result in the changes of repeat length; however, such a mechanism seems less likely for the AT-rich (GAA)_n·(TTC)_n repeats. Using two-dimensional agarose gels, chemical probing and atomic force microscopy, we characterized the formation of non-B-DNA structures in the Friedreich ataxia-associated (GAA)_n·(TTC)_n repeats from the FRDA gene that were cloned with flanking genomic sequences into plasmids. For the normal genomic repeat length (n = 9) our data are consistent with the formation of a very stable protonated intramolecular triplex (H-DNA). Its stability at pH 7.4 is likely due to the high proportion of the T·A·T triads which form within the repeats as well as in the immediately adjacent AT-rich sequences with a homopurine· homopyrimidine bias. At the long normal repeat length (n = 23), a family of H-DNAs of slightly different sizes has been detected. At the premutation repeat length (n = 42) and higher negative supercoiling, the formation of a single H-DNA structure becomes less favorable and the data are consistent with the formation of a bi-triplex structure.     

Accuracy of genotype imputation in Nelore cattle

Background

Genotype imputation from low-density (LD) to high-density single nucleotide polymorphism (SNP) chips is an important step before applying genomic selection, since denser chips tend to provide more reliable genomic predictions. Imputation methods rely partially on linkage disequilibrium between markers to infer unobserved genotypes. Bos indicus cattle (e.g. Nelore breed) are characterized, in general, by lower levels of linkage disequilibrium between genetic markers at short distances, compared to taurine breeds. Thus, it is important to evaluate the accuracy of imputation to better define which imputation method and chip are most appropriate for genomic applications in indicine breeds.

Methods

Accuracy of genotype imputation in Nelore cattle was evaluated using different LD chips, imputation software and sets of animals. Twelve commercial and customized LD chips with densities ranging from 7 K to 75 K were tested. Customized LD chips were virtually designed taking into account minor allele frequency, linkage disequilibrium and distance between markers. Software programs FImpute and BEAGLE were applied to impute genotypes. From 995 bulls and 1247 cows that were genotyped with the Illumina® BovineHD chip (HD), 793 sires composed the reference set, and the remaining 202 younger sires and all the cows composed two separate validation sets for which genotypes were masked except for the SNPs of the LD chip that were to be tested.

Results

Imputation accuracy increased with the SNP density of the LD chip. However, the gain in accuracy with LD chips with more than 15 K SNPs was relatively small because accuracy was already high at this density. Commercial and customized LD chips with equivalent densities presented similar results. FImpute outperformed BEAGLE for all LD chips and validation sets. Regardless of the imputation software used, accuracy tended to increase as the relatedness between imputed and reference animals increased, especially for the 7 K chip.

Conclusions

If the Illumina® BovineHD is considered as the target chip for genomic applications in the Nelore breed, cost-effectiveness can be improved by genotyping part of the animals with a chip containing around 15 K useful SNPs and imputing their high-density missing genotypes with FImpute.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-014-0069-1) contains supplementary material, which is available to authorized users.     

Thermal performance of the Chagas disease vector,Triatoma infestans,under thermal variability

Vector-borne diseases (VBD) are particularly susceptible to climate change because most of the diseases’ vectors are ectotherms, which themselves are susceptible to thermal changes. The Chagas disease is one neglected tropical disease caused by the protozoan parasite, Trypanosoma cruzi. One of the main vectors of the Chagas disease in South America is Triatoma infestans, a species traditionally considered to be restricted to domestic or peridomestic habitats, but sylvatic foci have also been described along its distribution. The infestation of wild individuals, together with the projections of environmental changes due to global warming, urge the need to understand the relationship between temperature and the vector’s performance. Here, we evaluated the impact of temperature variability on the thermal response of T. infestans. We acclimated individuals to six thermal treatments for five weeks to then estimate their thermal performance curves (TPCs) by measuring the walking speed of the individuals. We found that the TPCs varied with thermal acclimation and body mass. Individuals acclimated to a low and variable ambient temperature (18°C ± 5°C) exhibited lower performances than those individuals acclimated to an optimal temperature (27°C ± 0°C); while those individuals acclimated to a low but constant temperature (18°C ± 0°C) did not differ in their maximal performance from those at an optimal temperature. Additionally, thermal variability (i.e., ± 5°C) at a high temperature (30°C) increased performance. These results evidenced the plastic response of T. infestans to thermal acclimation. This plastic response and the non-linear effect of thermal variability on the performance of T. infestans posit challenges when predicting changes in the vector’s distribution range under climate change.     

Low-temperature leaf photosynthesis of a Miscanthus germplasm collection correlates positively to shoot growth rate and specific leaf area

Background and Aims The C₄ perennial grass miscanthus has been found to be less sensitive to cold than most other C₄ species, but still emerges later in spring than C₃ species. Genotypic differences in miscanthus were investigated to identify genotypes with a high cold tolerance at low temperatures and quick recovery upon rising temperatures to enable them to exploit the early growing season in maritime cold climates. Suitable methods for field screening of cold tolerance in miscanthus were also identified.Methods Fourteen genotypes of M. sacchariflorus, M. sinensis, M. tinctorius and M. × giganteus were selected and grown under warm (24 °C) and cold (14 °C) conditions in a controlled environment. Dark-adapted chlorophyll fluorescence, specific leaf area (SLA) and net photosynthetic rate at a photosynthetically active radiation (PAR) of 1000 μmol m^–2 s^–1 (A₁₀₀₀) were measured. Photosynthetic light and CO₂ response curves were obtained from 11 of the genotypes, and shoot growth rate was measured under field conditions.Key Results A positive linear relationship was found between SLA and light-saturated photosynthesis (A_sat) across genotypes, and also between shoot growth rate under cool field conditions and A₁₀₀₀ at 14 °C in a climate chamber. When lowering the temperature from 24 to 14 °C, one M. sacchariflorus exhibited significantly higher A_sat and maximum photosynthetic rate in the CO₂ response curve (V_max) than other genotypes at 14 °C, except M. × giganteus ‘Hornum’. Several genotypes returned to their pre-chilling A₁₀₀₀ values when the temperature was increased to 24 °C after 24 d growth at 14 °C.Conclusions One M. sacchariflorus genotype had similar or higher photosynthetic capacity than M. × giganteus, and may be used for cultivation together with M. × giganteus or for breeding new interspecies hybrids with improved traits for temperate climates. Two easily measured variables, SLA and shoot growth rate, may be useful for genotype screening of productivity and cold tolerance.     

Kinetics of Phosphomevalonate Kinase from Saccharomyces cerevisiae

The mevalonate-based isoprenoid biosynthetic pathway is responsible for producing cholesterol in humans and is used commercially to produce drugs, chemicals, and fuels. Heterologous expression of this pathway in Escherichia coli has enabled high-level production of the antimalarial drug artemisinin and the proposed biofuel bisabolane. Understanding the kinetics of the enzymes in the biosynthetic pathway is critical to optimize the pathway for high flux. We have characterized the kinetic parameters of phosphomevalonate kinase (PMK, EC 2.7.4.2) from Saccharomyces cerevisiae, a previously unstudied enzyme. An E. coli codon-optimized version of the S. cerevisiae gene was cloned into pET-52b+, then the C-terminal 6X His-tagged protein was expressed in E. coli BL21(DE3) and purified on a Ni²⁺ column. The K_M of the ATP binding site was determined to be 98.3 µM at 30°C, the optimal growth temperature for S. cerevisiae, and 74.3 µM at 37°C, the optimal growth temperature for E. coli. The K_M of the mevalonate-5-phosphate binding site was determined to be 885 µM at 30°C and 880 µM at 37°C. The V_max was determined to be 4.51 µmol/min/mg enzyme at 30°C and 5.33 µmol/min/mg enzyme at 37°C. PMK is Mg²⁺ dependent, with maximal activity achieved at concentrations of 10 mM or greater. Maximum activity was observed at pH = 7.2. PMK was not found to be substrate inhibited, nor feedback inhibited by FPP at concentrations up to 10 µM FPP.     

Mitochondrial DNA D‐loop sequence analysis reveals high variation and multiple maternal origins of indigenous Tanzanian goat populations

The Small East African (SEA) goat are widely distributed in different agro‐ecological zones of Tanzania. We report the genetic diversity, maternal origin, and phylogenetic relationship among the 12 Tanzanian indigenous goat populations, namely Fipa, Songwe, Tanga, Pwani, Iringa, Newala, Lindi, Gogo, Pare, Maasai, Sukuma, and Ujiji, based on the mitochondrial DNA (mtDNA) D‐loop. High haplotype (H _d = 0.9619–0.9945) and nucleotide (π = 0.0120–0.0162) diversities were observed from a total of 389 haplotypes. The majority of the haplotypes (n = 334) belonged to Haplogroup A which was consistent with the global scenario on the genetic pattern of maternal origin of all goat breeds in the world. Haplogroup G comprised of 45 haplotypes drawn from all populations except the Ujiji goat population while Haplogroup B with 10 haplotypes was dominated by Ujiji goats (41%). Tanzanian goats shared four haplotypes with the Kenyan goats and two with goats from South Africa, Namibia, and Mozambique. There was no sharing of haplotypes observed between individuals from Tanzanian goat populations with individuals from North or West Africa. The indigenous goats in Tanzania have high genetic diversity defined by 389 haplotypes and multiple maternal origins of haplogroup A, B, and G. There is a lot of intermixing and high genetic variation within populations which represent an abundant resource for selective breeding in the different agro‐ecological regions of the country.     

Structural analysis of DNA complexation with cationic lipids

Complexes of cationic liposomes with DNA are promising tools to deliver genetic information into cells for gene therapy and vaccines. Electrostatic interaction is thought to be the major force in lipid–DNA interaction, while lipid-base binding and the stability of cationic lipid–DNA complexes have been the subject of more debate in recent years. The aim of this study was to examine the complexation of calf-thymus DNA with cholesterol (Chol), 1,2-dioleoyl-3-trimethylammonium-propane (DOTAP), dioctadecyldimethylammoniumbromide (DDAB) and dioleoylphosphatidylethanolamine (DOPE), at physiological condition, using constant DNA concentration and various lipid contents. Fourier transform infrared (FTIR), UV-visible, circular dichroism spectroscopic methods and atomic force microscopy were used to analyse lipid-binding site, the binding constant and the effects of lipid interaction on DNA stability and conformation. Structural analysis showed a strong lipid–DNA interaction via major and minor grooves and the backbone phosphate group with overall binding constants of K_Chol = 1.4 (±0.5) × 10⁴ M⁻¹, K_DDAB = 2.4 (±0.80) × 10⁴ M⁻¹, K_DOTAP = 3.1 (±0.90) × 10⁴ M⁻¹ and K_DOPE = 1.45 (± 0.60) × 10⁴ M⁻¹. The order of stability of lipid–DNA complexation is DOTAP>DDAB>DOPE>Chol. Hydrophobic interactions between lipid aliphatic tails and DNA were observed. Chol and DOPE induced a partial B to A-DNA conformational transition, while a partial B to C-DNA alteration occurred for DDAB and DOTAP at high lipid concentrations. DNA aggregation was observed at high lipid content.     

Cloning of the xynB Gene from Dictyoglomus thermophilum Rt46B.1 and Action of the Gene Product on Kraft Pulp

A two-step PCR protocol was used to identify and sequence a family 11 xylanase gene from Dictyoglomus thermophilum Rt46B.1. Family 11 xylanase consensus fragments (GXCFs) were amplified from Rt46B.1 genomic DNA by using different sets of consensus PCR primers that exhibited broad specificity for conserved motifs within fungal and/or bacterial family 11 xylanase genes. On the basis of the sequences of a representative sample of the GXCFs a single family 11 xylanase gene (xynB) was identified. The entire gene sequence was obtained in the second step by using genomic walking PCR to amplify Rt46B.1 genomic DNA fragments upstream and downstream of the xynB GXCF region. The putative XynB peptide (M_r, 39,800) encoded by the Rt46B.1 xynB open reading frame was a multidomain enzyme comprising an N-terminal catalytic domain (M_r, 22,000) and a possible C-terminal substrate-binding domain (M_r, 13,000) that were separated by a short serine-glycine-rich 23-amino-acid linker peptide. Seven xylanases which differed at their N and C termini were produced from different xynB expression plasmids. All seven xylanases exhibited optimum activity at pH 6.5. However, the temperature optima of the XynB xylanases varied from 70 to 85°C. Pretreatment of Pinus radiata and eucalypt kraft-oxygen pulps with XynB resulted in moderate xylan solubilization and a substantial improvement in the bleachability of these pulps.     

The physical properties of inulin solutions

1. Carefully recrystallized samples of inulin have been analysed chromatographically for low-molecular-weight contaminants. The analysis of the samples revealed amounts of low-molecular-weight polyfructosans amounting to 5–10% of the inulin. 2. The possible polymorphism of inulin, its solubility behaviour and suggestions of association in solution are discussed. 3. Physical techniques suggest considerable polydispersity. Ultracentrifugation gives a weight-average molecular weight (M_w) 7250, whereas osmotic-pressure measurements suggest a number-average molecular weight (M_n) 5600. 4. A probable steric model is suggested of a helix repeating every four residues, of diameter 12Å and axial ratio 7·5–10. 5. Fractionation of the inulin on columns of Sephadex G-25 show a pronounced polydispersity. 6. Some doubt is cast on the reliability of inulin as an indicator of physiological volumes of distribution in body fluids.