Development and utility of SSR markers based on Brassica sp. whole-genome in triangle of U

Introduction

Simple sequence repeats (SSR), also known as microsatellites, are crucial molecular markers in both animals and plants. Despite extensive previous research on SSRs, the development of microsatellite markers in Brassica crops remains limited and inefficient.

Methods

Krait software was used to identify microsatellites by genome-wide and marker development based on three recently sequenced basic species of Brassica crops in the triangle of U (Brassica rapa, B. nigra and B. oleracea), as well as three allotetraploids (B. juncea, B. napus and B. carinata) using public databases. Subsequently, the primers and the characteristics of microsatellites for most of them were accordingly designed on each chromosome of each of the six Brassica species, and their physical locations were identified,and the cross-transferability of primers have been carried out. In addition, a B-genome specific SSR marker was screened out.

Results

A total of 79341, 92089, 125443, 173964, 173604, and 222160 SSR loci have been identified from the whole genome sequences of Brassica crops within the triangle of U crops, B. rapa (AA), B. nigra (BB), B. oleracea (CC), B. napus (AACC), B. juncea (AABB) and B. carinata (BBCC), respectively. Comparing the number distribution of the three allotetraploid SSR loci in the three subgenomes AA, BB and CC, results indicate that the allotetraploid species have significant reduction in the number of SSR loci in the genome compared with their basic diploid counterparts. Moreover, we compared the basic species with their corresponding varieties, and found that the microsatellite characters between the allotetraploids and their corresponding basic species were very similar or almost identical. Subsequently, each of the 40 SSR primers was employed to investigate the polymorphism potential of B. rapa (85.27%), B. nigra (81.33%) and B. oleracea (73.45%), and B. rapa was found to have a higher cross-transfer rate among the basic species in the triangle of U. Meanwhile, a B-genome specific SSR marker, BniSSR23228 possessing the (AAGGA)3 sequence characteristics was obtained, and it located in chromosome B3 with a total length of 97 bp.

Discussion

In this study, results suggest that the pattern of distribution may be highly conserved during the differentiation of basic Brassica species and their allotetraploid counterparts. Our data indicated that the allotetraploidization process resulted in a significant reduction in SSR loci in the three subgenomes AA, BB and CC. The reasons may be partial gene dominated chromosomal homologous recombination and rearrangement during the evolution of basic diploid species into allotetraploids. This study provides a basis for future genomics and genetic research on the relatedness of Brassica species.

Genome-wide copy number variant screening of Saudi schizophrenia patients reveals larger deletions in cases versus controls

Introduction

Genome-wide association studies have discovered common polymorphisms in regions associated with schizophrenia. No genome-wide analyses have been performed in Saudi schizophrenia subjects.

Methods

Genome-wide genotyping data from 136 Saudi schizophrenia cases and 97 Saudi controls in addition to 4,625 American were examined for copy number variants (CNVs). A hidden Markov model approach was used to call CNVs.

Results

CNVs in schizophrenia cases were twice as large on average than CNVs in controls (p = 0.04). The analyses focused on extremely large >250 kilobases CNVs or homozygous deletions of any size. One extremely large deletion was noted in a single case (16.5 megabases on chromosome 10). Two cases had an 814 kb duplication of chromosome 7 spanning a cluster of genes, including circadian-related loci, and two other cases had 277 kb deletions of chromosome 9 encompassing an olfactory receptors gene family. CNVs were also seen in loci previously associated with schizophrenia, namely a 16p11 proximal duplication and two 22q11.2 deletions.

Discussion

Runs of homozygosity (ROHs) were analyzed across the genome to investigate correlation with schizophrenia risk. While rates and sizes of these ROHs were similar in cases and controls, we identified 10 regions where multiple cases had ROHs and controls did not.

Cytochrome P450 Superfamily: Evolutionary and Functional Divergence in Sorghum (Sorghum bicolor) Stress Resistance

Cytochrome P450 (CYP) genes encode enzymes that catalyze various growth-, development-, and stress-related reactions. Sorghum (Sorghum bicolor) is a type of C4 plant and an important cash crop. However, systematic identification and analysis of functional differentiation and evolution of CYP genes have not been carried out in this species. In the present study, we revealed that the sorghum genome contains 351 CYP genes, which can be divided into nine classes. These genes are from ancestors and repeated segments, rather than tandem repeats. Based on collinearity results, a large number of CYPs were extended before cotyledon differentiation, during the emergence of Gramineae, suggesting that genomewide duplication events and stress adaptation processes were important for the expansion of CYP genes. Their gene structure and motifs contain conserved regions and include various changes and loci. The expression characteristics and functional annotation of CYP genes indicated tissue specificity and selective expression. Overall, we identified all CYP genes in the sorghum genome and preliminarily explored their naming, structure, evolution, expression, and functional differentiation. The results advanced our understanding of plant gene family evolution and functional differentiation.

MITF controls the TCA cycle to modulate the melanoma hypoxia response

In response to the dynamic intra-tumor microenvironment, melanoma cells adopt distinct phenotypic states associated with differential expression of the microphthalmia-associated transcription factor (MITF). The response to hypoxia is driven by hypoxia-inducible transcription factors (HIFs) that reprogram metabolism and promote angiogenesis. HIF1α indirectly represses MITF that can activate HIF1α expression. Although HIF and MITF share a highly related DNA-binding specificity, it is unclear whether they co-regulate subset of target genes. Moreover, the genomewide impact of hypoxia on melanoma and whether melanoma cell lines representing different phenotypic states exhibit distinct hypoxic responses is unknown. Here we show that three different melanoma cell lines exhibit widely different hypoxia responses with only a core 23 genes regulated in common after 12 hr in hypoxia. Surprisingly, under hypoxia MITF is transiently up-regulated by HIF1α and co-regulates a subset of HIF targets including VEGFA. Significantly, we also show that MITF represses itself and also regulates SDHB to control the TCA cycle and suppress pseudo-hypoxia. Our results reveal a previously unsuspected role for MITF in metabolism and the network of factors underpinning the hypoxic response in melanoma.

Characterization and Expression Patterns of Auxin Response Factors in Wheat

Auxin response factors (ARFs) are important transcription factors involved in both the auxin signaling pathway and the regulatory development of various plant organs. In this study, 23 TaARF members encoded by a total of 68 homeoalleles were isolated from 18 wheat chromosomes (excluding chromosome 4). The TaARFs, including their conserved domains, exon/intron structures, related microRNAs, and alternative splicing (AS) variants, were then characterized. Phylogenetic analysis revealed that members of the TaARF family share close homology with ARFs in other grass species. qRT-PCR analyses revealed that 20 TaARF members were expressed in different organs and tissues and that the expression of some members significantly differed in the roots, stems, and leaves of wheat seedlings in response to exogenous auxin treatment. Moreover, protein network analyses and co-expression results showed that TaTIR1-TaARF15/18/19-TaIAA13 may interact at both the protein and genetic levels. The results of subsequent evolutionary analyses showed that three transcripts of TaARF15 in the A subgenome of wheat exhibited high evolutionary rate and underwent positive selection. Transgenic analyses indicated that TaARF15-A.1 promoted the growth of roots and leaves of Arabidopsis thaliana and was upregulated in the overexpression plants after auxin treatment. Our results will provide reference information for subsequent research and utilization of the TaARF gene family.

Anxiety genetics: Dispatches from the frontier

Anxiety disorders are the most common class of psychiatric disorders and incur an enormous burden in terms of economic costs, disability and personal suffering. Despite their public health importance and documented heritability, genetic research aimed at identifying the genetic contributions to these disorders has had limited success, particularly in comparison to recent advances in the genetics of other major psychiatric disorders. The major factor contributing to this lagging progress has been a dearth of well-powered genomic studies. As a result, the genetic basis of anxiety disorders remains a largely uncharted frontier. In this Special Issue, we bring together dispatches from this frontier, reflecting a range of strategies being pursued including genomewide common variant association studies, exome sequencing, genetic prediction of treatment response, cross-disorder genetic epidemiology, and the prospects for dissecting underlying anxiety domains using the RDoC framework. Overall, these efforts underscore the genetic and phenotypic complexity of pathologic anxiety and the urgent need for more powerful and comprehensive analyses if we are to make real progress in advancing this important frontier. © 2017 Wiley Periodicals, Inc.

Clinical validation of a noninvasive prenatal test for genomewide detection of fetal copy number variants

BACKGROUND

Current cell-free DNA assessment of fetal chromosomes does not analyze and report on all chromosomes. Hence, a significant proportion of fetal chromosomal abnormalities are not detectable by current noninvasive methods. Here we report the clinical validation of a novel noninvasive prenatal test (NIPT) designed to detect genomewide gains and losses of chromosomal material ≥7 Mb and losses associated with specific deletions <7 Mb.

OBJECTIVE

The objective of this study is to provide a clinical validation of the sensitivity and specificity of a novel NIPT for detection of genomewide abnormalities.

STUDY DESIGN

This retrospective, blinded study included maternal plasma collected from 1222 study subjects with pregnancies at increased risk for fetal chromosomal abnormalities that were assessed for trisomy 21 (T21), trisomy 18 (T18), trisomy 13 (T13), sex chromosome aneuploidies (SCAs), fetal sex, genomewide copy number variants (CNVs) ≥7 Mb, and select deletions <7 Mb. Performance was assessed by comparing test results with findings from G-band karyotyping, microarray data, or high coverage sequencing.

RESULTS

Clinical sensitivity within this study was determined to be 100% for T21 (95% confidence interval [CI], 94.6-100%), T18 (95% CI, 84.4-100%), T13 (95% CI, 74.7-100%), and SCAs (95% CI, 84-100%), and 97.7% for genomewide CNVs (95% CI, 86.2-99.9%). Clinical specificity within this study was determined to be 100% for T21 (95% CI, 99.6-100%), T18 (95% CI, 99.6-100%), and T13 (95% CI, 99.6-100%), and 99.9% for SCAs and CNVs (95% CI, 99.4-100% for both). Fetal sex classification had an accuracy of 99.6% (95% CI, 98.9-99.8%).

CONCLUSION

This study has demonstrated that genomewide NIPT for fetal chromosomal abnormalities can provide high resolution, sensitive, and specific detection of a wide range of subchromosomal and whole chromosomal abnormalities that were previously only detectable by invasive karyotype analysis. In some instances, this NIPT also provided additional clarification about the origin of genetic material that had not been identified by invasive karyotype analysis.

Genomewide introgressive hybridization patterns in wild Atlantic salmon influenced by inadvertent gene flow from hatchery releases

Many salmonid fish populations are threatened by genetic homogenization, primarily due to introgressive hybridization with hatchery-reared conspecifics. By applying genomewide analysis using two molecular marker types (1986 SNPs and 17 microsatellites), we assessed the genetic impacts of inadvertent gene flow via straying from hatchery releases on wild populations of Atlantic salmon in the Gulf of Finland, Baltic Sea, over 16 years (1996-2012). Both microsatellites and SNPs revealed congruent population genetic structuring, indicating that introgression changed the genetic make-up of wild populations by increasing genetic diversity and reducing genetic divergence. However, the degree of genetic introgression varied among studied populations, being higher in the eastern part and lower in the western part of Estonia, which most likely reflects the history of past stocking activities. Using kernel smoothing and permutation testing, we detected considerable heterogeneity in introgression patterns across the genome, with a large number of regions exhibiting nonrandom introgression widely dispersed across the genome. We also observed substantial variation in nonrandom introgression patterns within populations, as the majority of genomic regions showing elevated or reduced introgression were not consistently detected among temporal samples. This suggests that recombination, selection and stochastic processes may contribute to complex nonrandom introgression patterns. Our results suggest that (i) some genomic regions in Atlantic salmon are more vulnerable to introgressive hybridization, while others show greater resistance to unidirectional gene flow; and (ii) the hybridization of previously separated populations leads to complex and dynamic nonrandom introgression patterns that most likely have functional consequences for indigenous populations.

European multiple sclerosis risk variants in the south Asian population

BACKGROUND

In less than a decade, genomewide association studies have identified over 100 single-nucleotide variants that are associated with increased risk of developing multiple sclerosis. However, since these studies have focused almost exclusively on European populations, it is unclear what role these variants might play in determining risk in other ethnic groups.

OBJECTIVE

To assess the effects of European multiple sclerosis-associated risk variants in the south Asian population.

METHODS

Using a combination of chip-based genotyping and next-generation sequencing, we have assessed 109 European-associated variants in a total of 270 cases and 555 controls from the south Asian population.

RESULTS

We found that two-thirds of the tested variants (72/109) showed over representation of the European risk allele in south Asian cases (p < 0.0003). In the rest of the Immunochip array, the most associated variant was rs7318477 which maps close to TNFSF13B, the gene for the B-cell-related protein BAFF.

CONCLUSION

Our data indicate substantial overlap in genetic risk architecture between Europeans and south Asians and suggest that the aetiology of the disease may be largely independent of ethnicity.

Aberrant sperm DNA methylation predicts male fertility status and embryo quality

OBJECTIVE

To evaluate whether male fertility status and/or embryo quality during in vitro fertilization (IVF) therapy can be predicted based on genomewide sperm deoxyribonucleic acid (DNA) methylation patterns.

DESIGN

Retrospective cohort study.

SETTING

University-based fertility center.

PATIENT(S)

Participants were 127 men undergoing IVF treatment (where any major female factor cause of infertility had been ruled out), and 54 normozoospermic, fertile men. The IVF patients were stratified into 2 groups: patients who had generally good embryogenesis and a positive pregnancy (n = 55), and patients with generally poor embryogenesis (n = 72; 42 positive and 30 negative pregnancies) after IVF.

INTERVENTION(S)

Genomewide sperm DNA methylation analysis was performed to measure methylation at >485,000 sites across the genome.

MAIN OUTCOME MEASURE(S)

A comparison was made of DNA methylation patterns of IVF patients vs. normozoospermic, fertile men.

RESULT(S)

Predictive models proved to be highly accurate in classifying male fertility status (fertile or infertile), with 82% sensitivity, and 99% positive predictive value. Hierarchic clustering identified clusters enriched for IVF patient samples and for poor-quality-embryo samples. Models built to identify samples within these groups, from neat samples, achieved positive predictive value ≥ 94% while identifying >one fifth of all IVF patient and poor-quality-embryo samples in each case. Using density gradient prepared samples, the same approach recovered 46% of poor-quality-embryo samples with no false positives.

CONCLUSION(S)

Sperm DNA methylation patterns differ significantly and consistently for infertile vs. fertile, normozoospermic men. In addition, DNA methylation patterns may be predictive of embryo quality during IVF.

Genome-wide linkage analysis of obsessive-compulsive disorder implicates chromosome 1p36

BACKGROUND

Obsessive-compulsive disorder (OCD) has a complex etiology involving both genetic and environmental factors. However, the genetic causes of OCD are largely unknown, despite the identification of several promising candidate genes and linkage regions.

METHODS

Our objective was to conduct genetic linkage studies of the type of OCD thought to have the strongest genetic etiology (i.e., childhood-onset OCD), in 33 Caucasian families with ≥2 childhood-onset OCD-affected individuals from the United States (n = 245 individuals with genotype data). Parametric and nonparametric genome-wide linkage analyses were conducted with Morgan and Merlin in these families using a selected panel of single nucleotide repeat polymorphisms from the Illumina 610-Quad Bead Chip. The initial analyses were followed by fine-mapping analyses in genomic regions with initial heterogeneity logarithm of odds (HLOD) scores of ≥2.0.

RESULTS

We identified five areas of interest (HLOD score ≥2) on chromosomes 1p36, 2p14, 5q13, 6p25, and 10p13. The strongest result was on chromosome 1p36.33-p36.32 (HLOD = 3.77, suggestive evidence for linkage after fine mapping). At this location, several of the families showed haplotypes co-segregating with OCD.

CONCLUSIONS

The results of this study represent the strongest linkage finding for OCD in a primary analysis to date and suggest that chromosome 1p36, and possibly several other genomic regions, may harbor susceptibility loci for OCD. Multiple brain-expressed genes lie under the primary linkage peak (approximately 4 megabases in size). Follow-up studies, including replication in additional samples and targeted sequencing of the areas of interest, are needed to confirm these findings and to identify specific OCD risk variants.

Further evidence for linkage of bipolar disorder to chromosomes 6 and 17 in a new independent pedigree series

OBJECTIVES

We have previously reported the results of a linkage analysis of bipolar disorder in an initial set of 20 pedigrees ascertained through collaboration among three sites. We now report the results of our genome-wide linkage analysis in an independent sample of 34 pedigrees segregating bipolar disorder.

METHODS

Families were ascertained through a bipolar I or II disorder proband for the presence of bipolar I disorder, bipolar II disorder, or recurrent major depression in at least two other family members. A total of 440 markers at an average spacing of 8 cM were genotyped in 229 family members using fluorescent methods.

RESULTS

Initial nonparametric analyses of chromosomes 6 and 17 provided evidence for a modest replication of linkage to these chromosomes previously reported in other studies. Additional analyses using multipoint parametric methods provided further evidence to support the 6q25 region with a heterogeneity logarithm of odds score of 3.28. Evidence from two-point parametric analyses also provides a modest replication of our previous findings of linkage to the 23 cM region of chromosome 22q13 in our original University of California, San Diego sample of 20 families and 57 families from the National Institute of Mental Health bipolar disorder sample.

CONCLUSIONS

Our results suggest replication of some reported linkage peaks, such as 6q25 and 17p12; however, other peaks from our own previous study, such as 5p15, 13q32, and 22q13, were either not replicated or were only modestly replicated in these analyses.