Selecting high-dimensional mixed graphical models using minimal AIC or BIC forests

Background  

Chow and Liu showed that the maximum likelihood tree for multivariate discrete distributions may be found using a maximum weight spanning tree algorithm, for example Kruskal's algorithm. The efficiency of the algorithm makes it tractable for high-dimensional problems.     

Enlarging a training set for genomic selection by imputation of un-genotyped animals in populations of varying genetic architecture

Background

The most common application of imputation is to infer genotypes of a high-density panel of markers on animals that are genotyped for a low-density panel. However, the increase in accuracy of genomic predictions resulting from an increase in the number of markers tends to reach a plateau beyond a certain density. Another application of imputation is to increase the size of the training set with un-genotyped animals. This strategy can be particularly successful when a set of closely related individuals are genotyped.

Methods

Imputation on completely un-genotyped dams was performed using known genotypes from the sire of each dam, one offspring and the offspring’s sire. Two methods were applied based on either allele or haplotype frequencies to infer genotypes at ambiguous loci. Results of these methods and of two available software packages were compared. Quality of imputation under different population structures was assessed. The impact of using imputed dams to enlarge training sets on the accuracy of genomic predictions was evaluated for different populations, heritabilities and sizes of training sets.

Results

Imputation accuracy ranged from 0.52 to 0.93 depending on the population structure and the method used. The method that used allele frequencies performed better than the method based on haplotype frequencies. Accuracy of imputation was higher for populations with higher levels of linkage disequilibrium and with larger proportions of markers with more extreme allele frequencies. Inclusion of imputed dams in the training set increased the accuracy of genomic predictions. Gains in accuracy ranged from close to zero to 37.14%, depending on the simulated scenario. Generally, the larger the accuracy already obtained with the genotyped training set, the lower the increase in accuracy achieved by adding imputed dams.

Conclusions

Whenever a reference population resembling the family configuration considered here is available, imputation can be used to achieve an extra increase in accuracy of genomic predictions by enlarging the training set with completely un-genotyped dams. This strategy was shown to be particularly useful for populations with lower levels of linkage disequilibrium, for genomic selection on traits with low heritability, and for species or breeds for which the size of the reference population is limited.     

Bacterial colonization of Hydra hatchlings follows a robust temporal pattern

Animals are colonized by complex bacterial communities. The processes controlling community membership and influencing the establishment of the microbial ecosystem during development are poorly understood. Here we aimed to explore the assembly of bacterial communities in Hydra with the broader goal of elucidating the general rules that determine the temporal progression of bacterial colonization of animal epithelia. We profiled the microbial communities in polyps at various time points after hatching in four replicates. The composition and temporal patterns of the bacterial communities were strikingly similar in all replicates. Distinct features included high diversity of community profiles in the first week, a remarkable but transient adult-like profile 2 weeks after hatching, followed by progressive emergence of a stable adult-like pattern characterized by low species diversity and the preponderance of the Betaproteobacterium Curvibacter. Intriguingly, this process displayed important parallels to the assembly of human fecal communities after birth. In addition, a mathematical modeling approach was used to uncover the organizational principles of this colonization process, suggesting that both, local environmental or host-derived factor(s) modulating the colonization rate, as well as frequency-dependent interactions of individual bacterial community members are important aspects in the emergence of a stable bacterial community at the end of development.     

Does small-sided-games’ court area influence metabolic,perceptual, and physical performance parameters of young elite basketball players?

The purpose of this study was to investigate the effect of court size on physiological responses and physical performance of young elite basketball players. Twelve male basketball players (18.6 ± 0.5 years; 88.8 ± 14.5 kg; 192.6 ± 6.5 cm) from an under-19 team performed two small-sided games (matches) with different court areas (28x15 m and 28x9 m; 28x15 and 28x9 protocols). The number of players (3x3) was kept the same in each protocol. The players performed a repeated-sprint ability (RSA) test before and after each match. Blood lactate concentration was collected before (pre) and after (post) the matches, and the session rating of perceived exertion (session-RPE) was determined 30 minutes after the match. Best and mean time in the RSA test were not different between the 28x15 and the 28x9 match protocols (p > 0.05). A significant difference was observed for lactate concentration from pre- to post-match (p < 0.05) in both protocols (28x15 and 28x9); however, there was no significant interaction between protocols. A similar session-RPE mean score (28x15: 7.2 ± 1.4 and 28x9: 6.6 ± 1.4) was detected for both protocols (p > 0.05, ES=0.41). In summary, the results of the current study suggest that the different court areas induced similar responses. Although there was no significant difference in effort perception, players tended to perceive a greater effort in the larger court size.     

Data partitioning and correction for ascertainment bias reduce the uncertainty of placental mammal divergence times inferred from the morphological clock

Bayesian estimates of divergence times based on the molecular clock yield uncertainty of parameter estimates measured by the width of posterior distributions of node ages. For the relaxed molecular clock, previous works have reported that some of the uncertainty inherent to the variation of rates among lineages may be reduced by partitioning data. Here we test this effect for the purely morphological clock, using placental mammals as a case study. We applied the uncorrelated lognormal relaxed clock to morphological data of 40 extant mammalian taxa and 4,533 characters, taken from the largest published matrix of discrete phenotypic characters. The morphologically derived timescale was compared to divergence times inferred from molecular and combined data. We show that partitioning data into anatomical units significantly reduced the uncertainty of divergence time estimates for morphological data. For the first time, we demonstrate that ascertainment bias has an impact on the precision of morphological clock estimates. While analyses including molecular data suggested most divergences between placental orders occurred near the K‐Pg boundary, the partitioned morphological clock recovered older interordinal splits and some younger intraordinal ones, including significantly later dates for the radiation of bats and rodents, which accord to the short‐fuse hypothesis.     

Timing the origin of New World monkeys

The origin of New World monkeys (Infraorder Platyrrhini) has been an extensively debated issue. In this study, we analyzed mitochondrial genomes from Cebus (Platyrrhini), Homo, Hylobates, Pan, Pongo (Hominoids), Macaca, Papio (Cercopithecoids), and Tarsius (outgroup) to investigate this matter. Two distinct methodologies were employed on mitochondrial genes to estimate divergence times: the traditional likelihood ratio test performed in ML analyses of individual and concatenated gene sequences and the recent multigene Bayesian approach. Using the Cercopithecoid-Hominoid split as calibration point (25 MYA), our results show consistently that Platyrrhines split from Catarrhines at around 35 MYA. Although the main focus of the study is New World monkey origins, we have also estimated other primate divergence times: Homo-Pan at 5-7 MYA; Pongo-(Homo/Pan) at 13-16 MYA; Hylobates-(Pongo/Homo/Pan) at 15-19 MYA; and Macaca-Papio at 10-12 MYA. Our estimate for the origin of New World monkeys is in agreement with the hypothesis of a transatlantic journey from Africa to South America, as suggested by the fossil record.     

Diversification of the Genus Anopheles and a Neotropical Clade from the Late Cretaceous

The Anopheles genus is a member of the Culicidae family and consists of approximately 460 recognized species. The genus is composed of 7 subgenera with diverse geographical distributions. Despite its huge medical importance, a consensus has not been reached on the phylogenetic relationships among Anopheles subgenera. We assembled a comprehensive dataset comprising the COI, COII and 5.8S rRNA genes and used maximum likelihood and Bayesian inference to estimate the phylogeny and divergence times of six out of the seven Anopheles subgenera. Our analysis reveals a monophyletic group composed of the three exclusively Neotropical subgenera, Stethomyia, Kerteszia and Nyssorhynchus, which began to diversify in the Late Cretaceous, at approximately 90 Ma. The inferred age of the last common ancestor of the Anopheles genus was ca. 110 Ma. The monophyly of all Anopheles subgenera was supported, although we failed to recover a significant level of statistical support for the monophyly of the Anopheles genus. The ages of the last common ancestors of the Neotropical clade and the Anopheles and Cellia subgenera were inferred to be at the Late Cretaceous (ca. 90 Ma). Our analysis failed to statistically support the monophyly of the Anopheles genus because of an unresolved polytomy between Bironella and A. squamifemur.     

Stringent and reproducible tetracycline-regulated transgene expression by site-specific insertion at chromosomal loci with pre-characterised induction characteristics

Background  

The ability to regulate transgene expression has many applications, mostly concerning the analysis of gene function. Desirable induction characteristics, such as low un-induced expression, high induced expression and limited cellular heterogeneity, can be seriously impaired by chromosomal position effects at the site of transgene integration. Many clones may therefore need to be screened before one with optimal induction characteristics is identified. Furthermore, such screens must be repeated for each new transgene investigated, and comparisons between clones with different transgenes is complicated by their different integration sites.     

MamMiBase: a mitochondrial genome database for mammalian phylogenetic studies

MamMiBase, the mammalian mitochondrial genome database, is a relational database of complete mitochondrial genome sequences of mammalian species. The database is useful for phylogenetic analysis, since it allows a ready retrieval of nucleotide and aminoacid individual alignments, in three different formats (NEXUS for PAUP program, for MEGA program and for PHYLIP program) of the 13 protein coding mitochondrial genes. The user may download the sequences that are useful for him/her based on their parameters values, such as sequence length, p-distances, base content, transition transversion ratio, gamma, which are also given by MamMiBase. A simple phylogenetic tree (neighbor-joining tree with Jukes Cantor distance) is also available for download, useful for parameter calculations and other simple tasks. AVAILABILITY: MamMiBase is available at http://www.mammibase.lncc.br     

An IGF-I promoter polymorphism modifies the relationships between birth weight and risk factors for cardiovascular disease and diabetes at age 36

BACKGROUND: Biochemical testing for pheochromocytoma by measurement of fractionated plasma metanephrines is limited by false positive rates of up to 18% in people without known genetic predisposition to the disease. The plasma normetanephrine fraction is responsible for most false positives and plasma normetanephrine increases with age. The objective of this study was to determine if we could improve the specificity of fractionated plasma measurements, by statistically adjusting for age. METHODS: An age-adjusted metanephrine score was derived using logistic regression from 343 subjects (including 33 people with pheochromocytoma) who underwent fractionated plasma metanephrine measurements as part of investigations for suspected pheochromocytoma at Mayo Clinic Rochester (derivation set). The performance of the age-adjusted score was validated in a dataset of 158 subjects (including patients 23 with pheochromocytoma) that underwent measurements of fractionated plasma metanephrines at Mayo Clinic the following year (validation dataset). None of the participants in the validation dataset had known genetic predisposition to pheochromocytoma. RESULTS: The sensitivity of the age-adjusted metanephrine score was the same as that of traditional interpretation of fractionated plasma metanephrine measurements, yielding a sensitivity of 100% (23/23, 95% confidence interval [CI] 85.7%, 100%). However, the false positive rate with traditional interpretation of fractionated plasma metanephrine measurements was 16.3% (22/135, 95% CI, 11.0%, 23.4%) and that of the age-adjusted score was significantly lower at 3.0% (4/135, 95% CI, 1.2%, 7.4%) (p < 0.001 using McNemar's test). CONCLUSION: An adjustment for age in the interpretation of results of fractionated plasma metanephrines may significantly decrease false positives when using this test to exclude sporadic pheochromocytoma. Such improvements in false positive rate may result in savings of expenditures related to confirmatory imaging.