Visualization and quality assessment of de novo genome assemblies.

Recent technological progress has greatly facilitated de novo genome sequencing. However, de novo assemblies consist in many pieces of contiguous sequence (contigs) arranged in thousands of scaffolds instead of small numbers of chromosomes. Confirming and improving the quality of such assemblies is critical for subsequent analysis. We present a method to evaluate genome scaffolding by aligning independently obtained transcriptome sequences to the genome and visually summarizing the alignments using the Cytoscape software. Applying this method to the genome of the red fire ant Solenopsis invicta allowed us to identify inconsistencies in 7%, confirm contig order in 20% and extend 16% of scaffolds.Scripts that generate tables for visualization in Cytoscape from FASTA sequence and scaffolding information files are publicly available at https://github.com/ksanao/TGNet.

The genome of the fire ant Solenopsis invicta.

Ants have evolved very complex societies and are key ecosystem members. Some ants, such as the fire ant Solenopsis invicta, are also major pests. Here, we present a draft genome of S. invicta, assembled from Roche 454 and Illumina sequencing reads obtained from a focal haploid male and his brothers. We used comparative genomic methods to obtain insight into the unique features of the S. invicta genome. For example, we found that this genome harbors four adjacent copies of vitellogenin. A phylogenetic analysis revealed that an ancestral vitellogenin gene first underwent a duplication that was followed by possibly independent duplications of each of the daughter vitellogenins. The vitellogenin genes have undergone subfunctionalization with queen- and worker-specific expression, possibly reflecting differential selection acting on the queen and worker castes. Additionally, we identified more than 400 putative olfactory receptors of which at least 297 are intact. This represents the largest repertoire reported so far in insects. S. invicta also harbors an expansion of a specific family of lipid-processing genes, two putative orthologs to the transformer/feminizer sex differentiation gene, a functional DNA methylation system, and a single putative telomerase ortholog. EST data indicate that this S. invicta telomerase ortholog has at least four spliceforms that differ in their use of two sets of mutually exclusive exons. Some of these and other unique aspects of the fire ant genome are likely linked to the complex social behavior of this species.

Trunk neuromuscular responses to a single whole-body vibration session in patients with chronic low back pain: a cross-sectional study.

OBJECTIVE: Whole-body vibration (WBV) exercise is progressively adopted as an alternative therapeutic modality for enhancing muscle force and muscle activity via neurogenic potentiation. So far, possible changes in the recruitment patterns of the trunk musculature after WBV remain undetermined. The main objective of this study was to evaluate the short-term effects of a single WBV session on trunk neuromuscular responses in patients with chronic low back pain (cLBP) and healthy participants. METHODS: Twenty patients with cLBP and 21 healthy participants performed 10 trunk flexion-extensions before and after a single WBV session consisting of five 1-minute vibration sets. Surface electromyography (EMG) of erector spinae at L2-L3 and L4-L5 and lumbopelvic kinematic variables were collected during the trials. Data were analyzed using 2-way mixed analysis of variance models. RESULTS: The WBV session led to increased lumbar EMG activity during the flexion and extension phases but yielded no change in the quiet standing and fully flexed phases. Kinematic data showed a decreased contribution to the movement of the lumbar region in the second extension quartile. These effects were not different between patients with cLBP and healthy participants. CONCLUSIONS: Increased lumbar EMG activity after a single WBV session most probably results from potentiation effects of WBV on lumbar muscles reflex responses. Decreased EMG activity in full trunk flexion, usually observed in healthy individuals, was still present after WBV, suggesting that the ability of the spine stabilizing mechanisms to transfer the extension torque from muscles to passive structures was not affected.

Genome-Wide Prediction Methods in Highly Diverse and Heterozygous Species: Proof-of-Concept through Simulation in Grapevine.

Nowadays, genome-wide association studies (GWAS) and genomic selection (GS) methods which use genome-wide marker data for phenotype prediction are of much potential interest in plant breeding. However, to our knowledge, no studies have been performed yet on the predictive ability of these methods for structured traits when using training populations with high levels of genetic diversity. Such an example of a highly heterozygous, perennial species is grapevine. The present study compares the accuracy of models based on GWAS or GS alone, or in combination, for predicting simple or complex traits, linked or not with population structure. In order to explore the relevance of these methods in this context, we performed simulations using approx 90,000 SNPs on a population of 3,000 individuals structured into three groups and corresponding to published diversity grapevine data. To estimate the parameters of the prediction models, we defined four training populations of 1,000 individuals, corresponding to these three groups and a core collection. Finally, to estimate the accuracy of the models, we also simulated four breeding populations of 200 individuals. Although prediction accuracy was low when breeding populations were too distant from the training populations, high accuracy levels were obtained using the sole core-collection as training population. The highest prediction accuracy was obtained (up to 0.9) using the combined GWAS-GS model. We thus recommend using the combined prediction model and a core-collection as training population for grapevine breeding or for other important economic crops with the same characteristics.

Genome-wide association study for coronary artery calcification with follow-up in myocardial infarction.

BACKGROUND: Coronary artery calcification (CAC) detected by computed tomography is a noninvasive measure of coronary atherosclerosis, which underlies most cases of myocardial infarction (MI). We sought to identify common genetic variants associated with CAC and further investigate their associations with MI. METHODS AND RESULTS: Computed tomography was used to assess quantity of CAC. A meta-analysis of genome-wide association studies for CAC was performed in 9961 men and women from 5 independent community-based cohorts, with replication in 3 additional independent cohorts (n=6032). We examined the top single-nucleotide polymorphisms (SNPs) associated with CAC quantity for association with MI in multiple large genome-wide association studies of MI. Genome-wide significant associations with CAC for SNPs on chromosome 9p21 near CDKN2A and CDKN2B (top SNP: rs1333049; P=7.58×10(-19)) and 6p24 (top SNP: rs9349379, within the PHACTR1 gene; P=2.65×10(-11)) replicated for CAC and for MI. Additionally, there is evidence for concordance of SNP associations with both CAC and MI at a number of other loci, including 3q22 (MRAS gene), 13q34 (COL4A1/COL4A2 genes), and 1p13 (SORT1 gene). CONCLUSIONS: SNPs in the 9p21 and PHACTR1 gene loci were strongly associated with CAC and MI, and there are suggestive associations with both CAC and MI of SNPs in additional loci. Multiple genetic loci are associated with development of both underlying coronary atherosclerosis and clinical events.

Whole brain radiotherapy with a conformational external beam radiation boost for lung cancer patients with 1-3 brain metastasis: a multi institutional study.

BACKGROUND: To determine the outcome of patients with brain metastasis (BM) from lung cancer treated with an external beam radiotherapy boost (RTB) after whole brain radiotherapy (WBRT). METHODS: A total of 53 BM patients with lung cancer were treated sequentially with WBRT and RTB between 1996 and 2008 according to our institutional protocol. Mean age was 58.8 years. The median KPS was 90. Median recursive partitioning analysis (RPA) and graded prognostic assessment (GPA) grouping were 2 and 2.5, respectively. Surgery was performed on 38 (71%) patients. The median number of BM was 1 (range, 1-3). Median WBRT and RTB combined dose was 39 Gy (range, 37.5-54). Median follow-up was 12.0 months. RESULTS: During the period of follow-up, 37 (70%) patients died. The median overall survival (OS) was 14.5 months. Only 13 patients failed in the brain. The majority of patients (n = 29) failed distantly. The 1-year OS, -local control, extracranial failure rates were 61.2%, 75.2% and 60.8%, respectively. On univariate analysis, improved OS was found to be significantly associated with total dose (< or = 39 Gy vs. > 39 Gy; p < 0.01), age < 65 (p < 0.01), absence of extracranial metastasis (p < 0.01), GPA > or = 2.5 (p = 0.01), KPS > or = 90 (p = 0.01), and RPA < 2 (p = 0.04). On multivariate analysis, total dose (p < 0.01) and the absence of extracranial metastasis (p = 0.03) retained statistical significance. CONCLUSIONS: The majority of lung cancer patients treated with WBRT and RTB progressed extracranially. There might be a subgroup of younger patients with good performance status and no extracranial disease who may benefit from dose escalation after WBRT to the metastatic site.

The genome of the leaf-cutting ant Acromyrmex echinatior suggests key adaptations to advanced social life and fungus farming.

We present a high-quality (>100× depth) Illumina genome sequence of the leaf-cutting ant Acromyrmex echinatior, a model species for symbiosis and reproductive conflict studies. We compare this genome with three previously sequenced genomes of ants from different subfamilies and focus our analyses on aspects of the genome likely to be associated with known evolutionary changes. The first is the specialized fungal diet of A. echinatior, where we find gene loss in the ant's arginine synthesis pathway, loss of detoxification genes, and expansion of a group of peptidase proteins. One of these is a unique ant-derived contribution to the fecal fluid, which otherwise consists of "garden manuring" fungal enzymes that are unaffected by ant digestion. The second is multiple mating of queens and ejaculate competition, which may be associated with a greatly expanded nardilysin-like peptidase gene family. The third is sex determination, where we could identify only a single homolog of the feminizer gene. As other ants and the honeybee have duplications of this gene, we hypothesize that this may partly explain the frequent production of diploid male larvae in A. echinatior. The fourth is the evolution of eusociality, where we find a highly conserved ant-specific profile of neuropeptide genes that may be related to caste determination. These first analyses of the A. echinatior genome indicate that considerable genetic changes are likely to have accompanied the transition from hunter-gathering to agricultural food production 50 million years ago, and the transition from single to multiple queen mating 10 million years ago.

The predictive value of trabecular bone score (TBS) on whole lumbar vertebrae mechanics: an ex vivo study.

We investigated the association of trabecular bone score (TBS) with microarchitecture and mechanical behavior of human lumbar vertebrae. We found that TBS reflects vertebral trabecular microarchitecture and is an independent predictor of vertebral mechanics. However, the addition of TBS to areal BMD (aBMD) did not significantly improve prediction of vertebral strength. INTRODUCTION: The trabecular bone score (TBS) is a gray-level measure of texture using a modified experimental variogram which can be extracted from dual-energy X-ray absorptiometry (DXA) images. The current study aimed to confirm whether TBS is associated with trabecular microarchitecture and mechanics of human lumbar vertebrae, and if its combination with BMD improves prediction of fracture risk. METHODS: Lumbar vertebrae (L3) were harvested fresh from 16 donors. The anteroposterior and lateral bone mineral content (BMC) and areal BMD (aBMD) of the vertebral body were measured using DXA; then, the TBS was extracted using TBS iNsight software (Medimaps SA, France). The trabecular bone volume (Tb.BV/tissue volume, TV), trabecular thickness (Tb.Th), degree of anisotropy, and structure model index (SMI) were measured using microcomputed tomography. Quasi-static uniaxial compressive testing was performed on L3 vertebral bodies to assess failure load and stiffness. RESULTS: The TBS was significantly correlated to Tb.BV/TV and SMI (râeuro0/00=âeuro0/000.58 and -0.62; pâeuro0/00=âeuro0/000.02, 0.01), but not related to BMC and BMD. TBS was significantly correlated with stiffness (râeuro0/00=âeuro0/000.64; pâeuro0/00=âeuro0/000.007), independently of bone mass. Using stepwise multiple regression models, we failed to demonstrate that the combination of BMD and TBS was better at explaining mechanical behavior than either variable alone. However, the combination TBS, Tb.Th, and BMC did perform better than each parameter alone, explaining 79 % of the variability in stiffness. CONCLUSIONS: In our study, TBS was associated with microarchitecture parameters and with vertebral mechanical behavior, but TBS did not improve prediction of vertebral biomechanical properties in addition to aBMD.

Use of 2D Sensitivity Encoding for Slow-Infusion Contrast-Enhanced Isotropic 3-T Whole-Heart Coronary MR Angiography.

OBJECTIVE. The purpose of this study was to improve the blood-pool signal-to-noise ratio (SNR) and blood-myocardium contrast-to-noise ratio (CNR) of slow-infusion 3-T whole-heart coronary MR angiography (MRA).SUBJECTS AND METHODS. In 2D sensitivity encoding (SENSE), the number of acquired k-space lines is reduced, allowing less radiofrequency excitation per cardiac cycle and a longer TR. The former can be exploited for signal enhancement with a higher radiofrequency excitation angle, and the latter leads to noise reduction due to lower data-sampling bandwidth. Both effects contribute to SNR gain in coronary MRA when spatial and temporal resolution and acquisition time remain identical. Numeric simulation was performed to select the optimal 2D SENSE pulse sequence parameters and predict the SNR gain. Eleven patients underwent conventional unenhanced and the proposed 2D SENSE contrast-enhanced coronary MRA acquisition. Blood-pool SNR, blood-myocardium CNR, visible vessel length, vessel sharpness, and number of side branches were evaluated.RESULTS. Consistent with the numeric simulation, using 2D SENSE in contrast-enhanced coronary MRA resulted in significant improvement in aortic blood-pool SNR (unenhanced vs contrast-enhanced, 37.5 +/- 14.7 vs 121.3 +/- 44.0; p < 0.05) and CNR (14.4 +/- 6.9 vs 101.5 +/- 40.8; p < 0.05) in the patient sample. A longer length of left anterior descending coronary artery was visualized, but vessel sharpness, coronary artery coverage, and image quality score were not improved with the proposed approach.CONCLUSION. In combination with contrast administration, 2D SENSE was found effective in improving SNR and CNR in 3-T whole-heart coronary MRA. Further investigation of cardiac motion compensation is necessary to exploit the SNR and CNR advantages and to achieve submillimeter spatial resolution.

A burst of segmental duplications in the genome of the African great ape ancestor

It is generally accepted that the extent of phenotypic change between human and great apes is dissonant with the rate of molecular change. Between these two groups, proteins are virtually identical, cytogenetically there are few rearrangements that distinguish ape-human chromosomes, and rates of single-base-pair change and retrotransposon activity have slowed particularly within hominid lineages when compared to rodents or monkeys. Studies of gene family evolution indicate that gene loss and gain are enriched within the primate lineage. Here, we perform a systematic analysis of duplication content of four primate genomes (macaque, orang-utan, chimpanzee and human) in an effort to understand the pattern and rates of genomic duplication during hominid evolution. We find that the ancestral branch leading to human and African great apes shows the most significant increase in duplication activity both in terms of base pairs and in terms of events. This duplication acceleration within the ancestral species is significant when compared to lineage-specific rate estimates even after accounting for copy-number polymorphism and homoplasy. We discover striking examples of recurrent and independent gene-containing duplications within the gorilla and chimpanzee that are absent in the human lineage. Our results suggest that the evolutionary properties of copy-number mutation differ significantly from other forms of genetic mutation and, in contrast to the hominid slowdown of single-base-pair mutations, there has been a genomic burst of duplication activity at this period during human evolution.

Development and cross-amplification of thirty microsatellite loci in five diploid and polyploid Central Asian species of Palearctic green toads (Bufo viridis subgroup)

We report 30 polymorphic microsatellite markers for five species of Palearctic green toads (Bufo viridis subgroup): 23 in the diploid B. latastii, 19 in diploid B. turanensis, 20 in diploid B. shaartusiensis, 27 in tetraploid B. pewzowi, and 30 in triploid B. baturae. Genetic diversity at these loci, measured for 10-18 individuals per species, ranged from 2 to 19 alleles. These microsatellite loci will be useful for conservation plans (genetic diversity, population structure, evolutionary units), inheritance patterns, and evolution of green toads.

Neuromuscular fatigue induced by whole-body vibration exercise

The aim of this study was to examine the magnitude and the origin of neuromuscular fatigue induced by half-squat static whole-body vibration (WBV) exercise, and to compare it to a non-WBV condition. Nine healthy volunteers completed two fatiguing protocols (WBV and non-WBV, randomly presented) consisting of five 1-min bouts of static half-squat exercise with a load corresponding to 50 % of their individual body mass. Neuromuscular fatigue of knee and ankle muscles was investigated before and immediately after each fatiguing protocol. The main outcomes were maximal voluntary contraction (MVC) torque, voluntary activation, and doublet peak torque. Knee extensor MVC torque decreased significantly (P < 0.01) and to the same extent after WBV (-23 %) and non-WBV (-25 %), while knee flexor, plantar flexor, and dorsiflexor MVC torque was not affected by the treatments. Voluntary activation of knee extensor and plantar flexor muscles was unaffected by the two fatiguing protocols. Doublet peak torque decreased significantly and to a similar extent following WBV and non-WBV exercise, for both knee extensors (-25 %; P < 0.01) and plantar flexors (-7 %; P < 0.05). WBV exercise with additional load did not accentuate fatigue and did not change its causative factors compared to non-WBV half-squat resistive exercise in recreationally active subjects.

Influence of neuroblastoma stage on serum-based detection of MYCN amplification.

BACKGROUND: MYCN oncogene amplification has been defined as the most important prognostic factor for neuroblastoma (NB), the most common solid extracranial neoplasm in children. High copy numbers are strongly associated with rapid tumor progression and poor outcome, independently of tumor stage or patient age, and this has become an important factor in treatment stratification. PROCEDURE: By real-time quantitative PCR analysis, we evaluated the clinical relevance of circulating MYCN DNA of 267 patients with locoregional or metastatic NB in children less than 18 months of age. RESULTS: For patients in this age group with INSS stage 4 or 4S NB and stage 3 patients, serum-based determination of MYCN DNA sequences had good sensitivity (85%, 83%, and 75% respectively) and high specificity (100%) when compared to direct tumor gene determination. In contrast, the approach showed low sensitivity patients with stages 1 and 2 disease. CONCLUSION: Our results show that the sensitivity of the serum-based MYCN DNA sequence determination depends on the stage of the disease. However, this simple, reproducible assay may represent a reasonably sensitive and very specific tool to assess tumor MYCN status in cases with stage 3 and metastatic disease for whom a wait and see strategy is often recommended.

Integration of genome-wide association studies with biological knowledge identifies six novel genes related to kidney function.

In conducting genome-wide association studies (GWAS), analytical approaches leveraging biological information may further understanding of the pathophysiology of clinical traits. To discover novel associations with estimated glomerular filtration rate (eGFR), a measure of kidney function, we developed a strategy for integrating prior biological knowledge into the existing GWAS data for eGFR from the CKDGen Consortium. Our strategy focuses on single nucleotide polymorphism (SNPs) in genes that are connected by functional evidence, determined by literature mining and gene ontology (GO) hierarchies, to genes near previously validated eGFR associations. It then requires association thresholds consistent with multiple testing, and finally evaluates novel candidates by independent replication. Among the samples of European ancestry, we identified a genome-wide significant SNP in FBXL20 (P = 5.6 × 10(-9)) in meta-analysis of all available data, and additional SNPs at the INHBC, LRP2, PLEKHA1, SLC3A2 and SLC7A6 genes meeting multiple-testing corrected significance for replication and overall P-values of 4.5 × 10(-4)-2.2 × 10(-7). Neither the novel PLEKHA1 nor FBXL20 associations, both further supported by association with eGFR among African Americans and with transcript abundance, would have been implicated by eGFR candidate gene approaches. LRP2, encoding the megalin receptor, was identified through connection with the previously known eGFR gene DAB2 and extends understanding of the megalin system in kidney function. These findings highlight integration of existing genome-wide association data with independent biological knowledge to uncover novel candidate eGFR associations, including candidates lacking known connections to kidney-specific pathways. The strategy may also be applicable to other clinical phenotypes, although more testing will be needed to assess its potential for discovery in general.