Haplotype variation and linkage disequilibrium in 313 human genes

DeepHapNet: a haplotype assembly method based on RetNet and deep spectral clustering

Gene polymorphism originates from single-nucleotide polymorphisms (SNPs), and the analysis and study of SNPs are of great significance in the field of biogenetics. The haplotype, which consists of the sequence of SNP loci, carries more genetic information than a single SNP. Haplotype assembly plays a significant role in understanding gene function, diagnosing complex diseases, and pinpointing species genes. We propose a novel method, DeepHapNet, for haplotype assembly through the clustering of reads and learning correlations between read pairs. We employ a sequence model called Retentive Network (RetNet), which utilizes a multiscale retention mechanism to extract read features and learn the global relationships among them. Based on the feature representation of reads learned from the RetNet model, the clustering process of reads is implemented using the SpectralNet model, and, finally, haplotypes are constructed based on the read clusters. Experiments with simulated and real datasets show that the method performs well in the haplotype assembly problem of diploid and polyploid based on either long or short reads. The code implementation of DeepHapNet and the processing scripts for experimental data are publicly available at https://github.com/wjj6666/DeepHapNet.

Genetic Variation and Population Structure of Clonorchis sinensis: An In Silico Analysis

Clonorchis sinensis is an important zoonotic parasite that is mainly prevalent in China, Korea, Vietnam and the Russian Far East. To explore the genetic variation and population structure of C. sinensis, an in silico analysis was conducted based on mitochondrial cytochrome c oxidase subunit 1 (COX1), ribosomal internal transcribed spacer 1 (ITS1) and ribosomal internal transcribed spacer 2 (ITS2) sequences. The sequences obtained from NCBI were truncated for further analyses, including haplotype network, phylogenetic, gene flow, diversity and neutrality analyses. The results showed that there were 20, 11 and 4 haplotypes for COX1, ITS1 and ITS2, respectively. The results of both the haplotype network and phylogenetic analyses indicated that the haplotypes for each type of sequence from the same country were not all clustered together. Haplotype diversity values were all lower than 0.5. Values of nucleotide diversity were higher than 0.005, except for ITS2. Tajima's D and Fu's Fs values were all negative, and p-values showed significant differences, indicating that the population of C. sinensis is growing. Fst values were all lower than 0.05. In conclusion, this study found that there are specific variations of C. sinensis in different countries, and the population of this parasite is growing with less genetic variation. The findings provide a crucial foundation for understanding the molecular epidemiology and population dynamics of C. sinensis.

Genetic variation and population structure of Taenia multiceps (Coenurus cerebralis) based on mitochondrial cox1 gene: A comprehensive global analysis

Taenia multiceps is a neglected parasite having veterinary and public health importance. The predilection sites of the parasite larva (Coenurus cerebralis) are brain (cerebral coenurosis) and subcutaneous (non-cerebral coenurosis). There is a dearth of data regarding molecular characterization of T. multiceps and even fewer population structure-based studies on T. multiceps. The current study was conducted to provide epidemiological information regarding the global population structure of the parasite. The NCBI GenBank database was accessed to download the sequences of cox1 gene, which were further subjected to PopArt software to construct median-joining networks. The DnaSp software was used to compute neutrality and diversity indices. Host and region-wise indices of neutrality and diversity were also computed. There were 166 gene sequences found in the NCBI database. Followed by removal of short gene sequences, 143 were considered to perform bioinformatic analyses. A total of 30 haplotypes with 46 mutations and 23 parsimony informative sites were found. High diversity (Hd = 0.889, π = 0.01186) and negative but statistically insignificant neutrality indices (Tajima's D = -1.57659, Fu's Fs = -10.552) were found. Region-wise results revealed highest haplotype diversities in isolates from KSA (Hd = 1.00) followed by Greece and Italy (Hd = 0.962), and China (Hd = 0.931). Host-wise data analysis showed an overall negative Tajima's D value and there exists highest haplotype diversity in cattle (Hd = 1.00) followed by dogs (Hd = 0.833), sheep (Hd = 0.795) and goats (Hd = 0.788). The findings of the study indicate that the population diversity of T. multiceps will increase worldwide as shown by high diversity and negative neutrality indices. The findings of the study significantly add-in to the existing bank of knowledge about population structure of T. multiceps. We recommend conducting more studies employing different genetic markers to better comprehend the epidemiology of the parasite.

Haplotype-based breeding: A new insight in crop improvement

Haplotype-based breeding (HBB) is one of the cutting-edge technologies in the realm of crop improvement due to the increasing availability of Single Nucleotide Polymorphisms identified by Next Generation Sequencing technologies. The complexity of the data can be decreased with fewer statistical tests and a lower probability of spurious associations by combining thousands of SNPs into a few hundred haplotype blocks. The presence of strong genomic regions in breeding lines of most crop species facilitates the use of haplotypes to improve the efficiency of genomic and marker-assisted selection. Haplotype-based breeding as a Genomic Assisted Breeding (GAB) approach harnesses the genome sequence data to pinpoint the allelic variation used to hasten the breeding cycle and circumvent the challenges associated with linkage drag. This review article demonstrates ways to identify candidate genes, superior haplotype identification, haplo-pheno analysis, and haplotype-based marker-assisted selection. The crop improvement strategies that utilize superior haplotypes will hasten the breeding progress to safeguard global food security.

A fine-scale genetic map of the Japanese population

Genetic maps are fundamental resources for linkage and association studies. A fine-scale genetic map can be constructed by inferring historical recombination events from the genome-wide structure of linkage disequilibrium-a non-random association of alleles among loci-by using population-scale sequencing data. We constructed a fine-scale genetic map and identified recombination hotspots from 10 092 551 bi-allelic high-quality autosomal markers segregating among 150 unrelated Japanese individuals whose genotypes were determined by high-coverage (30×) whole-genome sequencing, and the genotype quality was carefully controlled by using their parents' and offspring's genotypes. The pedigree information was also utilized for haplotype phasing. The resulting genome-wide recombination rate profiles were concordant with those of the worldwide population on a broad scale, and the resolution was much improved. We identified 9487 recombination hotspots and confirmed the enrichment of previously known motifs in the hotspots. Moreover, we demonstrated that the Japanese genetic map improved the haplotype phasing and genotype imputation accuracy for the Japanese population. The construction of a population-specific genetic map will help make genetics research more accurate.

Unravelling the potential of TIM-3 gene polymorphism in allogeneic hematopoietic stem cell transplantation - a preliminary study

BACKGROUND

T-cell immunoglobulin and mucin-domain containing-3 (TIM-3) molecule is a key regulator of the immune response by exerting an inhibitory effect on various types of immune cells. Understanding the role of TIM-3 in hematopoietic stem cell transplantation (HSCT) may improve transplant outcomes. Our study evaluated the potential association between TIM-3 polymorphisms, namely rs1036199 (A > C) or rs10515746 (C > A), changes which are located in exon 3 and the promoter region of the TIM-3 gene, and post-HSCT outcomes.

METHODS

One-hundred and twenty allogeneic HSCT patients and their respective donors were enrolled and genotyped for TIM-3 single nucleotide polymorphisms (SNPs) using real-time PCR with TaqMan assays.

RESULTS

We found that the presence of the rare alleles and heterozygous genotypes of studied SNP in recipients tended to protect against or increase the risk for acute graft-versus-host disease (aGvHD). For the rs1036199 polymorphism, recipients with the AC heterozygous genotype (p = 0.0287) or carrying the rarer C allele (p = 0.0334) showed a lower frequency of aGvHD development along all I-IV grades. A similar association was detected for the rs10515746 polymorphism as recipients with the CA genotype (p = 0.0095) or the recessive A allele (p = 0.0117) less frequently developed aGvHD. Furthermore, the rarer A allele of rs10515746 SNP was also associated with a prolonged aGvHD-free survival (p = 0.0424). Cytomegalovirus (CMV) infection was more common in patients transplanted with TIM-3 rs10515746 mismatched donors (p = 0.0229) and this association was also found to be independent of HLA incompatibility and pre-transplant CMV-IgG status. Multivariate analyses confirmed the role of these recessive alleles and donor-recipient TIM-3 incompatibility as an independent factor in aGvHD and CMV development.

CONCLUSIONS

Polymorphism of TIM-3 molecule may affect the immune response in HSCT patients. The recessive alleles of rs1036199 and rs10515746 SNPs decreased the risk of developing aGvHD. TIM-3 donor-recipient genetic matching may also affect the risk of post-transplant CMV infection, indicating the potential value of genetic profiling in optimizing transplant strategies.

Genome-Wide Association Analyses and Population Verification Highlight the Potential Genetic Basis of Horned Morphology during Polled Selection in Tibetan Sheep

The types and morphology of sheep horns have been extensively researched, yet the genetic foundation underlying the emergence of diverse horn characteristics during the breeding of polled Tibetan sheep has remained elusive. Genome-wide association analysis (GWAS) was performed on 103 subtypes (normal large horn, scurs, and polled) differentiated from G2 (offspring (G2) of parent (G1) of polled) of the polled core herd. Six single nucleotide polymorphisms (SNPs) located on chromosome 10 of the relaxin family peptide receptor 2 (RXFP2) gene exhibited positive correlations with horn length, horn base circumference, and horn base interval. Furthermore, in genotyping 382 G2 individuals, significant variations were observed for each specific horn type. Three additional mutations were identified near the target SNP upstream of the amplification product. Finally, the RXFP2-specific haplotype associated with the horned trait effectively maintained horn length, horn base circumference, and horn base interval in Tibetan sheep, as confirmed by population validation of nine loci in a sample size of 1125 individuals. The present study offers novel insights into the genetic differentiation of the horned type during improvement breeding and evolution, thereby establishing a robust theoretical foundation for polled Tibetan sheep breeding and providing valuable guidance for practical production.

Polymorphisms of ITGA9 Gene and Their Correlation with Milk Quality Traits in Yak (Bos grunniens)

A single-nucleotide polymorphism (SNP) is a genome-level trait that arises from a variation in a single nucleotide, leading to diversity in DNA sequences. SNP screening is commonly used to provide candidate genes for yak breeding efforts. Integrin Subunit Alpha 9 (ITGA9) is an integrin protein. It plays an important role in cell adhesion, signalling, and other processes. The aim of this study was to discuss the association between genetic polymorphisms in the ITGA9 gene and milk quality traits and to identify potential molecular marker loci for yak breeding quality. We genotyped 162 yaks using an Illumina Yak cGPS 7K liquid chip and identified the presence of polymorphisms at nine SNP loci in the ITGA9 gene of yaks. The results showed that the mutant genotypes in the loci g.285,808T>A, g.306,600T>C, and g.315,413C>T were positively correlated with the contents of casein, protein, total solids (TS), and solid nonfat (SNF) in yak milk. In other loci, heterozygous genotypes had a positive correlation with nutrient content in yak milk. Then, two ITGA9 haplotype blocks were constructed based on linkage disequilibrium, which facilitated a more accurate screening of ITGA9 as a candidate gene for yak milk quality improvement. In conclusion, we identified SNPs and haplotype blocks related to yak milk quality traits and provided genetic resources for marker-assisted selection in yak breeding.

Global scenario of genetic diversity in cox1 and nad1 genes of Moniezia expansa

Monieziasis is a parasite-borne production-limiting disease of livestock. Moniezia expansa is the most important species having cosmopolitan distribution. Despite of numerous prevalence reports, very little information is available about the evolutionary biology and population genetics of M. expansa. To close this research gap, this study was undertaken to recognize and inspect the genetic variation of M. expansa populations around the world using the cox1 and nad1 genes and deduce phylogenetic relationships with M. expansa populations. The cox1 and nad1 gene sequences were downloaded from the NCBI GenBank database. Followed by sequence alignment, median-joining networks were constructed using PopArt software. Diversity and neutrality indices were computed through DnaSp software while MEGA software was used to draw the maximum-likelihood phylogenetic tree. Thirty-two cox1 sequences, from five different countries, and 9 nad1 sequences from three different countries, were among the sequences used in this study. The cox1 and nad1 gene sequences had mutations in 97 and 36 different places, respectively. Twenty and 7 unique haplotypes were discovered for the cox1 and nad1 gene sequences, respectively. Comparable haplotype diversities were observed for both the genes under study (cox1 = 0.950; nad1 = 0.944). Negative Tajima's D and Fu Fs were found for the cox1 gene while these indices were positive for the nad1 gene. Phylogenetic analysis also showed the existence of unique haplotypes for both the cox1 and nad1 genes. The results of this study indicate that there is the existence of a huge genetic diversity in M. expansa isolates. For future studies, it is recommended that longer gene sequences should be used to describe genetic variation among M. expansa isolates as the length of the gene under study affects the genetic variation. Moreover, additional mitochondrial markers should also be investigated because the assertive strength of a group of gene targets is superior to defining genetic diversity.

A computational analysis to evaluate deleterious SNPs of GSK3β, a multifunctional and regulatory protein, for metabolism, wound healing, and migratory processes

This study focused on GSK-3β, a critical serine/threonine kinase with diverse cellular functions. However, there is limited understanding of the impact of non-synonymous single nucleotide polymorphisms (nsSNPs) on its structure and function. Through an exhaustive in-silico investigation 12 harmful nsSNPs were predicted from a pool of 172 acquired from the NCBI dbSNP database using 12 established tools that detects deleterious SNPs. Consistently, these nsSNPs were discovered in locations with high levels of conservation. Notably, the three harmful nsSNPs F67C, A83T, and T138I were situated in the active/binding site of GSK-3β, which may affect the protein's capacity to bind to substrates and other proteins. Molecular dynamics simulations revealed that the F67C and T138I mutants had stable structures, indicating rigidness, whereas the A83T mutant was unstable. Analysis of secondary structures revealed different modifications in all mutant forms, which may affect the stability, functioning, and interactions of the protein. These mutations appear to alter the structural dynamics of GSK-3β, which may have functional ramifications, such as the formation of novel secondary structures and variations in coil-to-helix transitions. In conclusion, this study illuminates the possible structural and functional ramifications of these GSK-3 nsSNPs, revealing how protein compactness, stiffness, and interactions may affect biological activities.

Insertion/deletions within the bovine FoxO1 gene and their association analysis with growth traits in three Chinese cattle breeds

FOXO1 (FKHR) gene, as a transcription factor, plays a vital role in animal growth and development, participating in many biological processes. The aim of this study was to ascertain Insertion/deletions (Indels) polymorphism within bovine FoxO1 gene in 679 Chinese adult cows and associate them with stature traits. Two Indels (named as Indel-3 and Indel-4, recorded as rs383545622 and rs525318770 in NCBI, respectively) were successfully genotyped by the Once PCR method, which was reliable, rapid and cost effective for simultaneous detection of two or more Indels. Indel-3 and Indel-4 were located at the second intron. All four different haplotypes (H1: D3D4, H2: I3D4, H3: D3I4, H4: I3I4) could be identified, and the D (del-) allele, DD (del-/del-) genotype and D3D4 haplotype retained the highest frequency. However, individuals with DI (D3I3, D4I4 or H1H4/H2H3 genotype) showed significantly better phenotypic traits than those with the other genotypes in Nanyang cattle, showing a hybrid vigor. The results implied that this DI genotype can be applied to early selective breeding to improve the productivity of Nanyang cattle. Our results suggested that these two Indels within the bovine FoxO1 gene might be used as genetic markers for marker-assisted selection (MAS) in cattle breeding and genetics.

Serum Phosphorus, Serum Bicarbonate, and Renal Function in Relation to Liver CYP1A2 Activity

The liver plays an important role in normal metabolism and physiological functions such as acid-base balance; however, limited epidemiologic studies have investigated how the liver contributes toward acid-base balance using non-invasive biomarkers. We determined associations between serum biomarkers related to acid-base balance and renal function with liver CYP1A2 activity. We used data from 1381 participants of the 2009-2010 National Health and Nutrition Examination Survey (NHANES) with measurements of serum phosphorus, serum bicarbonate, caffeine intake, caffeine metabolites, and estimated glomerular filtration rate (eGFR). Liver CYP1A2 activity was estimated using urine caffeine metabolite indices, which were calculated as the ratio of one of the urine caffeine metabolites (i.e., paraxanthine and 1-methyluric acid) to caffeine intake. We analyzed associations in the whole data set and in different strata of hepatic steatosis index (HSI) based on different cut-points. We found that serum bicarbonate was positively associated with CYP1A2 activity in the whole data set when comparing persons with bicarbonate at Q4 to Q1 (β = 0.18, p = 0.10 for paraxanthine; β = 0.20, p = 0.02 for 1-methyluric acid). Furthermore, serum phosphorus was positively associated with CYP1A2 activity only in the stratum of 30 ≤ HSI < 36. Lastly, low eGFR was significantly associated with lower CYP1A2 activity measured with paraxanthine in the whole dataset and in all the strata with HSI < 42; when comparing eGFR < 60 to eGFR > 90, β estimates ranged from -0.41 to -1.38, p-values ranged from 0.0018 to 0.004. We observed an opposite trend in the highest stratum (HSI ≥ 42). Non-invasive measurements of serum bicarbonate, serum phosphorus, and eGFR have dynamic associations with CYP1A2 activity. These associations depend on the extent of liver damage and the caffeine metabolite used to assess CYP1A2 activity.

Genetic variation and population structure of Fasciola hepatica: an in silico analysis

Fasciola hepatica is a trematode leading to heavy economic setbacks to the livestock sector globally. The population's genetic information and intimate kinship level are frequently assessed using analysis of mitochondrial DNA. In this analysis, we retrieved cox1 (n = 247) and nad1 (n = 357) sequences of F. hepatica from the NCBI GenBank database and aligned the sequences with the respective reference sequences using MEGA software. The median joining network was drawn using PopArt software while neutrality and diversity indices were estimated with the help of DnaSp software. Neighbor-joining phylogenetic tree was constructed using the MEGA software package. A total of 46 and 98 distinctive haplotypes were observed for cox1 and nad1 genes, respectively. Diversity indices indicated high haplotype and nucleotide diversities in both genes. Positive Tajima's D and Fu's Fs values were found for the entire population of both the genes under study. The cox1 and nad1 gene segments in this study showed high Tajima's D values, suggesting a low likelihood of future population growth. The Tajima's D value of the nad1 gene sequence is lower (2.14910) than that of the cox1 gene sequence (3.40314), which suggests that the former is growing at a slower rate. However, the region-wise analysis revealed that both the cox1 and nad1 genes showed deviation from neutrality suggesting a recent population expansion as a result of an excess of low-frequency polymorphism. Furthermore, the overall host-wise analysis showed positive and significant Tajima's D values for the cox1 and nad1 gene sequences. To the best of our knowledge, this is the first attempt to provide insights into genetic variations and population structure of F. hepatica at a global scale using cox1 and nad1 genes. Our findings suggest the existence of specific variants of F. hepatica in different parts of the world and provide information on the molecular ecology of F. hepatica. The results of this study also mark a critical development in upcoming epidemiological investigations on F. hepatica and will also contribute to understanding the global molecular epidemiology and population structure of F. hepatica.

A dual gene-specific mutator system installs all transition mutations at similar frequencies in vivo

Targeted in vivo hypermutation accelerates directed evolution of proteins through concurrent DNA diversification and selection. Although systems employing a fusion protein of a nucleobase deaminase and T7 RNA polymerase present gene-specific targeting, their mutational spectra have been limited to exclusive or dominant C:G→T:A mutations. Here we describe eMutaT7transition, a new gene-specific hypermutation system, that installs all transition mutations (C:G→T:A and A:T→G:C) at comparable frequencies. By using two mutator proteins in which two efficient deaminases, PmCDA1 and TadA-8e, are separately fused to T7 RNA polymerase, we obtained similar numbers of C:G→T:A and A:T→G:C substitutions at a sufficiently high frequency (∼6.7 substitutions in 1.3 kb gene during 80-h in vivo mutagenesis). Through eMutaT7transition-mediated TEM-1 evolution for antibiotic resistance, we generated many mutations found in clinical isolates. Overall, with a high mutation frequency and wider mutational spectrum, eMutaT7transition is a potential first-line method for gene-specific in vivo hypermutation.

Mitochondrial DNA Deficiency and Supplementation in Sus scrofa Oocytes Influence Transcriptome Profiles in Oocytes and Blastocysts

Mitochondrial DNA (mtDNA) deficiency correlates with poor oocyte quality and fertilisation failure. However, the supplementation of mtDNA deficient oocytes with extra copies of mtDNA improves fertilisation rates and embryo development. The molecular mechanisms associated with oocyte developmental incompetence, and the effects of mtDNA supplementation on embryo development are largely unknown. We investigated the association between the developmental competence of Sus scrofa oocytes, assessed with Brilliant Cresyl Blue, and transcriptome profiles. We also analysed the effects of mtDNA supplementation on the developmental transition from the oocyte to the blastocyst by longitudinal transcriptome analysis. mtDNA deficient oocytes revealed downregulation of genes associated with RNA metabolism and oxidative phosphorylation, including 56 small nucleolar RNA genes and 13 mtDNA protein coding genes. We also identified the downregulation of a large subset of genes for meiotic and mitotic cell cycle process, suggesting that developmental competence affects the completion of meiosis II and first embryonic cell division. The supplementation of oocytes with mtDNA in combination with fertilisation improves the maintenance of the expression of several key developmental genes and the patterns of parental allele-specific imprinting gene expression in blastocysts. These results suggest associations between mtDNA deficiency and meiotic cell cycle and the developmental effects of mtDNA supplementation on Sus scrofa blastocysts.

In Silico Evaluation of the Haplotype Diversity, Phylogenetic Variation and Population Structure of Human E. granulosus sensu stricto (G1 Genotype) Sequences

Echinococcus granulosus sensu lato is the causative agent of cystic echinococcosis (CE), which is a neglected zoonotic disease with an important role in human morbidity. In this study, we aimed to investigate the haplotype diversity, genetic variation, population structure and phylogeny of human E. granulosus sensu stricto (s.s.) (G1 genotype) isolates submitted to GenBank from different parts of the world by sequencing the mitochondrial CO1 and ND1 genes. The sequences of the mt-CO1 (401 bp; n = 133) and mt-ND1 (407 bp; n = 140) genes were used to analyze the haplotype, polymorphism and phylogenetic of 273 E. granulosus s.s. (G1 genotype) isolates. Mutations were observed at 31 different points in the mt-CO1 gene sequences and at 100 different points in the mt-ND1 gene sequences. Furthermore, 34 haplotypes of the mt-CO1 sequences and 37 haplotypes of the mt-ND1 sequences were identified. Tajima's D, Fu's Fs, and Fu's LD values showed high negative values in both mt-CO1 and mt-ND1 gene fragments. The haplotype diversities in the sequences retrieved from GenBank in this study indicate that the genetic variation in human isolates of E. granulosus s.s. in western countries is higher than in eastern countries. This may be due to demographic expansions due to animal trades and natural selections.

Regional association and transcriptome analysis revealed candidate genes controlling plant height in Brassica napus

Plant height is a key morphological trait in rapeseed, which not only plays an important role in determining plant architecture, but is also an important characteristic related to yield. Presently, the improvement of plant architecture is a major challenge in rapeseed breeding. This work was carried out to identify genetic loci related to plant height in rapeseed. In this study, a genome-wide association study (GWAS) of plant height was performed using a Brassica 60 K Illumina Infinium SNP array and 203 Brassica napus accessions. Eleven haplotypes containing important candidate genes were detected and significantly associated with plant height on chromosomes A02, A03, A05, A07, A08, C03, C06, and C09. Moreover, regional association analysis of 50 resequenced rapeseed inbred lines was used to further analyze these eleven haplotypes and revealed nucleotide variation in the BnFBR12-A08 and BnCCR1-C03 gene regions related to the phenotypic variation in plant height. Furthermore, coexpression network analysis showed that BnFBR12-A08 and BnCCR1-C03 were directly connected with hormone genes and transcription factors and formed a potential network regulating the plant height of rapeseed. Our results will aid in the development of haplotype functional markers to further improve plant height in rapeseed.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-022-01337-1.

TaNAC020 homoeologous genes are associated with higher thousand kernel weight and kernel length in Chinese wheat

NAC proteins constitute one of the largest plant-specific transcription factor (TF) families and play significant roles in plant growth and development. In the present study, three TaNAC020 homoeologous genes located on chromosomes 7A, 7B, and 7D were isolated from wheat (Triticum aestivum L.). TaNAC020s were predominantly expressed in developing grains. The developed transgenic rice lines for TaNAC020-B showed higher starch density and lower amylose contents than those of the wild type (WT). Sequence polymorphism studies showed seven and eight SNPs in TaNAC020-A/B, making three and two haplotypes, respectively. No sequence polymorphism was identified in TaNAC020-D. Association analysis revealed that HAP-2 of TaNAC020-A and TaNAC020-B was the favored haplotype for higher thousand kernel weight and length. Geographic distribution and allelic frequency showed that our favored haplotype experienced strong selection in China, and likewise, diversity increased in TaNAC020s during wheat polyploidization. The results obtained in this study demonstrate that TaNAC020s positively influence starch synthesis and accumulation and are one of the key regulators of the kernel (seed) size and kernel number and have the potential for utilization in wheat breeding to improve grain yield. Molecular markers developed in this study stand instrumental in marker-assisted selection for genetic improvement and germplasm enhancement in wheat.

A Study of Combined Genotype Effects of SHCBP1 on Wool Quality Traits in Chinese Merino

SHCBP1 (Shc SH2-domain binding protein 1) is a member of the Src and collagen homolog (Shc) protein family and is closely associated with multiple signaling pathways that play important roles during hair follicle induction, morphogenesis, and cycling. The purpose of this study was to investigate SHCBP1 gene expression, polymorphisms, and the association between SHCBP1 and wool quality traits in Chinese Merino sheep. The SHCBP1 gene was shown, by qPCR, to be ubiquitously expressed in sheep tissues and differentially expressed in the adult skin of Chinese Merino and Suffolk sheep. Four SNPs (termed SHCBP1SNPs 1-4) were identified by Sanger sequencing and were located in exon 2, intron 9, intron 12, and exon 13 of the sheep SHCBP1 gene, respectively. SHCBP1SNPs 3 and 4 (rs411176240 and rs160910635) were significantly associated with wool crimp (P < 0.05). The combined polymorphism (SHCBP1SNP3-SHCBP1SNP4) was significantly associated with wool crimp (P < 0.05). Bioinformatics analysis showed that the SNPs associated with wool crimp (SHCBP1SNPs 3 and 4) might affect the pre-mRNA splicing by creating binding sites for serine-arginine-rich proteins and that SHCBP1SNP4 might alter the SHCBP1 mRNA and protein secondary structure. Our results suggest that SHCBP1 influences wool crimp and SHCBP1SNPs 3 and 4 might be useful markers for marker-assisted selection and sheep breeding.

Determining the Effects of Chronic Kidney Disease on Organic Anion Transporter1/3 Activity Through Physiologically Based Pharmacokinetic Modeling

BACKGROUND AND OBJECTIVE

The renal excretion of drugs via organic anion transporters 1 and 3 (OAT1/3) is significantly decreased in patients with renal impairment. This study uses physiologically based pharmacokinetic models to quantify the reduction in OAT1/3-mediated secretion of drugs throughout varying stages of chronic kidney disease.

METHODS

Physiologically based pharmacokinetic models were constructed for four OAT1/3 substrates in healthy individuals: acyclovir, meropenem, furosemide, and ciprofloxacin. Observed data from drug-drug interaction studies with probenecid, a potent OAT1/3 inhibitor, were used to parameterize the contribution of OAT1/3 to the renal elimination of each drug. The models were then translated to patients with chronic kidney disease by accounting for changes in glomerular filtration rate, kidney volume, renal blood flow, plasma protein binding, and hematocrit. Additionally, a relationship was derived between the estimated glomerular filtration rate and the reduction in OAT1/3-mediated secretion of drugs based on the renal extraction ratios of ƿ-aminohippuric acid in patients with varying degrees of renal impairment. The relationship was evaluated in silico by evaluating the predictive performance of each final model in describing the pharmacokinetics (PK) of drugs across stages of chronic kidney disease.

RESULTS

OAT1/3-mediated renal excretion of drugs was found to be decreased by 27-49%, 50-68%, and 70-96% in stage 3, stage 4, and stage 5 of chronic kidney disease, respectively. In support of the parameterization, physiologically based pharmacokinetic models of four OAT1/3 substrates were able to adequately characterize the PK in patients with different degrees of renal impairment. Total exposure after intravenous administration was predicted within a 1.5-fold error and 85% of the observed data points fell within a 1.5-fold prediction error. The models modestly under-predicted plasma concentrations in patients with end-stage renal disease undergoing intermittent hemodialysis. However, results should be interpreted with caution because of the limited number of molecules analyzed and the sparse sampling in observed chronic kidney disease pharmacokinetic studies.

CONCLUSIONS

A quantitative understanding of the reduction in OAT1/3-mediated excretion of drugs in differing stages of renal impairment will contribute to better predictive accuracy for physiologically based pharmacokinetic models in drug development, assisting with clinical trial planning and potentially sparing this population from unnecessary toxic exposures.

Genetic Polymorphisms of IGF1 and IGF1R Genes and Their Effects on Growth Traits in Hulun Buir Sheep

The identification of candidate genes and genetic variations associated with growth traits is important for sheep breeding. Insulin like growth factor 1 (IGF1) and insulin like growth factor 1 receptor (IGF1R) are well-accepted candidate genes that affect animal growth and development. The current study attempted to assess the association between IGF1 and IGF1R genetic polymorphisms and growth traits in Hulun Buir sheep. To achieve this goal, we first identified three and ten single nucleotide polymorphisms (SNPs) in exons of IGF1 and IGF1R in Hulun Buir sheep and then constructed six haplotypes of IGF1R based on linkage disequilibrium, respectively. Association studies were performed between SNPs and haplotypes of IGF1 and IGF1R with twelve growth traits in a population encompassing 229 Hulun Buir sheep using a general linear model. Our result indicated three SNPs in IGF1 were significantly associated with four growth traits (p < 0.05). In IGF1R, three SNPs and two haplotype blocks were significantly associated with twelve growth traits (p < 0.05). The combined haplotype H5H5 and H5H6 in IGF1R showed the strong association with 12 superior growth traits in Hulun Buir sheep (p < 0.05). In conclusion, we identified SNPs and haplotype combinations associated with the growth traits, which provided genetic resources for marker-assisted selection (MAS) in Hulun Buir sheep breeding.

Kam Sweet Rice (Oryza sativa L.) Is a Special Ecotypic Rice in Southeast Guizhou, China as Revealed by Genetic Diversity Analysis

Kam Sweet Rice (KSR) is a special kind of rice landrace that has been domesticated for thousands of years by the local Dong people in southeast Guizhou province, China. KSR has many distinguishing characteristics including strong fragrance; high resistance to diseases, pests, and adverse abiotic conditions; difficulty of threshing; and glutinous texture. There is a lack of systematic research on its genetic diversity. In this study, we analyzed the levels and patterns of genetic diversity and nucleotide variation in 1,481 rice germplasm using simple sequence repeat (SSR) markers and single nucleotide polymorphism (SNP) haplotype analysis of six unlinked nuclear loci. The accessions included 315 KSR resources from southeast Guizhou, 578 rice landraces from six rice-growing ecological zones in Guizhou, 546 rice landraces from nine provinces around Guizhou, and 42 wild rice sources. Genetic diversity and heterozygosity of KSR were both low, and thus KSR might be close to a pure rice line. Population structure analysis showed that KSR was isolated into a single type of rice, which had a large genetic distance and a unique genetic background compared to the local varieties in Guizhou province, indicating that KSR is a special rice ecotype. Haplotype analysis of the target genes showed that the population of KSR was rich in haplotypes for resistance to bacterial blight (Xa23) and rice blast (Pid3), and identified unique haplotypes that were different from those of the six rice ecotypes in Guizhou. This study shows that KSR is an excellent rice germplasm resource, provides important information for the improvement and utilization of rice landraces, and serves as a reference for formulating effective rice conservation measures.

The recombination landscapes of spiny lizards (genus Sceloporus)

Despite playing a critical role in evolutionary processes and outcomes, relatively little is known about rates of recombination in the vast majority of species, including squamate reptiles-the second largest order of extant vertebrates, many species of which serve as important model organisms in evolutionary and ecological studies. This paucity of data has resulted in limited resolution on questions related to the causes and consequences of rate variation between species and populations, the determinants of within-genome rate variation, as well as the general tempo of recombination rate evolution on this branch of the tree of life. In order to address these questions, it is thus necessary to begin broadening our phylogenetic sampling. We here provide the first fine-scale recombination maps for two species of spiny lizards, Sceloporus jarrovii and Sceloporus megalepidurus, which diverged at least 12 Mya. As might be expected from similarities in karyotype, population-scaled recombination landscapes are largely conserved on the broad-scale. At the same time, considerable variation exists at the fine-scale, highlighting the importance of incorporating species-specific recombination maps in future population genomic studies.

Mendelian randomization in the multivariate general linear model framework

Mendelian randomization (MR) is an application of instrumental variable (IV) methods to observational data in which the IV is a genetic variant. MR methods applicable to the general exponential family of distributions are currently not well characterized. We adapt a general linear model framework to the IV setting and propose a general MR method applicable to any full-rank distribution from the exponential family. Empirical bias and coverage are estimated via simulations. The proposed method is compared to several existing MR methods. Real data analyses are performed using data from the REGARDS study to estimate the potential causal effect of smoking frequency on stroke risk in African Americans. In simulations with binary variates and very weak instruments the proposed method had the lowest median [Q₁ , Q₃ ] bias (0.10 [-3.68 to 3.62]); compared with 2SPS (0.27 [-3.74 to 4.26]) and the Wald method (-0.69 [-1.72 to 0.35]). Low bias was observed throughout other simulation scenarios; as well as more than 90% coverage for the proposed method. In simulations with count variates, the proposed method performed comparably to 2SPS; the Wald method maintained the most consistent low bias; and 2SRI was biased towards the null. Real data analyses find no evidence for a causal effect of smoking frequency on stroke risk. The proposed MR method has low bias and acceptable coverage across a wide range of distributional scenarios and instrument strengths; and provides a more parsimonious framework for asymptotic hypothesis testing compared to existing two-stage procedures.

The Correlation Between FGB Promoter Polymorphism and Clotting Function in Patients With Idiopathic Lower Extremity Deep Venous Thrombosis

To explore the possible single nucleotide polymorphisms (SNPs) sites in the promoter region of fibrinogen B β (FGB), and construct logistic regression model and haplotype model, so as to reveal the influence of FGB promoter SNPs on susceptibility, hemodynamics and coagulation function of lower extremity deep venous thrombosis (LEDVT) in the genetic background. LEDVT patients (120) and healthy people (120) were taken as case and control objects, respectively. SNPs and their genotypes of FGB promoter were detected by promoter sequencing and PCR-RFLP. The parameters of coagulation system were evaluated. There were 6 SNPs in FGB promoter, which were β-148C/T, β-249C/T, β-455G/A, β-854G/A, β-993C/T and β-1420G/A. The genotype and allele frequency of β-1420 G/A, β-455G/A, β-249c/T and β-148C/T were significantly different between the LEDVT group and the control group, but not β-993C/T and β-854G/A. In addition, we found that the higher the content of Fibrinogen (FG), the higher the risk of LEDVT. The risk of LEDVT increased by 4.579 times for every unit increase of fibrinogen. We also found that FG, PT and APTT in LEDVT group were higher than those in control group, while TT was lower than those in control group; Furthermore, there was no significant difference in all coagulation indexes among 6 SNP genotypes in LEDVT group, while a significant difference was found between the 2 genotypes of β-993C/T in the control group. β-993C/T may indirectly affect the susceptibility of LEDVT by improving the basic level of plasma FG.

Evidence for opposing selective forces operating on human-specific duplicated TCAF genes in Neanderthals and humans

TRP channel-associated factor 1/2 (TCAF1/TCAF2) proteins antagonistically regulate the cold-sensor protein TRPM8 in multiple human tissues. Understanding their significance has been complicated given the locus spans a gap-ridden region with complex segmental duplications in GRCh38. Using long-read sequencing, we sequence-resolve the locus, annotate full-length TCAF models in primate genomes, and show substantial human-specific TCAF copy number variation. We identify two human super haplogroups, H4 and H5, and establish that TCAF duplications originated ~1.7 million years ago but diversified only in Homo sapiens by recurrent structural mutations. Conversely, in all archaic-hominin samples the fixation for a specific H4 haplotype without duplication is likely due to positive selection. Here, our results of TCAF copy number expansion, selection signals in hominins, and differential TCAF2 expression between haplogroups and high TCAF2 and TRPM8 expression in liver and prostate in modern-day humans imply TCAF diversification among hominins potentially in response to cold or dietary adaptations.

Duplications of gene segments can allow novel physiological adaptations to evolve. A detailed analysis of the TCAF gene family in primates and archaic humans suggest rapid duplication and diversification in this gene family is associated with cold or dietary adaptations.

Genomic regions associated with heat stress tolerance in tropical maize (Zea mays L.)

With progressive climate change and the associated increase in mean temperature, heat stress tolerance has emerged as one of the key traits in the product profile of the maize breeding pipeline for lowland tropics. The present study aims to identify the genomic regions associated with heat stress tolerance in tropical maize. An association mapping panel, called the heat tolerant association mapping (HTAM) panel, was constituted by involving a total of 543 tropical maize inbred lines from diverse genetic backgrounds, test-crossed and phenotyped across nine locations in South Asia under natural heat stress. The panel was genotyped using a genotyping-by-sequencing (GBS) platform. Considering the large variations in vapor pressure deficit (VPD) at high temperature (Tmax) across different phenotyping locations, genome-wide association study (GWAS) was conducted separately for each location. The individual location GWAS identified a total of 269 novel significant single nucleotide polymorphisms (SNPs) for grain yield under heat stress at a p value of < 10^–5. A total of 175 SNPs were found in 140 unique gene models implicated in various biological pathway responses to different abiotic stresses. Haplotype trend regression (HTR) analysis of the significant SNPs identified 26 haplotype blocks and 96 single SNP variants significant across one to five locations. The genomic regions identified based on GWAS and HTR analysis considering genomic region x environment interactions are useful for breeding efforts aimed at developing heat stress resilient maize cultivars for current and future climatic conditions through marker-assisted introgression into elite genetic backgrounds and/or genome-wide selection.

A novel Mendelian randomization method with binary risk factor and outcome

BACKGROUND

Mendelian randomization (MR) applies instrumental variable (IV) methods to observational data using a genetic variant as an IV. Several Monte-Carlo studies investigate the performance of MR methods with binary outcomes, but few consider them in conjunction with binary risk factors.

OBJECTIVE

To develop a novel MR estimator for scenarios with a binary risk factor and outcome; and compare to existing MR estimators via simulations and real data analysis.

METHODS

A bivariate Bernoulli distribution is adapted to the IV setting. Empirical bias and asymptotic coverage probabilities are estimated via simulations. The proposed method is compared to the Wald method, two-stage predictor substitution (2SPS), two-stage residual inclusion (2SRI), and the generalized method of moments (GMM). An analysis is performed using existing data from the CLEAR study to estimate the potential causal effect of smoking on rheumatoid arthritis risk in African Americans.

RESULTS

Bias was low for the proposed method and comparable to 2SPS. The Wald method was often biased towards the null. Coverage was adequate for the proposed method, 2SPS, and 2SRI. Coverage for the Wald and GMM methods was poor in several scenarios. The causal effect of ever smoking on rheumatoid arthritis risk was not statistically significant using a variety of genetic instruments.

CONCLUSIONS

Simulations suggest the proposed MR method is sound with binary risk factors and outcomes, and comparable to 2SPS and 2SRI in terms of bias. The proposed method also provides more natural framework for hypothesis testing compared to 2SPS or 2SRI, which require ad-hoc variance adjustments.

Soybean (Glycine max) Haplotype Map (GmHapMap): a universal resource for soybean translational and functional genomics

Here, we describe a worldwide haplotype map for soybean (GmHapMap) constructed using whole-genome sequence data for 1007 Glycine max accessions and yielding 14.9 million variants as well as 4.3 M tag single-nucleotide polymorphisms (SNPs). When sampling random subsets of these accessions, the number of variants and tag SNPs plateaued beyond approximately 800 and 600 accessions, respectively. This suggests extensive coverage of diversity within the cultivated soybean. GmHapMap variants were imputed onto 21 618 previously genotyped accessions with up to 96% success for common alleles. A local association analysis was performed with the imputed data using markers located in a 1-Mb region known to contribute to seed oil content and enabled us to identify a candidate causal SNP residing in the NPC1 gene. We determined gene-centric haplotypes (407 867 GCHs) for the 55 589 genes and showed that such haplotypes can help to identify alleles that differ in the resulting phenotype. Finally, we predicted 18 031 putative loss-of-function (LOF) mutations in 10 662 genes and illustrated how such a resource can be used to explore gene function. The GmHapMap provides a unique worldwide resource for applied soybean genomics and breeding.

Soybean (Glycine max) Haplotype Map (GmHapMap): a universal resource for soybean translational and functional genomics

Here, we describe a worldwide haplotype map for soybean (GmHapMap) constructed using whole-genome sequence data for 1007 Glycine max accessions and yielding 14.9 million variants as well as 4.3 M tag single-nucleotide polymorphisms (SNPs). When sampling random subsets of these accessions, the number of variants and tag SNPs plateaued beyond approximately 800 and 600 accessions, respectively. This suggests extensive coverage of diversity within the cultivated soybean. GmHapMap variants were imputed onto 21 618 previously genotyped accessions with up to 96% success for common alleles. A local association analysis was performed with the imputed data using markers located in a 1-Mb region known to contribute to seed oil content and enabled us to identify a candidate causal SNP residing in the NPC1 gene. We determined gene-centric haplotypes (407 867 GCHs) for the 55 589 genes and showed that such haplotypes can help to identify alleles that differ in the resulting phenotype. Finally, we predicted 18 031 putative loss-of-function (LOF) mutations in 10 662 genes and illustrated how such a resource can be used to explore gene function. The GmHapMap provides a unique worldwide resource for applied soybean genomics and breeding.

Association of UCP1 polymorphisms with type 2 diabetes mellitus and their interaction with physical activity and sedentary behavior

BACKGROUND

Uncoupling protein 1 (UCP1) has been reported to be associated with type 2 diabetes mellitus (T2DM) in different populations, however, little is reported in Chinese population. The present study aimed to explore the association between some polymorphisms of UCP1 with T2DM and the interactions between UCP1 and physical activity/sedentary behavior (PA/SB) lifestyle in Chinese population.

METHODS

Three polymorphisms (rs1472268, rs3811790 and rs3811791) were genotyped in 929 T2DM patients and 1044 nondiabetic controls. The data of PA and SB were acquired. Logistic regression and linear regression were conducted to assess the association of UCP1 and T2DM and related traits.

RESULTS

The CC genotype of rs3811791 was significantly associated with an increased risk of T2DM [odds ratio (OR) = 1.42, P = 0.042] and a higher level of triglyceride (TG) (β = 0.048, P = 0.034). This association still existed in the group of SB ≥ 3 h/d (OR = 1.66, P = 0.009) and the group of PA ≥ 150 min/week and SB ≥ 3 h/d (OR = 1.60, P = 0.034). In the group of PA < 150 min/week and SB < 3h/d, CC genotype was associated with a higher level of homeostatic model assessment of insulin resistance (HOMA-IR) index, and in the group of PA < 150 min/week and SB ≥ 3 h/d, CC genotype was associated with increased level of TG and decreased high-density lipoprotein cholesterol (HDL-C).

CONCLUSION

This study suggests that rs3811791 of UCP1 may be associated with T2DM and TG. Moreover, we demonstrate that SB interacted with rs3811791 of UCP1 was associated with T2DM, and PA interacted with rs3811791 of UCP1 was associated with the level of HOMA-IR, HDL-C, and TG.

Integrative phenotypic and gene expression data identify myostatin as a muscle growth inhibitor in Chinese shrimp Fenneropenaeus chinensis

Growth traits, largely determined by muscle growth, are the most critical economic traits in shrimp breeding. Myostatin (Mstn) is a conserved inhibitor of muscle growth in vertebrates, but until now solid evidence supporting a similar function of Mstn in invertebrates has been lacking. In the present study, we examined the Mstn expression along with growth trait data in a Fenneropenaeus chinensis population, to establish a potential correlation between Mstn and growth. The heritabilities of FcMstn expression, body weight at 190 days of culture, body weight and length at 230 days of culture, and average daily gain were estimated using 773 individuals and a thirteen-generation pedigree. The results showed FcMstn expression was negatively correlated with the growth traits, and the mean FcMstn expression in females was significantly lower than that of males, indicating Mstn negatively regulates muscle growth in shrimp, and its lower expression may underscore the faster growth of females. Low heritabilities were detected for FcMstn expression, suggesting that the expression of Mstn might be heritable in shrimp. These results provide strong support for a growth inhibitory function of Mstn in F. chinensis, and suggest a potential method for selective breeding of this species without substantial experimental resources and labor force.

Counting Oceanians of Non-European, Non-Asian Descent (ONENA) in the South Pacific to Make Them Count in Global Health

Several diseases and vulnerabilities associated with genetic or microbial factors are more frequent among populations of Oceanian, Non-European, Non-Asian descent (ONENA). ONENA are specific and have long been isolated geographically. To our knowledge, there are no published official, quantitative, aggregated data on the populations impacted by these excess vulnerabilities in Oceania. We searched official census reports for updated estimates of the total population for each of the Pacific Island Countries and Territories (including Australia) and the US State of Hawaii, privileging local official statistical or censual sources. We multiplied the most recent total population estimate by the cumulative percentage of the ONENA population as determined in official reports. Including Australia and the US State of Hawaii, Oceania counts 27 countries and territories, populated in 2016 by approximately 41 M inhabitants (17 M not counting Australia) among which approximately 12.5 M (11.6 M not counting Australia) consider themselves of entire or partial ONENA ancestry. Specific genetic and microbiome traits of ONENA may be unique and need further investigation to adjust risk estimates, risk prevention, diagnostic and therapeutic strategies, to the benefit of populations in the Pacific and beyond.

Association of TLR4 and TLR9 gene polymorphisms and haplotypes with cervicitis susceptibility

BACKGROUND

Cervicitis is one of the major health problems amongst women caused by infection of various pathogens including Chlamydia trachomatis (CT), Neisseria gonorrhoeae (NG), Trichomonas vaginalis (TV) as well as human papillomavirus (HPV), and persistent cervical inflammation is one of the etiologic agents of cervical cancer. Toll-like receptors (TLRs) play an important role in the recognition and subsequent elimination of these pathogens. Variations in the Toll-like receptor genes influence susceptibility to pathogens as well as disease progression independently.

METHODS

Ten single nucleotide polymorphisms, five each of TLR4 and TLR9 genes were analyzed among 130 cervicitis patients and 150 controls either using polymerase chain reaction-restriction fragment length polymorphism or allele specific-PCR.

RESULTS

T. vaginalis infection was found at the highest frequency (30.7%) as compared to C. trachomatis (1.5%), N. gonorrhoeae (2.3%) and HPV (4.6%) infections in cervicitis patients. TLR4 rs11536889 CC (age-adjusted OR, 2.469 [95% CI, 1.499 to 4.065]; p < 0.001) and TLR9 rs187084 TC (age-adjusted OR, 2.165 [95% CI, 1.267-3.699]; p = 0.005) genotypes showed the higher distribution in cervicitis patients compared to controls. In addition, TLR4 rs11536889 C allele was shown to increase the risk of cervicitis (age-adjusted OR, 1.632 [95% CI, 1.132 to 2.352]; p = 0.009) compared to controls. The TLR4 haplotype GCA (OR, 0.6 [95% CI, 0.38-0.95]; p = 0.0272) and TLR9 haplotype GTA (OR, 1.99 [95% CI, 1.14-3.48]; p = 0.014) were found to be associated with decreased and increased risk of cervicitis respectively.

CONCLUSIONS

TLR4 and TLR9 polymorphisms, as well as haplotypes were shown to modulate the cervicitis risk.

Association of TLR4 and TLR9 polymorphisms and haplotypes with cervical cancer susceptibility

Single nucleotide polymorphisms (SNPs) in TLR genes may serve as a crucial marker for early susceptibility of various cancers including cervical cancer. The present study was therefore designed to ascertain the role of TLR4 and TLR9 SNPs and haplotypes to hrHPV infection and cervical cancer susceptibility. The study included 110 cervical cancer biopsies and 141 cervical smears from age-matched healthy controls of Gujarati ethnicity of Western India. hrHPV 16 and 18 were detected using Real-time PCR. Eight SNPs, four each in TLR4 and TLR9 were analyzed using Polymerase Chain Reaction-Restriction Fragment Length Polymorphism and Allele-Specific PCR. HPV 16 and 18 were detected in 68% cervical cancer cases. TLR4 rs4986790, rs1927911 and TLR9 rs187084 showed association with HPV 16/18 infection. CC and CT genotypes of TLR4 rs11536889 and rs1927911 respectively, and TC, CC genotypes of TLR9 rs187084, as well as minor alleles of TLR4 rs4986790 and TLR9 rs187084, were associated with the increased risk of cervical cancer. Stage-wise analysis revealed TLR9 rs187084 and rs352140 to be associated with early-stage cancer. TLR4 haplotype GTAC and TLR9 haplotype GATC were associated with the increased risk of cervical cancer while TLR4 haplotype GCAG was associated with the decreased risk. TLR4 haplotype GCAG and TLR9 haplotype GATC showed association with increased susceptibility to hrHPV infection. In conclusion, the present study revealed association of TLR4 and TLR9 polymorphisms and haplotypes with hrHPV infection and cervical cancer risk. Further evaluation of a larger sample size covering diverse ethnic populations globally is warranted.

Genome-Wide Analyses Reveal Footprints of Divergent Selection and Drought Adaptive Traits in Synthetic-Derived Wheats

Crop-wild introgressions have long been exploited without knowing the favorable recombination points. Synthetic hexaploid wheats are one of the most exploited genetic resources for bread wheat improvement. However, despite some QTL with major effects, much less is known about genome-wide patterns of introgressions and their effects on phenotypes. We used two genome-wide association approaches: SNP-GWAS and haplotype-GWAS to identify SNPs and haplotypes associated with productivity under water-limited conditions in a synthetic-derived wheat (SYN-DER) population. Haplotype-GWAS further enriched and identified 20 more genomic regions associated with drought adaptability that did not overlap with SNP-GWAS. Since GWAS is biased to the phenotypes in the study and may fail to detect important genetic diversity during breeding, we used five complementary analytical approaches (t-test, Tajima’s D, nucleotide diversity (π), Fst, and EigenGWAS) to identify divergent selections in SYN-DER compared to modern bread wheat. These approaches consistently pinpointed 89 ‘selective sweeps’, out of which 30 selection loci were identified on D-genome. These key selections co-localized with important functional genes of adaptive traits such as TaElf3-D1 (1D) for earliness per se (Eps), TaCKX-D1 (3D), TaGS1a (6D) and TaGS-D1 (7D) for grain size, weight and morphology, TaCwi-D1 (5D) influencing drought tolerance, and Vrn-D3 (7D) for vernalization. Furthermore, 55 SNPs and 23 haplotypes of agronomic and physiological importance such as grain yield, relative water content and thousand grain weight in SYN-DER, were among the top 5% of divergent selections contributed by synthetic hexaploid wheats. These divergent selections associated with improved agronomic performance carry new alleles that have been introduced to wheat. Our results demonstrated that GWAS and selection sweep analyses are powerful approaches for investigating favorable introgressions under strong selection pressure and the use of crop-wild hybridization to assist the improvement of wheat yield and productivity under moisture limiting environments.

Hair Proteome Variation at Different Body Locations on Genetically Variant Peptide Detection for Protein-Based Human Identification

Human hair contains minimal intact nuclear DNA for human identification in forensic and archaeological applications. In contrast, proteins offer a pathway to exploit hair evidence for human identification owing to their persistence, abundance, and derivation from DNA. Individualizing single nucleotide polymorphisms (SNPs) are often conserved as single amino acid polymorphisms in genetically variant peptides (GVPs). Detection of GVP markers in the hair proteome via high-resolution tandem mass spectrometry permits inference of SNPs with known statistical probabilities. To adopt this approach for forensic investigations, hair proteomic variation and its effects on GVP identification must first be characterized. This research aimed to assess variation in single-inch head, arm, and pubic hair, and discover body location-invariant GVP markers to distinguish individuals. Comparison of protein profiles revealed greater body location-specific variation in keratin-associated proteins and intracellular proteins, allowing body location differentiation. However, robust GVP markers derive primarily from keratins that do not exhibit body location-specific differential expression, supporting GVP identification independence from hair proteomic variation at the various body locations. Further, pairwise comparisons of GVP profiles with 8 SNPs demonstrated greatest interindividual variation and high intraindividual consistency, enabling similar differentiative potential of individuals using single hairs irrespective of body location origin.

EcoTILLING Reveals Natural Allelic Variations in Starch Synthesis Key Gene TaSSIV and Its Haplotypes Associated with Higher Thousand Grain Weight

Wheat is a staple food commodity grown worldwide, and wheat starch is a valuable source of energy and carbon that constitutes 80% of the grain weight. Manipulation of genes involved in starch synthesis significantly affects wheat grain weight and yield. TaSSIV plays an important role in starch synthesis and its main function is granule formation. To mine and stack more favorable alleles, single nucleotide polymorphisms (SNPs) of TaSSIV-A, B, and D were investigated across 362 wheat accessions by Ecotype-Targeting Induced Local Lesions IN Genome (EcoTILLING). As a result, a total of 38 SNPs in the amplified regions of three TaSSIV genes were identified, of which 10, 15, and 13 were in TaSSIV-A, B, and D, respectively. These 38 SNPs were evaluated by using KASP and six SNPs showed an allele frequency >5% whereas the rest were <5%, i.e., considered to be minor alleles. In the Chinese mini core collection, three haplotypes were detected for TaSSIV–A and three for TaSSIV–B. The results of an association study in the Chinese mini core collection with thousand grain weight (TGW) and spike length (SPL) showed that Hap-2-1A was significantly associated with TGW and Hap-3-1B with SPL. Allelic frequency and geographic distribution indicated that the favored haplotype (Hap-2-1A) has been positively selected in Chinese wheat breeding. These results suggested that the Kompetitive Allele Specific PCR (KASP) markers can be applied in starch improvement to ultimately improve wheat yield by marker assisted selection in wheat breeding.

Apple whole genome sequences: recent advances and new prospects

In 2010, a major scientific milestone was achieved for tree fruit crops: publication of the first draft whole genome sequence (WGS) for apple (Malus domestica). This WGS, v1.0, was valuable as the initial reference for sequence information, fine mapping, gene discovery, variant discovery, and tool development. A new, high quality apple WGS, GDDH13 v1.1, was released in 2017 and now serves as the reference genome for apple. Over the past decade, these apple WGSs have had an enormous impact on our understanding of apple biological functioning, trait physiology and inheritance, leading to practical applications for improving this highly valued crop. Causal gene identities for phenotypes of fundamental and practical interest can today be discovered much more rapidly. Genome-wide polymorphisms at high genetic resolution are screened efficiently over hundreds to thousands of individuals with new insights into genetic relationships and pedigrees. High-density genetic maps are constructed efficiently and quantitative trait loci for valuable traits are readily associated with positional candidate genes and/or converted into diagnostic tests for breeders. We understand the species, geographical, and genomic origins of domesticated apple more precisely, as well as its relationship to wild relatives. The WGS has turbo-charged application of these classical research steps to crop improvement and drives innovative methods to achieve more durable, environmentally sound, productive, and consumer-desirable apple production. This review includes examples of basic and practical breakthroughs and challenges in using the apple WGSs. Recommendations for "what's next" focus on necessary upgrades to the genome sequence data pool, as well as for use of the data, to reach new frontiers in genomics-based scientific understanding of apple.

From markers to genome-based breeding in wheat

Recent technological advances in wheat genomics provide new opportunities to uncover genetic variation in traits of breeding interest and enable genome-based breeding to deliver wheat cultivars for the projected food requirements for 2050. There has been tremendous progress in development of whole-genome sequencing resources in wheat and its progenitor species during the last 5 years. High-throughput genotyping is now possible in wheat not only for routine gene introgression but also for high-density genome-wide genotyping. This is a major transition phase to enable genome-based breeding to achieve progressive genetic gains to parallel to projected wheat production demands. These advances have intrigued wheat researchers to practice less pursued analytical approaches which were not practiced due to the short history of genome sequence availability. Such approaches have been successful in gene discovery and breeding applications in other crops and animals for which genome sequences have been available for much longer. These strategies include, (i) environmental genome-wide association studies in wheat genetic resources stored in genbanks to identify genes for local adaptation by using agroclimatic traits as phenotypes, (ii) haplotype-based analyses to improve the statistical power and resolution of genomic selection and gene mapping experiments, (iii) new breeding strategies for genome-based prediction of heterosis patterns in wheat, and (iv) ultimate use of genomics information to develop more efficient and robust genome-wide genotyping platforms to precisely predict higher yield potential and stability with greater precision. Genome-based breeding has potential to achieve the ultimate objective of ensuring sustainable wheat production through developing high yielding, climate-resilient wheat cultivars with high nutritional quality.

Association of Genetic Variants of NLRP4 with Exacerbation of Asthma: The Effect of Smoking

Asthma exacerbation is induced by the interaction of genes and environmental factors such as cigarette smoke. NLRP4 counteracts the activity of the inflammasome, which is responsible for asthma exacerbation. In this study, we analyzed the association of single-nucleotide polymorphisms of NLRP4 with the annual rate of exacerbation and evaluated the additive effect of smoking in 1454 asthmatics. Asthmatics possessing the minor allele of rs1696718G > A had more frequent exacerbation episodes than those homozygous for the common allele (0.59 vs. 0.36/year) and the association was present only in current and ex-smokers. There was a significant interaction between the amount smoked and rs16986718 genotypes (p = 0.014) and a positive correlation between the number of annual exacerbation episodes and amount smoked only in rs16986718G > A AA homozygotes. The prevalence of frequent exacerbators (≥2 exacerbation episodes/year) was 2.5 times higher in rs16986718G > A minor allele homozygotes than in common allele homozygotes (12.0% vs. 5.9%). Furthermore, the prevalence was 6 times higher in rs16986718G > A minor allele homozygotes who were current and ex-smokers than in nonsmokers (25.6% vs. 4.1%). The minor allele of rs16986718G > A in NLRP4 may be a genetic marker that predicts asthma exacerbation in adult asthmatics who smoke.

Microhaplotypes in forensic genetics

Microhaplotype loci (microhaps, MHs) are a novel type of molecular marker of less than 300 nucleotides, defined by two or more closely linked SNPs associated in multiple allelic combinations. The value of these markers is enhanced by massively parallel sequencing (MPS), which allows the sequencing of both parental haplotypes at each of the many multiplexed loci. This review describes the features of these multi-SNP markers and documents their value in forensic genetics, focusing on individualization, biogeographic ancestry inference, and mixture deconvolution. Foreseeable applications also include missing person identification, relationship testing, and medical diagnostic applications. The technique is not restricted to humans.