Improving Illumina assemblies with Hi-C and long reads : An example with the North African dromedary

Researchers have assembled thousands of eukaryotic genomes using Illumina reads, but traditional mate-pair libraries cannot span all repetitive elements, resulting in highly fragmented assemblies. However, both chromosome conformation capture techniques, such as Hi-C and Dovetail Genomics Chicago libraries and long-read sequencing, such as Pacific Biosciences and Oxford Nanopore, help span and resolve repetitive regions and therefore improve genome assemblies. One important livestock species of arid regions that does not have a high-quality contiguous reference genome is the dromedary (Camelus dromedarius). Draft genomes exist but are highly fragmented, and a high-quality reference genome is needed to understand adaptation to desert environments and artificial selection during domestication. Dromedaries are among the last livestock species to have been domesticated, and together with wild and domestic Bactrian camels, they are the only representatives of the Camelini tribe, which highlights their evolutionary significance. Here we describe our efforts to improve the North African dromedary genome. We used Chicago and Hi-C sequencing libraries from Dovetail Genomics to resolve the order of previously assembled contigs, producing almost chromosome-level scaffolds. Remaining gaps were filled with Pacific Biosciences long reads, and then scaffolds were comparatively mapped to chromosomes. Long reads added 99.32 Mbp to the total length of the new assembly. Dovetail Chicago and Hi-C libraries increased the longest scaffold over 12-fold, from 9.71 Mbp to 124.99 Mbp and the scaffold N50 over 50-fold, from 1.48 Mbp to 75.02 Mbp. We demonstrate that Illumina de novo assemblies can be substantially upgraded by combining chromosome conformation capture and long-read sequencing.Peer reviewe

Assessing antibody decline after chemotherapy of early chronic Chagas disease patients

Background: Chagas disease remains a significant public health problem in Latin America. There are only two chemotherapy drugs, nifurtimox and benznidazole, and both may have severe side effects. After complete chemotherapy of acute cases, seropositive diagnosis may revert to negative. However, there are no definitive parasitological or serological biomarkers of cure. Methods: Following a pilot study with seven Bolivian migrants to Spain, we tested 71 serum samples from chronic patients (mean age 12.6 years) inhabiting the Argentine Chaco region. Benznidazole chemotherapy (5–8 mg/kg day, twice daily for 60 days) was administered during 2011–2016. Subsequently, pre-and post-chemotherapy serum samples were analysed in pairs by IgG1 and IgG ELISA using two different antigens and Chagas Sero K-SeT rapid diagnostic tests (RDT). Molecular diagnosis by kDNA-PCR was applied to post-treatment samples. Results: Pilot data demonstrated IgG1 antibody decline in three of seven patients from Bolivia 1 year post-treatment. All Argentine patients in 2017 (averaging 5 years post-treatment), except one, were positive by conventional serology. All were kDNA-PCR-negative. Most (91.5%) pre-treatment samples were positive by the Chagas Sero K-SeT RDT, confirming the predominance of TcII/V/VI. IgG1 and IgG of Argentine patients showed significant decline in antibody titres post-chemotherapy, with either lysate (IgG, P = 0.0001, IgG1, P = 0.0001) or TcII/V/VI peptide antigen (IgG, P = 0.0001, IgG1, P = 0.0001). IgG1 decline was more discriminative than IgG. Antibody decline after treatment was also detected by the RDT. Incomplete treatment was associated with high IgG1 post-treatment titres against lysate (P = 0.013), as were IgG post-treatment titres to TcII/V/VI peptide (P = 0.0001). High pre-treatment IgG1 with lysate was associated with Qom ethnicity (P = 0.045). No associations were found between gender, age, body mass index and pre- or post-treatment antibody titres. Conclusions: We show that following chemotherapy of early chronic Chagas disease, significant decline in IgG1 antibody suggests cure, whereas sustained or increased IgG1 is a potential indicator of treatment failure. Due to restricted sensitivity, IgG1 should not be used as a diagnostic marker but has promise, with further development, as a biomarker of cure. Graphical abstract: We show that following chemotherapy of early chronic Chagas disease, a significant decline in IgG1 antibody suggests cure, whereas sustained or increased IgG1 is a potential indicator of treatment failure. Due to restricted sensitivity, IgG1 should not be used as a diagnostic marker but has promise, with further development, as a biomarker of cure.Fil: Murphy, Niamh. University of London; Reino UnidoFil: Cardinal, Marta Victoria. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Departamento de Ecología, Genética y Evolución. Laboratorio de Eco-Epidemiología; Argentina. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Ecología, Genética y Evolución de Buenos Aires. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Ecología, Genética y Evolución de Buenos Aires; ArgentinaFil: Bhattacharyya, Tapan. University of London; Reino UnidoFil: Enriquez, Gustavo Fabián. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Ecología, Genética y Evolución de Buenos Aires. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Ecología, Genética y Evolución de Buenos Aires; ArgentinaFil: Macchiaverna, Natalia Paula. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Ecología, Genética y Evolución de Buenos Aires. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Ecología, Genética y Evolución de Buenos Aires; ArgentinaFil: Alvedro, Alejandra. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Ecología, Genética y Evolución de Buenos Aires. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Ecología, Genética y Evolución de Buenos Aires; ArgentinaFil: Freilij, Héctor. Gobierno de la Ciudad de Buenos Aires. Hospital General de Niños "Ricardo Gutiérrez"; ArgentinaFil: Martinez de Salazar, Pablo. Barcelona Institute For Global Health; EspañaFil: Molina, Israel. Barcelona Institute For Global Health; EspañaFil: Mertens, Pascal. No especifíca;Fil: Gilleman, Quentin. No especifíca;Fil: Gurtler, Ricardo Esteban. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Ecología, Genética y Evolución de Buenos Aires. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Ecología, Genética y Evolución de Buenos Aires; ArgentinaFil: Miles, Michael A.. University of London; Reino Unid

Lineage-specific rapid diagnostic tests can resolve Trypanosoma cruzi TcII/V/VI ecological and epidemiological associations in the Argentine Chaco.

BACKGROUND: Trypanosoma cruzi, the protozoan agent of Chagas disease, is comprised of at least 6 genetic lineages (TcI-TcVI). Their geographical distribution, clinical associations and reservoir hosts are not fully elucidated, as genotyping is hampered due to the difficulty in isolating representative populations of organisms. Lineage-specific serological techniques may address these issues. METHODS: Trypanosoma cruzi lineage-specific serological assays were performed on human, canine, feline and armadillo sera from the Gran Chaco in northern Argentina, a region of ongoing transmission. Synthetic peptides representing lineage-specific epitopes of the trypomastigote small surface antigen (TSSA) were used in ELISA, and the TcII/V/VI shared epitope peptide (TSSApep-II/V/VI) was used in the Chagas Sero K-SeT rapid diagnostic test (RDT). RESULTS: Chagas Sero K-SeT RDT, using Protein G to detect human and canine IgG, was at least as sensitive as TSSApep-II/V/VI ELISA using specific secondary antibodies. For sera from humans TSSApep-II/V/VI seroprevalence by Chagas Sero K-SeT was 273/393 (69.5%), for dogs 48/73 (65.8%) and for armadillos 1/7 (14.3%); by ELISA for cats 5/19 (26.3%). The seroprevalence for humans was similar to that for Bolivian patients, amongst whom we previously observed an association of TSSApep-II/V/VI seropositivity with severity of cardiomyopathy. In humans, prevalence of TSSApep-II/V/VI recognition was associated with locality, and with increasing and decreasing age within the Qom and Creole populations, respectively. For dogs TSSApep-II/V/VI recognition was associated with being born before community-wide insecticide spraying (P = 0.05) and with Qom household (P < 0.001). CONCLUSIONS: We show here that Chagas Sero K-SeT RDT can replace ELISA for TSSApep-II/V/VI serology of humans and dogs; for humans there were statistically significant associations between a positive Chagas Sero K-SeT RDT and being resident in Area IV, and for dogs association with Qom household or with being born before the mass spraying campaign; we also show that with cats the TcII/V/VI epitope can be detected by ELISA. We assessed the lineage distribution in an unprecedented 83% of the human T. cruzi-seropositive population. These results form the basis for more detailed studies, enabling rapid in-the-field surveillance of the distribution and clustering of these lineages among humans and mammalian reservoirs of T. cruzi infection

A robust methodology to subclassify pseudokinases based on their nucleotide-binding properties

Protein kinase-like domains that lack conserved residues known to catalyse phosphoryl transfer, termed pseudokinases, have emerged as important signalling domains across all kingdoms of life. Although predicted to function principally as catalysis-independent protein-interaction modules, several pseudokinase domains have been attributed unexpected catalytic functions, often amid controversy. We established a thermal-shift assay as a benchmark technique to define the nucleotide-binding properties of kinase-like domains. Unlike in vitro kinase assays, this assay is insensitive to the presence of minor quantities of contaminating kinases that may otherwise lead to incorrect attribution of catalytic functions to pseudokinases. We demonstrated the utility of this method by classifying 31 diverse pseudokinase domains into four groups: devoid of detectable nucleotide or cation binding; cation-independent nucleotide binding; cation binding; and nucleotide binding enhanced by cations. Whereas nine pseudokinases bound ATP in a divalent cation-dependent manner, over half of those examined did not detectably bind nucleotides, illustrating that pseudokinase domains predominantly function as non-catalytic protein-interaction modules within signalling networks and that only a small subset is potentially catalytically active. We propose that henceforth the thermal-shift assay be adopted as the standard technique for establishing the nucleotide-binding and catalytic potential of kinase-like domains

Initial Sequence and Comparative Analysis of the Cat Genome

The genome sequence (1.9-fold coverage) of an inbred Abyssinian domestic cat was assembled, mapped, and annotated with a comparative approach that involved cross-reference to annotated genome assemblies of six mammals (human, chimpanzee, mouse, rat, dog, and cow). The results resolved chromosomal positions for 663,480 contigs, 20,285 putative feline gene orthologs, and 133,499 conserved sequence blocks (CSBs). Additional annotated features include repetitive elements, endogenous retroviral sequences, nuclear mitochondrial (numt) sequences, micro-RNAs, and evolutionary breakpoints that suggest historic balancing of translocation and inversion incidences in distinct mammalian lineages. Large numbers of single nucleotide polymorphisms (SNPs), deletion insertion polymorphisms (DIPs), and short tandem repeats (STRs), suitable for linkage or association studies were characterized in the context of long stretches of chromosome homozygosity. In spite of the light coverage capturing ∼65% of euchromatin sequence from the cat genome, these comparative insights shed new light on the tempo and mode of gene/genome evolution in mammals, promise several research applications for the cat, and also illustrate that a comparative approach using more deeply covered mammals provides an informative, preliminary annotation of a light (1.9-fold) coverage mammal genome sequence

Parent-of-origin-specific allelic associations among 106 genomic loci for age at menarche.

Age at menarche is a marker of timing of puberty in females. It varies widely between individuals, is a heritable trait and is associated with risks for obesity, type 2 diabetes, cardiovascular disease, breast cancer and all-cause mortality. Studies of rare human disorders of puberty and animal models point to a complex hypothalamic-pituitary-hormonal regulation, but the mechanisms that determine pubertal timing and underlie its links to disease risk remain unclear. Here, using genome-wide and custom-genotyping arrays in up to 182,416 women of European descent from 57 studies, we found robust evidence (P < 5 × 10(-8)) for 123 signals at 106 genomic loci associated with age at menarche. Many loci were associated with other pubertal traits in both sexes, and there was substantial overlap with genes implicated in body mass index and various diseases, including rare disorders of puberty. Menarche signals were enriched in imprinted regions, with three loci (DLK1-WDR25, MKRN3-MAGEL2 and KCNK9) demonstrating parent-of-origin-specific associations concordant with known parental expression patterns. Pathway analyses implicated nuclear hormone receptors, particularly retinoic acid and γ-aminobutyric acid-B2 receptor signalling, among novel mechanisms that regulate pubertal timing in humans. Our findings suggest a genetic architecture involving at least hundreds of common variants in the coordinated timing of the pubertal transition

Development and validation of DNA Methylation scores in two European cohorts augment 10-year risk prediction of type 2 diabetes

This is the author accepted manuscriptAvailability of Data and Material: According to the terms of consent for Generation Scotland participants, access to data must be reviewed by the Generation Scotland Access Committee. Applications should be made to [email protected]. All code is available with open access at the following Gitlab repository: https://github.com/marioni-group MethylPipeR (version 1.0.0) is available at: https://github.com/marioni-group/MethylPipeR MethylPipeR-UI is available at: https://github.com/marioni-group/MethylPipeR-UI. The informed consents given by KORA study participants do not cover data posting in public databases. However, data are available upon request from KORA Project Application Self Service Tool (https://epi.helmholtz-muenchen.de/). Data requests can be submitted online and are subject to approval by the KORA Board.Type 2 diabetes mellitus (T2D) presents a major health and economic burden that could be alleviated with improved early prediction and intervention. While standard risk factors have shown good predictive performance, we show that the use of blood-based DNA methylation information leads to a significant improvement in the prediction of 10-year T2D incidence risk. Previous studies have been largely constrained by linear assumptions, the use of CpGs one-at43 a-time, and binary outcomes. We present a flexible approach (via an R package, MethylPipeR) based on a range of linear and tree-ensemble models that incorporate time-to-event data for prediction. Using the Generation Scotland cohort (training set ncases=374, ncontrols=9,461; test set ncases=252, ncontrols=4,526) our best-performing model (Area Under the Curve (AUC)=0.872, Precision Recall AUC (PRAUC)=0.302) showed notable improvement in 10-year onset prediction beyond standard risk factors (AUC=0.839, PRAUC=0.227). Replication was observed in the German-based KORA study (n=1,451, ncases = 142, p=1.6x10-5 49 ).Wellcome TrustChief Scientist Office of the Scottish Government Health DirectoratesScottish Funding Counci

Social factors influencing Russian male alcohol use over the life course: a qualitative study investigating age based social norms, masculinity, and workplace context

The massive fluctuations occurring in Russian alcohol-related mortality since the mid-1980s cannot be seen outside of the context of great social and economic change. There is a dearth of qualitative studies about Russian male drinking and especially needed are those that address social processes and individual changes in drinking. Conducted as part of a longitudinal study on men’s alcohol consumption in Izhevsk, this qualitative study uses 25 semi-structured biographical interviews with men aged 33–60 years to explore life course variation in drinking. The dominant pattern was decreasing binge and frequent drinking as men reached middle age which was precipitated by family building, reductions in drinking with work colleagues, and health concerns. A minority of men described chaotic drinking histories with periods of abstinence and heavy drinking. The results highlight the importance of the blue-collar work environment for conditioning male heavy drinking in young adulthood through a variety of social, normative and structural mechanisms. Post-Soviet changes had a structural influence on the propensity for workplace drinking but the important social function of male drinking sessions remained. Bonding with workmates through heavy drinking was seen as an unavoidable and essential part of young men’s social life. With age peer pressure to drink decreased and the need to perform the role of responsible breadwinner put different behavioural demands on men. For some resisting social pressure to drink became an important site of self-determination and a mark of masculine maturity. Over the lifetime the place where masculine identity was asserted shifted from the workplace to the home, which commonly resulted in a reduction in drinking. We contribute to existing theories of Russian male drinking by showing that the performance of age-related social roles influences Russian men’s drinking patterns, drinking contexts and their attitudes. Further research should be conducted investigating drinking trajectories in Russian men

Where do stars explode in the ISM? -- The distribution of dense gas around massive stars and supernova remnants in M33

Star formation in galaxies is regulated by turbulence, outflows, gas heating and cloud dispersal -- processes which depend sensitively on the properties of the interstellar medium (ISM) into which supernovae (SNe) explode. Unfortunately, direct measurements of ISM environments around SNe remain scarce, as SNe are rare and often distant. Here we demonstrate a new approach: mapping the ISM around the massive stars that are soon to explode. This provides a much larger census of explosion sites than possible with only SNe, and allows comparison with sensitive, high-resolution maps of the atomic and molecular gas from the Jansky VLA and ALMA. In the well-resolved Local Group spiral M33, we specifically observe the environments of red supergiants (RSGs, progenitors of Type II SNe), Wolf-Rayet stars (WRs, tracing stars $>$30 M$_{\odot}$, and possibly future stripped-envelope SNe), and supernova remnants (SNRs, locations where SNe have exploded). We find that massive stars evolve not only in dense, molecular-dominated gas (with younger stars in denser gas), but also a substantial fraction ($\sim$45\% of WRs; higher for RSGs) evolve in lower-density, atomic-gas-dominated, inter-cloud media. We show that these measurements are consistent with expectations from different stellar-age tracer maps, and can be useful for validating SN feedback models in numerical simulations of galaxies. Along with the discovery of a 20-pc diameter molecular gas cavity around a WR, these findings re-emphasize the importance of pre-SN/correlated-SN feedback evacuating the dense gas around massive stars before explosion, and the need for high-resolution (down to pc-scale) surveys of the multi-phase ISM in nearby galaxies.Comment: 34 pages, 14 figures. Submitted to ApJ. Comments welcome! The density distributions will be made publicly available after journal acceptance of manuscript. Please feel free to contact us in the meantime if you would like to use the

Brain Circuitries Involved in Semantic Interference by Demands of Emotional and Non-Emotional Distractors

BACKGROUND: Previous studies have indicated that the processes leading to the resolution of emotional and non-emotional interference conflicts are unrelated, involving separate networks. It is also known that conflict resolution itself suggests a considerable overlap of the networks. Our study is an attempt to examine how these findings may be related. METHODOLOGY/PRINCIPAL FINDINGS: We used functional magnetic resonance imaging (fMRI) to study neural responses of 24 healthy subjects to emotional and non-emotional conflict paradigms involving the presentation of congruent and incongruent word-face pairs based on semantic incompatibility between targets and distractors. In the emotional task, the behavioral interference conflict was greater (compared to the non-emotional task) and was paralleled by involvement of the extrastriate visual and posterodorsal medial frontal cortices. In both tasks, we also observed a common network including the dorsal anterior cingulate, the supplemental motor area, the anterior insula and the inferior prefrontal cortex, indicating that these brain structures are markers of experienced conflict. However, the emotional task involved conflict-triggered networks to a considerably higher degree. CONCLUSIONS/SIGNIFICANCE: Our findings indicate that responses to emotional and non-emotional distractors involve the same systems, which are capable of flexible adjustments based on conflict demands. The function of systems related to conflict resolution is likely to be adjusted on the basis of an evaluation process that primarily involves the extrastriate visual cortex, with target playing a significant role