γ-aminobutyric acid measurement in the human brain at 7 T: Short echo-time or Mescher-Garwood editing

The purposes of the current study were to introduce a Mescher-Garwood (MEGA) semi-adiabatic spin-echo full-intensity localization (MEGA-sSPECIAL) sequence with macromolecule (MM) subtraction and to compare the test-retest reproducibility of γ-aminobutyric acid (GABA) measurements at 7 T using the sSPECIAL and MEGA-sSPECIAL sequences. The MEGA-sSPECIAL editing scheme using asymmetric adiabatic and highly selective Gaussian pulses was used to compare its GABA measurement reproducibility with that of short echo-time (TE) sSPECIAL. Proton magnetic resonance spectra were acquired in the motor cortex (M1) and medial prefrontal cortex (mPFC) using the sSPECIAL (TR/TE = 4000/16 ms) and MEGA-sSPECIAL sequences (TR/TE = 4000/80 ms). The metabolites were quantified using LCModel with unsuppressed water spectra. The concentrations are reported in institutional units. The test-retest reproducibility was evaluated by scanning each subject twice. Between-session reproducibility was assessed using coefficients of variation (CVs), Pearson's r correlation coefficients, and intraclass correlation coefficients (ICCs). Intersequence agreement was evaluated using Pearson's r correlation coefficients and Bland-Altman plots. Regarding GABA measurements by sSPECIAL, the GABA concentrations were 0.92 ± 0.31 (IU) in the M1 and 1.56 ± 0.49 (IU) in the mPFC. This demonstrated strong between-session correlation across both regions (r = 0.81, p < 0.01; ICC = 0.82). The CVs between the two scans were 21.8% in the M1 and 10.2% in the mPFC. On the other hand, the GABA measurements by MEGA-sSPECIAL were 0.52 ± 0.04 (IU) in the M1 and 1.04 ± 0.24 (IU) in the mPFC. MEGA-sSPECIAL demonstrated strong between-session correlation across the two regions (r = 0.98, p < 0.001; ICC = 0.98) and lower CVs than sSPECIAL, providing 4.1% in the M1 and 5.8% in the mPFC. The MEGA-editing method showed better reproducibility of GABA measurements in both brain regions compared with the short-TE sSPECIAL method. Thus it is a more sensitive method with which to detect small changes in areas with low GABA concentrations. In GABA-rich brain regions, GABA measurements can be achieved reproducibly using both methods.

Quantifying the excitatory-inhibitory balance: A comparison of SemiLASER and MEGA-SemiLASER for simultaneously measuring GABA and glutamate at 7T

The importance of the excitatory-inhibitory (E/I) balance in a wide range of cognitive and behavioral processes has prompted a commensurate interest in methods for reliably quantifying it. Proton Magnetic Resonance Spectroscopy (¹H-MRS) remains the only method capable of safely and non-invasively measuring the concentrations of the brain's major excitatory (glutamate) and inhibitory (γ-aminobutyric-acid, GABA) neurotransmitters in-vivo. MRS relies on spectral Mescher-Garwood (MEGA) editing techniques at 3T to distinguish GABA from its overlapping resonances. However, with the increased spectral resolution at ultrahigh field strengths of 7T and above, non-edited spectroscopic techniques become potential viable alternatives to MEGA based approaches, and also address some of their shortcomings, such as signal loss, sensitivity to transmitter inhomogeneities and temporal resolution. We present a comprehensive comparison of both edited and non-edited strategies at 7T for simultaneously quantifying glutamate and GABA from the dorsal anterior cingulate cortex (dACC), and evaluate their reproducibility and relative bias. The combined root-mean-square test-retest reproducibility of Glu and GABA (CV_E/I) was as low as 13.3% for unedited MRS at TE=80 ms using SemiLASER localization, while edited MRS at TE=80 ms yielded CV_E/I=20% and 21% for asymmetric and symmetric MEGA editing, respectively. An unedited SemiLASER acquisition using a shorter echo time of TE=42 ms yielded CV_E/I as low as 24.9%. Our results show that non-edited sequences at an echo time of 80 ms provide better reproducibility than either edited sequences at the same TE, or non-edited sequences at a shorter TE of 42 ms. This is supported by numerical simulations and is driven in part by a pseudo-singlet appearance of the GABA multiplets at TE=80 ms, and the excellent spectral resolution at 7T. Our results uphold a transition to non-edited MRS for monitoring the E/I balance at ultrahigh fields, and stress the importance of using a properly-optimized echo time.

How to Build a Super Predator: From Genotype to Phenotype

We present a drawing discovery lab that crosscuts multiple disciplines in biology and links concepts in genetics and evolutionary thinking to enhance understanding of the genotype-to-phenotype transformation. These combined concepts are also linked to ecological frameworks in nature through the model of biological plasticity. Students and teachers explore drawing skills to flesh out the future of a predator while engaging with the computational software MEGA, which introduces students and teachers to nucleotide changes, mutations, variation, phylogenetics, and molecular evolution.

Flexible MEGA editing scheme with asymmetric adiabatic pulses applied for T 2 measurement of lactate in human brain

PURPOSE

A flexible MEGA editing scheme which decouples the editing efficiency from TE is proposed and the utility of asymmetric adiabatic pulses for this new technique is explored. It is demonstrated that the method enables robust T 2 measurement of lactate in healthy human brain.

METHODS

The proposed variation of the MEGA scheme applies editing pulses in both acquired spectra, ensuring that the difference in J-evolution of the target resonance leads to maximal signal yield in the difference spectrum for arbitrary TE. A MEGA-sLASER sequence is augmented with asymmetric adiabatic editing pulses for enhanced flexibility and immunity to B 1 + miscalibration and inhomogeneities. The technique is validated and optimized for flexible lactate editing via a simple analytical model, numerical simulations and in vitro experiments. The T 2 relaxation constant of lactate is determined in vivo via multiple-TE measurements with the proposed method and a dedicated postprocessing and quantification approach.

RESULTS

Asymmetric adiabatic editing pulses improve robustness and facilitate efficient J-editing in sequences or protocols with strong timing constraints. Single voxel measurements using the proposed MEGA scheme in the occipital cortex of six healthy subjects yield a relaxation constant of T 2 = 171 ± 19  ms for the methyl resonance of lactate at a field strength of 3T.

CONCLUSIONS

The proposed MEGA editing scheme allows for novel kinds of J-editing experiments and promises to be an asset to robust T 2 measurement of lactate and potentially other J-coupled metabolites in vivo.

Discovery of rare coding variants in OGDHL and BRCA2 in relation to breast cancer risk in Chinese women

The missing heritability of breast cancer could be partially attributed to rare variants (MAF < 0.5%). To identify breast cancer-associated rare coding variants, we conducted whole-exome sequencing (~50×) in genomic DNA samples obtained from 831 breast cancer cases and 839 controls of Chinese females. Using burden tests for each gene that included rare missense or predicted deleterious variants, we identified 29 genes showing promising associations with breast cancer risk. We replicated the association for two genes, OGDHL and BRCA2, at a Bonferroni-corrected p < 0.05, by genotyping an independent set of samples from 1,628 breast cancer cases and 1,943 controls. The association for OGDHL was primarily driven by three predicted deleterious variants (p.Val827Met, p.Pro839Leu, p.Phe836Ser; p < 0.01 for all). For BRCA2, we characterized a total of 27 disruptive variants, including 18 nonsense, six frameshift and three splicing variants, whereas they were only detected in cases, but none of the controls. All of these variants were either very rare (AF < 0.1%) or not detected in >4,500 East Asian women from the genome Aggregation database (gnomAD), providing additional support to our findings. Our study revealed a potential novel gene and multiple disruptive variants of BRCA2 for breast cancer risk, which may identify high-risk women in Chinese populations.

VCF2PopTree: a client-side software to construct population phylogeny from genome-wide SNPs

In the past decades a number of software programs have been developed to infer phylogenetic relationships between populations. However, most of these programs typically use alignments of sequences from genes to build phylogeny. Recently, many standalone or web applications have been developed to handle large-scale whole genome data, but they are either computationally intensive, dependent on third party software or required significant time and resource of a web server. In the post-genomic era, researchers are able to obtain bioinformatically processed high-quality publication-ready whole genome data for many individuals in a population from next generation sequencing companies due to the reduction in the cost of sequencing and analysis. Such genotype data is typically presented in the Variant Call Format (VCF) and there is no simple software available that directly uses this data format to construct the phylogeny of populations in a short time. To address this limitation, we have developed a user-friendly software, VCF2PopTree that uses genome-wide SNPs to construct and display phylogenetic trees in seconds to minutes. For example, it reads a VCF file containing 4 million SNPs and draws a tree in less than 30 seconds. VCF2PopTree accepts genotype data from a local machine, constructs a tree using UPGMA and Neighbour-Joining algorithms and displays it on a web-browser. It also produces pairwise-diversity matrix in MEGA and PHYLIP file formats as well as trees in the Newick format which could be directly used by other popular phylogenetic software programs. The software including the source code, a test VCF file and a documentation are available at: https://github.com/sansubs/vcf2pop.

In Silico Analysis of Hepatitis B Virus Genotype D Subgenotype D1 Circulating in Pakistan, China, and India

The focus of this study was the computational analysis of hepatitis B virus (HBV) genotype D subgenotype D1 in Pakistan, China, and India. In total, 54 complete genome sequences of HBV genotype D subgenotype D1 were downloaded from National Center for Biotechnology Information (NCBI). Of these, 6 complete genome sequences were from Pakistan, 14 were from China, and 34 were from India. Sequence alignment showed less than 4% divergence in these sequences. C and X genes showed divergence of less than 3%. Comparison over the S gene showed more than 97% similarity among the nucleotide sequences of genotype D subgenotype D1. The identity and similarity matrix of 54 nucleotide sequences of HBV genotype D subgenotype D1 from Pakistan, China, and India revealed more than 93% identity and 93% similarity. Phylogenetic analysis highlighted that complete genome isolates of HBV circulating in Pakistan had the closest evolutionary relationship with its neighboring countries China and India. China’s (HQ833466) and Pakistan’s (AB583680.1) isolates shared the same ancestor. Gene structure analysis showed that “P” gene exons were the longest, about three-fourth of the genome size, whereas gene “S” had the second longest coding regions with 2 exons and 1 intron. However, “C” and “X” genes had 1 smallest exon. X proteins had proven role in spreading of the HBV infection diseases. For HBx analysis, 1 X protein sequence of HBV genotype D subgenotype D1 belonging to each country was obtained. Homology models of the 3 X proteins generated using SWISS-MODEL revealed GMQE (Global Model Quality Estimation) = 0.1. Global and local quality estimate scores including Z-scores for Qualitative Model Energy Analysis (QMEAN) C-beta, all-atom, solvation, and torsion energy scores were similar indicating good quality, accuracy, and reliability of the predicted models. Three-dimensional (3D) visualization showed similar structures and Ramachandran plots showed a high percentage of protein residues into the favorable region for X protein models.

Scientific Evolution of Clinical and Risk Assessment of Sexually Abusive Youth: A Comprehensive Review of Empirical Tools

Risk assessment of sexually abusive youth is a specialty, which ideally includes a clinical assessment (e.g., psychological evaluation) utilizing a validated risk assessment tool. This article reviews tools for assessing sexual improprieties and/or sexually abusive behaviors in youth: JSORRAT-II, J-SOAP-II, J-RAS, AIM2, and MEGA^♪,1 and a clinical assessment tool, MIDSA. Untested, structured, clinical checklists with face validity, J-RAT-4 and PROFESOR, are included. The authors also highlight clinical dilemmas using "utility tools" (i.e., polygraph and plethysmograph) with sexually abusive youth. The comprehensive review sustains a twenty-first-century New Paradigm of inclusive ecologically based, developmentally and gender-sensitive assessment tools that definitively and accurately assess risk and protective factors of sexually abusive youth.

Genome-wide analysis and identification of stress-responsive genes of the NAM-ATAF1,2-CUC2 transcription factor family in apple

NAC (NAM, ATAF1,2, and CUC2) proteins constitute one of the largest families of plant-specific transcription factors. To date, little is known about the NAC genes in the apple (Malus domestica). In this study, a total of 180 NAC genes were identified in the apple genome and were phylogenetically clustered into six groups (I-VI) with the NAC genes from Arabidopsis and rice. The predicted apple NAC genes were distributed across all of 17 chromosomes at various densities. Additionally, the gene structure and motif compositions of the apple NAC genes were analyzed. Moreover, the expression of 29 selected apple NAC genes was analyzed in different tissues and under different abiotic stress conditions. All of the selected genes, with the exception of four genes, were expressed in at least one of the tissues tested, which indicates that the NAC genes are involved in various aspects of the physiological and developmental processes of the apple. Encouragingly, 17 of the selected genes were found to respond to one or more of the abiotic stress treatments, and these 17 genes included not only the expected 7 genes that were clustered with the well-known stress-related marker genes in group IV but also 10 genes located in other subgroups, none of which contains members that have been reported to be stress-related. To the best of our knowledge, this report describes the first genome-wide analysis of the apple NAC gene family, and the results should provide valuable information for understanding the classification and putative functions of this family.

A functional (E)-4-hydroxy-3-methylbut-2-enyl diphosphate reductase exhibits diurnal regulation of expression in Stevia rebaudiana (Bertoni)

The leaves of stevia [Stevia rebaudiana (Bertoni)] are a rich source of steviol glycosides that are used as non-calorific sweetener in many countries around the world. Steviol moiety of steviol glycosides is synthesized via plastidial 2C-methyl-D-erythritol 4-phosphate pathway, where (E)-4-hydroxy-3-methylbut-2-enyl diphosphate reductase (HDR) is the key enzyme. HDR catalyzes the simultaneous conversion of (E)-4-hydroxy-3-methylbut-2-enyl diphosphate into five carbon isoprenoid units, isopentenyl diphosphate and dimethylallyl diphosphate. Stevia HDR (SrHDR) successfully rescued HDR lethal mutant strain MG1655 ara<>ispH upon genetic complementation, suggesting SrHDR to encode a functional protein. The gene exhibited diurnal variation in expression. To identify the possible regulatory elements, upstream region of the gene was cloned and putative cis-acting elements were detected by in silico analysis. Electrophoretic mobility shift assay, using a putative light responsive element GATA showed the binding of nuclear proteins (NP) isolated from leaves during light period of the day, but not with the NP from leaves during the dark period. Data suggested the involvement of GATA box in light mediated gene regulation of SrHDR in stevia.

Abundant sequence divergence in the native Japanese cattle Mishima-Ushi (Bos taurus) detected using whole-genome sequencing

The native Japanese cattle Mishima-Ushi, a designated national natural treasure, are bred on a remote island, which has resulted in the conservation of their genealogy. We examined the genetic characteristics of 8 Mishima-Ushi individuals by using single nucleotide polymorphisms (SNPs), insertions, and deletions obtained by whole-genome sequencing. Mapping analysis with various criteria showed that predicted heterozygous SNPs were more prevalent than predicted homozygous SNPs in the exonic region, especially non-synonymous SNPs. From the identified 6.54 million polymorphisms, we found 400 non-synonymous SNPs in 313 genes specific to each of the 8 Mishima-Ushi individuals. Additionally, 3,170,833 polymorphisms were found between the 8 Mishima-Ushi individuals. Phylogenetic analysis confirmed that the Mishima-Ushi population diverged from another strain of Japanese cattle. This study provides a framework for further genetic studies of Mishima-Ushi and research on the function of SNP-containing genes as well as understanding the genetic relationship between the domestic and native Japanese cattle breeds.

Genome-wide identification and expression analysis of MAPK and MAPKK gene family in Malus domestica

MAPK signal transduction modules play crucial roles in regulating many biological processes in plants, which are composed of three classes of hierarchically organized protein kinases, namely MAPKKKs, MAPKKs, and MAPKs. Although genome-wide analysis of this family has been carried out in some species, little is known about MAPK and MAPKK genes in apple (Malus domestica). In this study, a total of 26 putative apple MAPK genes (MdMPKs) and 9 putative apple MAPKK genes (MdMKKs) have been identified and located within the apple genome. Phylogenetic analysis revealed that MdMAPKs and MdMAPKKs could be divided into 4 subfamilies (groups A, B, C and D), respectively. The predicted MdMAPKs and MdMAPKKs were distributed across 13 out of 17 chromosomes with different densities. In addition, analysis of exon-intron junctions and of intron phase inside the predicted coding region of each candidate gene has revealed high levels of conservation within and between phylogenetic groups. According to the microarray and expressed sequence tag (EST) analysis, the different expression patterns indicate that they may play different roles during fruit development and rootstock-scion interaction process. Moreover, MAPK and MAPKK genes were performed expression profile analyses in different tissues (root, stem, leaf, flower and fruit), and all of the selected genes were expressed in at least one of the tissues tested, indicating that the MAPKs and MAPKKs are involved in various aspects of physiological and developmental processes of apple. To our knowledge, this is the first report of a genome-wide analysis of the apple MAPK and MAPKK gene family. This study provides valuable information for understanding the classification and putative functions of the MAPK signal in apple.

Modeling, docking and dynamics simulations of a non-specific lipid transfer protein from Peganum harmala L

Non-specific lipid transfer proteins (ns-LTPs), ubiquitously found in various types of plants, have been well-known to transfer amphiphilic lipids and promote the lipid exchange between mitochondria and microbody. In this study, an in silico analysis was proposed to study ns-LTP in Peganum harmala L., which may belong to ns-LTP1 family, aiming at constructing its three-dimensional structure. Moreover, we adopted MEGA to analyze ns-LTPs and other species phylogenetically, which brought out an initial sequence alignment of ns-LTPs. In addition, we used molecular docking and molecular dynamics simulations to further investigate the affinities and stabilities of ns-LTP with several ligands complexes. Taken together, our results about ns-LTPs and their ligand-binding activities can provide a better understanding of the lipid-protein interactions, indicating some future applications of ns-LTP-mediated transport.

Novel polymorphisms in UTR and coding region of inducible heat shock protein 70.1 gene in tropically adapted Indian zebu cattle (Bos indicus) and riverine buffalo (Bubalus bubalis)

Due to evolutionary divergence, cattle (taurine, and indicine) and buffalo are speculated to have different responses to heat stress condition. Variation in candidate genes associated with a heat-shock response may provide an insight into the dissimilarity and suggest targets for intervention. The present work was undertaken to characterize one of the inducible heat shock protein genes promoter and coding regions in diverse breeds of Indian zebu cattle and buffaloes. The genomic DNA from a panel of 117 unrelated animals representing 14 diversified native cattle breeds and 6 buffalo breeds were utilized to determine the complete sequence and gene diversity of HSP70.1 gene. The coding region of HSP70.1 gene in Indian zebu cattle, Bos taurus and buffalo was similar in length (1,926 bp) encoding a HSP70 protein of 641 amino acids with a calculated molecular weight (Mw) of 70.26 kDa. However buffalo had a longer 5' and 3' untranslated region (UTR) of 204 and 293 nucleotides respectively, in comparison to Indian zebu cattle and Bos taurus wherein length of 5' and 3'-UTR was 172 and 286 nucleotides, respectively. The increased length of buffalo HSP70.1 gene compared to indicine and taurine gene was due to two insertions each in 5' and 3'-UTR. Comparative sequence analysis of cattle (taurine and indicine) and buffalo HSP70.1 gene revealed a total of 54 gene variations (50 SNPs and 4 INDELs) among the three species in the HSP70.1 gene. The minor allele frequencies of these nucleotide variations varied from 0.03 to 0.5 with an average of 0.26. Among the 14 B. indicus cattle breeds studied, a total of 19 polymorphic sites were identified: 4 in the 5'-UTR and 15 in the coding region (of these 2 were non-synonymous). Analysis among buffalo breeds revealed 15 SNPs throughout the gene: 6 at the 5' flanking region and 9 in the coding region. In bubaline 5'-UTR, 2 additional putative transcription factor binding sites (Elk-1 and C-Re1) were identified, other than three common sites (CP2, HSE and Pax-4) observed across all the analyzed animals. No polymorphism was found within the 3'-UTR of Indian cattle or buffalo as it was found to be monomorphic. The promoter sequences generated in 117 individuals showed a rich array of sequence elements known to be involved in transcription regulation. A total of 11 nucleotide changes were observed in the promoter sequence across the analyzed species, 3 of these changes were located within the potential transcription factor binding domains. We also identified 4 microsatellite markers within the buffalo HSP70.1 gene and 3 microsatellites within bovine HSP70.1. The present study identified several distinct changes across indicine, taurine and bubaline HSP70.1 genes that could further be evaluated as molecular markers for thermotolerance.

Phylogenetic analysis of chloroplast matK gene from Zingiberaceae for plant DNA barcoding

MaturaseK gene (MatK) of chloroplast is highly conserved in plant systematics which is involved in Group II intron splicing. The size of the gene is 1500 bp in length, located with in the intron of trnK. In the present study, matK gene from Zingiberaceae was taken for the analysis of variants, parsimony site, patterns, transition/tranversion rates and phylogeny. The family of Zingiberaceae comprises 47 genera with medicinal values. The matK gene sequence have been obtained from genbank and used for the analysis. The sequence alignments were performed by Clustal X, transition/transversion rates were predicted by MEGA and phylogenetic analyses were carried out by PHYLIP package. The result indicates that the Zingiberaceae genus Afromonum, Alpinia, Globba, Curcuma and Zingiber shows polyphylogeny. The overall variants between the species are 24% and transition/transversion rate is 1.54. Phylogenetic tree was designed to identify the ideal regions that could be used for defining the inter and intera-generic relationships. From this study it could be concluded that the matK gene is a good candidate for DNA barcoding of plant family Zingiberaceae.