Copy Number Variations in Tilapia Genomes

Discovering the nature and pattern of genome variation is fundamental in understanding phenotypic diversity among populations. Although several millions of single nucleotide polymorphisms (SNPs) have been discovered in tilapia, the genome-wide characterization of larger structural variants, such as copy number variation (CNV) regions has not been carried out yet. We conducted a genome-wide scan for CNVs in 47 individuals from three tilapia populations. Based on 254 Gb of high-quality paired-end sequencing reads, we identified 4642 distinct high-confidence CNVs. These CNVs account for 1.9% (12.411 Mb) of the used Nile tilapia reference genome. A total of 1100 predicted CNVs were found overlapping with exon regions of protein genes. Further association analysis based on linear model regression found 85 CNVs ranging between 300 and 27,000 base pairs significantly associated to population types (R ² > 0.9 and P > 0.001). Our study sheds first insights on genome-wide CNVs in tilapia. These CNVs among and within tilapia populations may have functional effects on phenotypes and specific adaptation to particular environments.     

Genome-wide copy number variation study and gene expression analysis identify ABI3BP as a susceptibility gene for Kashin–Beck disease

Kashin–Beck disease (KBD) is a chronic osteochondropathy. In this study, we conducted the first genome-wide copy number variation study (GCNVS) of KBD totally involving 2,743 Chinese Han adults. GCNVS was first performed using Affymetrix Human SNP6.0 Arrays. The identified copy number variations (CNVs) were then replicated in an independent Chinese Han sample containing 1,026 subjects. SNP genotyping, CNV identification and quality control were implemented by Birdsuite. STRUCTURE and EIGENSTRAT were applied for controlling potential population stratification in the GCNVS. Association analysis was conducted using PLINK. Microarray and qRT-PCR were also conducted to compare the expression levels of the genes overlapping with identified CNVs between KBD patients and healthy controls. GCNVS found that CNV452 (P value = 7.78 × 10^?5) overlapping with ABI3BP gene was significantly associated with KBD. Replication association study observed that rs9850273 (P value = 0.008) and rs7613610 (P value = 0.021) in ABI3BP gene were significantly associated with KBD. Gene expression analysis also found that ABI3BP was up-regulated in KBD patients compared to healthy controls. Our results suggest that ABI3BP was a novel susceptibility gene for KBD.     

Identification of BACH2 as a susceptibility gene for Graves’ disease in the Chinese Han population based on a three-stage genome-wide association study

The BACH2 gene regulates B cell differentiation and function and has been reported to be a shared susceptibility gene for several autoimmune diseases. Our previous genome-wide association study (GWAS) indicated that several single nucleotide polymorphisms (SNPs) in the BACH2 gene are associated with Graves’ disease (GD) in the Chinese Han population; however, the association did not achieve genome-wide significance levels. Recently, this association of BACH2 with GD was confirmed in Caucasians in the UK population, but fine mapping in this region has not yet been reported. Here, we provide a refined analysis of a 331-kb region in the BACH2 gene, which harbors 359 SNPs, using GWAS data from 1,442 GD patients and 1,468 controls. The SNPs rs2474619 and rs9344996 were implied as the independent variants associated with GD by forward and two-locus logistic regression analysis. We genotyped eight out of 10 tagSNPs with P < 1 × 10^?3 in 3,508 GD patients and 3,209 controls, the results also showed that rs2474619 was independently associated with GD in the combined population from GWAS and the second stage (P = 1.81 × 10^?5). The rs2474619 and rs9344996 were further genotyped in the third stage cohorts, and rs2474619 showed evidence of association with GD at genome-wide significance levels in the combined population (P = 3.28 × 10^?8, odds ratio = 1.13). The association of rs9344996 with GD can be explained by its linkage to rs2474619 in the combined population. Our study clearly demonstrated that BACH2 is a susceptibility gene for GD in the Chinese Han population and further supported rs2474619, in intron 2 of BACH2, is the best association signal with GD. However, the mechanism by which BACH2 confers increased risk of GD requires further study.     

A GWA study reveals genetic loci for body conformation traits in Chinese Laiwu pigs and its implications for human BMI

Pigs share numerous physiological and phenotypic similarities with human and thus have been considered as a good model in nonrodent mammals for the study of genetic basis of human obesity. Researches on candidate genes for obesity traits have successfully identified some common genes between humans and pigs. However, few studies have assessed how many similarities exist between the genetic architecture of obesity in pigs and humans by large-scale comparative genomics. Here, we performed a genome-wide association study (GWAS) using the porcine 60 K SNP Beadchip for BMI and other four conformation traits at three different ages in a Chinese Laiwu pig population, which shows a large variability in fat deposition. In total, 35 SNPs were found to be significant at Bonferroni-corrected 5 % chromosome-wise level (P = 2.13 × 10^?5) and 88 SNPs had suggestive (P < 10^?4) association with the conformation traits. Some SNPs showed age-dependent association. Intriguingly, out of 32 regions associated with BMI in pigs, 18 were homologous with the loci for BMI in humans. Furthermore, five closest genes to GWAS peaks including HIF1AN, SMYD3, COX10, SLMAP, and GBE1 have been already associated with BMI in humans, which makes them very promising candidates for these QTLs. The result of GO analysis provided strong support to the fact that mitochondria and synapse play important roles in obesity susceptibility, which is consistent with previous findings on human obesity, and it also implicated new gene sets related to chromatin modification and Ig-like C2-type 5 domain. Therefore, these results not only provide new insights into the genetic architecture of BMI in pigs but also highlight that humans and pigs share the significant overlap of obesity-related genes.     

Genome-wide association study identified the human leukocyte antigen region as a novel locus for plasma beta-2 microglobulin

Beta-2 microglobulin (B2M) is a component of the major histocompatibility complex (MHC) class I molecule and has been studied as a biomarker of kidney function, cardiovascular diseases and mortality. Little is known about the genes influencing its levels directly or through glomerular filtration rate (GFR). We conducted a genome-wide association study of plasma B2M levels in 6738 European Americans from the Atherosclerosis Risk in Communities study to identify novel loci for B2M and assessed its association with known estimated GFR (eGFR) loci. We identified 2 genome-wide significant loci. One was in the human leukocyte antigen (HLA) region on chromosome 6 (lowest p value = 1.8 × 10^?23 for rs9264638). At this locus, 6 index SNPs accounted for 3.2 % of log(B2M) variance, and their association with B2M could largely be explained by imputed classical alleles of the MHC class I genes: HLA-A, HLA-B, or HLA-C. The index SNPs at this locus were not associated with eGFR based on serum creatinine (eGFRcr). The other locus of B2M was on chromosome 12 (rs3184504 at SH2B3, beta = 0.02, p value = 3.1 × 10^?8), which was previously implicated as an eGFR locus. In conclusion, although B2M is known to be a component of MHC class I molecule, the association between HLA class I alleles and plasma B2M levels in a community-based population is novel. The identification of the two novel loci for B2M extends our understanding of its metabolism and informs its use as a kidney filtration biomarker.     

Association between genome-wide association studies reported SNPs and pediatric-onset Crohn’s disease in Canadian children

A recent pediatric-focused genome-wide association study has implicated three novel susceptibility loci for Crohn’ disease (CD).We aimed to investigate whether the three recently reported and other previously reported genes/loci were also associated with CD in Canadian children. A case–control design was implemented at three pediatric gastroenterology clinics in Canada. Children <19 years of age with a confirmed diagnosis of CD were recruited along with controls. Single nucleotide polymorphisms (SNPs) in 19 reported genes/loci were genotyped. Associations between individual SNPs and CD were examined. A total of 563 cases and 553 controls were studied. The mean (±SD) age of the cases was 12.3 (±3.2) years. Most cases were male (56.0%), had ileo-colonic disease (L3 ± L4, 48.8%) and inflammatory behavior (B1 ± p, 87.9%) at diagnosis. Allelic association analysis (two-tailed) showed that 8 of the 19 targeted SNPs were significantly associated with overall susceptibility for CD. Associations with one additional SNP was borderline non-significant. Significantly associated SNPs included SNPs rs1250550 (p = 0.026) and rs8049439 (p = 0.04), recently reported to be specifically associated with pediatric-onset CD.Based on the results, we confirmed associations between two of the three novel pediatric-CD loci and other regions reported for associations with either pediatric and/or adult-onset CD.     

Metabolic heritability at birth: implications for chronic disease research

Recent genome-wide association studies of the adult human metabolome have identified genetic variants associated with relative levels of several acylcarnitines, which are important clinical correlates for chronic conditions such as type 2 diabetes and obesity. We have previously shown that these same metabolite levels are highly heritable at birth; however, no studies to our knowledge have examined genetic associations with these metabolites measured at birth. Here, we examine, in 743 newborns, 58 single nucleotide polymorphisms (SNPs) in 11 candidate genes previously associated with differing relative levels of short-chain acylcarnitines in adults. Six SNPs (rs2066938, rs3916, rs3794215, rs555404, rs558314, rs1799958) in the short-chain acyl-CoA dehydrogenase gene (ACADS) were associated with neonatal C4 levels. Most significant was the G allele of rs2066938, which was associated with significantly higher levels of C4 (P = 1.5 × 10^?29). This SNP explains 25 % of the variation in neonatal C4 levels, which is similar to the variation previously reported in adult C4 levels. There were also significant (P < 1 × 10^?4) associations between neonatal levels of C5-OH and SNPs in the solute carrier family 22 genes (SLC22A4 and SLC22A5) and the 3-methylcrotonyl-CoA carboxylase 1 gene (MCCC1). We have replicated, in newborns, SNP associations between metabolic traits and the ACADS and SLC22A4 genes observed in adults. This research has important implications not only for the identification of rare inborn errors of metabolism but also for personalized medicine and early detection of later life risks for chronic conditions.     

Increased paternal age and the influence on burden of genomic copy number variation in the general population

Genomic copy number variations (CNVs) and increased parental age are both associated with the risk to develop a variety of clinical neuropsychiatric disorders such as autism, schizophrenia and bipolar disorder. At the same time, it has been shown that the rate of transmitted de novo single nucleotide mutations is increased with paternal age. To address whether paternal age also affects the burden of structural genomic deletions and duplications, we examined various types of CNV burden in a large population sample from the Netherlands. Healthy participants with parental age information (n = 6,773) were collected at different University Medical Centers. CNVs were called with the PennCNV algorithm using Illumina genome-wide SNP array data. We observed no evidence in support of a paternal age effect on CNV load in the offspring. Our results were negative for global measures as well as several proxies for de novo CNV events in this unique sample. While recent studies suggest de novo single nucleotide mutation rate to be dominated by the age of the father at conception, our results strongly suggest that at the level of global CNV burden there is no influence of increased paternal age. While it remains possible that local genomic effects may exist for specific phenotypes, this study indicates that global CNV burden and increased father’s age may be independent disease risk factors.     

Candidate gene association mapping of Sclerotinia stalk rot resistance in sunflower (Helianthus annuus L.) uncovers the importance of COI1 homologs

Key message

Functional markers for Sclerotinia basal stalk rot resistance in sunflower were obtained using gene-level information from the model species Arabidopsis thaliana.

Abstract

Sclerotinia stalk rot, caused by Sclerotinia sclerotiorum, is one of the most destructive diseases of sunflower (Helianthus annuus L.) worldwide. Markers for genes controlling resistance to S. sclerotiorum will enable efficient marker-assisted selection (MAS). We sequenced eight candidate genes homologous to Arabidopsis thaliana defense genes known to be associated with Sclerotinia disease resistance in a sunflower association mapping population evaluated for Sclerotinia stalk rot resistance. The total candidate gene sequence regions covered a concatenated length of 3,791 bp per individual. A total of 187 polymorphic sites were detected for all candidate gene sequences, 149 of which were single nucleotide polymorphisms (SNPs) and 38 were insertions/deletions. Eight SNPs in the coding regions led to changes in amino acid codons. Linkage disequilibrium decay throughout the candidate gene regions declined on average to an r ² = 0.2 for genetic intervals of 120 bp, but extended up to 350 bp with r ² = 0.1. A general linear model with modification to account for population structure was found the best fitting model for this population and was used for association mapping. Both HaCOI1-1 and HaCOI1-2 were found to be strongly associated with Sclerotinia stalk rot resistance and explained 7.4 % of phenotypic variation in this population. These SNP markers associated with Sclerotinia stalk rot resistance can potentially be applied to the selection of favorable genotypes, which will significantly improve the efficiency of MAS during the development of stalk rot resistant cultivars.     

Genome-wide association analysis of host genotype and plastic wing morphological variation of an endoparasitoid wasp <Emphasis Type="Italic">Asobara japonica</Emphasis> (Hymenoptera: Braconidae)

Accumulating evidence suggests that genotype of host insects influences the development of koinobiont endoparasitoids. Although there are many potential genetic variations that lead to the internal body environmental variations of host insects, association between the host genotype and the parasitoid development has not been examined in a genome-wide manner. In the present study, we used highly inbred whole genome sequenced strains of Drosophila melanogaster to associate single nucleotide polymorphisms (SNPs) of host flies with morphological traits of Asobara japonica, a larval-pupal parasitoid wasp that infected those hosts. We quantified the outline shape of the forewings of A. japonica with two major principal components (PC1 and PC2) calculated from Fourier coefficients obtained from elliptic Fourier analysis. We also quantified wing size and estimated wasp survival. We then examined the association between the PC scores, wing size and 1,798,561 SNPs and  the association between the estimated wasp survival and 1,790,544 SNPs. As a result, we obtained 22, 24 and 14 SNPs for PC1, PC2 and wing size and four SNPs for the estimated survival with P values smaller than 10^?5. Based on the location of the SNPs, 12, 17, 11 and five protein coding genes were identified as potential candidates for PC1, PC2, wing size and the estimated survival, respectively. Based on the function of the candidate genes, it is suggested that the host genetic variation associated with the cell growth and morphogenesis may influence the wasp’s morphogenetic variation.     

Ambient Temperature is A Strong Selective Factor Influencing Human Development and Immunity

Solar energy, which is essential for the origin and evolution of all life forms on Earth, can be objectively recorded through attributes such as climatic ambient temperature (CAT), ultraviolet radiation (UVR), and sunlight duration (SD). These attributes have specific geographical variations and may cause different adaptation traits. However, the adaptation profile of each attribute and the selective role of solar energy as a whole during human evolution remain elusive. Here, we performed a genome-wide adaptation study with respect to CAT, UVR, and SD using the Human Genome Diversity Project-Centre Etude Polymorphism Humain (HGDP-CEPH) panel data. We singled out CAT as the most important driving force with the highest number of adaptive loci (6 SNPs at the genome-wide 1 × 10⁻⁷ level; 401 at the suggestive 1 × 10⁻⁵ level). Five of the six genome-wide significant adaptation SNPs were successfully replicated in an independent Chinese population (N = 1395). The corresponding 316 CAT adaptation genes were mostly involved in development and immunity. In addition, 265 (84%) genes were related to at least one genome-wide association study (GWAS)-mapped human trait, being significantly enriched in anthropometric loci such as those associated with body mass index (χ²; P < 0.005), immunity, metabolic syndrome, and cancer (χ²; P < 0.05). For these adaptive SNPs, balancing selection was evident in Euro-Asians, whereas obvious positive and/or purifying selection was observed in Africans. Taken together, our study indicates that CAT is the most important attribute of solar energy that has driven genetic adaptation in development and immunity among global human populations. It also supports the non-neutral hypothesis for the origin of disease-predisposition alleles in common diseases.     

Genome-wide association study for soybean cyst nematode resistance in Chinese elite soybean cultivars

Soybean cyst nematode (SCN) (Heterodera glycines Ichinohe) is a highly recalcitrant endoparasite of soybean roots, causing more yield loss than any other pest. To identify quantitative trait loci (QTL) controlling resistance to SCN (HG type 2.5.7, race 1), a genome-wide association study (GWAS) was performed. The association panel, consisting of 120 Chinese soybean cultivars, was genotyped with 7189 single nucleotide polymorphism (SNPs). A total of 6204 SNPs with minor allele frequency >0.05 were used to estimate linkage disequilibrium (LD) and population structure. The mean level of LD measured by r ² declined very rapidly to half its maximum value (0.51) at 220 kb. The overall population structure was approximately coincident with geographic origin. The GWAS results identified 13 SNPs in 7 different genomic regions significantly associated with SCN resistance. Of these, three SNPs were localized in previously mapped QTL intervals, including rhg1 and Rhg4. The GWAS results also detected 10 SNPs in 5 different genomic regions associated with SCN resistance. The identified loci explained an average of 95.5% of the phenotypic variance. The proportion of phenotypic variance was due to additive genetic variance of the validated SNPs. The present study identified multiple new loci and refined chromosomal regions of known loci associated with SCN resistance. The loci and trait-associated SNPs identified in this study can be used for developing soybean cultivars with durable resistance against SCN.     

Genome-wide association study identifies loci and candidate genes for meat quality traits in Simmental beef cattle

Improving meat quality is the best way to enhance profitability and strengthen competitiveness in beef industry. Identification of genetic variants that control beef quality traits can help breeders design optimal breeding programs to achieve this goal. We carried out a genome-wide association study for meat quality traits in 1141 Simmental cattle using the Illumina Bovine HD 770K SNP array to identify the candidate genes and genomic regions associated with meat quality traits for beef cattle, including fat color, meat color, marbling score, longissimus muscle area, and shear force. In our study, we identified twenty significant single-nucleotide polymorphisms (SNPs) (p < 1.47 × 10^?6) associated with these five meat quality traits. Notably, we observed several SNPs were in or near eleven genes which have been reported previously, including TMEM236, SORL1, TRDN, S100A10, AP2S1, KCTD16, LOC506594, DHX15, LAMA4, PREX1, and BRINP3. We identified a haplotype block on BTA13 containing five significant SNPs associated with fat color trait. We also found one of 19 SNPs was associated with multiple traits (shear force and longissimus muscle area) on BTA7. Our results offer valuable insights to further explore the potential mechanism of meat quality traits in Simmental beef cattle.     

Genome-wide scan revealed that polymorphisms in the PNPLA3, SAMM50, and PARVB genes are associated with development and progression of nonalcoholic fatty liver disease in Japan

We examined the genetic background of nonalcoholic fatty liver disease (NAFLD) in the Japanese population, by performing a genome-wide association study (GWAS). For GWAS, 392 Japanese NAFLD subjects and 934 control individuals were analyzed. For replication studies, 172 NAFLD and 1,012 control subjects were monitored. After quality control, 261,540 single-nucleotide polymorphisms (SNPs) in autosomal chromosomes were analyzed using a trend test. Association analysis was also performed using multiple logistic regression analysis using genotypes, age, gender and body mass index (BMI) as independent variables. Multiple linear regression analyses were performed to evaluate allelic effect of significant SNPs on biochemical traits and histological parameters adjusted by age, gender, and BMI. Rs738409 in the PNPLA3 gene was most strongly associated with NAFLD after adjustment (P = 6.8 × 10^?14, OR = 2.05). Rs2896019, and rs381062 in the PNPLA3 gene, rs738491, rs3761472, and rs2143571 in the SAMM50 gene, rs6006473, rs5764455, and rs6006611 in the PARVB gene had also significant P values (<2.0 × 10^?10) and high odds ratios (1.84–2.02). These SNPs were found to be in the same linkage disequilibrium block and were associated with decreased serum triglycerides and increased aspartate aminotransferase (AST) and alanine aminotransferase (ALT) in NAFLD patients. These SNPs were associated with steatosis grade and NAFLD activity score (NAS). Rs738409, rs2896019, rs738491, rs6006473, rs5764455, and rs6006611 were associated with fibrosis. Polymorphisms in the SAMM50 and PARVB genes in addition to those in the PNPLA3 gene were observed to be associated with the development and progression of NAFLD.     

Construction of a genetic map based on high-throughput SNP genotyping and genetic mapping of a TuMV resistance locus in Brassica rapa

Brassica rapa is a member of the Brassicaceae family and includes vegetables and oil crops that are cultivated worldwide. The introduction of durable resistance against turnip mosaic virus (TuMV) into agronomically important cultivars has been a significant challenge for genetic and horticultural breeding studies of B. rapa. Based on our previous genome-wide analysis of DNA polymorphisms between the TuMV-resistant doubled haploid (DH) line VC40 and the TuMV-susceptible DH line SR5, we constructed a core genetic map of the VCS-13M DH population, which is composed of 83 individuals derived from microspore cultures of a F₁ cross between VC40 and SR5, by analyzing the segregation of 314 sequence-characterized genetic markers. The genetic markers correspond to 221 SNPs and 31 InDels of genes as well as 62 SSRs, covering 1,115.9 cM with an average distance of 3.6 cM between the adjacent marker loci. The alignment and orientation of the constructed map showed good agreement with the draft genome sequence of Chiifu, thus providing an efficient strategy to map genic sequences. Using the genetic map, a novel dominant TuMV resistance locus (TuMV-R) in the VCS-13M DH population was identified as a 0.34 Mb region in the short arm of chromosome A6 in which four CC–NBS–LRR resistance genes and two pathogenesis-related-1 genes reside. The genetic map developed in this study can play an important role in the genetic study of TuMV resistance and the molecular breeding of B. rapa.     

Novel SNPs in the bovine ADIPOQ and PPARGC1A genes are associated with carcass traits in Hanwoo (Korean cattle)

Adiponectin (ADIPOQ) modulates several biological processes including energy homeostasis, glucose and lipid metabolism. The bovine ADIPOQ gene was located near the QTL affecting marbling, ribeye muscle area and fat thickness on BTA1. The gene encoding peroxisome proliferator-activated receptor-γ coactivator-1α (PPARGC1A) was located within the QTL region of the traits on BTA6. Moreover, its protein product has various biological functions such as cellular energy homeostasis, including adaptive thermogenesis, adipogenesis and gluconeogenesis. Therefore, the ADIPOQ and PPARGC1A genes are a positional and functional candidate gene for carcass traits in beef cattle. The objectives of this study were to identify polymorphisms in the bovine ADIPOQ and PPARGC1A genes, to evaluate their associations with carcass traits in Hanwoo (Korean cattle) population. We identified nine SNPs in the ADIPOQ gene. Two SNPs (DQ156119: g.1436T > C and DQ156119: g.1454A > G) in the promoter region were recognized as new SNPs identified in Hanwoo. Association analysis indicated that the g.1454A > G SNP genotype was significantly associated with effects on LMA (P = 0.004) and BF (P = 0.021). The ADIPOQ haplotype was also found to have significant effect on the LMA. In the PPARGC1A gene, we identified 11 SNPs in the two unexplored regions (intron 3 and 5). Among them, seven SNPs were located in intron 3 and four SNPs were located in intron 5. Of these 11 putative novel SNPs, two SNPs (AY839822: g.292C > T and AY839823: g.1064C > T) with minor allele frequency (MAF) > 0.20 were examined for associations with carcass traits. The association analysis revealed that both SNPs in PPARGC1A gene were significantly associated with LMA (P < 0.05). These findings suggest that the SNPs of bovine ADIPOQ and PPARGC1A genes may be a useful molecular marker for selection of carcass traits in Hanwoo.     

Single nucleotide polymorphism markers for rapid detection of the <Emphasis Type="Italic">Rsv4</Emphasis> locus for soybean mosaic virus resistance in diverse germplasm

Soybean mosaic virus (SMV) causes a substantial decrease in soybean yield and reduction of seed quality. The most effective management strategy to control the virus is the deployment of host resistance. Seven SMV strains and three independent multi-allelic loci for SMV resistance have been identified previously. The goal of this research was to detect single nucleotide polymorphisms (SNPs) associated with SMV resistance at the Rsv4 locus. Ten soybean accessions, with confirmed resistance genes, were used for sequencing the candidate gene Glyma.02g121400. Alignment of these sequences revealed three SNPs displaying 100% consistency for genotypes carrying the Rsv4 gene. These SNPs were applied for a rapid screen of diverse soybean germplasm using the Sequenom iPLEX Gold platform, phenotyped with SMV-G1 and G7 strains to determine phenotype and classified into several groups carrying the proposed R-gene. The population of V94-5152 (Rsv4) × Lee 68 (rsv) was screened using novel SNPs to create a genetic map with improved resolution to determine the location of the Rsv4. To observe the recombination frequencies within the population, three additional SNPs on both sides of the Glyma.02g121400 gene were added. A linkage map revealed a distance of 3.6 cM between the Rsv4 locus and the closest SNP, thus shifting the putative Rsv4 region downstream on chromosome 2. With this region, five candidate genes have been proposed. The genomic position of the discovered SNPs, linked to the Rsv4, could increase screening precision and accelerate breeding efforts to develop multi-strain-resistant crops.     

Novel quantitative trait loci for central corneal thickness identified by candidate gene analysis of osteogenesis imperfecta genes

Osteogenesis imperfecta (OI) is a rare connective tissue disorder caused by mutations in the type I collagen genes, COL1A1 and COL1A2, and is characterised by low bone mass and bone fragility. In this study, we explored the relationship between type 1 collagen genes and the quantitative trait central corneal thickness (CCT). CCT was measured in a cohort of 28 Australian type I OI patients and mean CCT was found to be significantly lower compared to a normal population (P < 0.001). We then investigated CCT and corneal collagen fibril diameter and density in a mouse model of OI with a col1a2 mutation. Mean CCT was significantly lower in mutant mice (P = 0.002), as was corneal collagen fibril diameter (P = 0.034), whilst collagen fibril density was significantly greater in mutants (P = 0.034). Finally, we conducted a genetic study to determine whether common single nucleotide polymorphisms (SNPs) in COL1A1 and COL1A2 are associated with CCT variation in the normal human population. Polymorphism rs2696297 (P = 0.003) in COL1A1 and a three SNP haplotype in COL1A2 (P = 0.007) were all significantly associated with normal CCT variation. These data implicate type 1 collagen in the determination of CCT in both OI patients and normal individuals. This provides the first evidence of quantitative trait loci that influence CCT in a normal population and has potential implications for investigating genes involved in glaucoma pathogenesis, a common eye disease in which the severity and progression is influenced by CCT.     

Genetic markers of ovarian follicle number and menopause in women of multiple ethnicities

Oocyte loss has a significant impact on fertility and somatic health. Yet, we know little about factors that impact this process. We sought to identify genetic variants associated with ovarian reserve (oocyte number as measured by antral follicle count, AFC). Based on recently published genome-wide scans that identified loci associated with age of menopause, we also sought to test our hypothesis that follicle number and menopausal age share underlying genetic associations. We analyzed menopause-related variants for association with follicle number in an independent population of approximately 450 reproductive-aged women of European and African ancestry; these women were assessed for AFC, anthropometric, clinical, and lifestyle factors. One SNP strongly associated with later menopausal age in Caucasian women (+1.07?±?0.11?years) in previous work was also associated with higher follicle counts in Caucasians (+2.79?±?1.67 follicles) in our study. This variant is within the Minichromosome Maintenance Complex Component 8 (MCM8) gene, which we found was expressed within oocytes in follicles of the human ovary. In genome-wide scans of AFC, we also identified one marginally genome-wide and several nominally significant SNPs within several other genes associated with follicle number in both ethnic groups. Further, there were overlapping variants associated with multiple ovarian reserve markers (AFC, serum hormone levels, menopausal age). This study provides the first evidence for direct genetic associations underlying both follicle number and menopause and identifies novel candidate genes. Genetic variants associated with ovarian reserve may facilitate discovery of genetic markers to predict reproductive health and lifespan in women.