An improved heuristic for haplotype inference

We cloned a gene, kexD, that provides a multidrug-resistant phenotype from multidrug-resistant Klebsiella pneumoniae MGH78578. The deduced amino acid sequence of KexD is similar to that of the inner membrane protein, RND-type multidrug efflux pump. Introduction of the kexD gene into Escherichia coli KAM32 resulted in a MIC that was higher for erythromycin, novobiocin, rhodamine 6G, tetraphenylphosphonium chloride, and ethidium bromide than that of the control. Intracellular ethidium bromide levels in E. coli cells carrying the kexD gene were lower than that in the control cells under energized conditions, suggesting that KexD is a component of an energy-dependent efflux pump. RND-type pumps typically consist of three components: an inner membrane protein, a periplasmic protein, and an outer membrane protein. We discovered that KexD functions with a periplasmic protein, AcrA, from E. coli and K. pneumoniae, but not with the periplasmic proteins KexA and KexG from K. pneumoniae. KexD was able to utilize either TolC of E. coli or KocC of K. pneumoniae as an outer membrane component. kexD mRNA was not detected in K. pneumoniae MGH78578 or ATCC10031. We isolated erythromycin-resistant mutants from K. pneumoniae ATCC10031, and some showed a multidrug-resistant phenotype similar to the drug resistance pattern of KexD. Two strains of multidrug-resistant mutants were investigated for kexD expression; kexD mRNA levels were increased in these strains. We conclude that changing kexD expression can contribute to the occurrence of multidrug-resistant K. pneumoniae.     

Ancestral modal Y-STR haplotype shared among Romani and South Indian populations

One of the primary unanswered questions regarding the dispersal of Romani populations concerns the geographical region and/or the Indian caste/tribe that gave rise to the proto-Romani group. To shed light on this matter, 161 Y-chromosomes from Roma, residing in two different provinces of Serbia, were analyzed. Our results indicate that the paternal gene pool of both groups is shaped by several strata, the most prominent of which, H1-M52, comprises almost half of each collection's patrilineages. The high frequency of M52 chromosomes in the two Roma populations examined may suggest that they descend from a single founder that has its origins in the Indian subcontinent. Moreover, when the Y-STR profiles of haplogroup H derived individuals in our Roma populations were compared to those typed in the South Indian emigrants from Malaysia and groups from Madras, Karnataka (Lingayat and Vokkaliga castes) and tribal Soligas, sharing of the two most common haplotypes was observed. These similarities suggest that South India may have been one of the contributors to the proto-Romanis. European genetic signatures (i.e., haplogroups E1b1b1a1b-V13, G2a-P15, I-M258, J2-M172 and R1-M173), on the other hand, were also detected in both groups, but at varying frequencies. The divergent European genetic signals in each collection are likely the result of differential gene flow and/or admixture with the European host populations but may also be attributed to dissimilar endogamous practices following the initial founder effect. Our data also support the notion that a number of haplogroups including G2a-P15, J2a3b-M67(xM92), I-M258 and E1b1b1-M35 were incorporated into the proto-Romani paternal lineages as migrants moved from northern India through Southwestern Asia, the Middle East and/or Anatolia into the Balkans.     

Hap-seqX: Expedite algorithm for haplotype phasing with imputation using sequence data

Haplotype phasing is one of the most important problems in population genetics as haplotypes can be used to estimate the relatedness of individuals and to impute genotype information which is a commonly performed analysis when searching for variants involved in disease. The problem of haplotype phasing has been well studied. Methodologies for haplotype inference from sequencing data either combine a set of reference haplotypes and collected genotypes using a Hidden Markov Model or assemble haplotypes by overlapping sequencing reads. A recent algorithm Hap-seq considers using both sequencing data and reference haplotypes and it is a hybrid of a dynamic programming algorithm and a Hidden Markov Model (HMM), which is shown to be optimal. However, the algorithm requires extremely large amount of memory which is not practical for whole genome datasets. The current algorithm requires saving intermediate results to disk and reads these results back when needed, which significantly affects the practicality of the algorithm. In this work, we proposed the expedited version of the algorithm Hap-seqX, which addressed the memory issue by using a posterior probability to select the records that should be saved in memory. We show that Hap-seqX can save all the intermediate results in memory and improves the execution time of the algorithm dramatically. Utilizing the strategy, Hap-seqX is able to predict haplotypes from whole genome sequencing data.     

Association of Interleukin-13 SNP rs20541 (Arg > Gln) to allergic rhinitis in an Asian population of ethnic Chinese in Singapore

In this study, we report our findings as we evaluated the association of IL-13 SNP rs20541 to AR as reported in the meta-analysis (Ying, X-J et al., April, Gene 2013) in an Asian population of 1322 ethnic Chinese recruited in Singapore well defined and characterized for allergic rhinitis. The SNP was significantly associated with AR at an odds ratio of 1.57 for the homozygous genotype comparisons indicating a strong risk for AR. This highlights the importance of this mutation in Asian populations for risk towards allergic phenotypes and also warrants further functional work to characterize the mechanism behind this genetic predisposition.     

An improved method for estimating sequence divergence of DNA using restriction endonuclease mappings

Summary The method proposed by Kaplan and Langley for estimating the extent of sequence divergence between related DNA's using restriction endonuclease maps is modified so that the estimates are easier to compute. In the two-species case, these modifications lead via a maximum likelihood approach to an estimate which is closely related to one recently suggested by Nei and Li (1979) and Gotoh et al. (1979). Simulation studies show that the modified estimates are comparable to those of Kaplan and Langley, providing that there is sufficient homology in the DNA segments of the related species. The M-species case, M 3, is also discussed.     

Influence of functional genetic polymorphism (-590C/T) in non-small cell lung cancer (NSCLC) development: the paradoxal role of IL-4

Lung cancer is the leading cause of death by cancer in the world, originating about 17.5% of total deaths from cancer (1.18 million). Inflammation plays an important role in the pathogenesis of lung cancer. IL-4 is an anti-inflammatory cytokine, which reduces the production of proinflammatory cytokines by monocytes and with direct antiproliferative effects in some tumors. The polymorphism -590C/T SNP is a C to T transition in the -590 position of the promoter region of the IL-4 gene. The aim of this study was to evaluate the influence of this polymorphism in the susceptibility to NSCLC. DNA was extracted from peripheral blood of 1060 individuals (391 patients diagnosed with NSCLC and a control group of 669 individuals without cancer). The characterization of IL-4 -590C/T genotypes was performed by PCR-RFLP (BsmFI). The -590C/T polymorphism genotypes were classified as low (CC) and high expression (TT). The frequencies obtained for the CC and TT genotypes were 90.1% and 9.9%, respectively, in the control group and 92.9% and 7.1%, respectively, in the case group. The analysis of the TT and CC genotype frequencies in the two groups showed a statistically significant difference in its distribution, indicating a protection of 80% for the development of NSCLC, type epidermoid in individuals with the TT genotype when compared with individuals with CC genotype (P=0.024, OR=0.221: 95% CI=0.053-0.928). We present for the first time that increased expression of IL-4 associated with the TT genotype may contribute to immune surveillance during NSCLC development.     

Screening of indigenous goats for prolificacy associated DNA markers of sheep

The present study was undertaken to explore the genetic basis of caprine prolificacy and to screen indigenous goats for prolificacy associated markers of sheep in BMPR1B, GDF9 and BMP15 genes. To detect the associated mutations and identify novel allelic variants in the candidate genes, representative samples were collected from the breeding tract of indigenous goat breeds varying in prolificacy and geographic distribution. DNA was extracted and PCR amplification was done using primers designed or available in literature for the coding DNA sequence of candidate genes. Direct sequencing was done to identify the genetic variations. Mutations in the candidate genes associated with fecundity in sheep were not detected in Indian goats. Three non-synonymous SNPs (C818T, A959C and G1189A) were identified in exon 2 of GDF9 gene out of which mutation A959C has been associated with prolificacy in exotic goats. Two novel SNPs (G735A and C808G) were observed in exon 2 of BMP15 gene.     

Analysis of the inheritance pattern of a Chinese family with phaeochromocytomas through whole exome sequencing

Phaeochromocytomas (PCCs) and paragangliomas (PGLs) are rare, catecholamine-producing tumors. Most familial PCC/PGLs have been detected to be autosomal dominantly inherited. However, this study was undertaken in a family with PCCs to determine candidate genes in a dominant or recessive inheritance pattern. After excluding mutations in ten PCC/PGL susceptibility genes by Sanger sequencing, we used whole exome sequencing for screening on the four family members to discover novel candidate genes associated with PCCs. Based on the inexistence of non-synonymous mutations or indels in the ten known genes and the structure of this pedigree, 3 damaging loci with dominant inheritance pattern, and 5 damaging loci with recessive homozygous inheritance pattern and 6 damaging genes with compound heterozygous inheritance pattern were narrowed down to indicate the association with PCCs. According to the Gene Ontology (GO) category analysis on the combined results, cell adhesion showed the most significant enrichment.     

Variable phenotype in 17q12 microdeletions: Clinical and molecular characterization of a new case

Microdeletions of 17q12 including the hepatocyte nuclear factor 1 beta (HNF1B) gene, as well as point mutations of this gene, are associated with the Renal Cysts and Diabetes syndrome (RCAD, OMIM 137920) and genitourinary alterations. Also, microdeletions encompassing HNF1B were identified as a cause of Mayer–Rokitansky–Küster–Hauser Syndrome (MRKH, OMIM 277000) in females and, recently, were associated with intellectual disability, autistic features, cerebral anomaly and facial dysmorphisms.     

DNA-maleimide: An improved maleimide compound for electrophoresis-based titration of reactive thiols in a specific protein

Background

Thiol-mediated redox regulation of proteins plays a key role in many cellular processes.

Methods

To understand the redox status of cysteinyl thiol groups of the desired proteins, we developed a new maleimide reagent: a maleimide-conjugated single strand DNA, DNA-maleimide (DNA-Mal).

Results

DNA-Mal labelled proteins run as a distinct band on SDS-PAGE, with a discrete 9.32 kDa mobility shift per label regardless of the protein species or electrophoretic conditions.

Conclusions

DNA-Mal labels free thiols like standard maleimide reagents, but possesses practical advantages in titration of the number and relative content of free thiols in a protein.

General significance

The versatility of DNA molecule enhances the application of DNA-Mal in a broader range of cysteine containing proteins.     

Genetic polymorphism of CYP4A11 and CYP4A22 genes and in silico insights from comparative 3D modelling in a French population

The CYP4A subfamily is known to ω-hydroxylate the endogenous arachidonic acid into 20-hydroxyeicosatetranoic acid, which has renovascular and tubular functions. The aim of this work was to report a comprehensive investigation of the CYP4A11 and CYP4A22 genetic polymorphisms in a French population. Using PCR-SSCP and sequencing strategies, a total of 26 sequence variations were identified comprising 3 missense mutations for CYP4A11 (Ser404Phe, Phe434Ser and Arg505His) and 7 missense mutations for CYP4A22 (Arg126Trp, Gly130Ser, Asn152Tyr, Val185Phe, Cys231Arg, Leu428Pro and Leu509Phe). In comparison with SNPs reported in the database (dbSNP) of the National Center for Biotechnology information (NCBI), 6 and 3 novel polymorphisms were identified in CYP4A11 and CYP4A22, respectively. The potential impact of the amino acid substitutions on the structure and/or catalytic activity of the enzymes has been estimated by the construction and validation of the CYP4A 3D models. These results could be helpful for further investigations of the potential role of CYP4A variants in the genetic susceptibility to cardiovascular diseases in humans such as arterial hypertension.     

Genomic diversity and affinities in population groups of North West India: An analysis of Alu insertion and a single nucleotide polymorphism

The North West region of India is extremely important to understand the peopling of India, as it acted as a corridor to the foreign invaders from Eurasia and Central Asia. A series of these invasions along with multiple migrations led to intermixture of variable populations, strongly contributing to genetic variations. The present investigation was designed to explore the genetic diversities and affinities among the five major ethnic groups from North West India; Brahmin, Jat Sikh, Bania, Rajput and Gujjar. A total of 327 individuals of the abovementioned ethnic groups were analyzed for 4 Alu insertion marker loci (ACE, PV92, APO and D1) and a Single Nucleotide Polymorphism (SNP) rs2234693 in the intronic region of the ESR1 gene. Statistical analysis was performed to interpret the genetic structure and diversity of the population groups. Genotypes for ACE, APO, ESR1 and PV92 loci were found to be in Hardy–Weinberg equilibrium in all the ethnic groups, while significant departures were observed at the D1 locus in every investigated population after Bonferroni's correction. The average heterozygosity for all the loci in these ethnic groups was fairly substantial ranging from 0.3927 ± 0.1877 to 0.4333 ± 0.1416. Inbreeding coefficient indicated an overall 10% decrease in heterozygosity in these North West Indian populations. The gene differentiation among the populations was observed to be of the order of 0.013. Genetic distance estimates revealed that Gujjars were close to Banias and Jat Sikhs were close to Rajputs. Overall the study favored the recent division of the populations of North West India into largely endogamous groups. It was observed that the populations of North West India represent a more or less homogenous genetic entity, owing to their common ancestral history as well as geographical proximity.     

Genetic variants in selenoprotein P plasma 1 gene (SEPP1) are associated with fasting insulin and first phase insulin response in Hispanics

Context

Insulin resistance is not fully explained on a molecular level, though several genes and proteins have been tied to this defect. Knockdowns of the SEPP1 gene, which encodes the selenoprotein P (SeP) protein, have been shown to increase insulin sensitivity in mice. SeP is a liver-derived plasma protein and a major supplier of selenium, which is a proposed insulin mimetic and antidiabetic agent.

Objective

SEPP1 single nucleotide polymorphisms (SNPs) were selected for analysis with glucometabolic measures.

Participants and measures

The study included1424 Hispanics from families in the Insulin Resistance Atherosclerosis Family Study (IRASFS). Additionally, the multi-ethnic Insulin Resistance Atherosclerosis Study was used. A frequently sampled intravenous glucose tolerance test was used to obtain precise measures of acute insulin response (AIR) and the insulin sensitivity index (S_I).

Design

21 SEPP1 SNPs (tagging SNPs (n = 12) from HapMap, 4 coding variants and 6 SNPs in the promoter region) were genotyped and analyzed for association.

Results

Two highly correlated (r² = 1) SNPs showed association with AIR (rs28919926; Cys368Arg; p = 0.0028 and rs146125471; Ile293Met; p = 0.0026) while rs16872779 (intronic) was associated with fasting insulin levels (p = 0.0097). In the smaller IRAS Hispanic cohort, few of the associations seen in the IRASFS were replicated, but meta-analysis of IRASFS and all 3 IRAS cohorts (N = 2446) supported association of rs28919926 and rs146125471 with AIR (p = 0.013 and 0.0047, respectively) as well as rs7579 with S_I (p = 0.047).

Conclusions

Overall, these results in a human sample are consistent with the literature suggesting a role for SEPP1 in insulin resistance.     

Association of adiponectin (AdipoQ) and sulphonylurea receptor (ABCC8) gene polymorphisms with Type 2 Diabetes in North Indian population of Punjab

In Type 2 Diabetes (T2D), adiponectin (AdipoQ) and sulphonylurea receptor genes (ABCC8) are important targets for candidate gene association studies. The single nucleotide polymorphisms (SNPs) in these genes have been associated with features of the metabolic syndrome across various populations. The present case–control study undertaken in the population of Punjab, evaluates the association of + 45T>G polymorphism in AdipoQ gene; and Exon16-3C>T as well as Exon18C>T polymorphisms in ABCC8 gene with T2D. These SNPs were genotyped in 200 T2D cases and 200 non-diabetic healthy controls using the PCR-RFLP method. The frequency of the minor G-allele for AdipoQ+ 45(T>G) polymorphism was significantly higher in T2D cases (29.0%) than in controls (21.5%) [P = 0.02, OR = 1.49 (1.07–2.04)]. The genetic model analysis revealed that the G-allele cumulatively provides nearly 1.59–1.78 fold increased risk to T2D under the additive (P = 0.009; OR = 1.59, 1.12–2.25 at 95% CI), dominant (TG/GG vs. TT) (P = 0.034, OR = 1.64, 1.04–2.56 at 95% CI) and codominant model (TG vs. TT/GG) (P = 0.014; OR = 1.78, 1.12–2.82 at 95% CI) after adjusting for confounding factors. However, no difference in the distribution of genotype and allele frequencies was observed for both the ABCC8 polymorphisms. The distribution of obesity profiles (BMI, WC and WHR) was also significantly different between cases and controls (P < 0.05). Higher BMI and central obesity were observed to increase the risk of T2D. G-allele of + 45(T>G) polymorphism in the adiponectin gene appears to be associated with increased risk of T2D, but the polymorphisms in sulphonylurea receptor gene do not seem to be associated with T2D in the population of Punjab.     

Polymorphism of BMP4 gene in Indian goat breeds differing in prolificacy

Bone morphogenetic proteins (BMPs) are members of the TGF-β (transforming growth factor-beta) superfamily, of which BMP4 is the most important due to its crucial role in follicular growth and differentiation, cumulus expansion and ovulation. Reproduction is a crucial trait in goat breeding and based on the important role of BMP4 gene in reproduction it was considered as a possible candidate gene for the prolificacy of goats. The objective of the present study was to detect polymorphism in intronic, exonic and 3′ un-translated regions of BMP4 gene in Indian goats. Nine different goat breeds (Barbari, Beetal, Black Bengal, Malabari, Jakhrana (Twinning > 40%), Osmanabadi, Sangamneri (Twinning 20–30%), Sirohi and Ganjam (Twinning < 10%)) differing in prolificacy and geographic distribution were employed for polymorphism scanning. Cattle sequence (AC_000167.1) was used to design primers for the amplification of a targeted region followed by direct DNA sequencing to identify the genetic variations. Single nucleotide polymorphisms (SNPs) were not detected in exon 3, the intronic region and the 3′ flanking region. A SNP (G1534A) was identified in exon 2. It was a non-synonymous mutation resulting in an arginine to lysine change in a corresponding protein sequence. G to A transition at the 1534 locus revealed two genotypes GG and GA in the nine investigated goat breeds. The GG genotype was predominant with a genotype frequency of 0.98. The GA genotype was present in the Black Bengal as well as Jakhrana breed with a genotype frequency of 0.02. A microsatellite was identified in the 3′ flanking region, only 20 nucleotides downstream from the termination site of the coding region, as a short sequence with more than nineteen continuous and repeated CA dinucleotides. Since the gene is highly evolutionarily conserved, identification of a non-synonymous SNP (G1534A) in the coding region gains further importance. To our knowledge, this is the first report of a mutation in the coding region of the caprine BMP4 gene. But whether the reproduction trait of goat is associated with the BMP4 polymorphism, needs to be further defined by association studies in more populations so as to delineate an effect on it.